Clonal dynamics limits detection of selection in tumour xenograft CRISPR/Cas9 screens

Transplantable in vivo CRISPR/Cas9 knockout screens, in which cells are edited in vitro and inoculated into mice to form tumours, allow evaluation of gene function in a cancer model that incorporates the multicellular interactions of the tumour microenvironment. To improve our understanding of the key parameters for success with this method, we investigated the choice of cell line, mouse host, tumour harvesting timepoint and guide RNA (gRNA) library size. We found that high gRNA (80-95%) representation was maintained in a HCT116 subline transduced with the GeCKOv2 whole-genome gRNA library and transplanted into NSG mice when tumours were harvested at early (14 d) but not late time points (38-43 d). The decreased representation in older tumours was accompanied by large increases in variance in gRNA read counts, with notable expansion of a small number of random clones in each sample. The variable clonal dynamics resulted in a high level of 'noise' that limited the detection of gRNA-based selection. Using simulated datasets derived from our experimental data, we show that considerable reductions in count variance would be achieved with smaller library sizes. Based on our findings, we suggest a pathway to rationally design adequately powered in vivo CRISPR screens for successful evaluation of gene function.

Dynamic ctDNA Mutational Complexity in Patients with Melanoma Receiving Immunotherapy

BACKGROUND

Circulating tumour DNA (ctDNA) analysis promises to improve the clinical care of people with cancer, address health inequities and guide translational research. This observational cohort study used ctDNA to follow 29 patients with advanced-stage cutaneous melanoma through multiple cycles of immunotherapy.

METHOD

A melanoma-specific ctDNA next-generation sequencing (NGS) panel, droplet digital polymerase chain reaction (ddPCR) and mass spectrometry analysis were used to identify ctDNA mutations in longitudinal blood plasma samples from Aotearoa New Zealand (NZ) patients receiving immunotherapy for melanoma. These technologies were used in conjunction to identify the breadth and complexity of tumour genomic information that ctDNA analysis can reliably report.

RESULTS

During the course of immunotherapy treatment, a high level of dynamic mutational complexity was identified in blood plasma, including multiple BRAF mutations in the same patient, clinically relevant BRAF mutations emerging through therapy and co-occurring sub-clonal BRAF and NRAS mutations. The technical validity of this ctDNA analysis was supported by high sample analysis-reanalysis concordance, as well as concordance between different ctDNA measurement technologies. In addition, we observed > 90% concordance in the detection of ctDNA when using cell-stabilising collection tubes followed by 7-day delayed processing, compared with standard EDTA blood collection protocols with rapid processing. We also found that the undetectability of ctDNA at a proportion of treatment cycles was associated with durable clinical benefit (DCB).

CONCLUSION

We found that multiple ctDNA processing and analysis methods consistently identified complex longitudinal patterns of clinically relevant mutations, adding support for expanded clinical trials of this technology in a variety of oncology settings.

Chromosomal Aberrations Accumulate during Metastasis of Virus-Negative Merkel Cell Carcinoma

Merkel cell carcinoma is a rare, aggressive skin tumor initiated by polyomavirus integration or UV light DNA damage. In New Zealand, there is a propensity toward the UV-driven form (31 of 107, 29% virus positive). Using archival formalin-fixed, paraffin-embedded tissues, we report targeted DNA sequencing covering 246 cancer genes on 71 tumor tissues and 38 nonmalignant tissues from 37 individuals, with 33 of 37 being negative for the virus. Somatic variants of New Zealand virus-negative Merkel cell carcinomas partially overlapped with those reported overseas, including TP53 variants in all tumors and RB1, LRP1B, NOTCH1, and EPHA3/7 variants each found in over half of the cohort. Variants in genes not analyzed or reported in previous studies were also found. Cataloging variants in TP53 and RB1 from published datasets revealed a broad distribution across these genes. Chr 1p gain and Chr 3p loss were identified in around 50% of New Zealand virus-negative Merkel cell carcinomas, and RB1 loss of heterozygosity was found in 90% of cases. Copy number variants accumulate in most metastases. Virus-negative Merkel cell carcinomas have complex combinations of somatic DNA-sequence variants and copy number variants. They likely carry the small genomic changes permissive for metastasis from early tumor development; however, chromosomal alterations may contribute to driving metastatic progression.

Complex patterns of genomic heterogeneity identified in 42 tumor samples and ctDNA of a pulmonary atypical carcinoid patient

Tumor evolution underlies many challenges facing precision oncology, and improving our understanding has the potential to improve clinical care. This study represents a rare opportunity to study tumor heterogeneity and evolution in a patient with an understudied cancer type. A patient with pulmonary atypical carcinoid, a neuroendocrine tumor, metastatic to 90 sites, requested and consented to donate tissues for research. 42 tumor samples collected at rapid autopsy from 14 anatomically distinct sites were analyzed through DNA whole-exome sequencing and RNA-Seq, and five analyzed through linked-read sequencing. Targeted DNA sequencing was completed on two clinical tissue biopsies and one blood plasma sample. Chromosomal alterations and gene variants accumulated over time, and specific chromosomal alterations preceded the single predicted gene driver variant (ARID1A). At the time of autopsy, all sites shared the gain of one copy of Chr 5, loss of one copy of Chr 6 and 21, chromothripsis of one copy of Chr 11, and 39 small variants. Two tumor clones (carrying additional variants) were detected at metastatic sites, and occasionally in different regions of the same organ (e.g., within the pancreas). Circulating tumor DNA (ctDNA) sequencing detected shared tumor variants in the blood plasma and captured marked genomic heterogeneity, including all metastatic clones but few private tumor variants. This study describes genomic tumor evolution and dissemination of a pulmonary atypical carcinoid donated by a single generous patient. It highlights the critical role of chromosomal alterations in tumor initiation and explores the potential of ctDNA analysis to represent genomically heterogeneous disease.

Clinically Relevant Biology of Hyaluronic Acid in the Desmoplastic Stroma of Pancreatic Ductal Adenocarcinoma

ABSTRACT

Pancreatic ductal adenocarcinoma (PDAC) is notorious for its poor outcome. The presence of a dense desmoplastic stroma is a hallmark of this malignancy, and abundant hyaluronic acid (HA) within this stroma is a common feature of PDAC. At the end of 2019, an HA-targeting drug, after initial promise, failed phase 3 clinical trials in PDAC. This failure in the face of such strong evidence for biological importance forces us to turn back to the research and seek a better understanding of HA biology in PDAC. Therefore, in this review, we reexamine what is known about HA biology, the methods used to detect and quantify HA, and the ability of the biological models in which HA has been investigated to recapitulate an HA-rich desmoplastic tumor stroma. The role of HA in PDAC relies on its complex interplay with a range of HA-associated molecules, which have not been as extensively investigated as HA itself. Therefore, using large genomic data sets, we cataloged the abundance and activity in PDAC of molecules that modulate HA synthesis, degradation, protein interactions, and receptor binding. Based on their association with clinical characteristics and individual patient outcomes, we suggest a small number of HA-associated molecules that warrant further investigation as biomarkers and drug targets.

Innate immune checkpoint inhibitor resistance is associated with melanoma sub-types exhibiting invasive and de-differentiated gene expression signatures

Melanoma is a highly aggressive skin cancer, which, although highly immunogenic, frequently escapes the body’s immune defences. Immune checkpoint inhibitors (ICI), such as anti-PD1, anti-PDL1, and anti-CTLA4 antibodies lead to reactivation of immune pathways, promoting rejection of melanoma. However, the benefits of ICI therapy remain limited to a relatively small proportion of patients who do not exhibit ICI resistance. Moreover, the precise mechanisms underlying innate and acquired ICI resistance remain unclear. Here, we have investigated differences in melanoma tissues in responder and non-responder patients to anti-PD1 therapy in terms of tumour and immune cell gene-associated signatures. We performed multi-omics investigations on melanoma tumour tissues, which were collected from patients before starting treatment with anti-PD1 immune checkpoint inhibitors. Patients were subsequently categorized into responders and non-responders to anti-PD1 therapy based on RECIST criteria. Multi-omics analyses included RNA-Seq and NanoString analysis. From RNA-Seq data we carried out HLA phenotyping as well as gene enrichment analysis, pathway enrichment analysis and immune cell deconvolution studies. Consistent with previous studies, our data showed that responders to anti-PD1 therapy had higher immune scores (median immune score for responders = 0.1335, median immune score for non-responders = 0.05426, p-value = 0.01, Mann-Whitney U two-tailed exact test) compared to the non-responders. Responder melanomas were more highly enriched with a combination of CD8+ T cells, dendritic cells (p-value = 0.03) and an M1 subtype of macrophages (p-value = 0.001). In addition, melanomas from responder patients exhibited a more differentiated gene expression pattern, with high proliferative- and low invasive-associated gene expression signatures, whereas tumours from non-responders exhibited high invasive- and frequently neural crest-like cell type gene expression signatures. Our findings suggest that non-responder melanomas to anti-PD1 therapy exhibit a de-differentiated gene expression signature, associated with poorer immune cell infiltration, which establishes a gene expression pattern characteristic of innate resistance to anti-PD1 therapy. Improved understanding of tumour-intrinsic gene expression patterns associated with response to anti-PD1 therapy will help to identify predictive biomarkers of ICI response and may help to identify new targets for anticancer treatment, especially with a capacity to function as adjuvants to improve ICI outcomes.

