Spatial proteomic profiling of tumor and stromal compartments in non-small-cell lung cancer identifies signatures associated with overall survival.

Objectives: Non-small-cell lung carcinoma (NSCLC) is the most prevalent and lethal form of lung cancer. The need for biomarker-informed stratification of targeted therapies has underpinned the need to uncover the underlying properties of the tumor microenvironment (TME) through high-plex quantitative assays. Methods: In this study, we profiled resected NSCLC tissues from 102 patients by targeted spatial proteomics of 78 proteins across tumor, immune activation, immune cell typing, immune-oncology, drug targets, cell death and PI3K/AKT modules to identify the tumor and stromal signatures associated with overall survival (OS). Results: Survival analysis revealed that stromal CD56 (HR = 0.384, P = 0.06) and tumoral TIM3 (HR = 0.703, P = 0.05) were associated with better survival in univariate Cox models. In contrast, after adjusting for stage, BCLXL (HR = 2.093, P = 0.02) and cleaved caspase 9 (HR = 1.575, P = 0.1) negatively influenced survival. Delta testing indicated the protective effect of TIM-3 (HR = 0.614, P = 0.04) on OS. In multivariate analysis, CD56 (HR = 0.172, P = 0.001) was associated with better survival in the stroma, while B7.H3 (HR = 1.72, P = 0.008) was linked to poorer survival in the tumor. Conclusions: Deciphering the TME using high-plex spatially resolved methods is giving us new insights into compartmentalised tumor and stromal protein signatures associated with clinical endpoints in NSCLC.

Spatial resolution of the head and neck cancer tumor microenvironment to identify tumor and stromal features associated with therapy response.

Head and neck cancer (HNC) is the seventh most common cancer globally, resulting in 440 000 deaths per year. While there have been advancements in chemoradiotherapy and surgery, relapse occurs in more than half of HNCs, and these patients have a median survival of 10 months and a 2-year survival of < 20%. Only a subset of patients displays durable benefits from immunotherapies in metastatic and recurrent HNC, making it critical to understand the tumor microenvironment (TME) underpinning therapy responses in HNC. To recognize biological differences within the TME that may be predictive of immunotherapy response, we applied cutting-edge geospatial whole-transcriptome profiling (NanoString GeoMx Digital Spatial Profiler) and spatial proteomics profiling (Akoya PhenoCycler-Fusion) on a tumor microarray consisting of 25 cores from 12 patients that included 4 immunotherapy-unresponsive (8 cores) and 2 immunotherapy-responsive patients (5 cores), as well as 6 immunotherapy naïve patients (12 cores). Through high-plex, regional-based transcriptomic mapping of the tumor and TME, pathways involved with the complement system and hypoxia were identified to be differentially expressed in patients who went on to experience a poor immunotherapy response. Single-cell, targeted proteomic analysis found that immune cell infiltration of the cancer cell mass and interactions of CD8 T cells with tumor and other immune cells were associated with positive immunotherapy response. The relative abundance of specific tumor phenotypes and their interactions with various immune cells was identified to be different between response groups. This study demonstrates how spatial transcriptomics and proteomics can resolve novel alterations in the TME of HNC that may contribute to therapy sensitivity and resistance.

In situ single-cell profiling sheds light on IFI27 localisation during SARS-CoV-2 infection.

The utilization of single-cell resolved spatial transcriptomics to delineate immune responses during SARS-CoV-2 infection was able to identify M1 macrophages to have elevated expression of IFI27 in areas of infection.

The in situ transcriptomic landscape of breast tumour-associated and normal adjacent endothelial cells.

