Leukemias exhibit a dysregulated developmental program mediated through both genetic and epigenetic mechanisms. Although IKZF2 is expressed in hematopoietic stem cells (HSCs), we found that it is dispensable for mouse and human HSC function. In contrast to its role as a tumor suppressor in hypodiploid B-acute lymphoblastic leukemia, we found that IKZF2 is required for myeloid leukemia. IKZF2 is highly expressed in leukemic stem cells (LSCs), and its deficiency results in defective LSC function. IKZF2 depletion in acute myeloid leukemia (AML) cells reduced colony formation, increased differentiation and apoptosis, and delayed leukemogenesis. Gene expression, chromatin accessibility, and direct IKZF2 binding in MLL-AF9 LSCs demonstrate that IKZF2 regulates a HOXA9 self-renewal gene expression program and inhibits a C/EBP-driven differentiation program. Ectopic HOXA9 expression and CEBPE depletion rescued the effects of IKZF2 depletion. Thus, our study shows that IKZF2 regulates the AML LSC program and provides a rationale to therapeutically target IKZF2 in myeloid leukemia.

Increasing evidence indicates CD4+ T cells can recognize cancer-specific antigens and control tumor growth. However, it remains difficult to predict the antigens that will be presented by human leukocyte antigen class II molecules (HLA-II), hindering efforts to optimally target them therapeutically. Obstacles include inaccurate peptide-binding prediction and unsolved complexities of the HLA-II pathway. To address these challenges, we developed an improved technology for discovering HLA-II binding motifs and conducted a comprehensive analysis of tumor ligandomes to learn processing rules relevant in the tumor microenvironment. We profiled >40 HLA-II alleles and showed that binding motifs were highly sensitive to HLA-DM, a peptide-loading chaperone. We also revealed that intratumoral HLA-II presentation was dominated by professional antigen-presenting cells (APCs) rather than cancer cells. Integrating these observations, we developed algorithms that accurately predicted APC ligandomes, including peptides from phagocytosed cancer cells. These tools and biological insights will enable improved HLA-II-directed cancer therapies.

Recent human decedent model studies1,2 and compassionate xenograft use3 have explored the promise of porcine organs for human transplantation. To proceed to human studies, a clinically ready porcine donor must be engineered and its xenograft successfully tested in nonhuman primates. Here we describe the design, creation and long-term life-supporting function of kidney grafts from a genetically engineered porcine donor transplanted into a cynomolgus monkey model. The porcine donor was engineered to carry 69 genomic edits, eliminating glycan antigens, overexpressing human transgenes and inactivating porcine endogenous retroviruses. In vitro functional analyses showed that the edited kidney endothelial cells modulated inflammation to an extent that was indistinguishable from that of human endothelial cells, suggesting that these edited cells acquired a high level of human immune compatibility. When transplanted into cynomolgus monkeys, the kidneys with three glycan antigen knockouts alone experienced poor graft survival, whereas those with glycan antigen knockouts and human transgene expression demonstrated significantly longer survival time, suggesting the benefit of human transgene expression in vivo. These results show that preclinical studies of renal xenotransplantation could be successfully conducted in nonhuman primates and bring us closer to clinical trials of genetically engineered porcine renal grafts.

