In vivo interrogation of regulatory genomes reveals extensive quasi-insufficiency in cancer evolution.

In contrast to mono- or biallelic loss of tumor-suppressor function, effects of discrete gene dysregulations, as caused by non-coding (epi)genome alterations, are poorly understood. Here, by perturbing the regulatory genome in mice, we uncover pervasive roles of subtle gene expression variation in cancer evolution. Genome-wide screens characterizing 1,450 tumors revealed that such quasi-insufficiency is extensive across entities and displays diverse context dependencies, such as distinct cell-of-origin associations in T-ALL subtypes. We compile catalogs of non-coding regions linked to quasi-insufficiency, show their enrichment with human cancer risk variants, and provide functional insights by engineering regulatory alterations in mice. As such, kilo-/megabase deletions in a Bcl11b-linked non-coding region triggered aggressive malignancies, with allele-specific tumor spectra reflecting gradual gene dysregulations through modular and cell-type-specific enhancer activities. Our study constitutes a first survey toward a systems-level understanding of quasi-insufficiency in cancer and gives multifaceted insights into tumor evolution and the tissue-specific effects of non-coding mutations.

Development and application of a next-generation sequencing protocol and bioinformatics pipeline for the comprehensive analysis of the canine immunoglobulin repertoire.

Profiling the adaptive immune repertoire using next generation sequencing (NGS) has become common in human medicine, showing promise in characterizing clonal expansion of B cell clones through analysis of B cell receptors (BCRs) in patients with lymphoid malignancies. In contrast, most work evaluating BCR repertoires in dogs has employed traditional PCR-based approaches analyzing the IGH locus only. The objectives of this study were to: (1) describe a novel NGS protocol to evaluate canine BCRs; (2) develop a bioinformatics pipeline for processing canine BCR sequencing data; and (3) apply these methods to derive insights into BCR repertoires of healthy dogs and dogs undergoing treatment for B-cell lymphoma. RNA from peripheral blood mononuclear cells of healthy dogs (n = 25) and dogs newly diagnosed with intermediate-to-large B-cell lymphoma (n = 18) with intent to pursue chemotherapy was isolated, converted into cDNA and sequenced by NGS. The BCR repertoires were identified and quantified using a novel analysis pipeline. The IGK repertoires of the healthy dogs were far less diverse compared to IGL which, as with IGH, was highly diverse. Strong biases at key positions within the CDR3 sequence were identified within the healthy dog BCR repertoire. For a subset of the dogs with B-cell lymphoma, clonal expansion of specific IGH sequences pre-treatment and reduction post-treatment was observed. The degree of expansion and reduction correlated with the clinical outcome in this subset. Future studies employing these techniques may improve disease monitoring, provide earlier recognition of disease progression, and ultimately lead to more targeted therapeutics.

A purine metabolic checkpoint that prevents autoimmunity and autoinflammation.

Still's disease, the paradigm of autoinflammation-cum-autoimmunity, predisposes for a cytokine storm with excessive T lymphocyte activation upon viral infection. Loss of function of the purine nucleoside enzyme FAMIN is the sole known cause for monogenic Still's disease. Here we discovered that a FAMIN-enabled purine metabolon in dendritic cells (DCs) restrains CD4+ and CD8+ T cell priming. DCs with absent FAMIN activity prime for enhanced antigen-specific cytotoxicity, IFNγ secretion, and T cell expansion, resulting in excessive influenza A virus-specific responses. Enhanced priming is already manifest with hypomorphic FAMIN-I254V, for which ∼6% of mankind is homozygous. FAMIN controls membrane trafficking and restrains antigen presentation in an NADH/NAD+-dependent manner by balancing flux through adenine-guanine nucleotide interconversion cycles. FAMIN additionally converts hypoxanthine into inosine, which DCs release to dampen T cell activation. Compromised FAMIN consequently enhances immunosurveillance of syngeneic tumors. FAMIN is a biochemical checkpoint that protects against excessive antiviral T cell responses, autoimmunity, and autoinflammation.

METTL1-mediated m7G modification of Arg-TCT tRNA drives oncogenic transformation.

