Distinct and predictive chromatin signatures of transcriptional promoters and enhancers in the human genome.

Eukaryotic gene transcription is accompanied by acetylation and methylation of nucleosomes near promoters, but the locations and roles of histone modifications elsewhere in the genome remain unclear. We determined the chromatin modification states in high resolution along 30 Mb of the human genome and found that active promoters are marked by trimethylation of Lys4 of histone H3 (H3K4), whereas enhancers are marked by monomethylation, but not trimethylation, of H3K4. We developed computational algorithms using these distinct chromatin signatures to identify new regulatory elements, predicting over 200 promoters and 400 enhancers within the 30-Mb region. This approach accurately predicted the location and function of independently identified regulatory elements with high sensitivity and specificity and uncovered a novel functional enhancer for the carnitine transporter SLC22A5 (OCTN2). Our results give insight into the connections between chromatin modifications and transcriptional regulatory activity and provide a new tool for the functional annotation of the human genome.

Histone modifications at human enhancers reflect global cell-type-specific gene expression.

The human body is composed of diverse cell types with distinct functions. Although it is known that lineage specification depends on cell-specific gene expression, which in turn is driven by promoters, enhancers, insulators and other cis-regulatory DNA sequences for each gene, the relative roles of these regulatory elements in this process are not clear. We have previously developed a chromatin-immunoprecipitation-based microarray method (ChIP-chip) to locate promoters, enhancers and insulators in the human genome. Here we use the same approach to identify these elements in multiple cell types and investigate their roles in cell-type-specific gene expression. We observed that the chromatin state at promoters and CTCF-binding at insulators is largely invariant across diverse cell types. In contrast, enhancers are marked with highly cell-type-specific histone modification patterns, strongly correlate to cell-type-specific gene expression programs on a global scale, and are functionally active in a cell-type-specific manner. Our results define over 55,000 potential transcriptional enhancers in the human genome, significantly expanding the current catalogue of human enhancers and highlighting the role of these elements in cell-type-specific gene expression.

SHP2 Inhibition Sensitizes Diverse Oncogene-Addicted Solid Tumors to Re-treatment with Targeted Therapy.

Rationally targeted therapies have transformed cancer treatment, but many patients develop resistance through bypass signaling pathway activation. PF-07284892 (ARRY-558) is an allosteric SHP2 inhibitor designed to overcome bypass-signaling-mediated resistance when combined with inhibitors of various oncogenic drivers. Activity in this setting was confirmed in diverse tumor models. Patients with ALK fusion-positive lung cancer, BRAFV600E-mutant colorectal cancer, KRASG12D-mutant ovarian cancer, and ROS1 fusion-positive pancreatic cancer who previously developed targeted therapy resistance were treated with PF-07284892 on the first dose level of a first-in-human clinical trial. After progression on PF-07284892 monotherapy, a novel study design allowed the addition of oncogene-directed targeted therapy that had previously failed. Combination therapy led to rapid tumor and circulating tumor DNA (ctDNA) responses and extended the duration of overall clinical benefit. SIGNIFICANCE: PF-07284892-targeted therapy combinations overcame bypass-signaling-mediated resistance in a clinical setting in which neither component was active on its own. This provides proof of concept of the utility of SHP2 inhibitors in overcoming resistance to diverse targeted therapies and provides a paradigm for accelerated testing of novel drug combinations early in clinical development. See related commentary by Hernando-Calvo and Garralda, p. 1762. This article is highlighted in the In This Issue feature, p. 1749.

Molecular characterization of renal cell carcinoma tumors from a phase III anti-angiogenic adjuvant therapy trial.

