Single-cell atlas of the human brain vasculature across development, adulthood and disease.

A broad range of brain pathologies critically relies on the vasculature, and cerebrovascular disease is a leading cause of death worldwide. However, the cellular and molecular architecture of the human brain vasculature remains incompletely understood1. Here we performed single-cell RNA sequencing analysis of 606,380 freshly isolated endothelial cells, perivascular cells and other tissue-derived cells from 117 samples, from 68 human fetuses and adult patients to construct a molecular atlas of the developing fetal, adult control and diseased human brain vasculature. We identify extensive molecular heterogeneity of the vasculature of healthy fetal and adult human brains and across five vascular-dependent central nervous system (CNS) pathologies, including brain tumours and brain vascular malformations. We identify alteration of arteriovenous differentiation and reactivated fetal as well as conserved dysregulated genes and pathways in the diseased vasculature. Pathological endothelial cells display a loss of CNS-specific properties and reveal an upregulation of MHC class II molecules, indicating atypical features of CNS endothelial cells. Cell-cell interaction analyses predict substantial endothelial-to-perivascular cell ligand-receptor cross-talk, including immune-related and angiogenic pathways, thereby revealing a central role for the endothelium within brain neurovascular unit signalling networks. Our single-cell brain atlas provides insights into the molecular architecture and heterogeneity of the developing, adult/control and diseased human brain vasculature and serves as a powerful reference for future studies.

Interaction domains of Sos1/Grb2 are finely tuned for cooperative control of embryonic stem cell fate.

Metazoan evolution involves increasing protein domain complexity, but how this relates to control of biological decisions remains uncertain. The Ras guanine nucleotide exchange factor (RasGEF) Sos1 and its adaptor Grb2 are multidomain proteins that couple fibroblast growth factor (FGF) signaling to activation of the Ras-Erk pathway during mammalian development and drive embryonic stem cells toward the primitive endoderm (PrE) lineage. We show that the ability of Sos1/Grb2 to appropriately regulate pluripotency and differentiation factors and to initiate PrE development requires collective binding of multiple Sos1/Grb2 domains to their protein and phospholipid ligands. This provides a cooperative system that only allows lineage commitment when all ligand-binding domains are occupied. Furthermore, our results indicate that the interaction domains of Sos1 and Grb2 have evolved so as to bind ligands not with maximal strength but with specificities and affinities that maintain cooperativity. This optimized system ensures that PrE lineage commitment occurs in a timely and selective manner during embryogenesis.

Obesity and the cerebral cortex: Underlying neurobiology in mice and humans.

Obesity is a major modifiable risk factor for Alzheimer's disease (AD), characterized by progressive atrophy of the cerebral cortex. The neurobiology of obesity contributions to AD is poorly understood. Here we show with in vivo MRI that diet-induced obesity decreases cortical volume in mice, and that higher body adiposity associates with lower cortical volume in humans. Single-nuclei transcriptomics of the mouse cortex reveals that dietary obesity promotes an array of neuron-adverse transcriptional dysregulations, which are mediated by an interplay of excitatory neurons and glial cells, and which involve microglial activation and lowered neuronal capacity for neuritogenesis and maintenance of membrane potential. The transcriptional dysregulations of microglia, more than of other cell types, are like those in AD, as assessed with single-nuclei cortical transcriptomics in a mouse model of AD and two sets of human donors with the disease. Serial two-photon tomography of microglia demonstrates microgliosis throughout the mouse cortex. The spatial pattern of adiposity-cortical volume associations in human cohorts interrogated together with in silico bulk and single-nucleus transcriptomic data from the human cortex implicated microglia (along with other glial cells and subtypes of excitatory neurons), and it correlated positively with the spatial profile of cortical atrophy in patients with mild cognitive impairment and AD. Thus, multi-cell neuron-adverse dysregulations likely contribute to the loss of cortical tissue in obesity. The dysregulations of microglia may be pivotal to the obesity-related risk of AD.

Mitochondrial carrier homolog 2 is necessary for AML survival.