BACKGROUND AND AIMS: Triple Negative Breast Cancer (TNBC) is associated with increased angiogenesis, which is known to aid tumour growth and metastasis. Anti-angiogenic therapies that have been developed to target this feature have mostly generated disappointing clinical results. Further research into targeted approaches is limited by a lack of understanding of the in situ molecular profile of tumour-associated vasculature. In this study, we aimed to understand the differences in the molecular profiles of tumour endothelial cells vs normal-adjacent endothelial cells in TNBC tissues. METHOD: We have applied unbiased whole transcriptome spatial profiling of in situ gene expressions of endothelial cells localized in full-face patient TNBC tissues (n = 4) and normal-adjacent regions of the same patient breast tissues. RESULTS: Our comparative analysis revealed that 2412 genes were differentially expressed (padj < 0.05) between the tumour endothelial cells and normal-adjacent endothelial cells. Pathway enrichment showed the enrichment of gene sets related to cell-cell, cell-ECM adhesion, chromatin organization and remodeling, and protein-DNA complex subunit organization. CONCLUSION: Overall, the results revealed unique molecular profiles and signalling pathways of tumour-associated vasculature, which is a critical step towards larger cohort studies investigating potential targets for TNBC prognosis and anti-angiogenic treatments.

Metabolomics at the tumor microenvironment interface: Decoding cellular conversations.

Cancer heterogeneity remains a significant challenge for effective cancer treatments. Altered energetics is one of the hallmarks of cancer and influences tumor growth and drug resistance. Studies have shown that heterogeneity exists within the metabolic profile of tumors, and personalized-combination therapy with relevant metabolic interventions could improve patient response. Metabolomic studies are identifying novel biomarkers and therapeutic targets that have improved treatment response. The spatial location of elements in the tumor microenvironment are becoming increasingly important for understanding disease progression. The evolution of spatial metabolomics analysis now allows scientists to deeply understand how metabolite distribution contributes to cancer biology. Recently, these techniques have spatially resolved metabolite distribution to a subcellular level. It has been proposed that metabolite mapping could improve patient outcomes by improving precision medicine, enabling earlier diagnosis and intraoperatively identifying tumor margins. This review will discuss how altered metabolic pathways contribute to cancer progression and drug resistance and will explore the current capabilities of spatial metabolomics technologies and how these could be integrated into clinical practice to improve patient outcomes.

A robust platform for integrative spatial multi-omics analysis to map immune responses to SARS-CoV-2 infection in lung tissues.

The SARS-CoV-2 (COVID-19) virus has caused a devastating global pandemic of respiratory illness. To understand viral pathogenesis, methods are available for studying dissociated cells in blood, nasal samples, bronchoalveolar lavage fluid and similar, but a robust platform for deep tissue characterization of molecular and cellular responses to virus infection in the lungs is still lacking. We developed an innovative spatial multi-omics platform to investigate COVID-19-infected lung tissues. Five tissue-profiling technologies were combined by a novel computational mapping methodology to comprehensively characterize and compare the transcriptome and targeted proteome of virus infected and uninfected tissues. By integrating spatial transcriptomics data (Visium, GeoMx and RNAScope) and proteomics data (CODEX and PhenoImager HT) at different cellular resolutions across lung tissues, we found strong evidence for macrophage infiltration and defined the broader microenvironment surrounding these cells. By comparing infected and uninfected samples, we found an increase in cytokine signalling and interferon responses at different sites in the lung and showed spatial heterogeneity in the expression level of these pathways. These data demonstrate that integrative spatial multi-omics platforms can be broadly applied to gain a deeper understanding of viral effects on cellular environments at the site of infection and to increase our understanding of the impact of SARS-CoV-2 on the lungs.

Portable NMR for quantification of breast density in vivo: Proof-of-concept measurements and comparison with quantitative MRI.