Mutations in the transcription factor FOXA1 define a unique subset of prostate cancers but the functional consequences of these mutations and whether they confer gain or loss of function is unknown1-9. Here, by annotating the landscape of FOXA1 mutations from 3,086 human prostate cancers, we define two hotspots in the forkhead domain: Wing2 (around 50% of all mutations) and the highly conserved DNA-contact residue R219 (around 5% of all mutations). Wing2 mutations are detected in adenocarcinomas at all stages, whereas R219 mutations are enriched in metastatic tumours with neuroendocrine histology. Interrogation of the biological properties of wild-type FOXA1 and fourteen FOXA1 mutants reveals gain of function in mouse prostate organoid proliferation assays. Twelve of these mutants, as well as wild-type FOXA1, promoted an exaggerated pro-luminal differentiation program, whereas two different R219 mutants blocked luminal differentiation and activated a mesenchymal and neuroendocrine transcriptional program. Assay for transposase-accessible chromatin using sequencing (ATAC-seq) of wild-type FOXA1 and representative Wing2 and R219 mutants revealed marked, mutant-specific changes in open chromatin at thousands of genomic loci and exposed sites of FOXA1 binding and associated increases in gene expression. Of note, ATAC-seq peaks in cells expressing R219 mutants lacked the canonical core FOXA1-binding motifs (GTAAAC/T) but were enriched for a related, non-canonical motif (GTAAAG/A), which was preferentially activated by R219-mutant FOXA1 in reporter assays. Thus, FOXA1 mutations alter its pioneering function and perturb normal luminal epithelial differentiation programs, providing further support for the role of lineage plasticity in cancer progression.

Genomic instability is a fundamental feature of human cancer often resulting from impaired genome maintenance. In prostate cancer, structural genomic rearrangements are a common mechanism driving tumorigenesis. However, somatic alterations predisposing to chromosomal rearrangements in prostate cancer remain largely undefined. Here, we show that SPOP, the most commonly mutated gene in primary prostate cancer modulates DNA double strand break (DSB) repair, and that SPOP mutation is associated with genomic instability. In vivo, SPOP mutation results in a transcriptional response consistent with BRCA1 inactivation resulting in impaired homology-directed repair (HDR) of DSB. Furthermore, we found that SPOP mutation sensitizes to DNA damaging therapeutic agents such as PARP inhibitors. These results implicate SPOP as a novel participant in DSB repair, suggest that SPOP mutation drives prostate tumorigenesis in part through genomic instability, and indicate that mutant SPOP may increase response to DNA-damaging therapeutics.

The initial next-generation sequencing technologies produced reads of 25 or 36 bp, and only from a single-end of the library sequence. Currently, it is possible to reliably produce 300 bp paired-end sequences for RNA expression analysis. While read lengths have consistently increased, people have assumed that longer reads are more informative and that paired-end reads produce better results than single-end reads. We used paired-end 101 bp reads and trimmed them to simulate different read lengths, and also separated the pairs to produce single-end reads. For each read length and paired status, we evaluated differential expression levels between two standard samples and compared the results to those obtained by qPCR.

The panoply of microorganisms and other species present in our environment influence human health and disease, especially in cities, but have not been profiled with metagenomics at a city-wide scale. We sequenced DNA from surfaces across the entire New York City (NYC) subway system, the Gowanus Canal, and public parks. Nearly half of the DNA (48%) does not match any known organism; identified organisms spanned 1,688 bacterial, viral, archaeal, and eukaryotic taxa, which were enriched for harmless genera associated with skin (e.g., Acinetobacter). Predicted ancestry of human DNA left on subway surfaces can recapitulate U.S. Census demographic data, and bacterial signatures can reveal a station's history, such as marine-associated bacteria in a hurricane-flooded station. Some evidence of pathogens was found (Bacillus anthracis), but a lack of reported cases in NYC suggests that the pathogens represent a normal, urban microbiome. This baseline metagenomic map of NYC could help long-term disease surveillance, bioterrorism threat mitigation, and health management in the built environment of cities.

Malignancy can be suppressed by the immune system in a process termed immunosurveillance. However, to what extent immunosurveillance occurs in spontaneous cancers and the composition of participating cell types remains obscure. Here, we show that cell transformation triggers a tissue-resident lymphocyte response in oncogene-induced murine cancer models. Non-circulating cytotoxic lymphocytes, derived from innate, T cell receptor (TCR)αβ, and TCRγδ lineages, expand in early tumors. Characterized by high expression of NK1.1, CD49a, and CD103, these cells share a gene-expression signature distinct from those of conventional NK cells, T cells, and invariant NKT cells. Generation of these lymphocytes is dependent on the cytokine IL-15, but not the transcription factor Nfil3 that is required for the differentiation of tumor-infiltrating NK cells, and IL-15 deficiency, but not Nfil3 deficiency, results in accelerated tumor growth. These findings reveal a tumor-elicited immunosurveillance mechanism that engages unconventional type-1-like innate lymphoid cells and type 1 innate-like T cells.

