Deciphering the Correlation between Breast Tumor Samples and Cell Lines by Integrating Copy Number Changes and Gene Expression Profiles

Identifying cell lines across pan-cancer to be used in preclinical research as a proxy for patient tumor samples

In pre-clinical trials of anti-cancer drugs, cell lines are utilized as a model for patient tumor samples to understand the response of drugs. However, in vitro culture of cell lines, in general, alters the biology of the cell lines and likely gives rise to systematic differences from the tumor samples' genomic profiles; hence the drug response of cell lines may deviate from actual patients' drug response. In this study, we computed a similarity score for the selection of cell lines depicting the close and far resemblance to patient tumor samples in twenty-two different cancer types at genetic, genomic, and epigenetic levels integrating multi-omics datasets. We also considered the presence of immune cells in tumor samples and cancer-related biological pathways in this score which aids personalized medicine research in cancer. We showed that based on these similarity scores, cell lines were able to recapitulate the drug response of patient tumor samples for several FDA-approved cancer drugs in multiple cancer types. Based on these scores, several of the high-rank cell lines were shown to have a close likeness to the corresponding tumor type in previously reported in vitro experiments.

Computational estimation of quality and clinical relevance of cancer cell lines

Immortal cancer cell lines (CCLs) are the most widely used system for investigating cancer biology and for the preclinical development of oncology therapies. Pharmacogenomic and genome-wide editing screenings have facilitated the discovery of clinically relevant gene-drug interactions and novel therapeutic targets via large panels of extensively characterised CCLs. However, tailoring pharmacological strategies in a precision medicine context requires bridging the existing gaps between tumours and in vitro models. Indeed, intrinsic limitations of CCLs such as misidentification, the absence of tumour microenvironment and genetic drift have highlighted the need to identify the most faithful CCLs for each primary tumour while addressing their heterogeneity, with the development of new models where necessary. Here, we discuss the most significant limitations of CCLs in representing patient features, and we review computational methods aiming at systematically evaluating the suitability of CCLs as tumour proxies and identifying the best patient representative in vitro models. Additionally, we provide an overview of the applications of these methods to more complex models and discuss future machine-learning-based directions that could resolve some of the arising discrepancies.

CELLector: Genomics-Guided Selection of Cancer In Vitro Models

Selecting appropriate cancer models is a key prerequisite for maximizing translational potential and clinical relevance of in vitro oncology studies. We developed CELLector: an R package and R Shiny application allowing researchers to select the most relevant cancer cell lines in a patient-genomic-guided fashion. CELLector leverages tumor genomics to identify recurrent subtypes with associated genomic signatures. It then evaluates these signatures in cancer cell lines to prioritize their selection. This enables users to choose appropriate in vitro models for inclusion or exclusion in retrospective analyses and future studies. Moreover, this allows bridging outcomes from cancer cell line screens to precisely defined sub-cohorts of primary tumors. Here, we demonstrate the usefulness and applicability of CELLector, showing how it can aid prioritization of in vitro models for future development and unveil patient-derived multivariate prognostic and therapeutic markers. CELLector is freely available at https://ot-cellector.shinyapps.io/CELLector_App/ (code at https://github.com/francescojm/CELLector and https://github.com/francescojm/CELLector_App).

Three-Dimensional Reconstructed Bone Marrow Matrix Culture Improves the Viability of Primary Myeloma Cells In-Vitro via a STAT3-Dependent Mechanism

Primary myeloma (PM) cells are short-lived in conventional culture, which limited their usefulness as a study model. Here, we evaluated if three-dimensional (3D) culture can significantly prolong the longevity of PM cells in-vitro. We employed a previously established 3D model for culture of bone marrow mononuclear cells isolated from 15 patients. We assessed the proportion of PM cells, viability and proliferation using CD38 staining, trypan blue exclusion assays and carboxy fluorescein succinimidyl ester (CFSE) staining, respectively. We observed significantly more CD38+ viable cells in 3D than in conventional culture (65% vs. 25%, p = 0.006) on day 3. CFSE staining showed no significant difference in cell proliferation between the two culture systems. Moreover, we found that PM cells in 3D culture are more STAT3 active by measure of pSTAT3 staining (66% vs. 10%, p = 0.008). Treatment of IL6, a STAT3 activator significantly increased CD38+ cell viability (41% to 68%, p = 0.021). In comparison, inhibition of STAT3 with Stattic significantly decreased PM cell viability in 3D culture (38% to 17% p = 0.010). Neither IL6 nor Stattic affected the PM cell viability in conventional culture. This study suggests that 3D culture can significantly improve the longevity of PM cells in-vitro, and STAT3 activation can further improve their viability.

A Clinical Genomics-Guided Prioritizing Strategy Enables Selecting Proper Cancer Cell Lines for Biomedical Research