Mammographic Density (MD) is the degree of radio-opacity of the breast in an X-ray mammogram. It is determined by the Fibroglandular: Adipose tissue ratio. MD has major implications in breast cancer risk and breast cancer chemoprevention. This study aimed to investigate the feasibility of accurate, low-cost quantification of MD in vivo without ionising radiation. We used single-sided portable nuclear magnetic resonance ("Portable NMR") due to its low cost and the absence of radiation-related safety concerns. Fifteen (N = 15) healthy female volunteers were selected for the study and underwent an imaging routine consisting of 2D X-ray mammography, quantitative breast 3T MRI (Dixon and T1-based 3D compositional breast imaging), and 1D compositional depth profiling of the right breast using Portable NMR. For each participant, all the measurements were made within 3-4 h of each other. MRI-determined tissue water content was used as the MD-equivalent quantity. Portable NMR depth profiles of tissue water were compared with the equivalent depth profiles reconstructed from Dixon and T1-based MR images, which were used as the MD-equivalent reference standard. The agreement between the depth profiles acquired using Portable NMR and the reconstructed reference-standard profiles was variable but overall encouraging. The agreement was somewhat inferior to that seen in breast tissue explant measurements conducted in vitro, where quantitative micro-CT was used as the reference standard. The lower agreement in vivo can be attributed to an uncertainty in the positioning of the Portable NMR sensor on the breast surface and breast compression in Portable NMR measurements. The degree of agreement between Portable NMR and quantitative MRI is encouraging. While the results call for further development of quantitative Portable NMR, they demonstrate the in-principle feasibility of Portable NMR-based quantitative compositional imaging in vivo and show promise for the development of safe and low-cost protocols for quantification of MD suitable for clinical applications.

Circulating Tumour Cells Indicate the Presence of Residual Disease Post-Castration in Prostate Cancer Patient-Derived Xenograft Models.

Castrate-resistant prostate cancer (CRPC) is the lethal form of prostate cancer. Epithelial mesenchymal plasticity (EMP) has been associated with disease progression to CRPC, and prostate cancer therapies targeting the androgen signalling axis, including androgen deprivation therapy (ADT), promote EMP. We explored effects of castration on EMP in the tumours and circulating tumour cells (CTCs) of patient-derived xenograft (PDX)-bearing castrated mice using human-specific RT-qPCR assays and immunocytochemistry. Expression of prostate epithelial cell marker KLK3 was below detection in most tumours from castrated mice (62%, 23/37 mice), consistent with its known up-regulation by androgens. Endpoint tumour size after castration varied significantly in a PDX model-specific pattern; while most tumours were castration-sensitive (BM18, LuCaP70), the majority of LuCaP105 tumours continued to grow following castration. By contrast, LuCaP96 PDX showed a mixed response to castration. CTCs were detected in 33% of LuCaP105, 43% of BM18, 47% of LuCaP70, and 54% of LuCaP96 castrated mice using RPL32 mRNA measurement in plasma. When present, CTC numbers estimated using human RPL32 expression ranged from 1 to 458 CTCs per ml blood, similar to our previous observations in non-castrated mice. In contrast to their non-castrated counterparts, there was no relationship between tumour size and CTC burden in castrated mice. Unsupervised hierarchical clustering of the gene expression profiles of CTCs collected from castrated and non-castrated mice revealed distinct CTC sub-groups within the pooled population that were classified as having mesenchymal, epithelial, or EMP hybrid gene expression profiles. The epithelial signature was only found in CTCs from non-castrated mice. Hybrid and mesenchymal signatures were detected in CTCs from both castrated and non-castrated mice, with an emphasis towards mesenchymal phenotypes in castrated mice. Post-castration serum PSA levels were either below detection or very low for all the CTC positive samples highlighting the potential usefulness of CTCs for disease monitoring after androgen ablation therapy. In summary, our study of castration effects on prostate cancer PDX CTCs showed that CTCs were often detected in the castrate setting, even in mice with no palpable tumours, and demonstrated the superior ability of CTCs to reveal residual disease over the conventional clinical biomarker serum PSA.

RASSF1A Suppression as a Potential Regulator of Mechano-Pathobiology Associated with Mammographic Density in BRCA Mutation Carriers.