The emerging "epitranscriptomics" field is providing insights into the biological and pathological roles of different RNA modifications. The RNA methyltransferase METTL1 catalyzes N7-methylguanosine (m7G) modification of tRNAs. Here we find METTL1 is frequently amplified and overexpressed in cancers and is associated with poor patient survival. METTL1 depletion causes decreased abundance of m7G-modified tRNAs and altered cell cycle and inhibits oncogenicity. Conversely, METTL1 overexpression induces oncogenic cell transformation and cancer. Mechanistically, we find increased abundance of m7G-modified tRNAs, in particular Arg-TCT-4-1, and increased translation of mRNAs, including cell cycle regulators that are enriched in the corresponding AGA codon. Accordingly, Arg-TCT expression is elevated in many tumor types and is associated with patient survival, and strikingly, overexpression of this individual tRNA induces oncogenic transformation. Thus, METTL1-mediated tRNA modification drives oncogenic transformation through a remodeling of the mRNA "translatome" to increase expression of growth-promoting proteins and represents a promising anti-cancer target.

Genetic Screens Identify a Context-Specific PI3K/p27Kip1 Node Driving Extrahepatic Biliary Cancer.

Biliary tract cancer ranks among the most lethal human malignancies, representing an unmet clinical need. Its abysmal prognosis is tied to an increasing incidence and a fundamental lack of mechanistic knowledge regarding the molecular basis of the disease. Here, we show that the Pdx1-positive extrahepatic biliary epithelium is highly susceptible toward transformation by activated PIK3CAH1047R but refractory to oncogenic KrasG12D. Using genome-wide transposon screens and genetic loss-of-function experiments, we discover context-dependent genetic interactions that drive extrahepatic cholangiocarcinoma (ECC) and show that PI3K signaling output strength and repression of the tumor suppressor p27Kip1 are critical context-specific determinants of tumor formation. This contrasts with the pancreas, where oncogenic Kras in concert with p53 loss is a key cancer driver. Notably, inactivation of p27Kip1 permits KrasG12D-driven ECC development. These studies provide a mechanistic link between PI3K signaling, tissue-specific tumor suppressor barriers, and ECC pathogenesis, and present a novel genetic model of autochthonous ECC and genes driving this highly lethal tumor subtype. SIGNIFICANCE: We used the first genetically engineered mouse model for extrahepatic bile duct carcinoma to identify cancer genes by genome-wide transposon-based mutagenesis screening. Thereby, we show that PI3K signaling output strength and p27Kip1 function are critical determinants for context-specific ECC formation. This article is highlighted in the In This Issue feature, p. 2945.

miR-200 deficiency promotes lung cancer metastasis by activating Notch signaling in cancer-associated fibroblasts.

Lung adenocarcinoma, the most prevalent lung cancer subtype, is characterized by its high propensity to metastasize. Despite the importance of metastasis in lung cancer mortality, its underlying cellular and molecular mechanisms remain largely elusive. Here, we identified miR-200 miRNAs as potent suppressors for lung adenocarcinoma metastasis. miR-200 expression is specifically repressed in mouse metastatic lung adenocarcinomas, and miR-200 decrease strongly correlates with poor patient survival. Consistently, deletion of mir-200c/141 in the KrasLSL-G12D/+ ; Trp53flox/flox lung adenocarcinoma mouse model significantly promoted metastasis, generating a desmoplastic tumor stroma highly reminiscent of metastatic human lung cancer. miR-200 deficiency in lung cancer cells promotes the proliferation and activation of adjacent cancer-associated fibroblasts (CAFs), which in turn elevates the metastatic potential of cancer cells. miR-200 regulates the functional interaction between cancer cells and CAFs, at least in part, by targeting Notch ligand Jagged1 and Jagged2 in cancer cells and inducing Notch activation in adjacent CAFs. Hence, the interaction between cancer cells and CAFs constitutes an essential mechanism to promote metastatic potential.

Gain-of-function factor H-related 5 protein impairs glomerular complement regulation resulting in kidney damage.