The identity of the RNA-binding proteins (RBPs) that govern cancer stem cells remains poorly characterized. The MSI2 RBP is a central regulator of translation of cancer stem cell programs. Through proteomic analysis of the MSI2-interacting RBP network and functional shRNA screening, we identified 24 genes required for in vivo leukemia. Syncrip was the most differentially required gene between normal and myeloid leukemia cells. SYNCRIP depletion increased apoptosis and differentiation while delaying leukemogenesis. Gene expression profiling of SYNCRIP-depleted cells demonstrated a loss of the MLL and HOXA9 leukemia stem cell program. SYNCRIP and MSI2 interact indirectly though shared mRNA targets. SYNCRIP maintains HOXA9 translation, and MSI2 or HOXA9 overexpression rescued the effects of SYNCRIP depletion. Altogether, our data identify SYNCRIP as a new RBP that controls the myeloid leukemia stem cell program. We propose that targeting these RBP complexes might provide a novel therapeutic strategy in leukemia.

Thyroid cancer is common, yet the sequence of alterations that promote tumor formation are incompletely understood. Here, we describe a novel model of thyroid carcinoma in zebrafish that reveals temporal changes due to BRAFV600E. Through the use of real-time in vivo imaging, we observe disruption in thyroid follicle structure that occurs early in thyroid development. Combinatorial treatment using BRAF and MEK inhibitors reversed the developmental effects induced by BRAFV600E. Adult zebrafish expressing BRAFV600E in thyrocytes developed invasive carcinoma. We identified a gene expression signature from zebrafish thyroid cancer that is predictive of disease-free survival in patients with papillary thyroid cancer. Gene expression studies nominated TWIST2 as a key effector downstream of BRAF. Using CRISPR/Cas9 to genetically inactivate a TWIST2 orthologue, we suppressed the effects of BRAFV600E and restored thyroid morphology and hormone synthesis. These data suggest that expression of TWIST2 plays a role in an early step of BRAFV600E-mediated transformation.

The Hippo pathway is an important regulator of organ size and tumorigenesis. It is unclear, however, how Hippo signalling provides the cellular building blocks required for rapid growth. Here, we demonstrate that transgenic zebrafish expressing an activated form of the Hippo pathway effector Yap1 (also known as YAP) develop enlarged livers and are prone to liver tumour formation. Transcriptomic and metabolomic profiling identify that Yap1 reprograms glutamine metabolism. Yap1 directly enhances glutamine synthetase (glul) expression and activity, elevating steady-state levels of glutamine and enhancing the relative isotopic enrichment of nitrogen during de novo purine and pyrimidine biosynthesis. Genetic or pharmacological inhibition of GLUL diminishes the isotopic enrichment of nitrogen into nucleotides, suppressing hepatomegaly and the growth of liver cancer cells. Consequently, Yap-driven liver growth is susceptible to nucleotide inhibition. Together, our findings demonstrate that Yap1 integrates the anabolic demands of tissue growth during development and tumorigenesis by reprogramming nitrogen metabolism to stimulate nucleotide biosynthesis.

The two major melanoma histologic subtypes, superficial spreading and nodular melanomas, differ in their speed of dermal invasion but converge biologically once they invade and metastasize. Herein, we tested the hypothesis that distinct molecular alterations arising in primary melanoma cells might persist as these tumors progress to invasion and metastasis. Ribosomal protein S6 kinase, 90 kDa, polypeptide 1 (RSK1; official name RPS6KA1) was significantly hyperactivated in human melanoma lines and metastatic tissues derived from nodular compared with superficial spreading melanoma. RSK1 was constitutively phosphorylated at Ser-380 in nodular but not superficial spreading melanoma and did not directly correlate with BRAF or MEK activation. Nodular melanoma cells were more sensitive to RSK1 inhibition using siRNA and the pharmacological inhibitor BI-D1870 compared with superficial spreading cells. Gene expression microarray analyses revealed that RSK1 orchestrated a program of gene expression that promoted cell motility and invasion. Differential overexpression of the prometastatic matrix metalloproteinase 8 and tissue inhibitor of metalloproteinases 1 in metastatic nodular compared with metastatic superficial spreading melanoma was observed. Finally, using an in vivo zebrafish model, constitutive RSK1 activation increased melanoma invasion. Together, these data reveal a novel role for activated RSK1 in the progression of nodular melanoma and suggest that melanoma originating from different histologic subtypes may be biologically distinct and that these differences are maintained as the tumors invade and metastasize.