High mammographic density (MD) increases breast cancer (BC) risk and creates a stiff tissue environment. BC risk is also increased in BRCA1/2 gene mutation carriers, which may be in part due to genetic disruption of the tumour suppressor gene Ras association domain family member 1 (RASSF1A), a gene that is also directly regulated by tissue stiffness. High MD combined with BRCA1/2 mutations further increase breast cancer risk, yet BRCA1/2 mutations alone or in combination do not increase MD. The molecular basis for this additive effect therefore remains unclear. We studied the interplay between MD, stiffness, and BRCA1/2 mutation status in human mammary tissue obtained after prophylactic mastectomy from women at risk of developing BC. Our results demonstrate that RASSF1A expression increased in MCF10DCIS.com cell cultures with matrix stiffness up until ranges corresponding with BiRADs 4 stiffnesses (~16 kPa), but decreased in higher stiffnesses approaching malignancy levels (>50 kPa). Similarly, higher RASSF1A protein was seen in these cells when co-cultivated with high MD tissue in murine biochambers. Conversely, local stiffness, as measured by collagen I versus III abundance, repressed RASSF1A protein expression in BRCA1, but not BRCA2 gene mutated tissues; regional density as measured radiographically repressed RASSF1A in both BRCA1/2 mutated tissues. The combinatory effect of high MD and BRCA mutations on breast cancer risk may be due to RASSF1A gene repression in regions of increased tissue stiffness.

Diversity of Epithelial-Mesenchymal Phenotypes in Circulating Tumour Cells from Prostate Cancer Patient-Derived Xenograft Models.

Metastasis is the leading cause of cancer-related deaths worldwide. The epithelial-mesenchymal plasticity (EMP) status of primary tumours has relevance to metastatic potential and therapy resistance. Circulating tumour cells (CTCs) provide a window into the metastatic process, and molecular characterisation of CTCs in comparison to their primary tumours could lead to a better understanding of the mechanisms involved in the metastatic cascade. In this study, paired blood and tumour samples were collected from four prostate cancer patient-derived xenograft (PDX) models (BM18, LuCaP70, LuCaP96, LuCaP105) and assessed using an EMP-focused, 42 gene human-specific, nested quantitative RT-PCR assay. CTC burden varied amongst the various xenograft models with LuCaP96 having the highest number of CTCs per mouse (mean: 704; median: 31) followed by BM18 (mean: 101; median: 21), LuCaP70 (mean: 73; median: 16) and LuCaP105 (mean: 57; median: 6). A significant relationship was observed between tumour size and CTC number (p = 0.0058). Decreased levels of kallikrein-related peptidase 3 (KLK3) mRNA (which encodes prostate-specific antigen; PSA) were observed in CTC samples from all four models compared to their primary tumours. Both epithelial- and mesenchymal-associated genes were commonly expressed at higher levels in CTCs compared to the bulk primary tumour, although some common EMT-associated genes (CDH1, VIM, EGFR, EPCAM) remained unchanged. Immunofluorescence co-staining for pan-cytokeratin (KRT) and vimentin (VIM) indicated variable proportions of CTCs across the full EMP axis, even in the same model. EMP hybrids predominated in the BM18 and LuCaP96 models, but were not detected in the LuCaP105 model, and variable numbers of KRT+ and human VIM+ cells were observed in each model. SERPINE1, which encodes plasminogen activator inhibitor-1 (PAI-1), was enriched at the RNA level in CTCs compared to primary tumours and was the most commonly expressed mesenchymal gene in the CTCs. Co-staining for SERPINE1 and KRT revealed SERPINE1+ cells in 7/11 samples, six of which had SERPINE+KRT+ CTCs. Cell size variation was observed in CTCs. The majority of samples (8/11) contained larger CTCs ranging from 15.3 to 37.8 µm, whilst smaller cells (10.7 ± 4.1 µm, similar in size to peripheral blood mononuclear cells (PBMCs)) were identified in 6 of 11 samples. CTC clusters were also identified in 9/11 samples, containing 2-100 CTCs per cluster. Where CTC heterogeneity was observed in the clusters, epithelial-like cells (KRT+VIM-) were located on the periphery of the cluster, forming a layer around hybrid (KRT+VIM+) or mesenchymal-like (KRT-VIM+) cells. The CTC heterogeneity observed in these models emphasises the complexity in CTC isolation and classification and supports the increasingly recognised importance of the epithelial-mesenchymal hybrid state in cancer progression and metastasis.

Overexpression of miRNA-9 enhances galectin-3 levels in oral cavity cancers.

