MAC-Seq: Coupling Low-Cost, High-Throughput RNA-Seq with Image-Based Phenotypic Screening in 2D and 3D Cell Models.

Transcriptomic profiling has fundamentally influenced our understanding of cancer pathophysiology and response to therapeutic intervention and has become a relatively routine approach. However, standard protocols are usually low-throughput, single-plex assays and costs are still quite prohibitive. With the evolving complexity of in vitro cell model systems, there is a need for resource-efficient high-throughput approaches that can support detailed time-course analytics, accommodate limited sample availability, and provide the capacity to correlate phenotype to genotype at scale. MAC-seq (multiplexed analysis of cells) is a low-cost, ultrahigh-throughput RNA-seq workflow in plate format to measure cell perturbations and is compatible with high-throughput imaging. Here we describe the steps to perform MAC-seq in 384-well format and apply it to 2D and 3D cell cultures. On average, our experimental conditions identified over ten thousand expressed genes per well when sequenced to a depth of one million reads. We discuss technical aspects, make suggestions on experimental design, and document critical operational procedures. Our protocol highlights the potential to couple MAC-seq with high-throughput screening applications including cell phenotyping using high-content cell imaging.

Tumor-Infiltrating Neutrophils after Neoadjuvant Therapy are Associated with Poor Prognosis in Esophageal Cancer.

BACKGROUND: In esophageal cancer (EC), there is a paucity of knowledge regarding the interplay between the tumor immune microenvironment and response to neoadjuvant treatment and, therefore, which factors may influence outcomes. Thus, our goal was to investigate the changes in the immune microenvironment with neoadjuvant treatment in EC by assessing the expression of immune related genes and their association with prognosis. METHODS: We examined the transcriptome of paired pre- and post-neoadjuvant treated EC specimens. Based on these findings, we validated the presence of tumor-infiltrating neutrophils using CD15+ immunohistochemistry in a discovery cohort of patients with residual pathologic disease. We developed a nomogram as a predictor of progression-free survival (PFS) incorporating the variables CD15+ cell count, tumor regression grade, and tumor grade. RESULTS: After neoadjuvant treatment, there was an increase in genes related to myeloid cell differentiation and a poor prognosis associated with high neutrophil (CD15+) counts. Our nomogram incorporating CD15+ cell count was predictive of PFS with a C-index of 0.80 (95% confidence interval [CI] 0.68-0.9) and a concordance probability estimate (CPE) of 0.77 (95% CI 0.69-0.86), which indicates high prognostic ability. The C-index and CPE of the validation cohort were 0.81 (95% CI 0.69-0.91) and 0.78 (95% CI 0.7-0.86), respectively. CONCLUSIONS: Our nomogram incorporating CD15+ cell count can potentially be used to identify patients at high risk of recurrent disease and thus stratify patients who will benefit most from adjuvant treatment.

Epigenetic modulators of B cell fate identified through coupled phenotype-transcriptome analysis.

High-throughput methodologies are the cornerstone of screening approaches to identify novel compounds that regulate immune cell function. To identify novel targeted therapeutics to treat immune disorders and haematological malignancies, there is a need to integrate functional cellular information with the molecular mechanisms that regulate changes in immune cell phenotype. We facilitate this goal by combining quantitative methods for dissecting complex simultaneous cell phenotypic effects with genomic analysis. This combination strategy we term Multiplexed Analysis of Cells sequencing (MAC-seq), a modified version of Digital RNA with perturbation of Genes (DRUGseq). We applied MAC-seq to screen compounds that target the epigenetic machinery of B cells and assess altered humoral immunity by measuring changes in proliferation, survival, differentiation and transcription. This approach revealed that polycomb repressive complex 2 (PRC2) inhibitors promote antibody secreting cell (ASC) differentiation in both murine and human B cells in vitro. This is further validated using T cell-dependent immunization in mice. Functional dissection of downstream effectors of PRC2 using arrayed CRISPR screening uncovered novel regulators of B cell differentiation, including Mybl1, Myof, Gas7 and Atoh8. Together, our findings demonstrate that integrated phenotype-transcriptome analyses can be effectively combined with drug screening approaches to uncover the molecular circuitry that drives lymphocyte fate decisions.

Epigenetic Activation of Plasmacytoid DCs Drives IFNAR-Dependent Therapeutic Differentiation of AML.