We present GuideScan software for the design of CRISPR guide RNA libraries that can be used to edit coding and noncoding genomic regions. GuideScan produces high-density sets of guide RNAs (gRNAs) for single- and paired-gRNA genome-wide screens. We also show that the trie data structure of GuideScan enables the design of gRNAs that are more specific than those designed by existing tools.

Aerobic glycolysis-the Warburg effect-converts glucose to lactate via the enzyme lactate dehydrogenase A (LDHA) and is a metabolic feature of effector T cells. Cells generate ATP through various mechanisms and Warburg metabolism is comparatively an energy-inefficient glucose catabolism pathway. Here, we examined the effect of ATP generated via aerobic glycolysis in antigen-driven T cell responses. Cd4CreLdhafl/fl mice were resistant to Th17-cell-mediated experimental autoimmune encephalomyelitis and exhibited defective T cell activation, migration, proliferation, and differentiation. LDHA deficiency crippled cellular redox balance and inhibited ATP production, diminishing PI3K-dependent activation of Akt kinase and thereby phosphorylation-mediated inhibition of Foxo1, a transcriptional repressor of T cell activation programs. Th17-cell-specific expression of an Akt-insensitive Foxo1 recapitulated the defects seen in Cd4CreLdhafl/fl mice. Induction of LDHA required PI3K signaling and LDHA deficiency impaired PI3K-catalyzed PIP3 generation. Thus, Warburg metabolism augments glycolytic ATP production, fueling a PI3K-centered positive feedback regulatory circuit that drives effector T cell responses.

The PI3K pathway integrates extracellular stimuli to phosphorylate effectors such as AKT and serum-and-glucocorticoid-regulated kinase (SGK1). We have previously reported that the PI3K pathway regulates estrogen receptor (ER)-dependent transcription in breast cancer through the phosphorylation of the lysine methyltransferase KMT2D by AKT. Here, we show that PI3Kα inhibition, via a negative-feedback loop, activates SGK1 to promote chromatin-based regulation of ER-dependent transcription. PI3K/AKT inhibitors activate ER, which promotes SGK1 transcription through direct binding to its promoter. Elevated SGK1, in turn, phosphorylates KMT2D, suppressing its function, leading to a loss of methylation of lysine 4 on histone H3 (H3K4) and a repressive chromatin state at ER loci to attenuate ER activity. Thus, SGK1 regulates the chromatin landscape and ER-dependent transcription via the direct phosphorylation of KMT2D. These findings reveal an ER-SGK1-KMT2D signaling circuit aimed to attenuate ER response through a role for SGK1 to program chromatin and ER transcriptional output.

Myelodysplastic syndrome (MDS) is a heterogeneous group of clonal hematopoietic disorders characterized by ineffective hematopoiesis, cytopenias, and dysplasia. The gene encoding ten-eleven translocation 2 (tet2), a dioxygenase enzyme that catalyzes the conversion of 5-methylcytosine (5mC) to 5-hydroxymethylcytosine, is a recurrently mutated tumor suppressor gene in MDS and other myeloid malignancies. Previously, we reported a stable zebrafish line with a loss-of-function mutation in the tet2 gene. The tet2m/m-mutant zebrafish developed a pre-MDS state with kidney marrow dysplasia, but normal circulating blood counts by 11 months of age and accompanying anemia, signifying the onset of MDS, by 24 months of age.