Tuberous sclerosis complex: a complex case

Tuberous sclerosis complex (TSC) is an inheritable disorder characterized by the formation of benign yet disorganized tumors in multiple organ systems. Germline mutations in the TSC1 (hamartin) or more frequently TSC2 (tuberin) genes are causative for TSC. The malignant manifestations of TSC, pulmonary lymphangioleiomyomatosis (LAM) and renal angiomyolipoma (AML), may also occur as independent sporadic perivascular epithelial cell tumor (PEComa) characterized by somatic TSC2 mutations. Thus, discerning TSC from the copresentation of sporadic LAM and sporadic AML may be obscured in TSC patients lacking additional features. In this report, we present a case study on a single patient initially reported to have sporadic LAM and a mucinous duodenal adenocarcinoma deficient in DNA mismatch repair proteins. Moreover, the patient had a history of Wilms' tumor, which was reclassified as AML following the LAM diagnosis. Therefore, we investigated the origins and relatedness of these tumors. Using germline whole-genome sequencing, we identified a premature truncation in one of the patient's TSC2 alleles. Using immunohistochemistry, loss of tuberin expression was observed in AML and LAM tissue. However, no evidence of a somatic loss of heterozygosity or DNA methylation epimutations was observed at the TSC2 locus, suggesting alternate mechanisms may contribute to loss of the tumor suppressor protein. In the mucinous duodenal adenocarcinoma, no causative mutations were found in the DNA mismatch repair genes MLH1, MSH2, MSH6, or PMS2 Rather, clonal deconvolution analyses were used to identify mutations contributing to pathogenesis. This report highlights both the utility of using multiple sequencing techniques and the complexity of interpreting the data in a clinical context.

Breast Cancer Patient Prognosis Is Determined by the Interplay between TP53 Mutation and Alternative Transcript Expression: Insights from TP53 Long Amplicon Digital PCR Assays

Simple Summary

The TP53 gene, the most commonly mutated gene in human cancers, is capable of producing multiple RNAs (transcripts). The aim of our study was to measure the abundance of each TP53 transcript, combined with TP53 gene mutation information, to determine the interplay between these in a cohort of breast tumors from New Zealand patients. To do this, we devised a new assay which then enabled the measurement of all known TP53 transcripts. We showed how TP53 gene mutations influenced the levels of specific TP53 transcripts in breast tumors. We evaluated whether a combination of TP53 tumor information, including TP53 mutation status and the levels of certain TP53 transcripts, with standard clinical and pathological information, was associated with breast cancer patient outcome. We recommend that a truly comprehensive analysis of TP53 needs to incorporate data about both TP53 DNA mutations and the expression of the alternative TP53 transcripts.

Abstract

The TP53 gene locus is capable of producing multiple RNA transcripts encoding the different p53 protein isoforms. We recently described multiplex long amplicon droplet digital PCR (ddPCR) assays to quantify seven of eight TP53 reference transcripts in human tumors. Here, we describe a new long amplicon ddPCR assay to quantify expression of the eighth TP53 reference transcript encoding ∆40p53α. We then applied these assays, alongside DNA sequencing of the TP53 gene locus, to tumors from a cohort of New Zealand (NZ) breast cancer patients. We found a high prevalence of mutations at TP53 splice sites in the NZ breast cancer cohort. Mutations at TP53 intron 4 splice sites were associated with overexpression of ∆133TP53 transcripts. Cox proportional hazards survival analysis showed that interplay between TP53 mutation status and expression of TP53 transcript variants was significantly associated with patient outcome, over and above standard clinical and pathological information. In particular, patients with no TP53 mutation and a low ratio of TP53 transcripts t2 to t1, which derive from alternative intron 1 acceptor splice sites, had a remarkably good outcome. We suggest that this type of analysis, integrating mutation and transcript expression, provides a step-change in our understanding of TP53 in cancer.

Specialized Cell-Free DNA Blood Collection Tubes Can Be Repurposed for Extracellular Vesicle Isolation: A Pilot Study

Background: Liquid biopsies offer a minimally invasive approach to patient disease diagnosis and monitoring. However, these are highly affected by preprocessing variables with many protocols designed for downstream analysis of a single molecular biomarker. Here we investigate whether specialized blood tubes could be repurposed for the analysis of an increasingly valuable biomarker, extracellular vesicles (EVs). Methods: Blood was collected from three donors into K3-EDTA, Roche, or Streck cell-free DNA (cfDNA) collection tubes and processed using sequential centrifugation either immediately or after storage for 3 days. MicroEV were collected from platelet-poor plasma by 10,000 g centrifugation and NanoEVs isolated using size exclusion chromatography. Particle size and counts were assessed by Nanoparticle Tracking Analysis, protein quantitation by bicinchoninic acid assay (BCA) assay, and dot blotting for blood cell surface proteins. Results: MicroEVs and NanoEVs could be isolated from plasma collected using all three tube types. Major variations were seen with delayed time to processing. Both MicroEV particle number and protein content increased with the processing delay. The NanoEV number did not change with the time-delay but their protein quantity increased. EV-associated proteins predominantly arose from platelets (CD61) and erythrocytes (CD235a). However, leukocyte marker CD45 was only increased in NanoEVs from ethylenediaminetetraacetic acid (EDTA) tubes, suggestive of stabilization of nucleated cells by the specialized blood tubes. Epithelial cell surface marker EpCAM, often used as a marker of cancer, remained the same across conditions in both MicroEV and NanoEV preparations indicating that these EVs were stable with time. Conclusions: Specialized cfDNA collection tubes can be repurposed for MicroEV and NanoEV analysis; however, simple counting or using protein quantity as a surrogate of EV number may be confounded by preanalytical processing. The EVs would be suitable for disease selective EV subtype analysis if the molecular target of interest is not present in blood cells.

Genomic and signalling pathway characterization of the NZM panel of melanoma cell lines: A valuable model for studying the impact of genetic diversity in melanoma

Melanoma is a disease associated with a very high mutation burden and thus the possibility of a diverse range of oncogenic mechanisms that allow it to evade therapeutic interventions and the immune system. Here, we describe the characterization of a panel of 102 cell lines from metastatic melanomas (the NZM lines), including using whole‐exome and RNA sequencing to analyse genetic variants and gene expression changes in a subset of this panel. Lines possessing all major melanoma genotypes were identified, and hierarchical clustering of gene expression profiles revealed four broad subgroups of cell lines. Immunogenotyping identified a range of HLA haplotypes as well as expression of neoantigens and cancer–testis antigens in the lines. Together, these characteristics make the NZM panel a valuable resource for cell‐based, immunological and xenograft studies to better understand the diversity of melanoma biology and the responses of melanoma to therapeutic interventions.

Accessing a New Dimension in TP53 Biology: Multiplex Long Amplicon Digital PCR to Specifically Detect and Quantitate Individual TP53 Transcripts

TP53, the most commonly-mutated gene in cancer, undergoes complex alternative splicing. Different TP53 transcripts play different biological roles, both in normal function and in the progression of diseases such as cancer. The study of TP53’s alternative RNA splice forms and their use as clinical biomarkers has been hampered by limited specificity and quantitative accuracy of current methods. TP53 RNA splice variants differ at both 5’ and 3’ ends, but because they have a common central region of 618 bp, the individual TP53 transcripts are impossible to specifically detect and precisely quantitate using standard PCR-based methods or short-read RNA sequencing. Therefore, we devised multiplex probe-based long amplicon droplet digital PCR (ddPCR) assays, which for the first time allow precise end-to-end quantitation of the seven major TP53 transcripts, with amplicons ranging from 0.85 to 1.85 kb. Multiple modifications to standard ddPCR assay procedures were required to enable specific co-amplification of these long transcripts and to overcome issues with secondary structure. Using these assays, we show that several TP53 transcripts are co-expressed in breast cancers, and illustrate the potential for this method to identify novel TP53 transcripts in tumour cells. This capability will facilitate a new level of biological and clinical understanding of the alternatively-spliced TP53 isoforms.