Oral cavity cancer (OCC) is the predominant subtype of head and neck cancer (HNC) and has up to 50% mortality. Genome-wide microRNA (miR) sequencing data indicates overexpression of miR-9-5p in HNC tumours, however, the biological role of miR-9-5p in OCC is complex; it can either act as a tumour suppressor or an oncomir, regulating many target genes at the post-transcriptional level. We have investigated the overexpression of miR-9-5p in three OCC cell lines. We have evaluated its expression levels and Galectin-3 as potential biomarkers in saliva samples collected from controls and OCC patients. We found that over expression of miR-9-5p in OCC cell lines resulted in a significant reduction in cell proliferation and migration, and an increase in apoptosis, which was paralleled by an increase in Galectin-3 secretion and export of Galectin-3 protein. Our data are consistent with miR-9-5p being a modulator of Galectin-3 via the AKT/γ-catenin pathway. In addition, the positive correlation between the levels of miR-9-5p expression and secreted Galectin-3 in saliva reflects a similar relationship in vivo, and supports the utility of their integrative evaluation in OCC. Our findings indicate that both miR-9-5p and Galectin-3 are critical biomolecules in the progression of OCC.

Corrigendum: Heparanase Promotes Syndecan-1 Expression to Mediate Fibrillar Collagen and Mammographic Density in Human Breast Tissue Cultured ex vivo.

[This corrects the article DOI: 10.3389/fcell.2020.00599.].

Integrin alpha-2 and beta-1 expression increases through multiple generations of the EDW01 patient-derived xenograft model of breast cancer-insight into their role in epithelial mesenchymal transition in vivo gained from an in vitro model system.

BACKGROUND: Breast cancers acquire aggressive capabilities via epithelial to mesenchymal transition (EMT), in which various integrins/integrin-linked kinase signalling are upregulated. METHODS: We investigated this in two patient-derived xenografts (PDXs) developed from breast-to-bone metastases, and its functional significance in a breast cancer cell line system. ED03 and EDW01 PDXs were grown subcutaneously in immunocompromised SCID mice through 11 passages and 7 passages, respectively. Tumour tissue was assessed using immunohistochemistry (IHC) for oestrogen receptor (ER)-alpha, E-cadherin, vimentin, Twist1, beta-catenin, P120-RasGAP, CD44, CD24 and Ki67, and RT-qPCR of EMT-related factors (CDH1, VIM, CD44, CD24), integrins beta 1 (ITGB1), alpha 2 (ITGA2) and ILK. Integrin and ILK expression in epidermal growth factor (EGF)-induced EMT of the PMC42-ET breast cancer cell line was assessed by RT-qPCR and Western blotting, as were the effects of their transient knockdown via small interfering RNA +/- EGF. Cell migration, changes in cell morphology and adhesion of siRNA-transfected PMC42-ET cells to various extracellular matrix (ECM) substrates was assessed. RESULTS: The ED03 (ER+/PR-/HER2-/lobular) and EDW01 (ER+/PR-/HER2-/ductal) PDXs were both classified as molecular subtype luminal A. ED03 xenografts exhibited mutated E-cadherin with minimal expression, but remained vimentin-negative across all passages. In EDW01, the hypoxic indicator gene CAIX and Twist1 were co-ordinately upregulated at passages 4-5, corresponding with a decrease in E-cadherin. At passages 6-7, VIM was upregulated along with ITGB1 and ITGA2, consistent with an increasing EMT. The ED03 PDX displayed minimal change over passages in mice, for all genes examined. ILK, ITGB1 and ITGA2 mRNAs were also increased in the EGF-induced EMT of PMC42-ET cells (in which CDH1 was downregulated) although siRNA against these targets revealed that this induction was not necessary for the observed EMT. However, their knockdown significantly reduced EMT-associated adhesion and Transwell migration. CONCLUSION: Our data suggest that despite an increase in ITGA2 and ITGB1 gene expression in the EMT exhibited by EDW01 PDX over multiple generations, this pathway may not necessarily drive the EMT process.

Heparanase Promotes Syndecan-1 Expression to Mediate Fibrillar Collagen and Mammographic Density in Human Breast Tissue Cultured ex vivo.