Through a clustered regularly insterspaced short palindromic repeats (CRISPR) screen to identify mitochondrial genes necessary for the growth of acute myeloid leukemia (AML) cells, we identified the mitochondrial outer membrane protein mitochondrial carrier homolog 2 (MTCH2). In AML, knockdown of MTCH2 decreased growth, reduced engraftment potential of stem cells, and induced differentiation. Inhibiting MTCH2 in AML cells increased nuclear pyruvate and pyruvate dehydrogenase (PDH), which induced histone acetylation and subsequently promoted the differentiation of AML cells. Thus, we have defined a new mechanism by which mitochondria and metabolism regulate AML stem cells and gene expression.

Functional CRISPR and shRNA Screens Identify Involvement of Mitochondrial Electron Transport in the Activation of Evofosfamide.

Evofosfamide (TH-302) is a hypoxia-activated DNA-crosslinking prodrug currently in clinical development for cancer therapy. Oxygen-sensitive activation of evofosfamide depends on one-electron reduction, yet the reductases that catalyze this process in tumors are unknown. We used RNA sequencing, whole-genome CRISPR knockout, and reductase-focused short hairpin RNA screens to interrogate modifiers of evofosfamide activation in cancer cell lines. Involvement of mitochondrial electron transport in the activation of evofosfamide and the related nitroaromatic compounds EF5 and FSL-61 was investigated using 143B ρ 0 (ρ zero) cells devoid of mitochondrial DNA and biochemical assays in UT-SCC-74B cells. The potency of evofosfamide in 30 genetically diverse cancer cell lines correlated with the expression of genes involved in mitochondrial electron transfer. A whole-genome CRISPR screen in KBM-7 cells identified the DNA damage-response factors SLX4IP, C10orf90 (FATS), and SLFN11, in addition to the key regulator of mitochondrial function, YME1L1, and several complex I constituents as modifiers of evofosfamide sensitivity. A reductase-focused shRNA screen in UT-SCC-74B cells similarly identified mitochondrial respiratory chain factors. Surprisingly, 143B ρ 0 cells showed enhanced evofosfamide activation and sensitivity but had global transcriptional changes, including increased expression of nonmitochondrial flavoreductases. In UT-SCC-74B cells, evofosfamide oxidized cytochromes a, b, and c and inhibited respiration at complexes I, II, and IV without quenching reactive oxygen species production. Our results suggest that the mitochondrial electron transport chain contributes to evofosfamide activation and that predicting evofosfamide sensitivity in patients by measuring the expression of canonical bioreductive enzymes such as cytochrome P450 oxidoreductase is likely to be futile.

Forward genetic screen in human podocytes identifies diphthamide biosynthesis genes as regulators of adhesion.

Podocyte function is tightly linked to the complex organization of its cytoskeleton and adhesion to the underlying glomerular basement membrane. Adhesion of cultured podocytes to a variety of substrates is reported to correlate with podocyte health. To identify novel genes that are important for podocyte function, we designed an in vitro genetic screen based on podocyte adhesion to plates coated with either fibronectin or soluble Fms-like tyrosine kinase-1 (sFLT1)/Fc. A genome-scale pooled RNA interference screen on immortalized human podocytes identified 77 genes that increased adhesion to fibronectin, 101 genes that increased adhesion to sFLT1/Fc, and 44 genes that increased adhesion to both substrates when knocked down. Multiple shRNAs against diphthamide biosynthesis protein 1-4 (DPH1-DPH4) were top hits for increased adhesion. Immortalized human podocyte cells stably expressing these hairpins displayed increased adhesion to both substrates. We then used CRISPR-Cas9 to generate podocyte knockout cells for DPH1, DPH2, or DPH3, which also displayed increased adhesion to both fibronectin and sFLT1/Fc, as well as a spreading defect. Finally, we showed that Drosophila nephrocyte-specific knockdown of Dph1, Dph2, and Dph4 resulted in altered nephrocyte function. In summary, we report here a novel high-throughput method to identify genes important for podocyte function. Given the central role of podocyte adhesion as a marker of podocyte health, these data are a rich source of candidate regulators of glomerular disease.

Functional Genomic Analysis of Candida albicans Adherence Reveals a Key Role for the Arp2/3 Complex in Cell Wall Remodelling and Biofilm Formation.