Pharmacologic inhibition of epigenetic enzymes can have therapeutic benefit against hematologic malignancies. In addition to affecting tumor cell growth and proliferation, these epigenetic agents may induce antitumor immunity. Here, we discovered a novel immunoregulatory mechanism through inhibition of histone deacetylases (HDAC). In models of acute myeloid leukemia (AML), leukemia cell differentiation and therapeutic benefit mediated by the HDAC inhibitor (HDACi) panobinostat required activation of the type I interferon (IFN) pathway. Plasmacytoid dendritic cells (pDC) produced type I IFN after panobinostat treatment, through transcriptional activation of IFN genes concomitant with increased H3K27 acetylation at these loci. Depletion of pDCs abrogated panobinostat-mediated induction of type I IFN signaling in leukemia cells and impaired therapeutic efficacy, whereas combined treatment with panobinostat and IFNα improved outcomes in preclinical models. These discoveries offer a new therapeutic approach for AML and demonstrate that epigenetic rewiring of pDCs enhances antitumor immunity, opening the possibility of exploiting this approach for immunotherapies. SIGNIFICANCE: We demonstrate that HDACis induce terminal differentiation of AML through epigenetic remodeling of pDCs, resulting in production of type I IFN that is important for the therapeutic effects of HDACis. The study demonstrates the important functional interplay between the immune system and leukemias in response to HDAC inhibition. This article is highlighted in the In This Issue feature, p. 1397.

Evolution of late-stage metastatic melanoma is dominated by aneuploidy and whole genome doubling.

Although melanoma is initiated by acquisition of point mutations and limited focal copy number alterations in melanocytes-of-origin, the nature of genetic changes that characterise lethal metastatic disease is poorly understood. Here, we analyze the evolution of human melanoma progressing from early to late disease in 13 patients by sampling their tumours at multiple sites and times. Whole exome and genome sequencing data from 88 tumour samples reveals only limited gain of point mutations generally, with net mutational loss in some metastases. In contrast, melanoma evolution is dominated by whole genome doubling and large-scale aneuploidy, in which widespread loss of heterozygosity sculpts the burden of point mutations, neoantigens and structural variants even in treatment-naïve and primary cutaneous melanomas in some patients. These results imply that dysregulation of genomic integrity is a key driver of selective clonal advantage during melanoma progression.

Biodosimetric transcriptional and proteomic changes are conserved in irradiated human tissue.

Transcriptional dosimetry is an emergent field of radiobiology aimed at developing robust methods for detecting and quantifying absorbed doses using radiation-induced fluctuations in gene expression. A combination of RNA sequencing, array-based and quantitative PCR transcriptomics in cellular, murine and various ex vivo human models has led to a comprehensive description of a fundamental set of genes with demonstrable dosimetric qualities. However, these are yet to be validated in human tissue due to the scarcity of in situ-irradiated source material. This represents a major hurdle to the continued development of transcriptional dosimetry. In this study, we present a novel evaluation of a previously reported set of dosimetric genes in human tissue exposed to a large therapeutic dose of radiation. To do this, we evaluated the quantitative changes of a set of dosimetric transcripts consisting of FDXR, BAX, BCL2, CDKN1A, DDB2, BBC3, GADD45A, GDF15, MDM2, SERPINE1, TNFRSF10B, PLK3, SESN2 and VWCE in guided pre- and post-radiation (2 weeks) prostate cancer biopsies from seven patients. We confirmed the prolonged dose-responsivity of most of these transcripts in in situ-irradiated tissue. BCL2, GDF15, and to some extent TNFRSF10B, were markedly unreliable single markers of radiation exposure. Nevertheless, as a full set, these genes reliably segregated non-irradiated and irradiated tissues and predicted radiation absorption on a patient-specific basis. We also confirmed changes in the translated protein product for a small subset of these dosimeters. This study provides the first confirmatory evidence of an existing dosimetric gene set in less-accessible tissues-ensuring peripheral responses reflect tissue-specific effects. Further work will be required to determine if these changes are conserved in different tissue types, post-radiation times and doses.

A Comprehensive Protocol Resource for Performing Pooled shRNA and CRISPR Screens.

This chapter details a compendium of protocols that collectively enable the reader to perform a pooled shRNA and/or CRISPR screen-with methods to identify and validate positive controls and subsequent hits; establish a viral titer in the cell line of choice; create and screen libraries, sequence strategies, and bioinformatics resources to analyze outcomes. Collectively, this provides an overarching resource from the start to finish of a screening project, making this technology possible in all laboratories.

The Transcriptional Landscape of Radiation-Treated Human Prostate Cancer: Analysis of a Prospective Tissue Cohort.

