CXCL8/CXCR2 signaling mediates bone marrow fibrosis and is a therapeutic target in myelofibrosis.

Proinflammatory signaling is a hallmark feature of human cancer, including in myeloproliferative neoplasms (MPNs), most notably myelofibrosis (MF). Dysregulated inflammatory signaling contributes to fibrotic progression in MF; however, the individual cytokine mediators elicited by malignant MPN cells to promote collagen-producing fibrosis and disease evolution are yet to be fully elucidated. Previously, we identified a critical role for combined constitutive JAK/STAT and aberrant NF-κB proinflammatory signaling in MF development. Using single-cell transcriptional and cytokine-secretion studies of primary cells from patients with MF and the human MPLW515L (hMPLW515L) murine model of MF, we extend our previous work and delineate the role of CXCL8/CXCR2 signaling in MF pathogenesis and bone marrow fibrosis progression. Hematopoietic stem/progenitor cells from patients with MF are enriched for a CXCL8/CXCR2 gene signature and display enhanced proliferation and fitness in response to an exogenous CXCL8 ligand in vitro. Genetic deletion of Cxcr2 in the hMPLW515L-adoptive transfer model abrogates fibrosis and extends overall survival, and pharmacologic inhibition of the CXCR1/2 pathway improves hematologic parameters, attenuates bone marrow fibrosis, and synergizes with JAK inhibitor therapy. Our mechanistic insights provide a rationale for therapeutic targeting of the CXCL8/CXCR2 pathway among patients with MF.

Predicting tumour content of liquid biopsies from cell-free DNA.

BACKGROUND: Liquid biopsy is a minimally-invasive method of sampling bodily fluids, capable of revealing evidence of cancer. The distribution of cell-free DNA (cfDNA) fragment lengths has been shown to differ between healthy subjects and cancer patients, whereby the distributional shift correlates with the sample's tumour content. These fragmentomic data have not yet been utilised to directly quantify the proportion of tumour-derived cfDNA in a liquid biopsy. RESULTS: We used statistical learning to predict tumour content from Fourier and wavelet transforms of cfDNA length distributions in samples from 118 cancer patients. The model was validated on an independent dilution series of patient plasma. CONCLUSIONS: This proof of concept suggests that our fragmentomic methodology could be useful for predicting tumour content in liquid biopsies.

Ubiquitously transcribed genes use alternative polyadenylation to achieve tissue-specific expression.

More than half of human genes use alternative cleavage and polyadenylation (ApA) to generate mRNA transcripts that differ in the lengths of their 3' untranslated regions (UTRs), thus altering the post-transcriptional fate of the message and likely the protein output. The extent of 3' UTR variation across tissues and the functional role of ApA remain poorly understood. We developed a sequencing method called 3'-seq to quantitatively map the 3' ends of the transcriptome of diverse human tissues and isogenic transformation systems. We found that cell type-specific gene expression is accomplished by two complementary programs. Tissue-restricted genes tend to have single 3' UTRs, whereas a majority of ubiquitously transcribed genes generate multiple 3' UTRs. During transformation and differentiation, single-UTR genes change their mRNA abundance levels, while multi-UTR genes mostly change 3' UTR isoform ratios to achieve tissue specificity. However, both regulation programs target genes that function in the same pathways and processes that characterize the new cell type. Instead of finding global shifts in 3' UTR length during transformation and differentiation, we identify tissue-specific groups of multi-UTR genes that change their 3' UTR ratios; these changes in 3' UTR length are largely independent from changes in mRNA abundance. Finally, tissue-specific usage of ApA sites appears to be a mechanism for changing the landscape targetable by ubiquitously expressed microRNAs.

Development of Resistance to Type II JAK2 Inhibitors in MPN Depends on AXL Kinase and Is Targetable.

PURPOSE: Myeloproliferative neoplasms (MPN) dysregulate JAK2 signaling. Because clinical JAK2 inhibitors have limited disease-modifying effects, type II JAK2 inhibitors such as CHZ868 stabilizing inactive JAK2 and reducing MPN clones, gain interest. We studied whether MPN cells escape from type ll inhibition. EXPERIMENTAL DESIGN: MPN cells were continuously exposed to CHZ868. We used phosphoproteomic analyses and ATAC/RNA sequencing to characterize acquired resistance to type II JAK2 inhibition, and targeted candidate mediators in MPN cells and mice. RESULTS: MPN cells showed increased IC50 and reduced apoptosis upon CHZ868 reflecting acquired resistance to JAK2 inhibition. Among >2,500 differential phospho-sites, MAPK pathway activation was most prominent, while JAK2-STAT3/5 remained suppressed. Altered histone occupancy promoting AP-1/GATA binding motif exposure associated with upregulated AXL kinase and enriched RAS target gene profiles. AXL knockdown resensitized MPN cells and combined JAK2/AXL inhibition using bemcentinib or gilteritinib reduced IC50 to levels of sensitive cells. While resistant cells induced tumor growth in NOD/SCID gamma mice despite JAK2 inhibition, JAK2/AXL inhibition largely prevented tumor progression. Because inhibitors of MAPK pathway kinases such as MEK are clinically used in other malignancies, we evaluated JAK2/MAPK inhibition with trametinib to interfere with AXL/MAPK-induced resistance. Tumor growth was halted similarly to JAK2/AXL inhibition and in a systemic cell line-derived mouse model, marrow infiltration was decreased supporting dependency on AXL/MAPK. CONCLUSIONS: We report on a novel mechanism of AXL/MAPK-driven escape from type II JAK2 inhibition, which is targetable at different nodes. This highlights AXL as mediator of acquired resistance warranting inhibition to enhance sustainability of JAK2 inhibition in MPN.