Genetic variation within the factor H-related (FHR) genes is associated with the complement-mediated kidney disease, C3 glomerulopathy (C3G). There is no definitive treatment for C3G, and a significant proportion of patients develop end-stage renal disease. The prototypical example is CFHR5 nephropathy, through which an internal duplication within a single CFHR5 gene generates a mutant FHR5 protein (FHR5mut) that leads to accumulation of complement C3 within glomeruli. To elucidate how abnormal FHR proteins cause C3G, we modeled CFHR5 nephropathy in mice. Animals lacking the murine factor H (FH) and FHR proteins, but coexpressing human FH and FHR5mut (hFH-FHR5mut), developed glomerular C3 deposition, whereas mice coexpressing human FH with the normal FHR5 protein (hFH-FHR5) did not. Like in patients, the FHR5mut had a dominant gain-of-function effect, and when administered in hFH-FHR5 mice, it triggered C3 deposition. Importantly, adeno-associated virus vector-delivered homodimeric mini-FH, a molecule with superior surface C3 binding compared to FH, reduced glomerular C3 deposition in the presence of the FHR5mut. Our data demonstrate that FHR5mut causes C3G by disrupting the homeostatic regulation of complement within the kidney and is directly pathogenic in C3G. These results support the use of FH-derived molecules with enhanced C3 binding for treating C3G associated with abnormal FHR proteins. They also suggest that targeting FHR5 represents a way to treat complement-mediated kidney injury.

CRISPR and transposon in vivo screens for cancer drivers and therapeutic targets.

Human cancers harbor substantial genetic, epigenetic, and transcriptional changes, only some of which drive oncogenesis at certain times during cancer evolution. Identifying the cancer-driver alterations amongst the vast swathes of "passenger" changes still remains a major challenge. Transposon and CRISPR screens in vivo provide complementary methods for achieving this, and each platform has its own advantages. Here, we review recent major technological breakthroughs made with these two approaches and highlight future directions. We discuss how each genetic screening platform can provide unique insight into cancer evolution, including intra-tumoral heterogeneity, metastasis, and immune evasion, presenting transformative opportunities for targeted therapeutic intervention.

Correction to: PiggyBac mutagenesis and exome sequencing identify genetic driver landscapes and potential therapeutic targets of EGFR-mutant gliomas.

An amendment to this paper has been published and can be accessed via the original article.

PiggyBac mutagenesis and exome sequencing identify genetic driver landscapes and potential therapeutic targets of EGFR-mutant gliomas.

BACKGROUND: Glioma is the most common intrinsic brain tumor and also occurs in the spinal cord. Activating EGFR mutations are common in IDH1 wild-type gliomas. However, the cooperative partners of EGFR driving gliomagenesis remain poorly understood. RESULTS: We explore EGFR-mutant glioma evolution in conditional mutant mice by whole-exome sequencing, transposon mutagenesis forward genetic screening, and transcriptomics. We show mutant EGFR is sufficient to initiate gliomagenesis in vivo, both in the brain and spinal cord. We identify significantly recurrent somatic alterations in these gliomas including mutant EGFR amplifications and Sub1, Trp53, and Tead2 loss-of-function mutations. Comprehensive functional characterization of 96 gliomas by genome-wide piggyBac insertional mutagenesis in vivo identifies 281 known and novel EGFR-cooperating driver genes, including Cdkn2a, Nf1, Spred1, and Nav3. Transcriptomics confirms transposon-mediated effects on expression of these genes. We validate the clinical relevance of new putative tumor suppressors by showing these are frequently altered in patients' gliomas, with prognostic implications. We discover shared and distinct driver mutations in brain and spinal gliomas and confirm in vivo differential tumor suppressive effects of Pten between these tumors. Functional validation with CRISPR-Cas9-induced mutations in novel genes Tead2, Spred1, and Nav3 demonstrates heightened EGFRvIII-glioma cell proliferation. Chemogenomic analysis of mutated glioma genes reveals potential drug targets, with several investigational drugs showing efficacy in vitro. CONCLUSION: Our work elucidates functional driver landscapes of EGFR-mutant gliomas, uncovering potential therapeutic strategies, and provides new tools for functional interrogation of gliomagenesis.

FAMIN Is a Multifunctional Purine Enzyme Enabling the Purine Nucleotide Cycle.