Selecting appropriate cell lines to represent a disease is crucial for the success of biomedical research, because the usage of less relevant cell lines could deliver misleading results. However, systematic guidance on cell line selection is unavailable. Here we developed a clinical Genomics-guided Prioritizing Strategy for Cancer Cell Lines (CCL-cGPS) and help to guide this process. Statistical analyses revealed CCL-cGPS selected cell lines were among the most appropriate models. Moreover, we observed a linear correlation between the drug response and CCL-cGPS score of cell lines for breast and thyroid cancers. Using RT4 cells selected by CCL-GPS, we identified mebendazole and digitoxin as candidate drugs against bladder cancer and validate their promising anticancer effect through in vitro and in vivo experiments. Additionally, a web tool was developed. In conclusion, CCL-cGPS bridges the gap between tumors and cell lines, presenting a helpful guide to select the most suitable cell line models.

•

Cell lines were ranked by the resemblance of transcriptional signatures to tumors

•

Among 44 tumor subtypes, CCL-cGPS provides proper cell lines for each subtype

•

CCL-cGPS was verified by the computational analysis, in vitro and in vivo assays

•

A web tool was developed to guide the selection of the most suitable cell lines

Personalized Network Modeling of the Pan-Cancer Patient and Cell Line Interactome

PURPOSE

Personalized network inference on diverse clinical and in vitro model systems across cancer types can be used to delineate specific regulatory mechanisms, uncover drug targets and pathways, and develop individualized predictive models in cancer.

METHODS

We developed TransPRECISE (personalized cancer-specific integrated network estimation model), a multiscale Bayesian network modeling framework, to analyze the pan-cancer patient and cell line interactome to identify differential and conserved intrapathway activities, to globally assess cell lines as representative models for patients, and to develop drug sensitivity prediction models. We assessed pan-cancer pathway activities for a large cohort of patient samples (> 7,700) from the Cancer Proteome Atlas across ≥ 30 tumor types, a set of 640 cancer cell lines from the MD Anderson Cell Lines Project spanning 16 lineages, and ≥ 250 cell lines' response to > 400 drugs.

RESULTS

TransPRECISE captured differential and conserved proteomic network topologies and pathway circuitry between multiple patient and cell line lineages: ovarian and kidney cancers shared high levels of connectivity in the hormone receptor and receptor tyrosine kinase pathways, respectively, between the two model systems. Our tumor stratification approach found distinct clinical subtypes of the patients represented by different sets of cell lines: patients with head and neck tumors were classified into two different subtypes that are represented by head and neck and esophagus cell lines and had different prognostic patterns (456 v 654 days of median overall survival; P = .02). High predictive accuracy was observed for drug sensitivities in cell lines across multiple drugs (median area under the receiver operating characteristic curve > 0.8) using Bayesian additive regression tree models with TransPRECISE pathway scores.

CONCLUSION

Our study provides a generalizable analytic framework to assess the translational potential of preclinical model systems and to guide pathway-based personalized medical decision making, integrating genomic and molecular data across model systems.

CellExpress: a comprehensive microarray-based cancer cell line and clinical sample gene expression analysis online system

Database URL

http://cellexpress.cgm.ntu.edu.tw/.

Multimodal Assessment of Estrogen Receptor mRNA Profiles to Quantify Estrogen Pathway Activity in Breast Tumors

BACKGROUND

Molecular markers have transformed our understanding of the heterogeneity of breast cancer and have allowed the identification of genomic profiles of estrogen receptor (ER)-α signaling. However, our understanding of the transcriptional profiles of ER signaling remains inadequate. Therefore, we sought to identify the genomic indicators of ER pathway activity that could supplement traditional immunohistochemical (IHC) assessments of ER status to better understand ER signaling in the breast tumors of individual patients.

MATERIALS AND METHODS

We reduced ESR1 (gene encoding the ER-α protein) mRNA levels using small interfering RNA in ER⁺ MCF7 breast cancer cells and assayed for transcriptional changes using Affymetrix HG U133 Plus 2.0 arrays. We also compared 1034 ER⁺ and ER^- breast tumors from publicly available microarray data. The principal components of ER activity generated from these analyses and from other published estrogen signatures were compared with ESR1 expression, ER-α IHC, and patient survival.

RESULTS

Genes differentially expressed in both analyses were associated with ER-α IHC and ESR1 mRNA expression. They were also significantly enriched for estrogen-driven molecular pathways associated with ESR1, cyclin D1 (CCND1), MYC (v-myc avian myelocytomatosis viral oncogene homolog), and NFKB (nuclear factor kappa B). Despite their differing constituent genes, the principal components generated from these new analyses and from previously published ER-associated gene lists were all associated with each other and with the survival of patients with breast cancer treated with endocrine therapies.

CONCLUSION

A biomarker of ER-α pathway activity, generated using ESR1-responsive mRNAs in MCF7 cells, when used alongside ER-α IHC and ESR1 mRNA expression, could provide a method for further stratification of patients and add insight into ER pathway activity in these patients.