Network-guided identification of cancer-selective combinatorial therapies in ovarian cancer.

Each patient's cancer consists of multiple cell subpopulations that are inherently heterogeneous and may develop differing phenotypes such as drug sensitivity or resistance. A personalized treatment regimen should therefore target multiple oncoproteins in the cancer cell populations that are driving the treatment resistance or disease progression in a given patient to provide maximal therapeutic effect, while avoiding severe co-inhibition of non-malignant cells that would lead to toxic side effects. To address the intra- and inter-tumoral heterogeneity when designing combinatorial treatment regimens for cancer patients, we have implemented a machine learning-based platform to guide identification of safe and effective combinatorial treatments that selectively inhibit cancer-related dysfunctions or resistance mechanisms in individual patients. In this case study, we show how the platform enables prediction of cancer-selective drug combinations for patients with high-grade serous ovarian cancer using single-cell imaging cytometry drug response assay, combined with genome-wide transcriptomic and genetic profiles. The platform makes use of drug-target interaction networks to prioritize those combinations that warrant further preclinical testing in scarce patient-derived primary cells. During the case study in ovarian cancer patients, we investigated (i) the relative performance of various ensemble learning algorithms for drug response prediction, (ii) the use of matched single-cell RNA-sequencing data to deconvolute cell population-specific transcriptome profiles from bulk RNA-seq data, (iii) and whether multi-patient or patient-specific predictive models lead to better predictive accuracy. The general platform and the comparison results are expected to become useful for future studies that use similar predictive approaches also in other cancer types.

Loss of DIAPH1 causes SCBMS, combined immunodeficiency, and mitochondrial dysfunction.

BACKGROUND: Homozygous loss of DIAPH1 results in seizures, cortical blindness, and microcephaly syndrome (SCBMS). We studied 5 Finnish and 2 Omani patients with loss of DIAPH1 presenting with SCBMS, mitochondrial dysfunction, and immunodeficiency. OBJECTIVE: We sought to further characterize phenotypes and disease mechanisms associated with loss of DIAPH1. METHODS: Exome sequencing, genotyping and haplotype analysis, B- and T-cell phenotyping, in vitro lymphocyte stimulation assays, analyses of mitochondrial function, immunofluorescence staining for cytoskeletal proteins and mitochondria, and CRISPR-Cas9 DIAPH1 knockout in heathy donor PBMCs were used. RESULTS: Genetic analyses found all Finnish patients homozygous for a rare DIAPH1 splice-variant (NM_005219:c.684+1G>A) enriched in the Finnish population, and Omani patients homozygous for a previously described pathogenic DIAPH1 frameshift-variant (NM_005219:c.2769delT;p.F923fs). In addition to microcephaly, epilepsy, and cortical blindness characteristic to SCBMS, the patients presented with infection susceptibility due to defective lymphocyte maturation and 3 patients developed B-cell lymphoma. Patients' immunophenotype was characterized by poor lymphocyte activation and proliferation, defective B-cell maturation, and lack of naive T cells. CRISPR-Cas9 knockout of DIAPH1 in PBMCs from healthy donors replicated the T-cell activation defect. Patient-derived peripheral blood T cells exhibited impaired adhesion and inefficient microtubule-organizing center repositioning to the immunologic synapse. The clinical symptoms and laboratory tests also suggested mitochondrial dysfunction. Experiments with immortalized, patient-derived fibroblasts indicated that DIAPH1 affects the amount of complex IV of the mitochondrial respiratory chain. CONCLUSIONS: Our data demonstrate that individuals with SCBMS can have combined immune deficiency and implicate defective cytoskeletal organization and mitochondrial dysfunction in SCBMS pathogenesis.

Improved detection of differentially represented DNA barcodes for high-throughput clonal phenomics.

