Novel stem cell technologies are powerful tools to understand the impact of human factors on Plasmodium falciparum malaria.

Plasmodium falciparum parasites have a complex life cycle, but the most clinically relevant stage of the disease is the invasion of erythrocytes and the proliferation of the parasite in the blood. The influence of human genetic traits on malaria has been known for a long time, however understanding the role of the proteins involved is hampered by the anuclear nature of erythrocytes that makes them inaccessible to genetic tools. Here we overcome this limitation using stem cells to generate erythroid cells with an in-vitro differentiation protocol and assess parasite invasion with an adaptation of flow cytometry to detect parasite hemozoin. We combine this strategy with reprogramming of patient cells to Induced Pluripotent Stem Cells and genome editing to understand the role of key genes and human traits in malaria infection. We show that deletion of basigin ablates invasion while deletion of ATP2B4 has a minor effect and that erythroid cells from reprogrammed patient-derived HbBart α-thalassemia samples poorly support infection. The possibility to obtain patient-secific and genetically modifed erythoid cells offers an unparalleled opportunity to study the role of human genes and polymorphisms in malaria allowing preservation of the genomic background to demonstrate their function and understand their mechanisms.

Genetic associations at regulatory phenotypes improve fine-mapping of causal variants for 12 immune-mediated diseases.

The resolution of causal genetic variants informs understanding of disease biology. We used regulatory quantitative trait loci (QTLs) from the BLUEPRINT, GTEx and eQTLGen projects to fine-map putative causal variants for 12 immune-mediated diseases. We identify 340 unique loci that colocalize with high posterior probability (≥98%) with regulatory QTLs and apply Bayesian frameworks to fine-map associations at each locus. We show that fine-mapping credible sets derived from regulatory QTLs are smaller compared to disease summary statistics. Further, they are enriched for more functionally interpretable candidate causal variants and for putatively causal insertion/deletion (INDEL) polymorphisms. Finally, we use massively parallel reporter assays to evaluate candidate causal variants at the ITGA4 locus associated with inflammatory bowel disease. Overall, our findings suggest that fine-mapping applied to disease-colocalizing regulatory QTLs can enhance the discovery of putative causal disease variants and enhance insights into the underlying causal genes and molecular mechanisms.

SETBP1 overexpression acts in the place of class-defining mutations to drive FLT3-ITD-mutant AML.

Advances in cancer genomics have revealed genomic classes of acute myeloid leukemia (AML) characterized by class-defining mutations, such as chimeric fusion genes or in genes such as NPM1, MLL, and CEBPA. These class-defining mutations frequently synergize with internal tandem duplications in FLT3 (FLT3-ITDs) to drive leukemogenesis. However, ∼20% of FLT3-ITD-positive AMLs bare no class-defining mutations, and mechanisms of leukemic transformation in these cases are unknown. To identify pathways that drive FLT3-ITD mutant AML in the absence of class-defining mutations, we performed an insertional mutagenesis (IM) screening in Flt3-ITD mice, using Sleeping Beauty transposons. All mice developed acute leukemia (predominantly AML) after a median of 73 days. Analysis of transposon insertions in 38 samples from Flt3-ITD/IM leukemic mice identified recurrent integrations at 22 loci, including Setbp1 (20/38), Ets1 (11/38), Ash1l (8/38), Notch1 (8/38), Erg (7/38), and Runx1 (5/38). Insertions at Setbp1 led exclusively to AML and activated a transcriptional program similar, but not identical, to those of NPM1-mutant and MLL-rearranged AMLs. Guide RNA targeting of Setbp1 was highly detrimental to Flt3ITD/+/Setbp1IM+, but not to Flt3ITD/+/Npm1cA/+, AMLs. Also, analysis of RNA-sequencing data from hundreds of human AMLs revealed that SETBP1 expression is significantly higher in FLT3-ITD AMLs lacking class-defining mutations. These findings propose that SETBP1 overexpression collaborates with FLT3-ITD to drive a subtype of human AML. To identify genetic vulnerabilities of these AMLs, we performed genome-wide CRISPR-Cas9 screening in Flt3ITD/+/Setbp1IM+ AMLs and identified potential therapeutic targets, including Kdm1a, Brd3, Ezh2, and Hmgcr. Our study gives new insights into epigenetic pathways that can drive AMLs lacking class-defining mutations and proposes therapeutic approaches against such cases.