Mutations in FAMIN cause arthritis and inflammatory bowel disease in early childhood, and a common genetic variant increases the risk for Crohn's disease and leprosy. We developed an unbiased liquid chromatography-mass spectrometry screen for enzymatic activity of this orphan protein. We report that FAMIN phosphorolytically cleaves adenosine into adenine and ribose-1-phosphate. Such activity was considered absent from eukaryotic metabolism. FAMIN and its prokaryotic orthologs additionally have adenosine deaminase, purine nucleoside phosphorylase, and S-methyl-5'-thioadenosine phosphorylase activity, hence, combine activities of the namesake enzymes of central purine metabolism. FAMIN enables in macrophages a purine nucleotide cycle (PNC) between adenosine and inosine monophosphate and adenylosuccinate, which consumes aspartate and releases fumarate in a manner involving fatty acid oxidation and ATP-citrate lyase activity. This macrophage PNC synchronizes mitochondrial activity with glycolysis by balancing electron transfer to mitochondria, thereby supporting glycolytic activity and promoting oxidative phosphorylation and mitochondrial H+ and phosphate recycling.

A Genome-Wide Knockout Screen in Human Macrophages Identified Host Factors Modulating Salmonella Infection.

A genome-scale CRISPR knockout library screen of THP-1 human macrophages was performed to identify loss-of-function mutations conferring resistance to Salmonella uptake. The screen identified 183 candidate genes, from which 14 representative genes involved in actin dynamics (ACTR3, ARPC4, CAPZB, TOR3A, CYFIP2, CTTN, and NHLRC2), glycosaminoglycan metabolism (B3GNT1), receptor signaling (PDGFB and CD27), lipid raft formation (CLTCL1), calcium transport (ATP2A2 and ITPR3), and cholesterol metabolism (HMGCR) were analyzed further. For some of these pathways, known chemical inhibitors could replicate the Salmonella resistance phenotype, indicating their potential as targets for host-directed therapy. The screen indicated a role for the relatively uncharacterized gene NHLRC2 in both Salmonella invasion and macrophage differentiation. Upon differentiation, NHLRC2 mutant macrophages were hyperinflammatory and did not exhibit characteristics typical of macrophages, including atypical morphology and inability to interact and phagocytose bacteria/particles. Immunoprecipitation confirmed an interaction of NHLRC2 with FRYL, EIF2AK2, and KLHL13.IMPORTANCESalmonella exploits macrophages to gain access to the lymphatic system and bloodstream to lead to local and potentially systemic infections. With an increasing number of antibiotic-resistant isolates identified in humans, Salmonella infections have become major threats to public health. Therefore, there is an urgent need to identify alternative approaches to anti-infective therapy, including host-directed therapies. In this study, we used a simple genome-wide screen to identify 183 candidate host factors in macrophages that can confer resistance to Salmonella infection. These factors may be potential therapeutic targets against Salmonella infections.

NPM1c alters FLT3-D835Y localization and signaling in acute myeloid leukemia.

Activating mutations in FMS-like tyrosine kinase receptor-3 (FLT3) and Nucleophosmin-1 (NPM1) are most frequent alterations in acute myeloid leukemia (AML), and are often coincidental. The mutational status of NPM1 has strong prognostic relevance to patients with point mutations of the FLT3 tyrosine kinase domain (TKD), but the biological mechanism underlying this effect remains unclear. In the present study, we investigated the effect of the coincidence of NPM1c and FLT3-TKD. Although expression of FLT3-TKD is not sufficient to induce a disease in mice, coexpression with NPM1c rapidly leads to an aggressive myeloproliferative disease in mice with a latency of 31.5 days. Mechanistically, we could show that FLT3-TKD is able to activate the downstream effector molecule signal transducer and activator of transcription 5 (STAT5) exclusively in the presence of mutated NPM1c. Moreover, NPM1c alters the cellular localization of FLT3-TKD from the cell surface to the endoplasmic reticulum, which might thereby lead to the aberrant STAT5 activation. Importantly, aberrant STAT5 activation occurs not only in primary murine cells but also in patients with AML with combined FLT3-TKD and NPM1c mutations. Thus, our data indicate a new mechanism, how NPM1c mislocalizes FLT3-TKD and changes its signal transduction ability.