Cellular DNA barcoding has become a popular approach to study heterogeneity of cell populations and to identify clones with differential response to cellular stimuli. However, there is a lack of reliable methods for statistical inference of differentially responding clones. Here, we used mixtures of DNA-barcoded cell pools to generate a realistic benchmark read count dataset for modelling a range of outcomes of clone-tracing experiments. By accounting for the statistical properties intrinsic to the DNA barcode read count data, we implemented an improved algorithm that results in a significantly lower false-positive rate, compared to current RNA-seq data analysis algorithms, especially when detecting differentially responding clones in experiments with strong selection pressure. Building on the reliable statistical methodology, we illustrate how multidimensional phenotypic profiling enables one to deconvolute phenotypically distinct clonal subpopulations within a cancer cell line. The mixture control dataset and our analysis results provide a foundation for benchmarking and improving algorithms for clone-tracing experiments.

High prevalence of GPRC5A germline mutations in BRCA1-mutant breast cancer patients.

In a search for new breast cancer (BC) predisposing genes, we performed a whole exome sequencing analysis using six patient samples of familial BC and identified a germline inactivating mutation c.183delG [p. Arg61fs] in an orphan G protein-coupled receptor GPRC5A. An extended case-control study revealed a tenfold enrichment for this mutation in BC patients carrying the 5382insC allele of BRCA1, the major founder mutation in the Russian population, compared to wild-type BRCA1 BC cases [6/117 (5.1%) vs. 8/1578 (0.5%), p = 0.0002]. In mammary tumors (n = 60), the mRNA expression of GPRC5A significantly correlated with that of BRCA1 (p = 0.00018). In addition, the amount of GPRC5A transcript was significantly lower in BC obtained from BRCA1 mutation carriers (n = 17) compared to noncarriers (n = 93) (p = 0.026). Accordingly, a siRNA-mediated knockdown of either BRCA1 or GPRC5A in the MDA-MB-231 human BC cell line reduced expression of GPRC5A or BRCA1, respectively. Knockdown of GPRC5A also attenuated radiation-induced BRCA1- and RAD51-containing nuclear DNA repair foci. Taken together, these data suggest that GPRC5A is a modifier of BC risk in BRCA1 mutation carriers and reveals a functional interaction of these genes.

A synthetic lethal dependency on casein kinase 2 in response to replication-perturbing therapeutics in RB1-deficient cancer cells.

Resistance to therapy commonly develops in patients with high-grade serous ovarian carcinoma (HGSC) and triple-negative breast cancer (TNBC), urging the search for improved therapeutic combinations and their predictive biomarkers. Starting from a CRISPR knockout screen, we identified that loss of RB1 in TNBC or HGSC cells generates a synthetic lethal dependency on casein kinase 2 (CK2) for surviving the treatment with replication-perturbing therapeutics such as carboplatin, gemcitabine, or PARP inhibitors. CK2 inhibition in RB1-deficient cells resulted in the degradation of another RB family cell cycle regulator, p130, which led to S phase accumulation, micronuclei formation, and accelerated PARP inhibition-induced aneuploidy and mitotic cell death. CK2 inhibition was also effective in primary patient-derived cells. It selectively prevented the regrowth of RB1-deficient patient HGSC organoids after treatment with carboplatin or niraparib. As about 25% of HGSCs and 40% of TNBCs have lost RB1 expression, CK2 inhibition is a promising approach to overcome resistance to standard therapeutics in large strata of patients.

Single-cell transcriptomes identify patient-tailored therapies for selective co-inhibition of cancer clones.

Intratumoral cellular heterogeneity necessitates multi-targeting therapies for improved clinical benefits in advanced malignancies. However, systematic identification of patient-specific treatments that selectively co-inhibit cancerous cell populations poses a combinatorial challenge, since the number of possible drug-dose combinations vastly exceeds what could be tested in patient cells. Here, we describe a machine learning approach, scTherapy, which leverages single-cell transcriptomic profiles to prioritize multi-targeting treatment options for individual patients with hematological cancers or solid tumors. Patient-specific treatments reveal a wide spectrum of co-inhibitors of multiple biological pathways predicted for primary cells from heterogenous cohorts of patients with acute myeloid leukemia and high-grade serous ovarian carcinoma, each with unique resistance patterns and synergy mechanisms. Experimental validations confirm that 96% of the multi-targeting treatments exhibit selective efficacy or synergy, and 83% demonstrate low toxicity to normal cells, highlighting their potential for therapeutic efficacy and safety. In a pan-cancer analysis across five cancer types, 25% of the predicted treatments are shared among the patients of the same tumor type, while 19% of the treatments are patient-specific. Our approach provides a widely-applicable strategy to identify personalized treatment regimens that selectively co-inhibit malignant cells and avoid inhibition of non-cancerous cells, thereby increasing their likelihood for clinical success.