Genetic perturbation of PU.1 binding and chromatin looping at neutrophil enhancers associates with autoimmune disease.

Neutrophils play fundamental roles in innate immune response, shape adaptive immunity, and are a potentially causal cell type underpinning genetic associations with immune system traits and diseases. Here, we profile the binding of myeloid master regulator PU.1 in primary neutrophils across nearly a hundred volunteers. We show that variants associated with differential PU.1 binding underlie genetically-driven differences in cell count and susceptibility to autoimmune and inflammatory diseases. We integrate these results with other multi-individual genomic readouts, revealing coordinated effects of PU.1 binding variants on the local chromatin state, enhancer-promoter contacts and downstream gene expression, and providing a functional interpretation for 27 genes underlying immune traits. Collectively, these results demonstrate the functional role of PU.1 and its target enhancers in neutrophil transcriptional control and immune disease susceptibility.

The Polygenic and Monogenic Basis of Blood Traits and Diseases.

Blood cells play essential roles in human health, underpinning physiological processes such as immunity, oxygen transport, and clotting, which when perturbed cause a significant global health burden. Here we integrate data from UK Biobank and a large-scale international collaborative effort, including data for 563,085 European ancestry participants, and discover 5,106 new genetic variants independently associated with 29 blood cell phenotypes covering a range of variation impacting hematopoiesis. We holistically characterize the genetic architecture of hematopoiesis, assess the relevance of the omnigenic model to blood cell phenotypes, delineate relevant hematopoietic cell states influenced by regulatory genetic variants and gene networks, identify novel splice-altering variants mediating the associations, and assess the polygenic prediction potential for blood traits and clinical disorders at the interface of complex and Mendelian genetics. These results show the power of large-scale blood cell trait GWAS to interrogate clinically meaningful variants across a wide allelic spectrum of human variation.

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.

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.

Genome-Scale CRISPRa Screen Identifies Novel Factors for Cellular Reprogramming.

Primed epiblast stem cells (EpiSCs) can be reverted to a pluripotent embryonic stem cell (ESC)-like state by expression of single reprogramming factor. We used CRISPR activation to perform a genome-scale, reprogramming screen in EpiSCs and identified 142 candidate genes. Our screen validated a total of 50 genes, previously not known to contribute to reprogramming, of which we chose Sall1 for further investigation. We show that Sall1 augments reprogramming of mouse EpiSCs and embryonic fibroblasts and that these induced pluripotent stem cells are indeed fully pluripotent including formation of chimeric mice. We also demonstrate that Sall1 synergizes with Nanog in reprogramming and that overexpression in ESCs delays their conversion back to EpiSCs. Lastly, using RNA sequencing, we identify and validate Klf5 and Fam189a2 as new downstream targets of Sall1 and Nanog. In summary, our work demonstrates the power of using CRISPR technology in understanding molecular mechanisms that mediate complex cellular processes such as reprogramming.

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.

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.

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.

Promoter-bound METTL3 maintains myeloid leukaemia by m6A-dependent translation control.