Fungal biofilms are complex, structured communities that can form on surfaces such as catheters and other indwelling medical devices. Biofilms are of particular concern with Candida albicans, one of the leading opportunistic fungal pathogens of humans. C. albicans biofilms include yeast and filamentous cells that are surrounded by an extracellular matrix, and they are intrinsically resistant to antifungal drugs such that resolving biofilm infections often requires surgery to remove the contaminated device. C. albicans biofilms form through a regulated process of adhesion to surfaces, filamentation, maturation, and ultimately dispersion. To uncover new strategies to block the initial stages of biofilm formation, we utilized a functional genomic approach to identify genes that modulate C. albicans adherence. We screened a library of 1,481 double barcoded doxycycline-repressible conditional gene expression strains covering ~25% of the C. albicans genome. We identified five genes for which transcriptional repression impaired adherence, including: ARC18, PMT1, MNN9, SPT7, and orf19.831. The most severe adherence defect was observed upon transcriptional repression of ARC18, which encodes a member of the Arp2/3 complex that is involved in regulation of the actin cytoskeleton and endocytosis. Depletion of components of the Arp2/3 complex not only impaired adherence, but also caused reduced biofilm formation, increased cell surface hydrophobicity, and increased exposure of cell wall chitin and ß-glucans. Reduced function of the Arp2/3 complex led to impaired cell wall integrity and activation of Rho1-mediated cell wall stress responses, thereby causing cell wall remodelling and reduced adherence. Thus, we identify important functional relationships between cell wall stress responses and a novel mechanism that controls adherence and biofilm formation, thereby illuminating novel strategies to cripple a leading fungal pathogen of humans.

Differential regulation of FGFR3 by PTPN1 and PTPN2.

Aberrant expression and activation of FGFR3 is associated with disease states including bone dysplasia and malignancies of bladder, cervix, and bone marrow. MS analysis of protein-phosphotyrosine in multiple myeloma cells revealed a prevalent phosphorylated motif, D/EYYR/K, derived from the kinase domain activation loops of tyrosine kinases including FGFR3 corresponding to a recognition sequence of protein-tyrosine phosphatase PTPN1. Knockdown of PTPN1 or the related enzyme PTPN2 by RNAi resulted in ligand-independent activation of FGFR3. Modulation of FGFR3 activation loop phosphorylation by both PTPN1 and PTPN2 was a function of receptor trafficking and phosphotyrosine phosphatase (PTP) compartmentalization. The FGFR3 activation loop motif DYYKK(650) is altered to DYYKE(650) in the oncogenic variant FGFR3(K650E) , and consequently it is constitutively fully activated and unaffected by activation loop phosphorylation. FGFR3(K650E) was nevertheless remarkably sensitive to negative regulation by PTPN1 and PTPN2. This suggests that in addition to modulating FGFR3 phosphorylation, PTPN1 and PTPN2 constrain the kinase domain by fostering an inactive-state. Loss of this constraint in response to ligand or impaired PTPN1/N2 may initiate FGFR3 activation. These results suggest a model wherein PTP expression levels may define conditions that select for ectopic FGFR3 expression and activation during tumorigenesis.

A negative genetic interaction map in isogenic cancer cell lines reveals cancer cell vulnerabilities.

Improved efforts are necessary to define the functional product of cancer mutations currently being revealed through large-scale sequencing efforts. Using genome-scale pooled shRNA screening technology, we mapped negative genetic interactions across a set of isogenic cancer cell lines and confirmed hundreds of these interactions in orthogonal co-culture competition assays to generate a high-confidence genetic interaction network of differentially essential or differential essentiality (DiE) genes. The network uncovered examples of conserved genetic interactions, densely connected functional modules derived from comparative genomics with model systems data, functions for uncharacterized genes in the human genome and targetable vulnerabilities. Finally, we demonstrate a general applicability of DiE gene signatures in determining genetic dependencies of other non-isogenic cancer cell lines. For example, the PTEN(-/-) DiE genes reveal a signature that can preferentially classify PTEN-dependent genotypes across a series of non-isogenic cell lines derived from the breast, pancreas and ovarian cancers. Our reference network suggests that many cancer vulnerabilities remain to be discovered through systematic derivation of a network of differentially essential genes in an isogenic cancer cell model.

COLT-Cancer: functional genetic screening resource for essential genes in human cancer cell lines.