Transcriptomic Features of T Cell-Barren Tumors Are Conserved Across Diverse Tumor Types

Background: Understanding how tumors subvert immune destruction is essential to the development of cancer immunotherapies. New evidence suggests that tumors limit anti-tumor immunity by exploiting transcriptional programs that regulate intratumoral trafficking and accumulation of effector cells. Here, we investigated the gene expression profiles that distinguish immunologically “cold” and “hot” tumors across diverse tumor types.

Methods: RNAseq profiles of tumors (n = 8,920) representing 23 solid tumor types were analyzed using immune gene signatures that quantify CD8+ T cell abundance. Genes and pathways associated with a low CD8+ T cell infiltration profile (CD8-Low) were identified by correlation, differential expression, and statistical ranking methods. Gene subsets were evaluated in immunotherapy treatment cohorts and functionally characterized in cell lines and mouse tumor models.

Results: Among different cancer types, we observed highly significant overlap of genes enriched in CD8-Low tumors, which included known immunomodulatory genes (e.g., BMP7, CMTM4, KDM5B, RCOR2) and exhibited significant associations with Wnt signaling, neurogenesis, cell-cell junctions, lipid biosynthesis, epidermal development, and cancer-testis antigens. Analysis of mutually exclusive gene clusters demonstrated that different transcriptional programs may converge on the T cell-cold phenotype as well as predict for response and survival of patients to Nivo treatment. Furthermore, we confirmed that a top-ranking candidate belonging to the TGF-β superfamily, BMP7, negatively regulates CD8+ T cell abundance in immunocompetent murine tumor models, with and without anti-PD-L1 treatment.

Conclusions: This study presents the first evidence that solid tumors of diverse anatomical origin acquire conserved transcriptional alterations that may be operative in the T cell-cold state. Our findings demonstrate the potential clinical utility of CD8-Low tumor-associated genes for predicting patient immunotherapy outcomes and point to novel mechanisms with potential for broad therapeutic exploitation.

Plasma miRNAs Display Limited Potential as Diagnostic Tools for Endometriosis

CONTEXT

Despite extensive searches for novel noninvasive diagnostics, laparoscopy remains the reference test for endometriosis. Circulating miRNAs are purported endometriosis biomarkers; however, the miRNA species and their diagnostic accuracy differ between studies and have not been validated in independent cohorts.

OBJECTIVE

Identify endometriosis-specific plasma miRNAs and determine their diagnostic test accuracy.

SETTING

Two university-based, public hospitals and a private gynecology practice in Australia.

DESIGN AND PARTICIPANTS

Four phases: (i) Explorative phase. Plasma miRNA menstrual cycle fluctuations were evaluated in women with endometriosis and asymptomatic controls (n = 16). (ii) Biomarker discovery. Endometriosis-specific plasma miRNAs were identified in (a) women with endometriosis and asymptomatic controls (n = 16) and (b) women with and without surgically defined endometriosis (n = 20). (iii) Biomarker selection. Plasma miRNAs with the best diagnostic potential for endometriosis were selected in a surgically defined selection cohort (n = 78). (iv) Biomarker validation. The diagnostic test accuracy of these miRNAs was calculated in an independent, surgically defined validation cohort (n = 119).

RESULTS

Forty-nine miRNAs were differentially expressed in women with endometriosis. Nine maintained dysregulation in the selection cohort, but only three (miR-155, miR574-3p and miR139-3p) did so in the validation cohort. Combined, these three miRNAs demonstrated a sensitivity and specificity of 83% and 51%, respectively.

CONCLUSION

Plasma miRNAs demonstrated modest sensitivity and specificity as diagnostic tests or triage tools for endometriosis. Other groups' findings were not replicated and accorded poorly with our results. Circulating miRNAs demonstrate diagnostic potential, but stringent, standardized methodological approaches are required for the development of a clinically applicable tool.

Functional CRISPR and shRNA Screens Identify Involvement of Mitochondrial Electron Transport in the Activation of Evofosfamide

Evofosfamide (TH-302) is a hypoxia-activated DNA-crosslinking prodrug currently in clinical development for cancer therapy. Oxygen-sensitive activation of evofosfamide depends on one-electron reduction, yet the reductases that catalyze this process in tumors are unknown. We used RNA sequencing, whole-genome CRISPR knockout, and reductase-focused short hairpin RNA screens to interrogate modifiers of evofosfamide activation in cancer cell lines. Involvement of mitochondrial electron transport in the activation of evofosfamide and the related nitroaromatic compounds EF5 and FSL-61 was investigated using 143B ρ ⁰ (ρ zero) cells devoid of mitochondrial DNA and biochemical assays in UT-SCC-74B cells. The potency of evofosfamide in 30 genetically diverse cancer cell lines correlated with the expression of genes involved in mitochondrial electron transfer. A whole-genome CRISPR screen in KBM-7 cells identified the DNA damage-response factors SLX4IP, C10orf90 (FATS), and SLFN11, in addition to the key regulator of mitochondrial function, YME1L1, and several complex I constituents as modifiers of evofosfamide sensitivity. A reductase-focused shRNA screen in UT-SCC-74B cells similarly identified mitochondrial respiratory chain factors. Surprisingly, 143B ρ ⁰ cells showed enhanced evofosfamide activation and sensitivity but had global transcriptional changes, including increased expression of nonmitochondrial flavoreductases. In UT-SCC-74B cells, evofosfamide oxidized cytochromes a, b, and c and inhibited respiration at complexes I, II, and IV without quenching reactive oxygen species production. Our results suggest that the mitochondrial electron transport chain contributes to evofosfamide activation and that predicting evofosfamide sensitivity in patients by measuring the expression of canonical bioreductive enzymes such as cytochrome P450 oxidoreductase is likely to be futile.

Imprinted and ancient gene: a potential mediator of cancer cell survival during tryptophan deprivation

BACKGROUND

Depletion of tryptophan and the accumulation of tryptophan metabolites mediated by the immunosuppressive enzyme indoleamine 2,3-dioxygenase 1 (IDO1), trigger immune cells to undergo apoptosis. However, cancer cells in the same microenvironment appear not to be affected. Mechanisms whereby cancer cells resist accelerated tryptophan degradation are not completely understood. We hypothesize that cancer cells co-opt IMPACT (the product of IMPrinted and AnCienT gene), to withstand periods of tryptophan deficiency.

METHODS

A range of bioinformatic techniques including correlation and gene set variation analyses was applied to genomic datasets of cancer (The Cancer Genome Atlas) and normal (Genotype Tissue Expression Project) tissues to investigate IMPACT's role in cancer. Survival of IMPACT-overexpressing GL261 glioma cells and their wild type counterparts cultured in low tryptophan media was assessed using fluorescence microscopy and MTT bio-reduction assay. Expression of the Integrated Stress Response proteins was measured using Western blotting.

RESULTS

We found IMPACT to be upregulated and frequently amplified in a broad range of clinical cancers relative to their non-malignant tissue counterparts. In a subset of clinical cancers, high IMPACT expression associated with decreased activity of pathways and genes involved in stress response and with increased activity of translational regulation such as the mTOR pathway. Experimental studies using the GL261 glioma line showed that cells engineered to overexpress IMPACT, gained a survival advantage over wild-type lines when cultured under limiting tryptophan concentrations. No significant difference in the expression of proteins in the Integrated Stress Response pathway was detected in tryptophan-deprived GL261 IMPACT-overexpressors compared to that in wild-type cells. IMPACT-overexpressing GL261 cells but not their wild-type counterparts, showed marked enlargement of their nuclei and cytoplasmic area when stressed by tryptophan deprivation.

CONCLUSIONS

The bioinformatics data together with our laboratory studies, support the hypothesis that IMPACT mediates a protective mechanism allowing cancer cells to overcome microenvironmental stresses such as tryptophan deficiency.

Evofosfamide for the treatment of human papillomavirus-negative head and neck squamous cell carcinoma

