Circulating tumor DNA is a sensitive marker for routine monitoring of treatment response in advanced colorectal cancer.

Accurate assessment of chemotherapy response provides the means to terminate ineffective treatment, trial alternative drug regimens or schedules and reduce dose to minimize toxicity. Here, we have compared circulating tumor DNA (ctDNA) with carcinoembryonic antigen (CEA) for the cycle by cycle assessment of chemotherapy response in 30 patients with metastatic colorectal cancer. CtDNA (quantified using individualized digital droplet PCR (ddPCR) assays) and CEA levels were determined immediately prior to each chemotherapy cycle over time periods ranging from 42-548 days (average of 10 time points/patient). Twenty-nine/thirty (97%) patients had detectable ctDNA compared with 83% whose tumors were CEA-positive (>5 ng/ml) during the monitoring course. Over the course of treatment, 20 disease progression events were detected by computed tomography; ctDNA predicted significantly more of these events than CEA (16 (80%) versus 6 (30%), respectively; P-value = 0.004). When progression was detected by both ctDNA and CEA, the rise in ctDNA occurred significantly earlier than CEA (P-value = 0.046). Partial responses to chemotherapy were also detected more frequently by ctDNA, although this was not significant (P-value = 0.07). In addition, another 28 colorectal cancer patients who underwent potentially curative surgery and showed no evidence of residual disease were monitored with ctDNA for up to 2 years. Clinical relapse was observed in 6/28 (21%) patients. Four out of 6 of these patients showed a significant increase in ctDNA at or prior to relapse. Overall, ctDNA analyses were able to be performed in a clinically relevant timeline and were a more sensitive and responsive measure of tumor burden than CEA.

Systems immunology reveals a linked IgG3-C4 response in patients with acute rheumatic fever.

Acute rheumatic fever (ARF) and chronic rheumatic heart disease (RHD) are autoimmune sequelae of a Group A streptococcal infection with significant global mortality and poorly understood pathogenesis. Immunoglobulin and complement deposition were observed in ARF/RHD valve tissue over 50 years ago, yet contemporary investigations have been lacking. This study applied systems immunology to investigate the relationships between the complement system and immunoglobulin in ARF. Patients were stratified by C-reactive protein (CRP) concentration into high (≥10 µg mL-1 ) and low (<10 µg mL-1 ) groups to distinguish those with clinically significant inflammatory processes from those with abating inflammation. The circulating concentrations of 17 complement factors and six immunoglobulin isotypes and subclasses were measured in ARF patients and highly matched healthy controls using multiplex bead-based immunoassays. An integrative statistical approach combining feature selection and principal component analysis revealed a linked IgG3-C4 response in ARF patients with high CRP that was absent in controls. Strikingly, both IgG3 and C4 were elevated above clinical reference ranges, suggesting these features are a marker of ARF-associated inflammation. Humoral immunity in response to M protein, an antigen implicated in ARF pathogenesis, was completely polarized to IgG3 in the patient group. Furthermore, the anti-M-protein IgG3 response was correlated with circulating IgG3 concentration, highlighting a potential role for this potent immunoglobulin subclass in disease. In conclusion, a linked IgG3-C4 response appears important in the initial, inflammatory stage of ARF and may have immediate utility as a clinical biomarker given the lack of specific diagnostic tests currently available.

Mapping a route to Indigenous engagement in cancer genomic research.