Jak2V617F Reversible Activation Shows Its Essential Requirement in Myeloproliferative Neoplasms.

Gain-of-function mutations activating JAK/STAT signaling are seen in the majority of patients with myeloproliferative neoplasms (MPN), most commonly JAK2V617F. Although clinically approved JAK inhibitors improve symptoms and outcomes in MPNs, remissions are rare, and mutant allele burden does not substantively change with chronic therapy. We hypothesized this is due to limitations of current JAK inhibitors to potently and specifically abrogate mutant JAK2 signaling. We therefore developed a conditionally inducible mouse model allowing for sequential activation, and then inactivation, of Jak2V617F from its endogenous locus using a combined Dre-rox/Cre-lox dual-recombinase system. Jak2V617F deletion abrogates MPN features, induces depletion of mutant-specific hematopoietic stem/progenitor cells, and extends overall survival to an extent not observed with pharmacologic JAK inhibition, including when cooccurring with somatic Tet2 loss. Our data suggest JAK2V617F represents the best therapeutic target in MPNs and demonstrate the therapeutic relevance of a dual-recombinase system to assess mutant-specific oncogenic dependencies in vivo. SIGNIFICANCE: Current JAK inhibitors to treat myeloproliferative neoplasms are ineffective at eradicating mutant cells. We developed an endogenously expressed Jak2V617F dual-recombinase knock-in/knock-out model to investigate Jak2V617F oncogenic reversion in vivo. Jak2V617F deletion abrogates MPN features and depletes disease-sustaining MPN stem cells, suggesting improved Jak2V617F targeting offers the potential for greater therapeutic efficacy. See related commentary by Celik and Challen, p. 701. This article is featured in Selected Articles from This Issue, p. 695.

A reversible epigenetic memory of inflammatory injury controls lineage plasticity and tumor initiation in the mouse pancreas.

Inflammation is essential to the disruption of tissue homeostasis and can destabilize the identity of lineage-committed epithelial cells. Here, we employ lineage-traced mouse models, single-cell transcriptomic and chromatin analyses, and CUT&TAG to identify an epigenetic memory of inflammatory injury in the pancreatic acinar cell compartment. Despite resolution of pancreatitis, our data show that acinar cells fail to return to their molecular baseline, with retention of elevated chromatin accessibility and H3K4me1 at metaplasia genes, such that memory represents an incomplete cell fate decision. In vivo, we find this epigenetic memory controls lineage plasticity, with diminished metaplasia in response to a second insult but increased tumorigenesis with an oncogenic Kras mutation. The lowered threshold for oncogenic transformation, in turn, can be restored by blockade of MAPK signaling. Together, we define the chromatin dynamics, molecular encoding, and recall of a prolonged epigenetic memory of inflammatory injury that impacts future responses but remains reversible.

Tumor fraction-guided cell-free DNA profiling in metastatic solid tumor patients.

BACKGROUND: Cell-free DNA (cfDNA) profiling is increasingly used to guide cancer care, yet mutations are not always identified. The ability to detect somatic mutations in plasma depends on both assay sensitivity and the fraction of circulating DNA in plasma that is tumor-derived (i.e., cfDNA tumor fraction). We hypothesized that cfDNA tumor fraction could inform the interpretation of negative cfDNA results and guide the choice of subsequent assays of greater genomic breadth or depth. METHODS: Plasma samples collected from 118 metastatic cancer patients were analyzed with cf-IMPACT, a modified version of the FDA-authorized MSK-IMPACT tumor test that can detect genomic alterations in 410 cancer-associated genes. Shallow whole genome sequencing (sWGS) was also performed in the same samples to estimate cfDNA tumor fraction based on genome-wide copy number alterations using z-score statistics. Plasma samples with no somatic alterations detected by cf-IMPACT were triaged based on sWGS-estimated tumor fraction for analysis with either a less comprehensive but more sensitive assay (MSK-ACCESS) or broader whole exome sequencing (WES). RESULTS: cfDNA profiling using cf-IMPACT identified somatic mutations in 55/76 (72%) patients for whom MSK-IMPACT tumor profiling data were available. A significantly higher concordance of mutational profiles and tumor mutational burden (TMB) was observed between plasma and tumor profiling for plasma samples with a high tumor fraction (z-score≥5). In the 42 patients from whom tumor data was not available, cf-IMPACT identified mutations in 16/42 (38%). In total, cf-IMPACT analysis of plasma revealed mutations in 71/118 (60%) patients, with clinically actionable alterations identified in 30 (25%), including therapeutic targets of FDA-approved drugs. Of the 47 samples without alterations detected and low tumor fraction (z-score<5), 29 had sufficient material to be re-analyzed using a less comprehensive but more sensitive assay, MSK-ACCESS, which revealed somatic mutations in 14/29 (48%). Conversely, 5 patients without alterations detected by cf-IMPACT and with high tumor fraction (z-score≥5) were analyzed by WES, which identified mutational signatures and alterations in potential oncogenic drivers not covered by the cf-IMPACT panel. Overall, we identified mutations in 90/118 (76%) patients in the entire cohort using the three complementary plasma profiling approaches. CONCLUSIONS: cfDNA tumor fraction can inform the interpretation of negative cfDNA results and guide the selection of subsequent sequencing platforms that are most likely to identify clinically-relevant genomic alterations.