Evofosfamide (TH-302) is a clinical-stage hypoxia-activated prodrug of a DNA-crosslinking nitrogen mustard that has potential utility for human papillomavirus (HPV) negative head and neck squamous cell carcinoma (HNSCC), in which tumor hypoxia limits treatment outcome. We report the preclinical efficacy, target engagement, preliminary predictive biomarkers and initial clinical activity of evofosfamide for HPV-negative HNSCC. Evofosfamide was assessed in 22 genomically characterized cell lines and 7 cell line-derived xenograft (CDX), patient-derived xenograft (PDX), orthotopic, and syngeneic tumor models. Biomarker analysis used RNA sequencing, whole-exome sequencing, and whole-genome CRISPR knockout screens. Five advanced/metastatic HNSCC patients received evofosfamide monotherapy (480 mg/m2 qw × 3 each month) in a phase 2 study. Evofosfamide was potent and highly selective for hypoxic HNSCC cells. Proliferative rate was a predominant evofosfamide sensitivity determinant and a proliferation metagene correlated with activity in CDX models. Evofosfamide showed efficacy as monotherapy and with radiotherapy in PDX models, augmented CTLA-4 blockade in syngeneic tumors, and reduced hypoxia in nodes disseminated from an orthotopic model. Of 5 advanced HNSCC patients treated with evofosfamide, 2 showed partial responses while 3 had stable disease. In conclusion, evofosfamide shows promising efficacy in aggressive HPV-negative HNSCC, with predictive biomarkers in development to support further clinical evaluation in this indication.

Regulation of the interferon-gamma (IFN-γ) pathway by p63 and Δ133p53 isoform in different breast cancer subtypes

The TP53 family consists of three sets of transcription factor genes, TP53, TP63 and TP73, each of which expresses multiple RNA variants and protein isoforms. Of these, TP53 is mutated in 25-30% of breast cancers. How TP53 mutations affect the interaction of TP53 family members and their isoforms in breast cancer is unknown. To investigate this, 3 independent breast cancer cohorts were stratified into 4 groups based on oestrogen receptor (ER) and TP53 mutation status. Using bioinformatic methodologies, principal signalling pathways associated with the expression of TP53 family members were identified. Results show an enrichment of IFN-γ signalling associated with TP63 RNA in wild type TP53 (wtTP53), ER negative (ER-) tumours and with Δ133TP53 RNA in mutant TP53 (mTP53) ER positive (ER+) tumours. Moreover, tumours with low IFN-γ signalling were associated with significantly poorer patient outcome. The predicted changes in expression of a subset of RNAs involved in IFN-γ signalling were confirmed in vitro. Our data show that different members of the TP53 family can drive transcription of genes involved in IFN-γ signalling in different breast cancer subgroups.

Gene expression profiling of breast tumours from New Zealand patients

AIMS

New Zealand has one of the highest rates of breast cancer incidence in the world. We investigated the gene expression profiles of breast tumours from New Zealand patients, compared them to gene expression profiles of international breast cancer cohorts and identified any associations between altered gene expression and the clinicopathological features of the tumours.

METHODS

Affymetrix microarrays were used to measure the gene expression profiles of 106 breast tumours from New Zealand patients. Gene expression data from six international breast cancer cohorts were collated, and all the gene expression data were analysed using standard bioinformatic and statistical tools.

RESULTS

Gene expression profiles associated with tumour ER and ERBB2 status, molecular subtype and selected gene expression signatures within the New Zealand cohort were consistent with those found in international cohorts. Significant differences in clinicopathological features such as tumour grade, tumour size and lymph node status were also observed between the New Zealand and international cohorts.

CONCLUSIONS

Gene expression profiles, which are a sensitive indicator of tumour biology, showed no clear difference between breast tumours from New Zealand patients and those from non-New Zealand patients. This suggests that other factors may contribute to the high and increasing breast cancer incidence in New Zealand compared to international populations.

Abstract 169: Preclinical efficacy and sensitivity determinants of evofosfamide in molecularly defined models of head and neck squamous cell carcinoma

Tumor hypoxia is prevalent in head and neck squamous cell carcinoma (HNSCC), where it limits radiotherapy outcomes. Hypoxia-activated prodrugs (HAPs) have been developed to target hypoxic regions of tumors. These agents undergo oxygen-sensitive reductive activation, thereby delivering cytotoxic species within hypoxic cells. This study investigated the efficacy and sensitivity determinants of the clinical-stage HAP evofosfamide (TH-302) using molecularly-characterized models of HNSCC. We deployed a collection of 27 HPV-negative HNSCC cell lines derived from lesions of varying TNM stages and primary, nodal or recurrent sites. The collection was characterized for gene expression by RNA-seq, from which somatic variants were also called. Their transcriptomic features were investigated in the context of pan-cancer TCGA data by hierarchical clustering. The potency and hypoxic selectivity of 3 HAPs - evofosfamide, PR-104A and SN30000 - were assessed by antiproliferative assay in 22 lines and compared to bromo-isophosphoramide mustard (Br-IPM), cisplatin and 5-FU. The antitumor activity of evofosfamide (50 mg/kg qdx5 for 2-3 cycles with or without a single 10 Gy dose of radiation on day 5 of cycle 1) was evaluated in HNSCC xenografts in addition to a PDX isolated from an SCC of the glottic larynx. The hypoxic fraction at baseline and after 5 days of treatment was quantified by pimonidazole staining. Genetic modifiers of sensitivity to evofosfamide and its cytotoxic metabolite Br-IPM were explored through whole-genome CRISPR-Cas9 screens using the GeCKO v2 library. High-throughput screens with a custom shRNA pool were performed in one HNSCC and two pancreatic ductal adenocarcinoma cell lines to identify reductases responsible for the activation of evofosfamide in hypoxic cells. Evofosfamide was more potent and more selective for hypoxic HNSCC cells in vitro than PR-104A or SN30000. Cell line sensitivity to evofosfamide was correlated with Br-IPM and cisplatin but not with PR-104A, SN30000 or 5-FU, indicating distinct sensitivity determinants. Evidence of antitumor activity with evofosfamide was observed in vivo. CRISPR screens identified potential evofosfamide sensitivity genes that were reproducibly enriched following drug exposure. Reductase-focused RNA interference screens defined a cluster of sensitivity genes that mapped to mitochondrial electron transport, whereas shRNA’s targeted against presumed activating enzymes such as POR were not enriched. Concentration-dependent oxidation of cytochrome a and decreased respiration was observed in cells exposed to evofosfamide, suggesting reduction by mitochondrial complexes. This study provides a rationale for the clinical evaluation of evofosfamide with radiotherapy in genetically defined subsets of HNSCC patients.

Citation Format: Francis W. Hunter, Avik Shome, Dan Li, Way W. Wong, Peter Tsai, Nooriyah Poonawala, Purvi M. Kakadiya, Troy M. Ketelä, Maria K. Kondratyev, Courtney R. Lynch, Tet-Woo Lee, Khanh B. Tran, Jules B. Devaux, Rachel Zussman, Cho R. Hong, Dennis Kee, Andrew M. Macann, Anthony J. Hickey, Stefan K. Bohlander, Cristin G. Print, William R. Wilson, Bradly G. Wouters, Stephen M. Jamieson. Preclinical efficacy and sensitivity determinants of evofosfamide in molecularly defined models of head and neck squamous cell carcinoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 169. doi:10.1158/1538-7445.AM2017-169

Sperm microRNA Content Is Altered in a Mouse Model of Male Obesity, but the Same Suite of microRNAs Are Not Altered in Offspring's Sperm

The prevalence of obesity is increasing worldwide and has tripled in men of reproductive age since the 1970s. Concerningly, obesity is not only comorbid with other chronic diseases, but there is mounting evidence that it increases the non-communicable disease load in their children (eg mortality, obesity, autism). Animal studies have demonstrated that paternal obesity increases the risk of metabolic (eg glucose metabolism defects, obesity) and reproductive disorders in offspring. Epigenetic changes within sperm are clear mechanistic candidates that are associated with both changes to the father’s environment and offspring phenotype. Specifically there is emerging evidence that a father’s sperm microRNA content both responds to paternal environmental cues and alters the gene expression profile and subsequent development of the early embryo. We used a mouse model of high fat diet (HFD) induced obesity to investigate whether male obesity could modulate sperm microRNA content. We also investigated whether this alteration to a father’s sperm microRNA content lead to a similar change in the sperm of male offspring. Our investigations were initially guided by a Taqman PCR array, which indicated the differential abundance of 28 sperm borne microRNAs in HFD mice. qPCR confirmation in a much larger cohort of founder males demonstrated that 13 of these microRNAs were differentially abundant (11 up-regulated; 2 down-regulated) due to HFD feeding. Despite metabolic and reproductive phenotypes also being observed in grand-offspring fathered via the male offspring lineage, there was no evidence that any of the 13 microRNAs were also dysregulated in male offspring sperm. This was presumably due to the variation seen within both groups of offspring and suggests other mechanisms might act between offspring and grand-offspring. Thus 13 sperm borne microRNAs are modulated by a father’s HFD and the presumed transfer of this altered microRNA payload to the embryo at fertilisation potentially acts to alter the embryonic molecular makeup post-fertilisation, altering its growth trajectory, ultimately affecting adult offspring phenotype and may contribute to paternal programming.

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.