Genome-wide pooled shRNA screens enable global identification of the genes essential for cancer cell survival and proliferation and provide a 'functional genetic' map of human cancer to complement genomic studies. Using a lentiviral shRNA library targeting approximately 16,000 human genes and a newly developed scoring approach, we identified essential gene profiles in more than 70 breast, pancreatic and ovarian cancer cell lines. We developed a web-accessible database system for capturing information from each step in our standardized screening pipeline and a gene-centric search tool for exploring shRNA activities within a given cell line or across multiple cell lines. The database consists of a laboratory information and management system for tracking each step of a pooled shRNA screen as well as a web interface for querying and visualization of shRNA and gene-level performance across multiple cancer cell lines. COLT-Cancer Version 1.0 is currently accessible at http://colt.ccbr.utoronto.ca/cancer.

Chemical-Genetic Approaches for Exploring Mode of Action of Antifungal Compounds in the Fungal Pathogen Candida albicans.

Candida albicans is a prevalent fungal pathogen of humans that can cause both superficial and life-threatening disease, primarily in immunocompromised populations. Currently, antifungal drug classes available to treat fungal infections remain limited and the emergence of drug-resistant strains threatens antifungal efficacy, necessitating the discovery and development of additional therapeutics. The construction of the C. albicans double-barcoded heterozygous deletion collection (DBC) enables the rapid and systematic assessment of haploinsufficiency phenotypes in a pooled format. Specifically, this functional genomics resource can be used to identify heterozygous deletion mutants that are hypersensitive to compounds in order to define putative cellular targets and/or other modifiers of compound activity. Here, we describe protocols to characterize the mode of action of small molecules using the C. albicans DBC, including how to prepare compound-treated cultures, isolate genomic DNA, amplify strain-specific barcodes, and prepare DNA libraries for high-throughput sequencing. This technique provides a powerful approach to elucidate the compound mechanism of action in order to bolster the antifungal pipeline.

Decreased left heart flow in fetal lambs causes left heart hypoplasia and pro-fibrotic tissue remodeling.

Low blood flow through the fetal left heart is often conjectured as an etiology for hypoplastic left heart syndrome (HLHS). To investigate if a decrease in left heart flow results in growth failure, we generate left ventricular inflow obstruction (LVIO) in mid-gestation fetal lambs by implanting coils in their left atrium using an ultrasound-guided percutaneous technique. Significant LVIO recapitulates important clinical features of HLHS: decreased antegrade aortic valve flow, compensatory retrograde perfusion of the brain and ascending aorta (AAo) from the arterial duct, severe left heart hypoplasia, a non-apex forming LV, and a thickened endocardial layer. The hypoplastic AAo have miRNA-gene pairs annotating to cell proliferation that are inversely differentially expressed by bulk RNA-seq. Single-nucleus RNA-seq of the hypoplastic LV myocardium shows an increase in fibroblasts with a reciprocal decrease in cardiomyocyte nuclei proportions. Fibroblasts, cardiomyocytes and endothelial cells from hypoplastic myocardium have increased expression of extracellular matrix component or fibrosis genes with dysregulated fibroblast growth factor signaling. Hence, a severe sustained ( ~ 1/3 gestation) reduction in fetal left heart flow is sufficient to cause left heart hypoplasia. This is accompanied by changes in cellular composition and gene expression consistent with a pro-fibrotic environment and aberrant induction of mesenchymal programs.

Identification of acquired Notch3 dependency in metastatic Head and Neck Cancer.

During cancer development, tumor cells acquire changes that enable them to invade surrounding tissues and seed metastasis at distant sites. These changes contribute to the aggressiveness of metastatic cancer and interfere with success of therapy. Our comprehensive analysis of "matched" pairs of HNSCC lines derived from primary tumors and corresponding metastatic sites identified several components of Notch3 signaling that are differentially expressed and/or altered in metastatic lines and confer a dependency on this pathway. These components were also shown to be differentially expressed between early and late stages of tumors in a TMA constructed from over 200 HNSCC patients. Finally, we show that suppression of Notch3 improves survival in mice in both subcutaneous and orthotopic models of metastatic HNSCC. Novel treatments targeting components of this pathway may prove effective in targeting metastatic HNSCC cells alone or in combination with conventional therapies.

Single-nucleus multiomic atlas of frontal cortex in amyotrophic lateral sclerosis with a deep learning-based decoding of alternative polyadenylation mechanisms.