Enhanced specificity of clinical high-sensitivity tumor mutation profiling in cell-free DNA via paired normal sequencing using MSK-ACCESS.

Circulating cell-free DNA from blood plasma of cancer patients can be used to non-invasively interrogate somatic tumor alterations. Here we develop MSK-ACCESS (Memorial Sloan Kettering - Analysis of Circulating cfDNA to Examine Somatic Status), an NGS assay for detection of very low frequency somatic alterations in 129 genes. Analytical validation demonstrated 92% sensitivity in de-novo mutation calling down to 0.5% allele frequency and 99% for a priori mutation profiling. To evaluate the performance of MSK-ACCESS, we report results from 681 prospective blood samples that underwent clinical analysis to guide patient management. Somatic alterations are detected in 73% of the samples, 56% of which have clinically actionable alterations. The utilization of matched normal sequencing allows retention of somatic alterations while removing over 10,000 germline and clonal hematopoiesis variants. Our experience illustrates the importance of analyzing matched normal samples when interpreting cfDNA results and highlights the importance of cfDNA as a genomic profiling source for cancer patients.

Cell-free DNA profiling in retinoblastoma patients with advanced intraocular disease: An MSKCC experience.

PURPOSE: The enucleation rate for retinoblastoma has dropped from over 95% to under 10% in the past 10 years as a result of improvements in therapy. This reduces access to tumor tissue for molecular profiling, especially in unilateral retinoblastoma, and hinders the confirmation of somatic RB1 mutations necessary for genetic counseling. Plasma cell-free DNA (cfDNA) has provided a platform for noninvasive molecular profiling in cancer, but its applicability in low tumor burden retinoblastoma has not been shown. We analyzed cfDNA collected from 10 patients with available tumor tissue to determine whether sufficient tumorderived cfDNA is shed in plasma from retinoblastoma tumors to enable noninvasive RB1 mutation detection. METHODS: Tumor tissue was collected from eye enucleations in 10 patients diagnosed with advanced intra-ocular unilateral retinoblastoma, three of which went on to develop metastatic disease. Tumor RB1 mutation status was determined using an FDA-cleared tumor sequencing assay, MSK-IMPACT. Plasma samples were collected before eye enucleation and analyzed with a customized panel targeting all exons of RB1. RESULTS: Tumor-guided genotyping detected 10 of the 13 expected somatic RB1 mutations in plasma cfDNA in 8 of 10 patients (average variant allele frequency 3.78%). Without referring to RB1 status in the tumor, de novo mutation calling identified 7 of the 13 expected RB1 mutations (in 6 of 10 patients) with high confidence. CONCLUSION: Plasma cfDNA can detect somatic RB1 mutations in patients with unilateral retinoblastoma. Since intraocular biopsies are avoided in these patients because of concern about spreading tumor, cfDNA can potentially offer a noninvasive platform to guide clinical decisions about treatment, follow-up schemes, and risk of metastasis.

Affinity regression predicts the recognition code of nucleic acid-binding proteins.

Predicting the affinity profiles of nucleic acid-binding proteins directly from the protein sequence is a challenging problem. We present a statistical approach for learning the recognition code of a family of transcription factors or RNA-binding proteins (RBPs) from high-throughput binding data. Our method, called affinity regression, trains on protein binding microarray (PBM) or RNAcompete data to learn an interaction model between proteins and nucleic acids using only protein domain and probe sequences as inputs. When trained on mouse homeodomain PBM profiles, our model correctly identifies residues that confer DNA-binding specificity and accurately predicts binding motifs for an independent set of divergent homeodomains. Similarly, when trained on RNAcompete profiles for diverse RBPs, our model correctly predicts the binding affinities of held-out proteins and identifies key RNA-binding residues, despite the high level of sequence divergence across RBPs. We expect that the method will be broadly applicable to modeling and predicting paired macromolecular interactions in settings where high-throughput affinity data are available.