PiggyBac transposon tools for recessive screening identify B-cell lymphoma drivers in mice.

B-cell lymphoma (BCL) is the most common hematologic malignancy. While sequencing studies gave insights into BCL genetics, identification of non-mutated cancer genes remains challenging. Here, we describe PiggyBac transposon tools and mouse models for recessive screening and show their application to study clonal B-cell lymphomagenesis. In a genome-wide screen, we discover BCL genes related to diverse molecular processes, including signaling, transcriptional regulation, chromatin regulation, or RNA metabolism. Cross-species analyses show the efficiency of the screen to pinpoint human cancer drivers altered by non-genetic mechanisms, including clinically relevant genes dysregulated epigenetically, transcriptionally, or post-transcriptionally in human BCL. We also describe a CRISPR/Cas9-based in vivo platform for BCL functional genomics, and validate discovered genes, such as Rfx7, a transcription factor, and Phip, a chromatin regulator, which suppress lymphomagenesis in mice. Our study gives comprehensive insights into the molecular landscapes of BCL and underlines the power of genome-scale screening to inform biology.

ATM orchestrates the DNA-damage response to counter toxic non-homologous end-joining at broken replication forks.

Mutations in the ATM tumor suppressor gene confer hypersensitivity to DNA-damaging chemotherapeutic agents. To explore genetic resistance mechanisms, we performed genome-wide CRISPR-Cas9 screens in cells treated with the DNA topoisomerase I inhibitor topotecan. Thus, we here establish that inactivating terminal components of the non-homologous end-joining (NHEJ) machinery or of the BRCA1-A complex specifically confer topotecan resistance to ATM-deficient cells. We show that hypersensitivity of ATM-mutant cells to topotecan or the poly-(ADP-ribose) polymerase (PARP) inhibitor olaparib reflects delayed engagement of homologous recombination at DNA-replication-fork associated single-ended double-strand breaks (DSBs), allowing some to be subject to toxic NHEJ. Preventing DSB ligation by NHEJ, or enhancing homologous recombination by BRCA1-A complex disruption, suppresses this toxicity, highlighting a crucial role for ATM in preventing toxic LIG4-mediated chromosome fusions. Notably, suppressor mutations in ATM-mutant backgrounds are different to those in BRCA1-mutant scenarios, suggesting new opportunities for patient stratification and additional therapeutic vulnerabilities for clinical exploitation.

PiggyBac Transposon-Based Insertional Mutagenesis in Mice.

While sequencing and array-based studies are creating catalogues of genetic alterations in cancer, discriminating cancer drivers among the large sets of epigenetically, transcriptionally or posttranslationally dysregulated genes remains a challenge. Transposon-based genetic screening in mice has proven to be a powerful approach to address this challenge. Insertional mutagenesis directly flags biologically relevant genes and, combined with the transposon's unique molecular fingerprint, facilitates the recovery of insertion sites. We have generated transgenic mouse lines harboring different versions of PiggyBac-based oncogenic transposons, which in conjunction with PiggyBac transposase mice can be used for whole-body or tissue-specific insertional mutagenesis screens. We have also developed QiSeq, a method for (semi-)quantitative transposon insertion site sequencing, which overcomes biasing limitations of previous library preparation methods. QiSeq can be used in multiplexed high-throughput formats for candidate cancer gene discovery and gives insights into the clonal distribution of insertions for the study of genetic tumor evolution.

SRPK1 maintains acute myeloid leukemia through effects on isoform usage of epigenetic regulators including BRD4.

We recently identified the splicing kinase gene SRPK1 as a genetic vulnerability of acute myeloid leukemia (AML). Here, we show that genetic or pharmacological inhibition of SRPK1 leads to cell cycle arrest, leukemic cell differentiation and prolonged survival of mice transplanted with MLL-rearranged AML. RNA-seq analysis demonstrates that SRPK1 inhibition leads to altered isoform levels of many genes including several with established roles in leukemogenesis such as MYB, BRD4 and MED24. We focus on BRD4 as its main isoforms have distinct molecular properties and find that SRPK1 inhibition produces a significant switch from the short to the long isoform at the mRNA and protein levels. This was associated with BRD4 eviction from genomic loci involved in leukemogenesis including BCL2 and MYC. We go on to show that this switch mediates at least part of the anti-leukemic effects of SRPK1 inhibition. Our findings reveal that SRPK1 represents a plausible new therapeutic target against AML.