Mammographic density (MD) is a strong and independent factor for breast cancer (BC) risk and is increasingly associated with BC progression. We have previously shown in mice that high MD, which is characterized by the preponderance of a fibrous stroma, facilitates BC xenograft growth and metastasis. This stroma is rich in extracellular matrix (ECM) factors, including heparan sulfate proteoglycans (HSPGs), such as the BC-associated syndecan-1 (SDC1). These proteoglycans tether growth factors, which are released by heparanase (HPSE). MD is positively associated with estrogen exposure and, in cell models, estrogen has been implicated in the upregulation of HPSE, the activity of which promotes SDC expression. Herein we describe a novel measurement approach (single-sided NMR) using a patient-derived explant (PDE) model of normal human (female) mammary tissue cultured ex vivo to investigate the role(s) of HPSE and SDC1 on MD. Relative HSPG gene and protein analyses determined in patient-paired high vs. low MD tissues identified SDC1 and SDC4 as potential mediators of MD. Using the PDE model we demonstrate that HPSE promotes SDC1 rather than SDC4 expression and cleavage, leading to increased MD. In this model system, synstatin (SSTN), an SDC1 inhibitory peptide designed to decouple SDC1-ITGαvß3 parallel collagen alignment, reduced the abundance of fibrillar collagen as assessed by picrosirius red viewed under polarized light, and reduced MD. Our results reveal a potential role for HPSE in maintaining MD via its direct regulation of SDC1, which in turn physically tethers collagen into aligned fibers characteristic of MD. We propose that inhibitors of HPSE and/or SDC1 may afford an opportunity to reduce MD in high BC risk individuals and reduce MD-associated BC progression in conjunction with established BC therapies.

Identifying Therapies to Combat Epithelial Mesenchymal Plasticity-Associated Chemoresistance to Conventional Breast Cancer Therapies Using An shRNA Library Screen.

BACKGROUND: Breast cancer (BC) is a heterogeneous disease for which the commonly used chemotherapeutic agents primarily include the anthracyclines (doxorubicin, epirubicin), microtubule inhibitors (paclitaxel, docetaxel, eribulin), and alkylating agents (cyclophosphamide). While these drugs can be highly effective, metastatic tumours are frequently refractory to treatment or become resistant upon tumour relapse. METHODS: We undertook a cell polarity/epithelial mesenchymal plasticity (EMP)-enriched short hairpin RNA (shRNA) screen in MDA-MB-468 breast cancer cells to identify factors underpinning heterogeneous responses to three chemotherapeutic agents used clinically in breast cancer: Doxorubicin, docetaxel, and eribulin. shRNA-transduced cells were treated for 6 weeks with the EC10 of each drug, and shRNA representation assessed by deep sequencing. We first identified candidate genes with depleted shRNA, implying that their silencing could promote a response. Using the Broad Institute's Connectivity Map (CMap), we identified partner inhibitors targeting the identified gene families that may induce cell death in combination with doxorubicin, and tested them with all three drug treatments. RESULTS: In total, 259 shRNAs were depleted with doxorubicin treatment (at p < 0.01), 66 with docetaxel, and 25 with eribulin. Twenty-four depleted hairpins overlapped between doxorubicin and docetaxel, and shRNAs for TGFB2, RUNX1, CCDC80, and HYOU1 were depleted across all the three drug treatments. Inhibitors of MDM/TP53, TGFBR, and FGFR were identified by CMap as the top pharmaceutical perturbagens and we validated the combinatorial benefits of the TGFBR inhibitor (SB525334) and MDM inhibitor (RITA) with doxorubicin treatment, and also observed synergy between the inhibitor SB525334 and eribulin in MDA-MB-468 cells. CONCLUSIONS: Taken together, a cell polarity/EMP-enriched shRNA library screen identified relevant gene products that could be targeted alongside current chemotherapeutic agents for the treatment of invasive BC.

Applications of RNA from circulating tumor cells.