The understanding of how different cell types contribute to amyotrophic lateral sclerosis (ALS) pathogenesis is limited. Here we generated a single-nucleus transcriptomic and epigenomic atlas of the frontal cortex of ALS cases with C9orf72 (C9) hexanucleotide repeat expansions and sporadic ALS (sALS). Our findings reveal shared pathways in C9-ALS and sALS, characterized by synaptic dysfunction in excitatory neurons and a disease-associated state in microglia. The disease subtypes diverge with loss of astrocyte homeostasis in C9-ALS, and a more substantial disturbance of inhibitory neurons in sALS. Leveraging high depth 3'-end sequencing, we found a widespread switch towards distal polyadenylation (PA) site usage across ALS subtypes relative to controls. To explore this differential alternative PA (APA), we developed APA-Net, a deep neural network model that uses transcript sequence and expression levels of RNA-binding proteins (RBPs) to predict cell-type specific APA usage and RBP interactions likely to regulate APA across disease subtypes.

Analysis of early C2C12 myogenesis identifies stably and differentially expressed transcriptional regulators whose knock-down inhibits myoblast differentiation.

Myogenesis is a tightly controlled process involving the transcriptional activation and repression of thousands of genes. Although many components of the transcriptional network regulating the later phases of myogenesis have been identified, relatively few studies have described the transcriptional landscape during the first 24 h, when myoblasts commit to differentiate. Through dense temporal profiling of differentiating C2C12 myoblasts, we identify 193 transcriptional regulators (TRs) whose expression is significantly altered within the first 24 h of myogenesis. A high-content shRNA screen of 77 TRs involving 427 stable lines identified 42 genes whose knockdown significantly inhibits differentiation of C2C12 myoblasts. Of the TRs that were differentially expressed within the first 24 h, over half inhibited differentiation when knocked down, including known regulators of myogenesis (Myod1, Myog, and Myf5), as well as 19 TRs not previously associated with this process. Surprisingly, a similar proportion (55%) of shRNAs targeting TRs whose expression did not change also inhibited C2C12 myogenesis. We further show that a subset of these TRs inhibits myogenesis by downregulating expression of known regulatory and structural proteins. Our findings clearly illustrate that several TRs critical for C2C12 myogenesis are not differentially regulated, suggesting that approaches that focus functional studies on differentially-expressed transcripts will fail to provide a comprehensive view of this complex process.

CD133 protein N-glycosylation processing contributes to cell surface recognition of the primitive cell marker AC133 epitope.

The AC133 epitope expressed on the CD133 glycoprotein has been widely used as a cell surface marker of numerous stem cell and cancer stem cell types. It has been recently proposed that posttranslational modification and regulation of CD133 may govern cell surface AC133 recognition. Therefore, we performed a large scale pooled RNA interference (RNAi) screen to identify genes involved in cell surface AC133 expression. Gene hits could be validated at a rate of 70.5% in a secondary assay using an orthogonal RNAi system, demonstrating that our primary RNAi screen served as a powerful genetic screening approach. Within the list of hits from the primary screen, genes involved in N-glycan biosynthesis were significantly enriched as determined by Ingenuity Canonical Pathway analyses. Indeed, inhibiting biosynthesis of the N-glycan precursor using the small molecule tunicamycin or inhibiting its transfer to CD133 by generating N-glycan-deficient CD133 mutants resulted in undetectable cell surface AC133. Among the screen hits involved in N-glycosylation were genes involved in complex N-glycan processing, including the poorly characterized MGAT4C, which we demonstrate to be a positive regulator of cell surface AC133 expression. Our study identifies a set of genes involved in CD133 N-glycosylation as a direct contributing factor to cell surface AC133 recognition and provides biochemical evidence for the function and structure of CD133 N-glycans.

A genome wide shRNA screen identifies α/ß hydrolase domain containing 4 (ABHD4) as a novel regulator of anoikis resistance.