Synthetic lethal interaction between WEE1 and PKMYT1 is a target for multiple low-dose treatment of high-grade serous ovarian carcinoma.

Ovarian cancer is driven by genetic alterations that necessitate protective DNA damage and replication stress responses through cell cycle control and genome maintenance. This creates specific vulnerabilities that may be exploited therapeutically. WEE1 kinase is a key cell cycle control kinase, and it has emerged as a promising cancer therapy target. However, adverse effects have limited its clinical progress, especially when tested in combination with chemotherapies. A strong genetic interaction between WEE1 and PKMYT1 led us to hypothesize that a multiple low-dose approach utilizing joint WEE1 and PKMYT1 inhibition would allow exploitation of the synthetic lethality. We found that the combination of WEE1 and PKMYT1 inhibition exhibited synergistic effects in eradicating ovarian cancer cells and organoid models at a low dose. The WEE1 and PKMYT1 inhibition synergistically promoted CDK activation. Furthermore, the combined treatment exacerbated DNA replication stress and replication catastrophe, leading to increase of the genomic instability and inflammatory STAT1 signalling activation. These findings suggest a new multiple low-dose approach to harness the potency of WEE1 inhibition through the synthetic lethal interaction with PKMYT1 that may contribute to the development of new treatments for ovarian cancer.

A platform for efficient establishment and drug-response profiling of high-grade serous ovarian cancer organoids.

The broad research use of organoids from high-grade serous ovarian cancer (HGSC) has been hampered by low culture success rates and limited availability of fresh tumor material. Here, we describe a method for generation and long-term expansion of HGSC organoids with efficacy markedly improved over previous reports (53% vs. 23%-38%). We established organoids from cryopreserved material, demonstrating the feasibility of using viably biobanked tissue for HGSC organoid derivation. Genomic, histologic, and single-cell transcriptomic analyses revealed that organoids recapitulated genetic and phenotypic features of original tumors. Organoid drug responses correlated with clinical treatment outcomes, although in a culture conditions-dependent manner and only in organoids maintained in human plasma-like medium (HPLM). Organoids from consenting patients are available to the research community through a public biobank and organoid genomic data are explorable through an interactive online tool. Taken together, this resource facilitates the application of HGSC organoids in basic and translational ovarian cancer research.

Selective drug combination vulnerabilities in STAT3- and TP53-mutant malignant NK cells.

Mature natural killer (NK) cell neoplasms are rare but very aggressive types of cancers. With currently available treatments, they have a very poor prognosis and, as such, are an example of group of cancers in which the development of effective precision therapies is needed. Using both short- and long-term drug sensitivity testing, we explored novel ways to target NK-cell neoplasms by combining the clinically approved JAK inhibitor ruxolitinib with other targeted agents. We profiled 7 malignant NK-cell lines in drug sensitivity screens and identified that these exhibit differential drug sensitivities based on their genetic background. In short-term assays, various classes of drugs combined with ruxolitinib seemed highly potent. Strikingly, resistance to most of these combinations emerged rapidly when explored in long-term assays. However, 4 combinations were identified that selectively eradicated the cancer cells and did not allow for development of resistance: ruxolitinib combined with the mouse double-minute 2 homolog (MDM2) inhibitor idasanutlin in STAT3-mutant, TP53 wild-type cell lines; ruxolitinib combined with the farnesyltransferase inhibitor tipifarnib in TP53-mutant cell lines; and ruxolitinib combined with either the glucocorticoid dexamethasone or the myeloid cell leukemia-1 (MCL-1) inhibitor S63845 but both without a clear link to underlying genetic features. In conclusion, using a new drug sensitivity screening approach, we identified drug combinations that selectively target mature NK-cell neoplasms and do not allow for development of resistance, some of which can be applied in a genetically stratified manner.