N6-methyladenosine (m6A) is an abundant internal RNA modification in both coding and non-coding RNAs that is catalysed by the METTL3-METTL14 methyltransferase complex. However, the specific role of these enzymes in cancer is still largely unknown. Here we define a pathway that is specific for METTL3 and is implicated in the maintenance of a leukaemic state. We identify METTL3 as an essential gene for growth of acute myeloid leukaemia cells in two distinct genetic screens. Downregulation of METTL3 results in cell cycle arrest, differentiation of leukaemic cells and failure to establish leukaemia in immunodeficient mice. We show that METTL3, independently of METTL14, associates with chromatin and localizes to the transcriptional start sites of active genes. The vast majority of these genes have the CAATT-box binding protein CEBPZ present at the transcriptional start site, and this is required for recruitment of METTL3 to chromatin. Promoter-bound METTL3 induces m6A modification within the coding region of the associated mRNA transcript, and enhances its translation by relieving ribosome stalling. We show that genes regulated by METTL3 in this way are necessary for acute myeloid leukaemia. Together, these data define METTL3 as a regulator of a chromatin-based pathway that is necessary for maintenance of the leukaemic state and identify this enzyme as a potential therapeutic target for acute myeloid leukaemia.

A single-copy Sleeping Beauty transposon mutagenesis screen identifies new PTEN-cooperating tumor suppressor genes.

The overwhelming number of genetic alterations identified through cancer genome sequencing requires complementary approaches to interpret their significance and interactions. Here we developed a novel whole-body insertional mutagenesis screen in mice, which was designed for the discovery of Pten-cooperating tumor suppressors. Toward this aim, we coupled mobilization of a single-copy inactivating Sleeping Beauty transposon to Pten disruption within the same genome. The analysis of 278 transposition-induced prostate, breast and skin tumors detected tissue-specific and shared data sets of known and candidate genes involved in cancer. We validated ZBTB20, CELF2, PARD3, AKAP13 and WAC, which were identified by our screens in multiple cancer types, as new tumor suppressor genes in prostate cancer. We demonstrated their synergy with PTEN in preventing invasion in vitro and confirmed their clinical relevance. Further characterization of Wac in vivo showed obligate haploinsufficiency for this gene (which encodes an autophagy-regulating factor) in a Pten-deficient context. Our study identified complex PTEN-cooperating tumor suppressor networks in different cancer types, with potential clinical implications.

Genome-wide transposon screening and quantitative insertion site sequencing for cancer gene discovery in mice.

Transposon-mediated forward genetics screening in mice has emerged as a powerful tool for cancer gene discovery. It pinpoints cancer drivers that are difficult to find with other approaches, thus complementing the sequencing-based census of human cancer genes. We describe here a large series of mouse lines for insertional mutagenesis that are compatible with two transposon systems, PiggyBac and Sleeping Beauty, and give guidance on the use of different engineered transposon variants for constitutive or tissue-specific cancer gene discovery screening. We also describe a method for semiquantitative transposon insertion site sequencing (QiSeq). The QiSeq library preparation protocol exploits acoustic DNA fragmentation to reduce bias inherent to widely used restriction-digestion-based approaches for ligation-mediated insertion site amplification. Extensive multiplexing in combination with next-generation sequencing allows affordable ultra-deep transposon insertion site recovery in high-throughput formats within 1 week. Finally, we describe principles of data analysis and interpretation for obtaining insights into cancer gene function and genetic tumor evolution.

A CRISPR Dropout Screen Identifies Genetic Vulnerabilities and Therapeutic Targets in Acute Myeloid Leukemia.

Acute myeloid leukemia (AML) is an aggressive cancer with a poor prognosis, for which mainstream treatments have not changed for decades. To identify additional therapeutic targets in AML, we optimize a genome-wide clustered regularly interspaced short palindromic repeats (CRISPR) screening platform and use it to identify genetic vulnerabilities in AML cells. We identify 492 AML-specific cell-essential genes, including several established therapeutic targets such as DOT1L, BCL2, and MEN1, and many other genes including clinically actionable candidates. We validate selected genes using genetic and pharmacological inhibition, and chose KAT2A as a candidate for downstream study. KAT2A inhibition demonstrated anti-AML activity by inducing myeloid differentiation and apoptosis, and suppressed the growth of primary human AMLs of diverse genotypes while sparing normal hemopoietic stem-progenitor cells. Our results propose that KAT2A inhibition should be investigated as a therapeutic strategy in AML and provide a large number of genetic vulnerabilities of this leukemia that can be pursued in downstream studies.