PURPOSE: The resistance of prostate cancer to radiation therapy (RT) is a significant clinical issue and still largely unable to be guided by patient-specific molecular characteristics. The present study describes the gene expression changes induced in response to RT in human prostate tissue obtained from a prospective tissue acquisition study designed for radiobiology research. METHODS AND MATERIALS: A prospective cohort of 5 men with intermediate-risk and clinically localized tumors were treated with high-dose-rate brachytherapy with 2 × 10-Gy fractions. Image-guided transperineal biopsy specimens were taken immediately before and 14 days after the first high-dose-rate brachytherapy fraction. Using genome-wide 3' RNA sequencing on total RNA extracted from 10 biopsy specimens, we obtained quantitative expression data for a median of 13,244 genes. We computed the fold-change information for each gene and extracted high-confidence lists of transcripts with either increased or decreased expression (≥1.5-fold) after radiation in ≥4 of the 5 patients. Several gene ontology analyses were then used to identify functionally enriched pathways. RESULTS: The predominant change in response to RT was elevation of the transcript levels, including that of DNA damage binding protein 2 and p21, and collagens, laminins, and integrins. We observed strong upregulation of the p53 pathway, without observable dysregulation of p53 itself. Interstitial remodeling, extracellular matrix proteins, and focal adhesion pathways were also strongly upregulated, as was inflammation. Functional network analysis showed clustering of the changes inherent in apoptosis and programmed cell death, extracellular matrix organization, and immune regulation. CONCLUSIONS: In the present prospective study of matched clinical tissues, we successfully recognized known radiation-sensitive transcriptional pathways and identified numerous other novel and significantly altered genes with no current association with RT. These data could be informative in the development of future personalized therapeutic agents.

Homologous Recombination DNA Repair Pathway Disruption and Retinoblastoma Protein Loss Are Associated with Exceptional Survival in High-Grade Serous Ovarian Cancer.

Purpose: Women with epithelial ovarian cancer generally have a poor prognosis; however, a subset of patients has an unexpected dramatic and durable response to treatment. We sought to identify clinical, pathological, and molecular determinants of exceptional survival in women with high-grade serous cancer (HGSC), a disease associated with the majority of ovarian cancer deaths.Experimental Design: We evaluated the histories of 2,283 ovarian cancer patients and, after applying stringent clinical and pathological selection criteria, identified 96 with HGSC that represented significant outliers in terms of treatment response and overall survival. Patient samples were characterized immunohistochemically and by genome sequencing.Results: Different patterns of clinical response were seen: long progression-free survival (Long-PFS), multiple objective responses to chemotherapy (Multiple Responder), and/or greater than 10-year overall survival (Long-Term Survivors). Pathogenic germline and somatic mutations in genes involved in homologous recombination (HR) repair were enriched in all three groups relative to a population-based series. However, 29% of 10-year survivors lacked an identifiable HR pathway alteration, and tumors from these patients had increased Ki-67 staining. CD8+ tumor-infiltrating lymphocytes were more commonly present in Long-Term Survivors. RB1 loss was associated with long progression-free and overall survival. HR deficiency and RB1 loss were correlated, and co-occurrence was significantly associated with prolonged survival.Conclusions: There was diversity in the clinical trajectory of exceptional survivors associated with multiple molecular determinants of exceptional outcome in HGSC patients. Concurrent HR deficiency and RB1 loss were associated with favorable outcomes, suggesting that co-occurrence of specific mutations might mediate durable responses in such patients. Clin Cancer Res; 24(3); 569-80. ©2017 AACRSee related commentary by Peng and Mills, p. 508.

Circulating tumour DNA reflects treatment response and clonal evolution in chronic lymphocytic leukaemia.

Several novel therapeutics are poised to change the natural history of chronic lymphocytic leukaemia (CLL) and the increasing use of these therapies has highlighted limitations of traditional disease monitoring methods. Here we demonstrate that circulating tumour DNA (ctDNA) is readily detectable in patients with CLL. Importantly, ctDNA does not simply mirror the genomic information contained within circulating malignant lymphocytes but instead parallels changes across different disease compartments following treatment with novel therapies. Serial ctDNA analysis allows clonal dynamics to be monitored over time and identifies the emergence of genomic changes associated with Richter's syndrome (RS). In addition to conventional disease monitoring, ctDNA provides a unique opportunity for non-invasive serial analysis of CLL for molecular disease monitoring.

Multiplexed transcriptome analysis to detect ALK, ROS1 and RET rearrangements in lung cancer.

ALK, ROS1 and RET gene fusions are important predictive biomarkers for tyrosine kinase inhibitors in lung cancer. Currently, the gold standard method for gene fusion detection is Fluorescence In Situ Hybridization (FISH) and while highly sensitive and specific, it is also labour intensive, subjective in analysis, and unable to screen a large numbers of gene fusions. Recent developments in high-throughput transcriptome-based methods may provide a suitable alternative to FISH as they are compatible with multiplexing and diagnostic workflows. However, the concordance between these different methods compared with FISH has not been evaluated. In this study we compared the results from three transcriptome-based platforms (Nanostring Elements, Agena LungFusion panel and ThermoFisher NGS fusion panel) to those obtained from ALK, ROS1 and RET FISH on 51 clinical specimens. Overall agreement of results ranged from 86-96% depending on the platform used. While all platforms were highly sensitive, both the Agena panel and Thermo Fisher NGS fusion panel reported minor fusions that were not detectable by FISH. Our proof-of-principle study illustrates that transcriptome-based analyses are sensitive and robust methods for detecting actionable gene fusions in lung cancer and could provide a robust alternative to FISH testing in the diagnostic setting.