Pan-RAF inhibition induces apoptosis in acute myeloid leukemia cells and synergizes with BCL2 inhibition.

Pan-RAF inhibitors have shown promise as antitumor agents in RAS and RAF mutated solid cancers. However, the efficacy of pan-RAF inhibitors in acute myeloid leukemia (AML) has not previously been explored. In AML, the RAS-RAF-MEK-ERK (MAPK) pathway is one of the most aberrantly activated oncogenic pathways, but previous targeting of this pathway by MEK inhibitors has not proven effective in clinical trials. Here we show that pan-RAF inhibition, but not MEK inhibition, induced cell death in 29% of AML samples while being nontoxic toward healthy bone marrow cells. Mechanistically, pan-RAF inhibition downregulated MCL1 protein synthesis and induced apoptosis in cells dependent on MCL1 for their survival. Furthermore, the combination of a pan-RAF and a BCL2 inhibitor overcame resistance to either compound alone in AML cell lines, as well as synergized and induced long-term responses ex vivo in AML patient samples relapsed or refractory to azacitidine + venetoclax treatment. Together, our results indicate that pan-RAF inhibition, alone or in combination with BCL2 inhibition, is a promising treatment strategy for AML.

Adaptive RSK-EphA2-GPRC5A signaling switch triggers chemotherapy resistance in ovarian cancer.

Metastatic cancers commonly activate adaptive chemotherapy resistance, attributed to both microenvironment-dependent phenotypic plasticity and genetic characteristics of cancer cells. However, the contribution of chemotherapy itself to the non-genetic resistance mechanisms was long neglected. Using high-grade serous ovarian cancer (HGSC) patient material and cell lines, we describe here an unexpectedly robust cisplatin and carboplatin chemotherapy-induced ERK1/2-RSK1/2-EphA2-GPRC5A signaling switch associated with cancer cell intrinsic and acquired chemoresistance. Mechanistically, pharmacological inhibition or knockdown of RSK1/2 prevented oncogenic EphA2-S897 phosphorylation and EphA2-GPRC5A co-regulation, thereby facilitating a signaling shift to the canonical tumor-suppressive tyrosine phosphorylation and consequent downregulation of EphA2. In combination with platinum, RSK inhibitors effectively sensitized even the most platinum-resistant EphA2high , GPRC5Ahigh cells to the therapy-induced apoptosis. In HGSC patient tumors, this orphan receptor GPRC5A was expressed exclusively in cancer cells and associated with chemotherapy resistance and poor survival. Our results reveal a kinase signaling pathway uniquely activated by platinum to elicit adaptive resistance. They further identify GPRC5A as a marker for abysmal HGSC outcome and putative vulnerability of the chemo-resistant cells to RSK1/2-EphA2-pS897 pathway inhibition.

Chemical, Physical and Biological Triggers of Evolutionary Conserved Bcl-xL-Mediated Apoptosis.

BACKGROUND: The evidence that pan-Bcl-2 or Bcl-xL-specific inhibitors prematurely kill virus-infected or RNA/DNA-transfected cells provides rationale for investigating these apoptotic inducers further. We hypothesized that not only invasive RNA or DNA (biological factors) but also DNA/RNA-damaging chemical or physical factors could trigger apoptosis that have been sensitized with pan-Bcl-2 or Bcl-xL-specific agents; Methods: We tested chemical and physical factors plus Bcl-xL-specific inhibitor A-1155463 in cells of various origins and the small roundworms (C. elegans); Results: We show that combination of a A-1155463 along with a DNA-damaging agent, 4-nitroquinoline-1-oxide (4NQO), prematurely kills cells of various origins as well as C. elegans. The synergistic effect is p53-dependent and associated with the release of Bad and Bax from Bcl-xL, which trigger mitochondrial outer membrane permeabilization. Furthermore, we found that combining Bcl-xL-specific inhibitors with various chemical compounds or physical insults also induced cell death; Conclusions: Thus, we were able to identify several biological, chemical and physical triggers of the evolutionarily conserved Bcl-xL-mediated apoptotic pathway, shedding light on strategies and targets for novel drug development.