The in vitro and in vivo effects of constitutive light expression on a bioluminescent strain of the mouse enteropathogen Citrobacter rodentium

Bioluminescent reporter genes, such as those from fireflies and bacteria, let researchers use light production as a non-invasive and non-destructive surrogate measure of microbial numbers in a wide variety of environments. As bioluminescence needs microbial metabolites, tagging microorganisms with luciferases means only live metabolically active cells are detected. Despite the wide use of bioluminescent reporter genes, very little is known about the impact of continuous (also called constitutive) light expression on tagged bacteria. We have previously made a bioluminescent strain of Citrobacter rodentium, a bacterium which infects laboratory mice in a similar way to how enteropathogenic Escherichia coli (EPEC) and enterohaemorrhagic E. coli (EHEC) infect humans. In this study, we compared the growth of the bioluminescent C. rodentium strain ICC180 with its non-bioluminescent parent (strain ICC169) in a wide variety of environments. To understand more about the metabolic burden of expressing light, we also compared the growth profiles of the two strains under approximately 2,000 different conditions. We found that constitutive light expression in ICC180 was near-neutral in almost every non-toxic environment tested. However, we also found that the non-bioluminescent parent strain has a competitive advantage over ICC180 during infection of adult mice, although this was not enough for ICC180 to be completely outcompeted. In conclusion, our data suggest that constitutive light expression is not metabolically costly to C. rodentium and supports the view that bioluminescent versions of microbes can be used as a substitute for their non-bioluminescent parents to study bacterial behaviour in a wide variety of environments.

A novel EGR-1 dependent mechanism for YB-1 modulation of paclitaxel response in a triple negative breast cancer cell line

Chemotherapy with taxanes such as paclitaxel (PTX) is a key component of triple negative breast cancer (TNBC) treatment. PTX is used in combination with other drugs in both the adjuvant setting and in advanced breast cancer. Because a proportion of patients respond poorly to PTX or relapse after its use, a greater understanding of the mechanisms conferring resistance to PTX is required. One protein shown to be involved in drug resistance is Y-box binding protein 1 (YB-1). High levels of YB-1 have previously been associated with resistance to PTX in TNBCs. In this study, we aimed to determine mechanisms by which YB-1 confers PTX resistance. We generated isogenic TNBC cell lines that differed by YB-1 levels and treated these with PTX. Using microarray analysis, we identified EGR1 as a potential target of YB-1. We found that low EGR1 mRNA levels are associated with poor breast cancer patient prognosis, and that EGR1 and YBX1 mRNA expression was inversely correlated in a TNBC line and in a proportion of TNBC tumours. Reducing the levels of EGR1 caused TNBC cells to become more resistant to PTX. Given that PTX targets cycling cells, we propose a model whereby high YB-1 levels in some TNBC cells can lead to reduced levels of EGR1, which in turn promotes slow cell cycling and resistance to PTX. Therefore YB-1 and EGR1 levels are biologically linked and may provide a biomarker for TNBC response to PTX.

Spatially transformed fluorescence image data for ERK-MAPK and selected proteins within human epidermis

Background

Phosphoprotein signalling pathways have been intensively studied in vitro, yet their role in regulating tissue homeostasis is not fully understood. In the skin, interfollicular keratinocytes differentiate over approximately 2 weeks as they traverse the epidermis. The extracellular signal-regulated kinase (ERK) branch of the mitogen-activated protein kinase (MAPK) pathway has been implicated in this process. Therefore, we examined ERK-MAPK activity within human epidermal keratinocytes in situ.

Findings

We used confocal microscopy and immunofluorescence labelling to measure the relative abundances of Raf-1, MEK1/2 and ERK1/2, and their phosphorylated (active) forms within three human skin samples. Additionally, we measured the abundance of selected proteins thought to modulate ERK-MAPK activity, including calmodulin, β1 integrin and stratifin (14-3-3σ); and of transcription factors known to act as effectors of ERK1/2, including the AP-1 components Jun-B, Fra2 and c-Fos. Imaging was performed with sufficient resolution to identify the plasma membrane, cytoplasm and nucleus as distinct domains within cells across the epidermis. The image field of view was also sufficiently large to capture the entire epidermis in cross-section, and thus the full range of keratinocyte differentiation in a single observation. Image processing methods were developed to quantify image data for mathematical and statistical analysis. Here, we provide raw image data and processed outputs.

Conclusions

These data indicate coordinated changes in ERK-MAPK signalling activity throughout the depth of the epidermis, with changes in relative phosphorylation-mediated signalling activity occurring along the gradient of cellular differentiation. We believe these data provide unique information about intracellular signalling as they are obtained from a homeostatic human tissue, and they might be useful for investigating intercellular heterogeneity.

Electronic supplementary material

The online version of this article (doi:10.1186/s13742-015-0102-5) contains supplementary material, which is available to authorized users.

Antiphospholipid antibodies bind syncytiotrophoblast mitochondria and alter the proteome of extruded syncytial nuclear aggregates

INTRODUCTION

Antiphospholipid antibodies (aPL) are autoantibodies that increase the risk of women developing the hypertensive disorder pre-eclampsia. aPL are internalised by the syncytiotrophoblast and increase extrusion of necrotic multinucleated syncytial nuclear aggregates (SNAs), which may trigger endothelial dysfunction in pre-eclampsia. The mechanisms by which aPL alter death processes in the syncytiotrophoblast leading to extrusion of SNAs are unknown.

METHODS

First trimester human placentae (n = 10) were dissected into explants and cultured either with aPL (50 μg/mL), isotype-matched control antibody (50 μg/mL), or media for 24 h. Harvested SNAs underwent iTRAQ proteomic analysis. Mitochondria in syncytiotrophoblast treated with aPL labelled with FluoroNanogold were visualised using transmission electron microscopy (TEM).

RESULTS

aPL altered the expression of 72 proteins in SNAs. Thirteen proteins were involved in mitochondrial function. TEM demonstrated that aPL bind to mitochondria in the syncytiotrophoblast and may cause mitochondrial swelling.

DISCUSSION

aPL disrupt mitochondria increasing the extrusion of SNAs with an altered proteome from the syncytiotrophoblast. These altered SNAs may trigger endothelial dysfunction and pre-eclampsia in these pregnancies.

Regulation of ERK-MAPK signaling in human epidermis

Background

The skin is largely comprised of keratinocytes within the interfollicular epidermis. Over approximately two weeks these cells differentiate and traverse the thickness of the skin. The stage of differentiation is therefore reflected in the positions of cells within the tissue, providing a convenient axis along which to study the signaling events that occur in situ during keratinocyte terminal differentiation, over this extended two-week timescale. The canonical ERK-MAPK signaling cascade (Raf-1, MEK-1/2 and ERK-1/2) has been implicated in controlling diverse cellular behaviors, including proliferation and differentiation. While the molecular interactions involved in signal transduction through this cascade have been well characterized in cell culture experiments, our understanding of how this sequence of events unfolds to determine cell fate within a homeostatic tissue environment has not been fully characterized.

Methods

We measured the abundance of total and phosphorylated ERK-MAPK signaling proteins within interfollicular keratinocytes in transverse cross-sections of human epidermis using immunofluorescence microscopy. To investigate these data we developed a mathematical model of the signaling cascade using a normalized-Hill differential equation formalism.

Results

These data show coordinated variation in the abundance of phosphorylated ERK-MAPK components across the epidermis. Statistical analysis of these data shows that associations between phosphorylated ERK-MAPK components which correspond to canonical molecular interactions are dependent upon spatial position within the epidermis. The model demonstrates that the spatial profile of activation for ERK-MAPK signaling components across the epidermis may be maintained in a cell-autonomous fashion by an underlying spatial gradient in calcium signaling.

Conclusions

Our data demonstrate an extended phospho-protein profile of ERK-MAPK signaling cascade components across the epidermis in situ, and statistical associations in these data indicate canonical ERK-MAPK interactions underlie this spatial profile of ERK-MAPK activation. Using mathematical modelling we have demonstrated that spatially varying calcium signaling components across the epidermis may be sufficient to maintain the spatial profile of ERK-MAPK signaling cascade components in a cell-autonomous manner. These findings may have significant implications for the wide range of cancer drugs which therapeutically target ERK-MAPK signaling components.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-015-0187-6) contains supplementary material, which is available to authorized users.