Circulating Tumor DNA Analysis and Functional Imaging Provide Complementary Approaches for Comprehensive Disease Monitoring in Metastatic Melanoma.

PURPOSE: Circulating tumor DNA (ctDNA) allows noninvasive disease monitoring across a range of malignancies. In metastatic melanoma, the extent to which ctDNA reflects changes in metabolic disease burden assessed by 18F-labeled fluorodeoxyglucose positron emission tomography (FDG-PET) is unknown. We assessed the role of ctDNA analysis in combination with FDG-PET to monitor tumor burden and genomic heterogeneity throughout treatment. PATIENTS AND METHODS: We performed a comprehensive analysis of serial ctDNA and FDG-PET in 52 patients who received systemic therapy for metastatic melanoma. Next-generation sequencing and digital polymerase chain reaction were used to analyze plasma samples from the cohort. RESULTS: ctDNA levels were monitored across patients with mutant BRAF, NRAS, and BRAF/NRAS wild type disease. Mutant BRAF and NRAS ctDNA levels correlated closely with changes in metabolic disease burden throughout treatment. TERT promoter mutant ctDNA levels also paralleled changes in tumor burden, which provide an alternative marker for disease monitoring. Of note, subcutaneous and cerebral disease sites were not well represented in plasma. Early changes in ctDNA and metabolic disease burden were important indicators of treatment response. Patients with an early decrease in ctDNA post-treatment had improved progression-free survival compared with patients in whom ctDNA levels remained unchanged or increased over time (hazard ratio, 2.6; P = .05). ctDNA analysis contributed key molecular information through the identification of putative resistance mechanisms to targeted therapy. A detailed comparison of the genomic architecture of plasma and multiregional tumor biopsy specimens at autopsy revealed the ability of ctDNA to comprehensively capture genomic heterogeneity across multiple disease sites. CONCLUSION: The findings highlight the powerful role of ctDNA in metastatic melanoma as a complementary modality to functional imaging that allows real-time monitoring of both tumor burden and genomic changes throughout therapy.

UV-Associated Mutations Underlie the Etiology of MCV-Negative Merkel Cell Carcinomas.

Merkel cell carcinoma (MCC) is an uncommon, but highly malignant, cutaneous tumor. Merkel cell polyoma virus (MCV) has been implicated in a majority of MCC tumors; however, viral-negative tumors have been reported to be more prevalent in some geographic regions subject to high sun exposure. While the impact of MCV and viral T-antigens on MCC development has been extensively investigated, little is known about the etiology of viral-negative tumors. We performed targeted capture and massively parallel DNA sequencing of 619 cancer genes to compare the gene mutations and copy number alterations in MCV-positive (n = 13) and -negative (n = 21) MCC tumors and cell lines. We found that MCV-positive tumors displayed very low mutation rates, but MCV-negative tumors exhibited a high mutation burden associated with a UV-induced DNA damage signature. All viral-negative tumors harbored mutations in RB1, TP53, and a high frequency of mutations in NOTCH1 and FAT1. Additional mutated or amplified cancer genes of potential clinical importance included PI3K (PIK3CA, AKT1, PIK3CG) and MAPK (HRAS, NF1) pathway members and the receptor tyrosine kinase FGFR2. Furthermore, looking ahead to potential therapeutic strategies encompassing immune checkpoint inhibitors such as anti-PD-L1, we also assessed the status of T-cell-infiltrating lymphocytes (TIL) and PD-L1 in MCC tumors. A subset of viral-negative tumors exhibited high TILs and PD-L1 expression, corresponding with the higher mutation load within these cancers. Taken together, this study provides new insights into the underlying biology of viral-negative MCC and paves the road for further investigation into new treatment opportunities.

Whole-genome characterization of chemoresistant ovarian cancer.

Patients with high-grade serous ovarian cancer (HGSC) have experienced little improvement in overall survival, and standard treatment has not advanced beyond platinum-based combination chemotherapy, during the past 30 years. To understand the drivers of clinical phenotypes better, here we use whole-genome sequencing of tumour and germline DNA samples from 92 patients with primary refractory, resistant, sensitive and matched acquired resistant disease. We show that gene breakage commonly inactivates the tumour suppressors RB1, NF1, RAD51B and PTEN in HGSC, and contributes to acquired chemotherapy resistance. CCNE1 amplification was common in primary resistant and refractory disease. We observed several molecular events associated with acquired resistance, including multiple independent reversions of germline BRCA1 or BRCA2 mutations in individual patients, loss of BRCA1 promoter methylation, an alteration in molecular subtype, and recurrent promoter fusion associated with overexpression of the drug efflux pump MDR1.