Acquisition of resistance to anchorage dependant cell death, a process termed anoikis, is a requirement for cancer cell metastasis. However, the molecular determinants of anoikis resistance and sensitivity are poorly understood. To better understand resistance to anoikis we conducted a genome wide lentiviral shRNA screen to identify genes whose knockdown render anoikis-sensitive RWPE-1 prostate cells resistant to anoikis. RWPE-1 cells were infected with a pooled lentiviral shRNA library with 54,021 shRNA targeting 11,255 genes. After infection, an anoikis-resistant cell population was selected and shRNA sequences were amplified and sequenced. Thirty-four shRNA sequences reproducibly protected RWPE-1 cells from anoikis after culture under suspension conditions including the top validated hit, α/ß hydrolase domain containing 4 (ABHD4). In validation studies, ABHD4 knockdown inhibited anoikis in RWPE-1 cells as well as anoikis sensitive NP69 nasopharyngeal and OVCAR3 ovarian cancer cells, while over-expression of the gene increased sensitivity. Induction of anoikis after ABHD4 knockdown was associated with cleavage of PARP and activation of caspases-3, but was independent in changes of FLIP, FAK and Src expression. Interestingly, induction of anoikis after ABHD4 knockdown was independent of the known role of ABHD4 in the anandamide synthesis pathway and the generation of glycerophospho-N-acyl ethanolamines. Thus, ABHD4 is a novel genetic regulator of anoikis sensitivity.

The ubiquitin-activating enzyme E1 as a therapeutic target for the treatment of leukemia and multiple myeloma.

The proteasomal pathway of protein degradation involves 2 discrete steps: ubiquitination and degradation. Here, we evaluated the effects of inhibiting the ubiquitination pathway at the level of the ubiquitin-activating enzyme UBA1 (E1). By immunoblotting, leukemia cell lines and primary patient samples had increased protein ubiquitination. Therefore, we examined the effects of genetic and chemical inhibition of the E1 enzyme. Knockdown of E1 decreased the abundance of ubiquitinated proteins in leukemia and myeloma cells and induced cell death. To further investigate effects of E1 inhibition in malignancy, we discovered a novel small molecule inhibitor, 3,5-dioxopyrazolidine compound, 1-(3-chloro-4-fluorophenyl)-4-[(5-nitro-2-furyl)methylene]-3,5-pyrazolidinedione (PYZD-4409). PYZD-4409 induced cell death in malignant cells and preferentially inhibited the clonogenic growth of primary acute myeloid leukemia cells compared with normal hematopoietic cells. Mechanistically, genetic or chemical inhibition of E1 increased expression of E1 stress markers. Moreover, BI-1 overexpression blocked cell death after E1 inhibition, suggesting ER stress is functionally important for cell death after E1 inhibition. Finally, in a mouse model of leukemia, intraperitoneal administration of PYZD-4409 decreased tumor weight and volume compared with control without untoward toxicity. Thus, our work highlights the E1 enzyme as a novel target for the treatment of hematologic malignancies.

SATB2 augments ΔNp63α in head and neck squamous cell carcinoma.

ΔNp63α is a critical pro-survival protein overexpressed in 80% of head and neck squamous cell carcinomas (HNSCCs) where it inhibits TAp73ß transcription of p53-family target genes, which is thought to increase HNSCC resistance to chemotherapy-induced cell death. However, the mechanisms governing ΔNp63α function are largely unknown. In this study, we identify special AT-rich-binding protein 2 (SATB2) as a new ΔNp63α-binding protein that is preferentially expressed in advanced-stage primary HNSCC and show that SATB2 promotes chemoresistance by enhancing ΔNp63α-mediated transrepression by augmenting ΔNp63α engagement to p53-family responsive elements. Furthermore, SATB2 expression positively correlates with HNSCC chemoresistance, and RNA interference-mediated knockdown of endogenous SATB2 re-sensitizes HNSCC cells to chemotherapy- and γ-irradiation-induced apoptosis, irrespective of p53 status. These findings unveil SATB2 as a pivotal modulator of ΔNp63α that governs HNSCC cell survival.

Genome-wide in vivo screen of circulating tumor cells identifies SLIT2 as a regulator of metastasis.

Circulating tumor cells (CTCs) break free from primary tumors and travel through the circulation system to seed metastatic tumors, which are the major cause of death from cancer. The identification of the major genetic factors that enhance production and persistence of CTCs in the bloodstream at a whole genome level would enable more comprehensive molecular mechanisms of metastasis to be elucidated and the identification of novel therapeutic targets, but this remains a challenging task due to the heterogeneity and extreme rarity of CTCs. Here, we describe an in vivo genome-wide CRISPR knockout screen using CTCs directly isolated from a mouse xenograft. This screen elucidated SLIT2-a gene encoding a secreted protein acting as a cellular migration cue-as the most significantly represented gene knockout in the CTC population. SLIT2 knockout cells are highly metastatic with hypermigratory and mesenchymal phenotype, resulting in enhanced cancer progression in xenograft models.