Transposon mutagenesis identifies genetic drivers of Braf(V600E) melanoma

Although nearly half of human melanomas harbor oncogenic BRAF(V600E) mutations, the genetic events that cooperate with these mutations to drive melanogenesis are still largely unknown. Here we show that Sleeping Beauty (SB) transposon-mediated mutagenesis drives melanoma progression in Braf(V600E) mutant mice and identify 1,232 recurrently mutated candidate cancer genes (CCGs) from 70 SB-driven melanomas. CCGs are enriched in Wnt, PI3K, MAPK and netrin signaling pathway components and are more highly connected to one another than predicted by chance, indicating that SB targets cooperative genetic networks in melanoma. Human orthologs of >500 CCGs are enriched for mutations in human melanoma or showed statistically significant clinical associations between RNA abundance and survival of patients with metastatic melanoma. We also functionally validate CEP350 as a new tumor-suppressor gene in human melanoma. SB mutagenesis has thus helped to catalog the cooperative molecular mechanisms driving BRAF(V600E) melanoma and discover new genes with potential clinical importance in human melanoma.

Cohesin modulates transcription of estrogen-responsive genes

The cohesin complex has essential roles in cell division, DNA damage repair and gene transcription. The transcriptional function of cohesin is thought to derive from its ability to connect distant regulatory elements with gene promoters. Genome-wide binding of cohesin in breast cancer cells frequently coincides with estrogen receptor alpha (ER), leading to the hypothesis that cohesin facilitates estrogen-dependent gene transcription. We found that cohesin modulates the expression of only a subset of genes in the ER transcription program, either activating or repressing transcription depending on the gene target. Estrogen-responsive genes most significantly influenced by cohesin were enriched in pathways associated with breast cancer progression such as PI3K and ErbB1. In MCF7 breast cancer cells, cohesin depletion enhanced transcription of TFF1 and TFF2, and was associated with increased ER binding and increased interaction between TFF1 and its distal enhancer situated within TMPRSS3. In contrast, cohesin depletion reduced c-MYC mRNA and was accompanied by reduced interaction between a distal enhancer of c-MYC and its promoters. Our data indicates that cohesin is not a universal facilitator of ER-induced transcription and can even restrict enhancer-promoter communication. We propose that cohesin modulates transcription of estrogen-dependent genes to achieve appropriate directionality and amplitude of expression.

Abstract A19: Characterizing the tumor stroma of B16 melanoma at different sites

The stromal cells of the tumor microenvironment play a critical role in tumor progression, and are therefore attractive targets for cancer therapy. We previously used the stromal targeting agent, 5,6-dimethylxanthenone-4-acetic acid (DMXAA), to demonstrate how tumor and stromal interactions can be manipulated to favor anti-tumor activity against human melanoma xenografts in immunodeficient mice, mediated by neutrophils and macrophages. An ideal stromal targeting agent would be effective against established metastases as well as primary tumors. In order to design such therapies, and to assess their potential, a better understanding of the tumor microenvironment at the various sites is required.

The aim of this study was to characterize the tumor stroma of established B16-F10 murine melanoma tumors at subcutaneous, pulmonary, and intracranial sites, in syngeneic C57BL/6 mice. Immunofluorescence-immunohistochemistry and flow cytometry were used to identify infiltrating leukocyte populations within B16-F10 tumors at these different sites. Furthermore, the antitumor activity of DMXAA against this model was investigated.

In subcutaneous tumors, monocytic myeloid cells (CD11b+ F4/80-/dim) and macrophages (CD11b+ F4/80+) were found predominantly around the tumor periphery. Lymphocyte infiltration was variable, with B-cells (CD19+ B220+) and T-cells (CD3+ CD4+ and CD3+ CD8+) rarely detected in some tumors, but scattered diffusely or formed dense lymphoid structures throughout necrotic regions in other tumors. Lung resident alveolar macrophages (CD11c+ CD11bdim F4/80+) appeared to be the most common leukocyte associated with metastatic foci, along with interstitial macrophages (CD11b+ F4/80+). In the presence of metastases, the total leukocyte abundance in the lungs was increased by nearly 70% (p < 0.05) compared to healthy lungs, owing largely to an increase in B-cells and T-cells scattered throughout the lung tissue; however these cells were not specifically associated with metastatic foci. In contrast, lymphocytes were scarcely detected in intracranial tumors or in brain tissue, where resident microglia (CD45dim CD11b+ CD11c+) were the most abundant. Microglia were found mainly at the invasive edge of intracranial tumors, while peripheral macrophages (CD45+ CD11b+ F4/80+) localised near vessels within the central core. Although DMXAA induced haemorrhagic necrosis in subcutaneous B16-F10 tumors, antitumor activity overall was modest, and no significant survival benefit was observed compared to untreated tumor-bearing mice.

Common at all sites was the presence of myeloid cells associated with tumor tissue, including resident and peripheral macrophages. DMXAA was recently shown to be an agonist for stimulator of interferon genes (STING) protein which triggers IFN-β signalling and innate immune gene expression in macrophages. DMXAA induces innate immunity against melanoma xenografts, where a similar presence of myeloid cells is observed before treatment. However, this activity was modest in aggressive B16-F10 melanomas. The precise reasons for this are unclear, but may involve immunosuppressive mechanisms. Our center is developing novel compounds to overcome immune-related resistance, such as inhibitors of indoleamine-2,3-dioxygenase, to be used with vaccines against melanoma. B16 metastatic melanoma is one challenging model that will be used to test these.

Citation Format: Kimiora Henare, Raymond Yung, P. Rod Dunbar, Cristin G. Print, Lai-Ming Ching. Characterizing the tumor stroma of B16 melanoma at different sites. [abstract]. In: Abstracts: AACR Special Conference on Cellular Heterogeneity in the Tumor Microenvironment; 2014 Feb 26-Mar 1; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2015;75(1 Suppl):Abstract nr A19. doi:10.1158/1538-7445.CHTME14-A19

NAIL, a software toolset for inferring, analyzing and visualizing regulatory networks

The wide variety of published approaches for the problem of regulatory network inference makes using multiple inference algorithms complex and time-consuming. Network Analysis and Inference Library (NAIL) is a set of software tools to simplify the range of computational activities involved in regulatory network inference. It uses a modular approach to connect different network inference algorithms to the same visualization and network-based analyses. NAIL is technology-independent and includes an interface layer to allow easy integration of components into other applications.

AVAILABILITY AND IMPLEMENTATION

NAIL is implemented in MATLAB, runs on Windows, Linux and OSX, and is available from SourceForge at https://sourceforge.net/projects/nailsystemsbiology/ for all researchers to use.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

In the secretory endometria of women, luminal epithelia exhibit gene and protein expressions that differ from those of glandular epithelia

OBJECTIVE

To characterize the transcriptome of luminal epithelia (LE) of fertile secretory endometria and compare the results with those from glandular epithelia (GE).

DESIGN

Endometrial samples were collected at 2 and 7 days after initial blood LH surge in separate menstrual cycles. LE were obtained with the use of laser microdissection. mRNA was amplified with the use of linear polymerase chain reaction and hybridized to Agilent 4×44 microarrays. Gene analysis was used to identify differentially expressed mRNAs. Immunohistochemistry was used to assess nine proteins.

SETTING

One IVF clinic.

PATIENT(S)

Seven Caucasian fertile cycling women.

INTERVENTION(S)

None.

MAIN OUTCOME MEASURE(S)

Cycle dating with the use of blood endocrinologic markers, microarrays of laser-microdissected LE, immunohistochemical analysis.

RESULT(S)

One hundred sixty-one (of 401) differentially expressed mRNAs in LE were identified from the metabolism pathway. Increased selective protein expression in LE at 7 days after initial LH surge was observed. LE mRNA expression was the converse of that in GE. The two cell types each had a different significant biologic pathway identified.

CONCLUSION(S)

Our results introduce a new concept that LE differentially expressed mRNAs are in the converse direction to that of GE, indicating different biologic processes despite the GE being continuous with the luminal monolayer. This probable distinction of biologic roles has not been noted previously. Further investigations must take cognizance of this observation.

The flavoprotein FOXRED2 reductively activates nitro-chloromethylbenzindolines and other hypoxia-targeting prodrugs