Multigene assays can provide insight into key biological processes and prognostic information to guide development and selection of adjuvant cancer therapy. We report a comprehensive genomic and transcriptomic analysis of tumor samples from 171 patients at high risk for recurrent renal cell carcinoma post nephrectomy from the S-TRAC trial (NCT00375674). We identify gene expression signatures, including STRAC11 (derived from the sunitinib-treated population). The overlap in key elements captured in these gene expression signatures, which include genes representative of the tumor stroma microenvironment, regulatory T cell, and myeloid cells, suggests they are likely to be both prognostic and predictive of the anti-angiogenic effect in the adjuvant setting. These signatures also point to the identification of potential therapeutic targets for development in adjuvant renal cell carcinoma, such as MERTK and TDO2. Finally, our findings suggest that while anti-angiogenic adjuvant therapy might be important, it may not be sufficient to prevent recurrence and that other factors such as immune response and tumor environment may be of greater importance.

Avelumab maintenance in advanced urothelial carcinoma: biomarker analysis of the phase 3 JAVELIN Bladder 100 trial.

In a recent phase 3 randomized trial of 700 patients with advanced urothelial cancer (JAVELIN Bladder 100; NCT02603432 ), avelumab/best supportive care (BSC) significantly prolonged overall survival relative to BSC alone as maintenance therapy after first-line chemotherapy. Exploratory biomarker analyses were performed to identify biological pathways that might affect survival benefit. Tumor molecular profiling by immunohistochemistry, whole-exome sequencing and whole-transcriptome sequencing revealed that avelumab survival benefit was positively associated with PD-L1 expression by tumor cells, tumor mutational burden, APOBEC mutation signatures, expression of genes underlying innate and adaptive immune activity and the number of alleles encoding high-affinity variants of activating Fcγ receptors. Pathways connected to tissue growth and angiogenesis might have been associated with reduced survival benefit. Individual biomarkers did not comprehensively identify patients who could benefit from therapy; however, multi-parameter models incorporating genomic alteration, immune responses and tumor growth showed promising predictive utility. These results characterize the complex biologic pathways underlying survival benefit from immune checkpoint inhibition in advanced urothelial cancer and suggest that multiple biomarkers might be needed to identify patients who would benefit from treatment.

An Engineered IL15 Cytokine Mutein Fused to an Anti-PD1 Improves Intratumoral T-cell Function and Antitumor Immunity.

The use of cytokines for immunotherapy shows clinical efficacy but is frequently accompanied by severe adverse events caused by excessive and systemic immune activation. Here, we set out to address these challenges by engineering a fusion protein of a single, potency-reduced, IL15 mutein and a PD1-specific antibody (anti-PD1-IL15m). This immunocytokine was designed to deliver PD1-mediated, avidity-driven IL2/15 receptor stimulation to PD1+ tumor-infiltrating lymphocytes (TIL) while minimally affecting circulating peripheral natural killer (NK) cells and T cells. Treatment of tumor-bearing mice with a mouse cross-reactive fusion, anti-mPD1-IL15m, demonstrated potent antitumor efficacy without exacerbating body weight loss in B16 and MC38 syngeneic tumor models. Moreover, anti-mPD1-IL15m was more efficacious than an IL15 superagonist, an anti-mPD-1, or the combination thereof in the B16 melanoma model. Mechanistically, anti-PD1-IL15m preferentially targeted CD8+ TILs and single-cell RNA-sequencing analyses revealed that anti-mPD1-IL15m treatment induced the expansion of an exhausted CD8+ TIL cluster with high proliferative capacity and effector-like signatures. Antitumor efficacy of anti-mPD1-IL15m was dependent on CD8+ T cells, as depletion of CD8+ cells resulted in the loss of antitumor activity, whereas depletion of NK cells had little impact on efficacy. The impact of anti-hPD1-IL15m on primary human TILs from patients with cancer was also evaluated. Anti-hPD1-IL15m robustly enhanced the proliferation, activation, and cytotoxicity of CD8+ and CD4+ TILs from human primary cancers in vitro, whereas tumor-derived regulatory T cells were largely unaffected. Taken together, our findings showed that anti-PD1-IL15m exhibits a high translational promise with improved efficacy and safety of IL15 for cancer immunotherapy via targeting PD1+ TILs.See related Spotlight by Felices and Miller, p. 1110.