Repair of double-strand breaks induced by CRISPR-Cas9 leads to large deletions and complex rearrangements.

CRISPR-Cas9 is poised to become the gene editing tool of choice in clinical contexts. Thus far, exploration of Cas9-induced genetic alterations has been limited to the immediate vicinity of the target site and distal off-target sequences, leading to the conclusion that CRISPR-Cas9 was reasonably specific. Here we report significant on-target mutagenesis, such as large deletions and more complex genomic rearrangements at the targeted sites in mouse embryonic stem cells, mouse hematopoietic progenitors and a human differentiated cell line. Using long-read sequencing and long-range PCR genotyping, we show that DNA breaks introduced by single-guide RNA/Cas9 frequently resolved into deletions extending over many kilobases. Furthermore, lesions distal to the cut site and crossover events were identified. The observed genomic damage in mitotically active cells caused by CRISPR-Cas9 editing may have pathogenic consequences.

Evolutionary routes and KRAS dosage define pancreatic cancer phenotypes.

The poor correlation of mutational landscapes with phenotypes limits our understanding of the pathogenesis and metastasis of pancreatic ductal adenocarcinoma (PDAC). Here we show that oncogenic dosage-variation has a critical role in PDAC biology and phenotypic diversification. We find an increase in gene dosage of mutant KRAS in human PDAC precursors, which drives both early tumorigenesis and metastasis and thus rationalizes early PDAC dissemination. To overcome the limitations posed to gene dosage studies by the stromal richness of PDAC, we have developed large cell culture resources of metastatic mouse PDAC. Integration of cell culture genomes, transcriptomes and tumour phenotypes with functional studies and human data reveals additional widespread effects of oncogenic dosage variation on cell morphology and plasticity, histopathology and clinical outcome, with the highest KrasMUT levels underlying aggressive undifferentiated phenotypes. We also identify alternative oncogenic gains (Myc, Yap1 or Nfkb2), which collaborate with heterozygous KrasMUT in driving tumorigenesis, but have lower metastatic potential. Mechanistically, different oncogenic gains and dosages evolve along distinct evolutionary routes, licensed by defined allelic states and/or combinations of hallmark tumour suppressor alterations (Cdkn2a, Trp53, Tgfß-pathway). Thus, evolutionary constraints and contingencies direct oncogenic dosage gain and variation along defined routes to drive the early progression of PDAC and shape its downstream biology. Our study uncovers universal principles of Ras-driven oncogenesis that have potential relevance beyond pancreatic cancer.

Comprehensive annotation and evolutionary insights into the canine (Canis lupus familiaris) antigen receptor loci.

Dogs are an excellent model for human disease. For example, the treatment of canine lymphoma has been predictive of the human response to that treatment. However, an incomplete picture of canine (Canis lupus familiaris) immunoglobulin (IG) and T cell receptor (TR)-or antigen receptor (AR)-gene loci has restricted their utility. This work advances the annotation of the canine AR loci and looks into breed-specific features of the loci. Bioinformatic analysis of unbiased RNA sequence data was used to complete the annotation of the canine AR genes. This annotation was used to query 107 whole genome sequences from 19 breeds and identified over 5500 alleles across the 550 genes of the seven AR loci: the IG heavy, kappa, and lambda loci; and the TR alpha, beta, gamma, and delta loci. Of note was the discovery that half of the IGK variable (V) genes were located downstream of, and inverted with respect to, the rest of the locus. Analysis of the germline sequences of all the AR V genes identified greater conservation between dog and human than mouse with either. This work brings our understanding of the genetic diversity and expression of AR in dogs to the same completeness as that of mice and men, making it the third species to have all AR loci comprehensively and accurately annotated. The large number of germline sequences serves as a reference for future studies, and has allowed statistically powerful conclusions to be drawn on the pressures that have shaped these loci.