Development and validation of a comprehensive genomic diagnostic tool for myeloid malignancies.

The diagnosis of hematologic malignancies relies on multidisciplinary workflows involving morphology, flow cytometry, cytogenetic, and molecular genetic analyses. Advances in cancer genomics have identified numerous recurrent mutations with clear prognostic and/or therapeutic significance to different cancers. In myeloid malignancies, there is a clinical imperative to test for such mutations in mainstream diagnosis; however, progress toward this has been slow and piecemeal. Here we describe Karyogene, an integrated targeted resequencing/analytical platform that detects nucleotide substitutions, insertions/deletions, chromosomal translocations, copy number abnormalities, and zygosity changes in a single assay. We validate the approach against 62 acute myeloid leukemia, 50 myelodysplastic syndrome, and 40 blood DNA samples from individuals without evidence of clonal blood disorders. We demonstrate robust detection of sequence changes in 49 genes, including difficult-to-detect mutations such as FLT3 internal-tandem and mixed-lineage leukemia (MLL) partial-tandem duplications, and clinically significant chromosomal rearrangements including MLL translocations to known and unknown partners, identifying the novel fusion gene MLL-DIAPH2 in the process. Additionally, we identify most significant chromosomal gains and losses, and several copy neutral loss-of-heterozygosity mutations at a genome-wide level, including previously unreported changes such as homozygosity for DNMT3A R882 mutations. Karyogene represents a dependable genomic diagnosis platform for translational research and for the clinical management of myeloid malignancies, which can be readily adapted for use in other cancers.

Leukemia-associated somatic mutations drive distinct patterns of age-related clonal hemopoiesis.

Clonal hemopoiesis driven by leukemia-associated gene mutations can occur without evidence of a blood disorder. To investigate this phenomenon, we interrogated 15 mutation hot spots in blood DNA from 4,219 individuals using ultra-deep sequencing. Using only the hot spots studied, we identified clonal hemopoiesis in 0.8% of individuals under 60, rising to 19.5% of those ≥90 years, thus predicting that clonal hemopoiesis is much more prevalent than previously realized. DNMT3A-R882 mutations were most common and, although their prevalence increased with age, were found in individuals as young as 25 years. By contrast, mutations affecting spliceosome genes SF3B1 and SRSF2, closely associated with the myelodysplastic syndromes, were identified only in those aged >70 years, with several individuals harboring more than one such mutation. This indicates that spliceosome gene mutations drive clonal expansion under selection pressures particular to the aging hemopoietic system and explains the high incidence of clonal disorders associated with these mutations in advanced old age.

A conditional piggyBac transposition system for genetic screening in mice identifies oncogenic networks in pancreatic cancer.

Here we describe a conditional piggyBac transposition system in mice and report the discovery of large sets of new cancer genes through a pancreatic insertional mutagenesis screen. We identify Foxp1 as an oncogenic transcription factor that drives pancreatic cancer invasion and spread in a mouse model and correlates with lymph node metastasis in human patients with pancreatic cancer. The propensity of piggyBac for open chromatin also enabled genome-wide screening for cancer-relevant noncoding DNA, which pinpointed a Cdkn2a cis-regulatory region. Histologically, we observed different tumor subentities and discovered associated genetic events, including Fign insertions in hepatoid pancreatic cancer. Our studies demonstrate the power of genetic screening to discover cancer drivers that are difficult to identify by other approaches to cancer genome analysis, such as downstream targets of commonly mutated human cancer genes. These piggyBac resources are universally applicable in any tissue context and provide unique experimental access to the genetic complexity of cancer.