Symmetric Arginine Dimethylation Is Selectively Required for mRNA Splicing and the Initiation of Type I and Type III Interferon Signaling.

Alternative splicing is well understood to enhance proteome diversity as cells respond to stimuli. However, mechanistic understanding for how the spliceosome processes precursor messenger RNA (mRNA) transcripts to achieve template diversification is incomplete. We use recently developed enzymatic inhibitors of protein arginine methyltransferase 5 (PRMT5) and human naive T lymphocyte activation as a model system to uncover a precise set of mRNA transcripts that require symmetric arginine dimethylation. This methylation-dependent splicing selectivity is associated with a limited set of signaling pathways that are affected when PRMT5 is inhibited. Specifically, we identify a conserved role for symmetric arginine dimethylation in the induction of antiviral type I and type III interferon signaling following T cell receptor and pattern recognition receptor stimulation in human T lymphocytes and undifferentiated human THP-1 monocytes. Altogether, these findings reveal a mechanism by which cells may be enabled to precisely modulate transcript heterogeneity to orchestrate specific functional outcomes.

First-in-Human Study of Utomilumab, a 4-1BB/CD137 Agonist, in Combination with Rituximab in Patients with Follicular and Other CD20+ Non-Hodgkin Lymphomas.

PURPOSE: In this phase I study (NCT01307267), we evaluated safety, pharmacokinetics, clinical activity, and pharmacodynamics of treatment with utomilumab plus rituximab in patients with relapsed/refractory follicular lymphoma (FL) and other CD20+ non-Hodgkin lymphomas (NHL). PATIENTS AND METHODS: Primary objectives were to assess treatment safety and tolerability for estimating the MTD, using a modified time-to-event continual reassessment method, and selecting the recommended phase II dose (RP2D). RESULTS: Sixty-seven patients received utomilumab (0.03-10.0 mg/kg every 4 weeks) and rituximab (375 mg/m2 weekly) in the dose-escalation groups or utomilumab (1.2 mg/kg every 4 weeks) plus rituximab in the dose-expansion cohort. No patient experienced dose-limiting toxicity. The MTD for utomilumab in combination with rituximab was not reached and estimated to be ≥10 mg/kg every 4 weeks. The majority of the utomilumab treatment-related adverse events (AE) were grade 1 to 2; the most common AE was fatigue (16.4%). The pharmacokinetics of utomilumab in combination with rituximab was linear in the 0.03 to 10 mg/kg dose range. A low incidence (1.5%) of treatment-induced antidrug antibodies against utomilumab was observed. The objective response rate was 21.2% (95% CI, 12.1%-33.0%) in all patients with NHL, including four complete and 10 partial responses. Analysis of paired biopsies from a relapsed/refractory FL patient with complete response showed increased T-cell infiltration and cytotoxic activity in tumors. Biomarker correlations with outcomes suggested that clinical benefit may be contingent on patient immune function. CONCLUSIONS: Utomilumab in combination with rituximab demonstrated clinical activity and a favorable safety profile in patients with CD20+ NHLs.

Avelumab plus axitinib versus sunitinib in advanced renal cell carcinoma: biomarker analysis of the phase 3 JAVELIN Renal 101 trial.

We report on molecular analyses of baseline tumor samples from the phase 3 JAVELIN Renal 101 trial (n = 886; NCT02684006 ), which demonstrated significantly prolonged progression-free survival (PFS) with first-line avelumab + axitinib versus sunitinib in advanced renal cell carcinoma (aRCC). We found that neither expression of the commonly assessed biomarker programmed cell death ligand 1 (PD-L1) nor tumor mutational burden differentiated PFS in either study arm. Similarly, the presence of FcɣR single nucleotide polymorphisms was unimpactful. We identified important biological features associated with differential PFS between the treatment arms, including new immunomodulatory and angiogenesis gene expression signatures (GESs), previously undescribed mutational profiles and their corresponding GESs, and several HLA types. These findings provide insight into the determinants of response to combined PD-1/PD-L1 and angiogenic pathway inhibition and may aid in the development of strategies for improved patient care in aRCC.