Precision oncology guided by genomic research has an increasingly important role in the care of people with cancer. However, substantial inequities remain in cancer outcomes of Indigenous peoples, including Indigenous Maori in Aotearoa New Zealand (New Zealand). These inequities will be perpetuated unless deliberate steps are taken to include Indigenous peoples in all parts of cancer research-as research participants, in research leadership, and in research governance. This approach is especially important when there have been historical breaches of trust that have discouraged their participation in health research. This Personal View describes a precision oncology research roadmap for neuroendocrine tumour research, which seeks to reflect the values of New Zealand's Indigenous Maori people. This roadmap includes facilitating ongoing dialogue, Maori leadership, reciprocity, agreed kawa (guiding principles), tikanga (cultural protocols), and honest monitoring of what is and what is not being achieved. We challenge cancer researchers worldwide to generate locally appropriate roadmaps that honestly assess their practices to benefit Indigenous people internationally.

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 ρ 0 (ρ 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 ρ 0 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.

Does the endometrial gene expression of fertile women vary within and between cycles?

STUDY QUESTION: Does gene expression of putative endometrial implantation markers vary in expression between menstrual cycles? SUMMARY ANSWER: In fertile women the expression of certain genes exhibits a pattern of stable regulation.which is not affected even when sampled twice in one cycle. WHAT IS KNOWN ALREADY: Successful implantation occurs in a minority of IVF embryo transfers. In contrast to knowledge regarding the ovulatory process, there is a sparse understanding of endometrial genes critical to implantation. This lack of knowledge hinders progress in this field. STUDY DESIGN, SIZE, DURATION: Endometrial pipelle samples were collected based on blood endocrinological markers at 2 and 7 days post initial LH surge. Five samples were collected over four cycles where the interval between collections ranged from sequential months to three years. PARTICIPANTS/MATERIALS, SETTING, METHODS: Six fertile women attending an IVF clinic for male factor infertility, had samples collected. Global gene expression profiles were obtained from laser-microdissected, endometrial glands and stroma. Nineteen potential proliferation, cytokine and adhesion markers based on previous validated reports were studied. MAIN RESULTS AND THE ROLE OF CHANCE: There was a significant modification between LH+2 and LH+7 of expression for 23 genes-11 in 8 in glands and stroma, 4 in stroma only and 3 in glands only suggesting stable, controlled regulation. Nevertheless, genes exhibited individual characteristics, e.g MKI67 exhibited lower expression at LH+7 than LH+2 and CCL4 higher, whereas TRO expressed limited difference in both cell types. Stability between cycles was demonstrated for gene expression at both LH+2-more than 60% of genes had <25% variation and at LH+7-60% had <30% variation. Further, effects of prior collection of an LH+2 sample on gene expression at LH+7 were not detected. The range of mRNA expression suggested that a clinical/diagnostic sample at LH+2 and LH+7 is likely to be a better index of endometrial function than a single sample. The possibility of redundancy suggests a panel would be more informative than a single marker. LARGE SCALE DATA: Raw and normalized microarray data have been deposited with the EMBL's European Genome-Phenome Archive for collaborative analysis, reference ega-box-815 (Lappalainen I, Almeida-King J, Kumanduri V, Senf A, Spalding JD, Ur-Rehman S, Saunders G, Kandasamy J, Caccamo M, Leinonen R et al. The European Genome-phenome Archive of human data consented for biomedical research. Nat Genet 2015;47:692-695.) [https://www.ebi.ac.uk/ega/home]. LIMITATIONS, REASONS FOR CAUTION: This type of research has difficulties of recruitment of fertile women for multiple blood testing and repeat endometrial biopsies. Therefore, these data had decreased statistical power due to the overall participant numbers. However, the inclusion of four cycles for each participant permitted the aim of obtaining information on intercycle and intracycle variability to be achieved. WIDER IMPLICATIONS OF THE FINDINGS: Our results support the feasibility of a clinical means of identification of a functional receptive endometrium. The robustness of data from individual women suggests that samples from one cycle can generally be applied to subsequent cycles. STUDY FUNDING/COMPETING INTEREST(S): Funding was granted from the Tertiary Education Commission of New Zealand, Contract I.D.:UOOX06007. There are no competing interests.

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.

Massive parallel sequencing of solid tumours - challenges and opportunities for pathologists.

The role of pathologists is to provide diagnostic, prognostic and predictive data to enable clinical colleagues to manage patients optimally. Current histo/anatomical pathology is predominantly morphology-based, with the addition of biomarkers, applied largely through immunohistochemistry, fluorescence in-situ hybridization (FISH) and a limited range of polymerase chain reaction (PCR)-based molecular tests. The desire to apply genomics to the clinical care of patients has been facilitated by the human genome project and subsequently by high-throughput technologies known collectively as massive parallel sequencing (MPS, also referred to as next-generation sequencing, NGS). The use of MPS to identify mutations/variants and tissue RNA expression profiles for diagnosis, prognostication and targeted therapy stratification is now a reality in many clinical specialities. If histopathologists are considered experts in solid tumour pathology, MPS potentially falls within their scope; however, it challenges our predominant morphology-based paradigm. This review summarizes and comments on the current and future state of play of MPS for the practising histopathologist. It will focus on somatic mutations in solid tumours and will challenge histopathologists to take further leadership roles in this area.

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.

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.

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

UNLABELLED: 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.

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.

YB-1: oncoprotein, prognostic marker and therapeutic target?

Hanahan and Weinberg have proposed the 'hallmarks of cancer' to cover the biological changes required for the development and persistence of tumours [Hanahan and Weinberg (2011) Cell 144, 646-674]. We have noted that many of these cancer hallmarks are facilitated by the multifunctional protein YB-1 (Y-box-binding protein 1). In the present review we evaluate the literature and show how YB-1 modulates/regulates cellular signalling pathways within each of these hallmarks. For example, we describe how YB-1 regulates multiple proliferation pathways, overrides cell-cycle check points, promotes replicative immortality and genomic instability, may regulate angiogenesis, has a role in invasion and metastasis, and promotes inflammation. We also argue that there is strong and sufficient evidence to suggest that YB-1 is an excellent molecular marker of cancer progression that could be used in the clinic, and that YB-1 could be a useful target for cancer therapy.

Gene network inference and visualization tools for biologists: application to new human transcriptome datasets.

Gene regulatory networks inferred from RNA abundance data have generated significant interest, but despite this, gene network approaches are used infrequently and often require input from bioinformaticians. We have assembled a suite of tools for analysing regulatory networks, and we illustrate their use with microarray datasets generated in human endothelial cells. We infer a range of regulatory networks, and based on this analysis discuss the strengths and limitations of network inference from RNA abundance data. We welcome contact from researchers interested in using our inference and visualization tools to answer biological questions.

A validated restriction enzyme ddPCR cg05575921 (AHRR) assay to accurately assess smoking exposure.

BACKGROUND & METHODS: In this study, a novel restriction enzyme (RE) digestion-based droplet digital polymerase chain reaction (ddPCR) assay was designed for cg005575921 within the AHRR gene body and compared with matching results obtained by bisulfite conversion (BIS) ddPCR and Illumina DNA methylation array. RESULTS: The RE ddPCR cg05575921 assay appeared concordant with BIS ddPCR (r2 = 0.94, P < 0.0001) and, when compared with the Illumina array, had significantly better smoking status classification performance for current versus never smoked (AUC 0.96 versus 0.93, P < 0.04) and current versus ex-smoker (AUC 0.88 versus 0.83, P < 0.04) comparisons. CONCLUSIONS: The RE ddPCR cg05575921 assay accurately predicts smoking status and could be a useful component of 'precision-medicine' chronic disease risk screening tools.

Combining TP53 mutation and isoform has the potential to improve clinical practice.

The clinical importance of assessing and combining data on TP53 mutations and isoforms is discussed in this article. It gives a succinct overview of the structural makeup and key biological roles of the isoforms. It then provides a comprehensive summary of the roles that p53 isoforms play in cancer development, therapy response and resistance. The review provides a summary of studies demonstrating the role of p53 isoforms as potential prognostic indicators. It further provides evidence on how the presence of TP53 mutations may affect one or more of these activities and the association of p53 isoforms with clinicopathological data in various tumour types. The review gives insight into the present diagnostic hurdles for identifying TP53 isoforms and makes recommendations to improve their evaluation. In conclusion, this review offers suggestions for enhancing the identification and integration of TP53 isoforms in conjunction with mutation data within the clinical context.

Biallelic ATP2B1 variants as a likely cause of a novel neurodevelopmental malformation syndrome with primary hypoparathyroidism.

ATP2B1 encodes plasma membrane calcium-transporting-ATPase1 and plays an essential role in maintaining intracellular calcium homeostasis that regulates diverse signaling pathways. Heterozygous de novo missense and truncating ATP2B1 variants are associated with a neurodevelopmental phenotype of variable expressivity. We describe a proband with distinctive craniofacial gestalt, Pierre-Robin sequence, neurodevelopmental and growth deficit, periventricular heterotopia, brachymesophalangy, cutaneous syndactyly, and persistent hypocalcemia from primary hypoparathyroidism. Proband-parent trio exome sequencing identified compound heterozygous ATP2B1 variants: a maternally inherited splice-site (c.3060+2 T > G) and paternally inherited missense c.2938 G > T; p.(Val980Leu). Reverse-transcription-PCR on the proband's fibroblast-derived mRNA showed aberrantly spliced ATP2B1 transcripts targeted for nonsense-mediated decay. All correctly-spliced ATP2B1 mRNA encoding p.(Val980Leu) functionally causes decreased cellular Ca2+ extrusion. Immunoblotting showed reduced fibroblast ATP2B1. We conclude that biallelic ATP2B1 variants are the likely cause of the proband's phenotype, strengthening the association of ATP2B1 as a neurodevelopmental gene and expanding the phenotypic characterization of a biallelic loss-of-function genotype.

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.

Vasohibin-1 is identified as a master-regulator of endothelial cell apoptosis using gene network analysis.

BACKGROUND: Apoptosis is a critical process in endothelial cell (EC) biology and pathology, which has been extensively studied at protein level. Numerous gene expression studies of EC apoptosis have also been performed, however few attempts have been made to use gene expression data to identify the molecular relationships and master regulators that underlie EC apoptosis. Therefore, we sought to understand these relationships by generating a Bayesian gene regulatory network (GRN) model. RESULTS: ECs were induced to undergo apoptosis using serum withdrawal and followed over a time course in triplicate, using microarrays. When generating the GRN, this EC time course data was supplemented by a library of microarray data from EC treated with siRNAs targeting over 350 signalling molecules.The GRN model proposed Vasohibin-1 (VASH1) as one of the candidate master-regulators of EC apoptosis with numerous downstream mRNAs. To evaluate the role played by VASH1 in EC, we used siRNA to reduce the expression of VASH1. Of 10 mRNAs downstream of VASH1 in the GRN that were examined, 7 were significantly up- or down-regulated in the direction predicted by the GRN.Further supporting an important biological role of VASH1 in EC, targeted reduction of VASH1 mRNA abundance conferred resistance to serum withdrawal-induced EC death. CONCLUSION: We have utilised Bayesian GRN modelling to identify a novel candidate master regulator of EC apoptosis. This study demonstrates how GRN technology can complement traditional methods to hypothesise the regulatory relationships that underlie important biological processes.

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.

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.