The nitro-chloromethylbenzindoline prodrug SN29428 has been rationally designed to target tumour hypoxia. SN29428 is metabolised to a DNA minor groove alkylator via oxygen-sensitive reductive activation initiated by unknown one-electron reductases. The present study sought to identify reductases capable of activating SN29428 in tumours. Expression of candidate reductases in cell lines was modulated using forced expression and, for P450 (cytochrome) oxidoreductase (POR), by zinc finger nuclease-mediated gene knockout. Affymetrix microarray mRNA expression of flavoreductases was correlated with SN29428 activation in a panel of 23 cancer cell lines. Reductive activation and cytotoxicity of prodrugs were measured using mass spectrometry and antiproliferative assays, respectively. SN29428 activation under hypoxia was strongly attenuated by the pan-flavoprotein inhibitor diphenyliodonium, but less so by knockout of POR suggesting other flavoreductases contribute. Forced expression of 5-methyltetrahydrofolate-homocysteine methyltransferase reductase (MTRR), as well as POR, increased activation of SN29428 in hypoxic HCT 116 cells. SN29428 activation strongly correlated with expression of POR and also FAD-dependent oxidoreductase domain containing 2 (FOXRED2), in cancer cell lines. This association persisted after removing the effect of POR enzyme activity using first-order partial correlation. Forced expression of FOXRED2 increased SN29428 activation and cytotoxicity in hypoxic HEK293 cells and also increased activation of hypoxia-targeted prodrugs PR-104A, tirapazamine and SN30000, and increased cytotoxicity of the clinical-stage prodrug TH-302. Thus this study has identified three flavoreductases capable of enzymatically activating SN29428, one of which (FOXRED2) has not previously been implicated in xenobiotic metabolism. These results will inform future development of biomarkers predictive of SN29428 sensitivity.

Integration of steady-state and temporal gene expression data for the inference of gene regulatory networks

We develop a new regression algorithm, cMIKANA, for inference of gene regulatory networks from combinations of steady-state and time-series gene expression data. Using simulated gene expression datasets to assess the accuracy of reconstructing gene regulatory networks, we show that steady-state and time-series data sets can successfully be combined to identify gene regulatory interactions using the new algorithm. Inferring gene networks from combined data sets was found to be advantageous when using noisy measurements collected with either lower sampling rates or a limited number of experimental replicates. We illustrate our method by applying it to a microarray gene expression dataset from human umbilical vein endothelial cells (HUVECs) which combines time series data from treatment with growth factor TNF and steady state data from siRNA knockdown treatments. Our results suggest that the combination of steady-state and time-series datasets may provide better prediction of RNA-to-RNA interactions, and may also reveal biological features that cannot be identified from dynamic or steady state information alone. Finally, we consider the experimental design of genomics experiments for gene regulatory network inference and show that network inference can be improved by incorporating steady-state measurements with time-series data.

BCL-2 hypermethylation is a potential biomarker of sensitivity to antimitotic chemotherapy in endocrine-resistant breast cancer

Overexpression of the antiapoptotic factor BCL-2 is a frequent feature of malignant disease and is commonly associated with poor prognosis and resistance to conventional chemotherapy. In breast cancer, however, high BCL-2 expression is associated with favorable prognosis, estrogen receptor (ER) positivity, and low tumor grade, whereas low expression is included in several molecular signatures associated with resistance to endocrine therapy. In the present study, we correlate BCL-2 expression and DNA methylation profiles in human breast cancer and in multiple cell models of acquired endocrine resistance to determine whether BCL-2 hypermethylation could provide a useful biomarker of response to cytotoxic therapy. In human disease, diminished expression of BCL-2 was associated with hypermethylation of the second exon, in a region that overlapped a CpG island and an ER-binding site. Hypermethylation of this region, which occurred in 10% of primary tumors, provided a stronger predictor of patient survival (P = 0.019) when compared with gene expression (n = 522). In multiple cell models of acquired endocrine resistance, BCL-2 expression was significantly reduced in parallel with increased DNA methylation of the exon 2 region. The reduction of BCL-2 expression in endocrine-resistant cells lowered their apoptotic threshold to antimitotic agents: nocodazole, paclitaxel, and the PLK1 inhibitor BI2536. This phenomenon could be reversed with ectopic expression of BCL-2, and rescued with the BCL-2 inhibitor ABT-737. Collectively, these data imply that BCL-2 hypermethylation provides a robust biomarker of response to current and next-generation cytotoxic agents in endocrine-resistant breast cancer, which may prove beneficial in directing therapeutic strategy for patients with nonresectable, metastatic disease.

MelanomaDB: A Web Tool for Integrative Analysis of Melanoma Genomic Information to Identify Disease-Associated Molecular Pathways

Despite on-going research, metastatic melanoma survival rates remain low and treatment options are limited. Researchers can now access a rapidly growing amount of molecular and clinical information about melanoma. This information is becoming difficult to assemble and interpret due to its dispersed nature, yet as it grows it becomes increasingly valuable for understanding melanoma. Integration of this information into a comprehensive resource to aid rational experimental design and patient stratification is needed. As an initial step in this direction, we have assembled a web-accessible melanoma database, MelanomaDB, which incorporates clinical and molecular data from publically available sources, which will be regularly updated as new information becomes available. This database allows complex links to be drawn between many different aspects of melanoma biology: genetic changes (e.g., mutations) in individual melanomas revealed by DNA sequencing, associations between gene expression and patient survival, data concerning drug targets, biomarkers, druggability, and clinical trials, as well as our own statistical analysis of relationships between molecular pathways and clinical parameters that have been produced using these data sets. The database is freely available at http://genesetdb.auckland.ac.nz/melanomadb/about.html. A subset of the information in the database can also be accessed through a freely available web application in the Illumina genomic cloud computing platform BaseSpace at http://www.biomatters.com/apps/melanoma-profiler-for-research. The MelanomaDB database illustrates dysregulation of specific signaling pathways across 310 exome-sequenced melanomas and in individual tumors and identifies the distribution of somatic variants in melanoma. We suggest that MelanomaDB can provide a context in which to interpret the tumor molecular profiles of individual melanoma patients relative to biological information and available drug therapies.

Paternal obesity initiates metabolic disturbances in two generations of mice with incomplete penetrance to the F2 generation and alters the transcriptional profile of testis and sperm microRNA content

Obesity is highly prevalent, and its incidence is increasing. The previous study showing a major effect of paternal obesity on metabolic health of offspring is confounded by comorbidity with diabetes. Therefore, we investigated the effect of diet-induced paternal obesity, in the absence of diabetes, on the metabolic health of two resultant generations and the molecular profiles of the testes and sperm. Founder (F0) male C57BL6 mice were fed either a high-fat diet (HFD) or a control diet (CD); n = 10/diet for a period of 10 wk. Testis expression of mRNA/microRNAs was analyzed by microarray and qPCR and sperm microRNA abundance by qPCR. Two subsequent generations were generated by mating F0 and then F1 mice to CD mice, and their metabolic health was investigated. All mice, other than F0 males, were maintained on a CD. HFD feeding induced paternal obesity with a 21% increase in adiposity, but not overt diabetes, and initiated intergenerational transmission of obesity and insulin resistance in two generations of offspring. This distinct phenotypic constellation is either partially or fully transmitted to both female and male F1 offspring and further transmitted through both parental lineages to the F2 generation, with a heightened effect on female F1 offspring (+67% in adiposity) and their F2 sons (+24% in adiposity). Founder male obesity altered the testes expression of 414 mRNAs by microarray and 11 microRNAs by qPCR, concomitant with alterations in sperm microRNA content and a 25% reduction in global methylation of germ cell DNA. Diet-induced paternal obesity modulates sperm microRNA content and germ cell methylation status, which are potential signals that program offspring health and initiate the transmission of obesity and impaired metabolic health to future generations. This study implicates paternal obesity in the transgenerational amplification of obesity and type 2 diabetes in humans.

Zinc finger nuclease mediated knockout of ADP-dependent glucokinase in cancer cell lines: effects on cell survival and mitochondrial oxidative metabolism