Chemotherapy induces dynamic immune responses in breast cancers that impact treatment outcome.

To elucidate the effects of neoadjuvant chemotherapy (NAC), we conduct whole transcriptome profiling coupled with histopathology analyses of a longitudinal breast cancer cohort of 146 patients including 110 pairs of serial tumor biopsies collected before treatment, after the first cycle of treatment and at the time of surgery. Here, we show that cytotoxic chemotherapies induce dynamic changes in the tumor immune microenvironment that vary by subtype and pathologic response. Just one cycle of treatment induces an immune stimulatory microenvironment harboring more tumor infiltrating lymphocytes (TILs) and up-regulation of inflammatory signatures predictive of response to anti-PD1 therapies while residual tumors are immune suppressed at end-of-treatment compared to the baseline. Increases in TILs and CD8+ T cell proportions in response to NAC are independently associated with pathologic complete response. Further, on-treatment immune response is more predictive of treatment outcome than immune features in paired baseline samples although these are strongly correlated.

Cell Index Database (CELLX): a web tool for cancer precision medicine.

The Cell Index Database, (CELLX) (http://cellx.sourceforge.net) provides a computational framework for integrating expression, copy number variation, mutation, compound activity, and meta data from cancer cells. CELLX provides the computational biologist a quick way to perform routine analyses as well as the means to rapidly integrate data for offline analysis. Data is accessible through a web interface which utilizes R to generate plots and perform clustering, correlations, and statistical tests for associations within and between data types for ~20,000 samples from TCGA, CCLE, Sanger, GSK, GEO, GTEx, and other public sources. We show how CELLX supports precision oncology through indications discovery, biomarker evaluation, and cell line screening analysis.

Molecular predictors of sensitivity to the insulin-like growth factor 1 receptor inhibitor Figitumumab (CP-751,871).

Figitumumab (CP-751,871), a potent and fully human monoclonal anti-insulin-like growth factor 1 receptor (IGF1R) antibody, has been investigated in clinical trials of several solid tumors. To identify biomarkers of sensitivity and resistance to figitumumab, its in vitro antiproliferative activity was analyzed in a panel of 93 cancer cell lines by combining in vitro screens with extensive molecular profiling of genomic aberrations. Overall response was bimodal and the majority of cell lines were resistant to figitumumab. Nine of 15 sensitive cell lines were derived from colon cancers. Correlations between genomic characteristics of cancer cell lines with figitumumab antiproliferative activity revealed that components of the IGF pathway, including IRS2 (insulin receptor substrate 2) and IGFBP5 (IGF-binding protein 5), played a pivotal role in determining the sensitivity of tumors to single-agent figitumumab. Tissue-specific differences among the top predictive genes highlight the need for tumor-specific patient selection strategies. For the first time, we report that alteration or expression of the MYB oncogene is associated with sensitivity to IGF1R inhibitors. MYB is dysregulated in hematologic and epithelial tumors, and IGF1R inhibition may represent a novel therapeutic opportunity. Although growth inhibitory activity with single-agent figitumumab was relatively rare, nine combinations comprising figitumumab plus chemotherapeutic agents or other targeted agents exhibited properties of synergy. Inhibitors of the ERBB family were frequently synergistic and potential biomarkers of drug synergy were identified. Several biomarkers of antiproliferative activity of figitumumab both alone and in combination with other therapies may inform the design of clinical trials evaluating IGF1R inhibitors.

Evaluation of the antitumor effects and mechanisms of PF00299804, a pan-HER inhibitor, alone or in combination with chemotherapy or targeted agents in gastric cancer.

Recently, HER2-directed treatment, such as trastuzumab, has shown clinical benefit in HER2-amplified gastric cancer. On the basis of recent studies about epidermal growth factor receptor (EGFR) or HER2-targeting agents (including gefitinib, lapatinib, and trastuzumab) in gastric cancer, the potent effects of pan-HER inhibitors targeting the HER family are anticipated. In this study, we evaluated the activity and mechanisms of PF00299804, an irreversible pan-HER inhibitor, in gastric cancer in vitro and in vivo models. PF00299804 showed significant growth-inhibitory effects in HER2-amplified gastric cancer cells (SNU216, N87), and it had lower 50% inhibitory concentration values compared with other EGFR tyrosine kinase inhibitors, including gefitinib, lapatinib, BIBW-2992, and CI-1033. PF00299804 induced apoptosis and G(1) arrest and inhibited phosphorylation of receptors in the HER family and downstream signaling pathways including STAT3, AKT, and extracellular signal-regulated kinases (ERK) in HER2-amplified gastric cancer cells. PF00299804 also blocked EGFR/HER2, HER2/HER3, and HER3/HER4 heterodimer formation as well as the association of HER3 with p85α in SNU216 cells. The combination of PF00299804 with clinically relevant chemotherapeutic agents or molecular-targeted agents including trastuzumab (an anti-HER2 monoclonal antibody), CP751871 (an IGF1R inhibitor), PD0325901 (an ERK1/2 inhibitor), and PF04691502 (a PI3K/mTOR inhibitor) produced synergistic effects. These findings indicate that PF00299804 can be used as a targeted therapy for the treatment of HER2-amplified gastric cancer through inhibition of HER family heterodimer formation and may augment antitumor efficacy of chemotherapeutic and/or molecular-targeted agents.

An integrated genomic approach to identify predictive biomarkers of response to the aurora kinase inhibitor PF-03814735.

PF-03814735 is a novel, reversible inhibitor of Aurora kinases A and B that finished a phase I clinical trial for the treatment of advanced solid tumors. To find predictive biomarkers of drug sensitivity, we screened a diverse panel of 87 cancer cell lines for growth inhibition upon PF-03814735 treatment. Small cell lung cancer (SCLC) and, to a lesser extent, colon cancer lines were very sensitive to PF-03814735. The status of the Myc gene family and retinoblastoma pathway members significantly correlated with the efficacy of PF-03814735. Whereas RB1 inactivation, intact CDKN2A/p16, and normal CCND1/Cyclin D1 status are hallmarks of SCLC, activation or amplification of any of the three Myc genes (MYC, MYCL1, and MYCN) clearly differentiated cell line sensitivity within the SCLC panel. By contrast, we found that expression of Aurora A and B were weak predictors of response. We observed a decrease in histone H3 phosphorylation and polyploidization of sensitive lines, consistent with the phenotype of Aurora B inhibition. In vivo experiments with two SCLC xenograft models confirmed the sensitivity of Myc gene-driven models to PF-03814735 and a possible schedule dependence of MYC/c-Myc-driven tumors. Altogether our results suggest that SCLC and other malignancies driven by the Myc family genes may be suitable indications for treatment by Aurora B kinase inhibitors.

A novel SND1-BRAF fusion confers resistance to c-Met inhibitor PF-04217903 in GTL16 cells through [corrected] MAPK activation.

Targeting cancers with amplified or abnormally activated c-Met (hepatocyte growth factor receptor) may have therapeutic benefit based on nonclinical and emerging clinical findings. However, the eventual emergence of drug resistant tumors motivates the pre-emptive identification of potential mechanisms of clinical resistance. We rendered a MET amplified gastric cancer cell line, GTL16, resistant to c-Met inhibition with prolonged exposure to a c-Met inhibitor, PF-04217903 (METi). Characterization of surviving cells identified an amplified chromosomal rearrangement between 7q32 and 7q34 which overexpresses a constitutively active SND1-BRAF fusion protein. In the resistant clones, hyperactivation of the downstream MAPK pathway via SND1-BRAF conferred resistance to c-Met receptor tyrosine kinase inhibition. Combination treatment with METi and a RAF inhibitor, PF-04880594 (RAFi) inhibited ERK activation and circumvented resistance to either single agent. Alternatively, treatment with a MEK inhibitor, PD-0325901 (MEKi) alone effectively blocked ERK phosphorylation and inhibited cell growth. Our results suggest that combination of a c-Met tyrosine kinase inhibitor with a BRAF or a MEK inhibitor may be effective in treating resistant tumors that use activated BRAF to escape suppression of c-Met signaling.

Dynamic chromatin states in human ES cells reveal potential regulatory sequences and genes involved in pluripotency.

Pluripotency, the ability of a cell to differentiate and give rise to all embryonic lineages, defines a small number of mammalian cell types such as embryonic stem (ES) cells. While it has been generally held that pluripotency is the product of a transcriptional regulatory network that activates and maintains the expression of key stem cell genes, accumulating evidence is pointing to a critical role for epigenetic processes in establishing and safeguarding the pluripotency of ES cells, as well as maintaining the identity of differentiated cell types. In order to better understand the role of epigenetic mechanisms in pluripotency, we have examined the dynamics of chromatin modifications genome-wide in human ES cells (hESCs) undergoing differentiation into a mesendodermal lineage. We found that chromatin modifications at promoters remain largely invariant during differentiation, except at a small number of promoters where a dynamic switch between acetylation and methylation at H3K27 marks the transition between activation and silencing of gene expression, suggesting a hierarchy in cell fate commitment over most differentially expressed genes. We also mapped over 50 000 potential enhancers, and observed much greater dynamics in chromatin modifications, especially H3K4me1 and H3K27ac, which correlate with expression of their potential target genes. Further analysis of these enhancers revealed potentially key transcriptional regulators of pluripotency and a chromatin signature indicative of a poised state that may confer developmental competence in hESCs. Our results provide new evidence supporting the role of chromatin modifications in defining enhancers and pluripotency.

Analysis of the vertebrate insulator protein CTCF-binding sites in the human genome.

Insulator elements affect gene expression by preventing the spread of heterochromatin and restricting transcriptional enhancers from activation of unrelated promoters. In vertebrates, insulator's function requires association with the CCCTC-binding factor (CTCF), a protein that recognizes long and diverse nucleotide sequences. While insulators are critical in gene regulation, only a few have been reported. Here, we describe 13,804 CTCF-binding sites in potential insulators of the human genome, discovered experimentally in primary human fibroblasts. Most of these sequences are located far from the transcriptional start sites, with their distribution strongly correlated with genes. The majority of them fit to a consensus motif highly conserved and suitable for predicting possible insulators driven by CTCF in other vertebrate genomes. In addition, CTCF localization is largely invariant across different cell types. Our results provide a resource for investigating insulator function and possible other general and evolutionarily conserved activities of CTCF sites.

Data and animal management software for large-scale phenotype screening.

The mouse N-ethyl-N-nitrosourea (ENU) mutagenesis program at the Genomics Institute of the Novartis Research Foundation (GNF) uses MouseTRACS to analyze phenotype screens and manage animal husbandry. MouseTRACS is a Web-based laboratory informatics system that electronically records and organizes mouse colony operations, prints cage cards, tracks inventory, manages requests, and reports Institutional Animal Care and Use Committee (IACUC) protocol usage. For efficient phenotype screening, MouseTRACS identifies mutants, visualizes data, and maps mutations. It displays and integrates phenotype and genotype data using likelihood odds ratio (LOD) plots of genetic linkage between genotype and phenotype. More detailed mapping intervals show individual single nucleotide polymorphism (SNP) markers in the context of phenotype. In addition, dynamically generated pedigree diagrams and inventory reports linked to screening results summarize the inheritance pattern and the degree of penetrance. MouseTRACS displays screening data in tables and uses standard charts such as box plots, histograms, scatter plots, and customized charts looking at clustered mice or cross pedigree comparisons. In summary, MouseTRACS enables the efficient screening, analysis, and management of thousands of animals to find mutant mice and identify novel gene functions. MouseTRACS is available under an open source license at http://www.mousetracs.sourceforge.net.