Orphan G protein-coupled receptor GPRC5A modulates integrin ß1-mediated epithelial cell adhesion.

G-Protein Coupled Receptor (GPCR), Class C, Group 5, Member A (GPRC5A) has been implicated in several malignancies. The underlying mechanisms, however, remain poorly understood. Using a panel of human cell lines, we demonstrate that CRISPR/Cas9-mediated knockout and RNAi-mediated depletion of GPRC5A impairs cell adhesion to integrin substrates: collagens I and IV, fibronectin, as well as to extracellular matrix proteins derived from the Engelbreth-Holm-Swarm (EHS) mouse sarcoma (Matrigel). Consistent with the phenotype, knock-out of GPRC5A correlated with a reduced integrin ß1 (ITGB1) protein expression, impaired phosphorylation of the focal adhesion kinase (FAK), and lower activity of small GTPases RhoA and Rac1. Furthermore, we provide the first evidence for a direct interaction between GPRC5A and a receptor tyrosine kinase EphA2, an upstream regulator of FAK, although its contribution to the observed adhesion phenotype is unclear. Our findings reveal an unprecedented role for GPRC5A in regulation of the ITGB1-mediated cell adhesion and it's downstream signaling, thus indicating a potential novel role for GPRC5A in human epithelial cancers.

Antiviral Properties of Chemical Inhibitors of Cellular Anti-Apoptotic Bcl-2 Proteins.

Viral diseases remain serious threats to public health because of the shortage of effective means of control. To combat the surge of viral diseases, new treatments are urgently needed. Here we show that small-molecules, which inhibit cellular anti-apoptotic Bcl-2 proteins (Bcl-2i), induced the premature death of cells infected with different RNA or DNA viruses, whereas, at the same concentrations, no toxicity was observed in mock-infected cells. Moreover, these compounds limited viral replication and spread. Surprisingly, Bcl-2i also induced the premature apoptosis of cells transfected with viral RNA or plasmid DNA but not of mock-transfected cells. These results suggest that Bcl-2i sensitizes cells containing foreign RNA or DNA to apoptosis. A comparison of the toxicity, antiviral activity, and side effects of six Bcl-2i allowed us to select A-1155463 as an antiviral lead candidate. Thus, our results pave the way for the further development of Bcl-2i for the prevention and treatment of viral diseases.

BRCA1-deficient breast cancer cell lines are resistant to MEK inhibitors and show distinct sensitivities to 6-thioguanine.

Germ-line or somatic inactivation of BRCA1 is a defining feature for a portion of human breast cancers. Here we evaluated the anti-proliferative activity of 198 FDA-approved and experimental drugs against four BRCA1-mutant (HCC1937, MDA-MB-436, SUM1315MO2, and SUM149PT) and four BRCA1-wild-type (MDA-MB-231, SUM229PE, MCF10A, and MCF7) breast cancer cell lines. We found that all BRCA1-mutant cell lines were insensitive to inhibitors of mitogen-activated protein kinase kinase 1 and 2 (MEK1/2) Selumetinib and Pimasertib in contrast to BRCA1-wildtype control cell lines. However, unexpectedly, only two BRCA1-mutant cell lines, HCC1937 and MDA-MB-436, were hypersensitive to a nucleotide analogue 6-thioguanine (6-TG). SUM149PT cells readily formed radiation-induced RAD51-positive nuclear foci indicating a functional homologous recombination, which may explain their resistance to 6-TG. However, the reason underlying 6-TG resistance of SUM1315MO2 cells remains unclear. Our data reveal a remarkable heterogeneity among BRCA1-mutant cell lines and provide a reference for future studies.