Zinc finger nucleases (ZFN) are powerful tools for editing genes in cells. Here we use ZFNs to interrogate the biological function of ADPGK, which encodes an ADP-dependent glucokinase (ADPGK), in human tumour cell lines. The hypothesis we tested is that ADPGK utilises ADP to phosphorylate glucose under conditions where ATP becomes limiting, such as hypoxia. We characterised two ZFN knockout clones in each of two lines (H460 and HCT116). All four clones had frameshift mutations in all alleles at the target site in exon 1 of ADPGK, and were ADPGK-null by immunoblotting. ADPGK knockout had little or no effect on cell proliferation, but compromised the ability of H460 cells to survive siRNA silencing of hexokinase-2 under oxic conditions, with clonogenic survival falling from 21±3% for the parental line to 6.4±0.8% (p = 0.002) and 4.3±0.8% (p = 0.001) for the two knockouts. A similar increased sensitivity to clonogenic cell killing was observed under anoxia. No such changes were found when ADPGK was knocked out in HCT116 cells, for which the parental line was less sensitive than H460 to anoxia and to hexokinase-2 silencing. While knockout of ADPGK in HCT116 cells caused few changes in global gene expression, knockout of ADPGK in H460 cells caused notable up-regulation of mRNAs encoding cell adhesion proteins. Surprisingly, we could discern no consistent effect on glycolysis as measured by glucose consumption or lactate formation under anoxia, or extracellular acidification rate (Seahorse XF analyser) under oxic conditions in a variety of media. However, oxygen consumption rates were generally lower in the ADPGK knockouts, in some cases markedly so. Collectively, the results demonstrate that ADPGK can contribute to tumour cell survival under conditions of high glycolytic dependence, but the phenotype resulting from knockout of ADPGK is cell line dependent and appears to be unrelated to priming of glycolysis in these lines.

Abstract 1116: Sleeping Beauty mutagenesis identifies genes and pathways that cooperate with BrafV600E in melanoma initiation and progression

Malignant melanoma is a devastating disease the prevalence of which is on the rise. In humans, somatic mutation of BRAF, to encode the BRAFV600E oncoprotein, is an early predisposing event in melanoma initiation, but is insufficient to drive malignant melanomagenesis in the absence of additional cooperating events. We describe a genome-wide, forward genetic screen in mice restricted to pigment producing melanocytes as an unbiased approach to identify new candidate melanoma progression genes. By combining Sleeping Beauty (SB) with a Cre-activated BRafCA allele we observed frank melanomas with average latency of 36 weeks. By contrast, melanocyte specific expression of BRAFV600E elicited benign melanocytic nevi that failed to progress to malignant melanoma by one year of age. These findings confirm that SB mobilization altered genes that cooperate with BRAFV600E to promote melanoma progression. A discovery set of 1,200 statistically-defined candidate melanoma genes was identified from 70 SB|BRAFV600E tumors. Importantly, our screen identified seven SB mutated genes with human orthologs that are well established human melanoma drivers including: Cdkn2a, Pten, Gnaq, Mitf, Nf1, Rac1 and Ppp6c. Many biological signaling pathways and cellular processes were enriched for candidate melanoma genes altered in SB|BRAFV600E melanomas, including Wnt/Beta-catenin, TGFβ, PI3’-kinase and MAPK signaling pathways and biological processes regulating ubiquitin mediated proteolysis, tight junctions, cell cycle and axonal guidance. Integrated comparative oncogenomic analysis with human melanoma data sets provided a data-driven method to filter and prioritize candidate melanoma genes. The genes identified using the SB|BRAFV600E model may represent new human melanoma candidate cancer genes with potential clinical relevance to the progression and maintenance of human BRAF mutated melanoma.

Citation Format: Michael B. Mann, Michael A. Black, Jerrold M. Ward, Marcus W. Bosenberg, Martin McMahon, Cristin G. Print, David J. Adams, Neal G. Copeland, Nancy A. Jenkins. Sleeping Beauty mutagenesis identifies genes and pathways that cooperate with BrafV600E in melanoma initiation and progression. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 1116. doi:10.1158/1538-7445.AM2013-1116

Biclustering reveals breast cancer tumour subgroups with common clinical features and improves prediction of disease recurrence

BACKGROUND

Many studies have revealed correlations between breast tumour phenotypes, variations in gene expression, and patient survival outcomes. The molecular heterogeneity between breast tumours revealed by these studies has allowed prediction of prognosis and has underpinned stratified therapy, where groups of patients with particular tumour types receive specific treatments. The molecular tests used to predict prognosis and stratify treatment usually utilise fixed sets of genomic biomarkers, with the same biomarker sets being used to test all patients. In this paper we suggest that instead of fixed sets of genomic biomarkers, it may be more effective to use a stratified biomarker approach, where optimal biomarker sets are automatically chosen for particular patient groups, analogous to the choice of optimal treatments for groups of similar patients in stratified therapy. We illustrate the effectiveness of a biclustering approach to select optimal gene sets for determining the prognosis of specific strata of patients, based on potentially overlapping, non-discrete molecular characteristics of tumours.

RESULTS

Biclustering identified tightly co-expressed gene sets in the tumours of restricted subgroups of breast cancer patients. The co-expressed genes in these biclusters were significantly enriched for particular biological annotations and gene regulatory modules associated with breast cancer biology. Tumours identified within the same bicluster were more likely to present with similar clinical features. Bicluster membership combined with clinical information could predict patient prognosis in conditional inference tree and ridge regression class prediction models.

CONCLUSIONS

The increasing clinical use of genomic profiling demands identification of more effective methods to segregate patients into prognostic and treatment groups. We have shown that biclustering can be used to select optimal gene sets for determining the prognosis of specific strata of patients.

Targeting p90 ribosomal S6 kinase eliminates tumor-initiating cells by inactivating Y-box binding protein-1 in triple-negative breast cancers

Y-box binding protein-1 (YB-1) is the first reported oncogenic transcription factor to induce the tumor-initiating cell (TIC) surface marker CD44 in triple-negative breast cancer (TNBC) cells. In order for CD44 to be induced, YB-1 must be phosphorylated at S102 by p90 ribosomal S6 kinase (RSK). We therefore questioned whether RSK might be a tractable molecular target to eliminate TICs. In support of this idea, injection of MDA-MB-231 cells expressing Flag-YB-1 into mice increased tumor growth as well as enhanced CD44 expression. Despite enrichment for TICs, these cells were sensitive to RSK inhibition when treated ex vivo with BI-D1870. Targeting RSK2 with small interfering RNA (siRNA) or small molecule RSK kinase inhibitors (SL0101 and BI-D1870) blocked TNBC monolayer cell growth by ∼100%. In a diverse panel of breast tumor cell line models RSK2 siRNA predominantly targeted models of TNBC. RSK2 inhibition decreased CD44 promoter activity, CD44 mRNA, protein expression, and mammosphere formation. CD44(+) cells had higher P-RSK(S221/227) , P-YB-1(S102) , and mitotic activity relative to CD44(-) cells. Importantly, RSK2 inhibition specifically suppressed the growth of TICs and triggered cell death. Moreover, silencing RSK2 delayed tumor initiation in mice. In patients, RSK2 mRNA was associated with poor disease-free survival in a cohort of 244 women with breast cancer that had not received adjuvant treatment, and its expression was highest in the basal-like breast cancer subtype. Taking this further, we report that P-RSK(S221/227) is present in primary TNBCs and correlates with P-YB-1(S102) as well as CD44. In conclusion, RSK2 inhibition provides a novel therapeutic avenue for TNBC and holds the promise of eliminating TICs.

The transcriptional targets of mutant FOXL2 in granulosa cell tumours

Background

Despite their distinct biology, granulosa cell tumours (GCTs) are treated the same as other ovarian tumours. Intriguingly, a recurring somatic mutation in the transcription factor Forkhead Box L2 (FOXL2) 402C>G has been found in nearly all GCTs examined. This investigation aims to identify the pathogenicity of mutant FOXL2 by studying its altered transcriptional targets.

Methods

The expression of mutant FOXL2 was reduced in the GCT cell line KGN, and wildtype and mutant FOXL2 were overexpressed in the GCT cell line COV434. Total RNA was hybridised to Affymetrix U133 Plus 2 microarrays. Comparisons were made between the transcriptomes of control cells and cells altered by FOXL2 knockdown and overexpression, to detect potential transcriptional targets of mutant FOXL2.

Results

The overexpression of wildtype and mutant FOXL2 in COV434, and the silencing of mutant FOXL2 expression in KGN, has shown that mutant FOXL2 is able to differentially regulate the expression of many genes, including two well known FOXL2 targets, StAR and CYP19A. We have shown that many of the genes regulated by mutant FOXL2 are clustered into functional annotations of cell death, proliferation, and tumourigenesis. Furthermore, TGF-β signalling was found to be enriched when using the gene annotation tools GATHER and GeneSetDB. This enrichment was still significant after performing a robust permutation analysis.

Conclusion

Given that many of the transcriptional targets of mutant FOXL2 are known TGF-β signalling genes, we suggest that deregulation of this key antiproliferative pathway is one way mutant FOXL2 contributes to the pathogenesis of adult-type GCTs. We believe this pathway should be a target for future therapeutic interventions, if outcomes for women with GCTs are to improve.