Circulating tumour cells (CTCs) are shed into the bloodstream from both primary and secondary tumours and provide a non-invasive means to study tumor progression and response to treatment. Assessment of ribonucleic acid (RNA) and monitoring dynamic changes in gene expression profiles of CTCs extends their clinical and prognostic power and establish their role in guiding treatment. Among these methods, droplet digital (RT-ddPCR) technique provides a high sensitivity and detectibility of CTCs. RNA-sequencing (RNAseq) is the most comprehensive method, that would allow the simultaneous measurement of a large number of genes and theoretically the whole transcriptome. Since CTCs are heterogeneous in nature, single cell RNAseq methods are very valuable in assessing population dynamics and functional states of CTCs. While RNA in situ hybridization (RNA-ISH) is used relatively less frequently, it also allows for the assessment of expression of multiple genes within individual CTCs. Epithelial to Mesenchymal Transition (EMT) is a major contributor to metastasis, providing a mechanism to allow cells to become migratory and invasive, and to survive in the bloodstream. Monitoring CTCs undergoing EMT may lead to improvement in their prognostic and predictive power. Here, we review various RNA analysis of CTCs and those that undergo EMT and their application in diagnosis, prognosis and treatment of cancers.

Multi-Omics Characterization of the Spontaneous Mesenchymal-Epithelial Transition in the PMC42 Breast Cancer Cell Lines.

Epithelial-mesenchymal plasticity (EMP), encompassing epithelial-mesenchymal transition (EMT) and mesenchymal-epithelial transition (MET), are considered critical events for cancer metastasis. We investigated chromosomal heterogeneity and chromosomal instability (CIN) profiles of two sister PMC42 breast cancer (BC) cell lines to assess the relationship between their karyotypes and EMP phenotypic plasticity. Karyotyping by GTG banding and exome sequencing were aligned with SWATH quantitative proteomics and existing RNA-sequencing data from the two PMC42 cell lines; the mesenchymal, parental PMC42-ET cell line and the spontaneously epithelially shifted PMC42-LA daughter cell line. These morphologically distinct PMC42 cell lines were also compared with five other BC cell lines (MDA-MB-231, SUM-159, T47D, MCF-7 and MDA-MB-468) for their expression of EMP and cell surface markers, and stemness and metabolic profiles. The findings suggest that the epithelially shifted cell line has a significantly altered ploidy of chromosomes 3 and 13, which is reflected in their transcriptomic and proteomic expression profiles. Loss of the TGFßR2 gene from chromosome 3 in the epithelial daughter cell line inhibits its EMT induction by TGF-ß stimulus. Thus, integrative 'omics' characterization established that the PMC42 system is a relevant MET model and provides insights into the regulation of phenotypic plasticity in breast cancer.

Interrogation of Phenotypic Plasticity between Epithelial and Mesenchymal States in Breast Cancer.

Dynamic interconversions between transitional epithelial and mesenchymal states underpin the epithelial mesenchymal plasticity (EMP) seen in some carcinoma cell systems. We have delineated epithelial and mesenchymal subpopulations existing within the PMC42-LA breast cancer cell line by their EpCAM expression. These purified but phenotypically plastic states, EpCAMHigh (epithelial) and EpCAMLow (mesenchymal), have the ability to regain the phenotypic equilibrium of the parental population (i.e., 80% epithelial and 20% mesenchymal) over time, although the rate of reversion in the mesenchymal direction (epithelial-mesenchymal transition; EMT) is higher than that in the epithelial direction (mesenchymal-epithelial transition; MET). Single-cell clonal propagation was implemented to delineate the molecular and cellular features of this intrinsic heterogeneity with respect to EMP flux. The dynamics of the phenotypic proportions of epithelial and mesenchymal states in single-cell generated clones revealed clonal diversity and intrinsic plasticity. Single cell-derived clonal progenies displayed differences in their functional attributes of proliferation, stemness marker (CD44/CD24), migration, invasion and chemo-sensitivity. Interrogation of genomic copy number variations (CNV) with whole exome sequencing (WES) in the context of chromosome count from metaphase spread indicated that chromosomal instability was not influential in driving intrinsic phenotypic plasticity. Overall, these findings reveal the stochastic nature of both the epithelial and mesenchymal subpopulations, and the single cell-derived clones for differential functional attributes.

Quantification of breast tissue density: Correlation between single-sided portable NMR and micro-CT measurements.

Mammographic density (MD) is a strong independent risk factor for breast cancer. Traditional screening for MD using X-ray mammography involves ionising radiation, which is not suitable for young women, those with previous radiation exposure, or those having undergone a partial mastectomy. Therefore, alternative approaches for MD screening that do not involve ionising radiation will be important as the clinical use of MD increases, and as more frequent MD testing becomes desirable for research purposes. We have previously demonstrated the potential utility of spin relaxation-based, single-sided portable-NMR measurements for the purpose of MD quantification. We present here a more refined analysis by quantifying breast tissue density in excised samples on a continuous scale (0% to 100% fibroglandular tissue content) using micro-CT (µCT), and comparing the results to spin-relaxation and diffusion portable-NMR measurements of the same samples. µCT analysis of mammary tissues containing high- and low-MD (HMD and LMD, respectively) regions had Hounsfield Unit (HU) histograms with a bimodal pattern, with HMD regions exhibiting significantly higher HU values than LMD regions. Quantitative MD (%HMD) values obtained using µCT exhibited an excellent correlation with portable-NMR results, namely longitudinal spin-relaxation time constants (T1) and the relative tissue water content obtained from portable-NMR diffusion measurements (R2 = 0.92, p < 0.0001 and R2 = 0.96, p < 0.0001, respectively). These findings are consistent with our previous results demonstrating relatively high water content in HMD breast tissue, consistent with the high proportion of fibroglandular tissue, FGT, which in turn contains more abundant water-carrying HSPG proteins. We observed an excellent correlation between the T1 values and diffusion NMR-measured relative tissue water content (R2 = 0.94, p < 0.0001). These findings demonstrate, for the first time, the ability of single-sided portable NMR to accurately quantify MD in vitro on a continuous scale. The results also indicate that portable-NMR analysis can assist in the identification of features underpinning MD, namely FGT and adipose tissue content. Future work will involve application of portable NMR to quantifying MD in vivo.

Human-specific RNA analysis shows uncoupled epithelial-mesenchymal plasticity in circulating and disseminated tumour cells from human breast cancer xenografts.

Blood samples, bone marrow, tumours and metastases where possible were collected from SCID mice bearing orthotopic xenografts of the triple-negative MDA-MB-468 cell line or a transplantable ER-positive patient derived xenograft (ED-03), and assessed using human-specific, tandem-nested RT-qPCR for markers relating to detection of circulating (CTCs) and disseminated tumour cells (DTCs), breast cancer clinicopathology, the 'cancer stem cell' phenotype, metabolism, hypoxia and epithelial-mesenchymal plasticity (EMP). Increased levels of SNAI1, ILK, NOTCH1, CK20, and PGR, and a decrease/loss of EPCAM in CTCs/DTCs were observed relative to the primary xenograft across both models. Decreased CD24 and EGFR was restricted to the MDA-MB-468 model, while increased TFF1 was seen in the ED-03 model. The major metabolic regulator PPARGC1A, and several hypoxia-related markers (HIF1A, APLN and BNIP3) were significantly elevated in both models. Increased expression of mesenchymal markers including SNAI1 was seen across both models, however CDH1 did not decrease concordantly, and several other epithelial markers were increased, suggesting an uncoupling of EMP to produce an EMP hybrid or partial-EMT. Single cell analysis of ED-03 CTCs, although limited, indicated uncoupling of the EMP axis in single hybrid cells, rather than distinct pools of epithelial or mesenchymal-enriched cells, however dynamic heterogeneity between CTCs/DTCs cannot be ruled out. Reduced CD24 expression was observed in the MDA-MB-468 CTCs, consistent with the 'breast cancer stem cell' phenotype, and metastatic deposits in this model mostly resembled the primary xenografts, consistent with the mesenchymal-epithelial transition paradigm.