Colorectal Cancer Cells Enter a Diapause-like DTP State to Survive Chemotherapy.

Cancer cells enter a reversible drug-tolerant persister (DTP) state to evade death from chemotherapy and targeted agents. It is increasingly appreciated that DTPs are important drivers of therapy failure and tumor relapse. We combined cellular barcoding and mathematical modeling in patient-derived colorectal cancer models to identify and characterize DTPs in response to chemotherapy. Barcode analysis revealed no loss of clonal complexity of tumors that entered the DTP state and recurred following treatment cessation. Our data fit a mathematical model where all cancer cells, and not a small subpopulation, possess an equipotent capacity to become DTPs. Mechanistically, we determined that DTPs display remarkable transcriptional and functional similarities to diapause, a reversible state of suspended embryonic development triggered by unfavorable environmental conditions. Our study provides insight into how cancer cells use a developmentally conserved mechanism to drive the DTP state, pointing to novel therapeutic opportunities to target DTPs.

Functional Genomic Landscape of Human Breast Cancer Drivers, Vulnerabilities, and Resistance.

Large-scale genomic studies have identified multiple somatic aberrations in breast cancer, including copy number alterations and point mutations. Still, identifying causal variants and emergent vulnerabilities that arise as a consequence of genetic alterations remain major challenges. We performed whole-genome small hairpin RNA (shRNA) "dropout screens" on 77 breast cancer cell lines. Using a hierarchical linear regression algorithm to score our screen results and integrate them with accompanying detailed genetic and proteomic information, we identify vulnerabilities in breast cancer, including candidate "drivers," and reveal general functional genomic properties of cancer cells. Comparisons of gene essentiality with drug sensitivity data suggest potential resistance mechanisms, effects of existing anti-cancer drugs, and opportunities for combination therapy. Finally, we demonstrate the utility of this large dataset by identifying BRD4 as a potential target in luminal breast cancer and PIK3CA mutations as a resistance determinant for BET-inhibitors.

High-Resolution CRISPR Screens Reveal Fitness Genes and Genotype-Specific Cancer Liabilities.

The ability to perturb genes in human cells is crucial for elucidating gene function and holds great potential for finding therapeutic targets for diseases such as cancer. To extend the catalog of human core and context-dependent fitness genes, we have developed a high-complexity second-generation genome-scale CRISPR-Cas9 gRNA library and applied it to fitness screens in five human cell lines. Using an improved Bayesian analytical approach, we consistently discover 5-fold more fitness genes than were previously observed. We present a list of 1,580 human core fitness genes and describe their general properties. Moreover, we demonstrate that context-dependent fitness genes accurately recapitulate pathway-specific genetic vulnerabilities induced by known oncogenes and reveal cell-type-specific dependencies for specific receptor tyrosine kinases, even in oncogenic KRAS backgrounds. Thus, rigorous identification of human cell line fitness genes using a high-complexity CRISPR-Cas9 library affords a high-resolution view of the genetic vulnerabilities of a cell.

Systematic exploration of dynamic splicing networks reveals conserved multistage regulators of neurogenesis.

Alternative splicing (AS) is a critical regulatory layer; yet, factors controlling functionally coordinated splicing programs during developmental transitions are poorly understood. Here, we employ a screening strategy to identify factors controlling dynamic splicing events important for mammalian neurogenesis. Among previously unknown regulators, Rbm38 acts widely to negatively control neural AS, in part through interactions mediated by the established repressor of splicing, Ptbp1. Puf60, a ubiquitous factor, is surprisingly found to promote neural splicing patterns. This activity requires a conserved, neural-differential exon that remodels Puf60 co-factor interactions. Ablation of this exon rewires distinct AS networks in embryonic stem cells and at different stages of mouse neurogenesis. Single-cell transcriptome analyses further reveal distinct roles for Rbm38 and Puf60 isoforms in establishing neuronal identity. Our results describe important roles for previously unknown regulators of neurogenesis and establish how an alternative exon in a widely expressed splicing factor orchestrates temporal control over cell differentiation.

Functional genomic landscape of cancer-intrinsic evasion of killing by T cells.

The genetic circuits that allow cancer cells to evade destruction by the host immune system remain poorly understood1-3. Here, to identify a phenotypically robust core set of genes and pathways that enable cancer cells to evade killing mediated by cytotoxic T lymphocytes (CTLs), we performed genome-wide CRISPR screens across a panel of genetically diverse mouse cancer cell lines that were cultured in the presence of CTLs. We identify a core set of 182 genes across these mouse cancer models, the individual perturbation of which increases either the sensitivity or the resistance of cancer cells to CTL-mediated toxicity. Systematic exploration of our dataset using genetic co-similarity reveals the hierarchical and coordinated manner in which genes and pathways act in cancer cells to orchestrate their evasion of CTLs, and shows that discrete functional modules that control the interferon response and tumour necrosis factor (TNF)-induced cytotoxicity are dominant sub-phenotypes. Our data establish a central role for genes that were previously identified as negative regulators of the type-II interferon response (for example, Ptpn2, Socs1 and Adar1) in mediating CTL evasion, and show that the lipid-droplet-related gene Fitm2 is required for maintaining cell fitness after exposure to interferon-γ (IFNγ). In addition, we identify the autophagy pathway as a conserved mediator of the evasion of CTLs by cancer cells, and show that this pathway is required to resist cytotoxicity induced by the cytokines IFNγ and TNF. Through the mapping of cytokine- and CTL-based genetic interactions, together with in vivo CRISPR screens, we show how the pleiotropic effects of autophagy control cancer-cell-intrinsic evasion of killing by CTLs and we highlight the importance of these effects within the tumour microenvironment. Collectively, these data expand our knowledge of the genetic circuits that are involved in the evasion of the immune system by cancer cells, and highlight genetic interactions that contribute to phenotypes associated with escape from killing by CTLs.

System-Wide Modulation of HECT E3 Ligases with Selective Ubiquitin Variant Probes.

HECT-family E3 ligases ubiquitinate protein substrates to control virtually every eukaryotic process and are misregulated in numerous diseases. Nonetheless, understanding of HECT E3s is limited by a paucity of selective and potent modulators. To overcome this challenge, we systematically developed ubiquitin variants (UbVs) that inhibit or activate HECT E3s. Structural analysis of 6 HECT-UbV complexes revealed UbV inhibitors hijacking the E2-binding site and activators occupying a ubiquitin-binding exosite. Furthermore, UbVs unearthed distinct regulation mechanisms among NEDD4 subfamily HECTs and proved useful for modulating therapeutically relevant targets of HECT E3s in cells and intestinal organoids, and in a genetic screen that identified a role for NEDD4L in regulating cell migration. Our work demonstrates versatility of UbVs for modulating activity across an E3 family, defines mechanisms and provides a toolkit for probing functions of HECT E3s, and establishes a general strategy for systematic development of modulators targeting families of signaling proteins.

A method for benchmarking genetic screens reveals a predominant mitochondrial bias.

We present FLEX (Functional evaluation of experimental perturbations), a pipeline that leverages several functional annotation resources to establish reference standards for benchmarking human genome-wide CRISPR screen data and methods for analyzing them. FLEX provides a quantitative measurement of the functional information captured by a given gene-pair dataset and a means to explore the diversity of functions captured by the input dataset. We apply FLEX to analyze data from the diverse cell line screens generated by the DepMap project. We identify a predominant mitochondria-associated signal within co-essentiality networks derived from these data and explore the basis of this signal. Our analysis and time-resolved CRISPR screens in a single cell line suggest that the variable phenotypes associated with mitochondria genes across cells may reflect screen dynamics and protein stability effects rather than genetic dependencies. We characterize this functional bias and demonstrate its relevance for interpreting differential hits in any CRISPR screening context. More generally, we demonstrate the utility of the FLEX pipeline for performing robust comparative evaluations of CRISPR screens or methods for processing them.

Functional genomic characterization of a synthetic anti-HER3 antibody reveals a role for ubiquitination by RNF41 in the anti-proliferative response.

Dysregulation of the ErbB family of receptor tyrosine kinases is involved in the progression of many cancers. Antibodies targeting the dimerization domains of family members EGFR and HER2 are approved cancer therapeutics, but efficacy is restricted to a subset of tumors and resistance often develops in response to treatment. A third family member, HER3, heterodimerizes with both EGFR and HER2 and has also been implicated in cancer. Consequently, there is strong interest in developing antibodies that target HER3, but to date, no therapeutics have been approved. To aid the development of anti-HER3 antibodies as cancer therapeutics, we combined antibody engineering and functional genomics screens to identify putative mechanisms of resistance or synthetic lethality with antibody-mediated anti-proliferative effects. We developed a synthetic antibody called IgG 95, which binds to HER3 and promotes ubiquitination, internalization, and receptor down-regulation. Using an shRNA library targeting enzymes in the ubiquitin proteasome system, we screened for genes that effect response to IgG 95 and uncovered the E3 ubiquitin ligase RNF41 as a driver of IgG 95 anti-proliferative activity. RNF41 has been shown previously to regulate HER3 levels under normal conditions and we now show that it is also responsible for down-regulation of HER3 upon treatment with IgG 95. Moreover, our findings suggest that down-regulation of RNF41 itself may be a mechanism for acquired resistance to treatment with IgG 95 and perhaps other anti-HER3 antibodies. Our work deepens our understanding of HER3 signaling by uncovering the mechanistic basis for the anti-proliferative effects of potential anti-HER3 antibody therapeutics.

CRISPR screens are feasible in TP53 wild-type cells.

A recent study by Haapaniemi et al (2018) reported that intact p53 signaling hampers CRISPR-based functional genomic screens. Brown et al report good performance of genome-scale screens in TP53 wild-type cells and reiterate best practices for CRISPR screening.

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.

Bmi1 regulates human glioblastoma stem cells through activation of differential gene networks in CD133+ brain tumor initiating cells.

PURPOSE: Glioblastoma (GBM) is the most aggressive adult brain cancer, with a 15 month median survivorship attributed to the existence of treatment-refractory brain tumor initiating cells (BTICs). In order to better understand the mechanisms regulating the tumorigenic properties of this population, we studied the role of the polycomb group member BMI1 in our patient-derived GBM BTICs and its relationship with CD133, a well-established marker of BTICs. METHODS: Using gain and loss-of-function studies for Bmi1 in neural stem cells (NSCs) and patient-derived GBM BTICs respectively, we assessed in vitro self-renewal and in vivo tumor formation in these two cell populations. We further explored the BMI1 transcriptional regulatory network through RNA sequencing of different GBM BTIC populations that were knocked down for Bmi1. RESULTS: There is a differential role of BMI1 in CD133-positive cells, notably involving cell metabolism. In addition, we identified pivotal targets downstream of BMI1 in CD133+ cells such as integrin alpha 2 (ITGA2), that may contribute to regulating GBM stem cell properties. CONCLUSIONS: Our work sheds light on the association of three genes with CD133-BMI1 circuitry, their importance as downstream effectors of the BMI1 signalling pathway, and their potential as future targets for tackling GBM treatment-resistant cell populations.

RNAi screen identifies essential regulators of human brain metastasis-initiating cells.

Brain metastases (BM) are the most common brain tumor in adults and are a leading cause of cancer mortality. Metastatic lesions contain subclones derived from their primary lesion, yet their functional characterization is limited by a paucity of preclinical models accurately recapitulating the metastatic cascade, emphasizing the need for a novel approach to BM and their treatment. We identified a unique subset of stem-like cells from primary human patient brain metastases, termed brain metastasis-initiating cells (BMICs). We now establish a BMIC patient-derived xenotransplantation (PDXT) model as an investigative tool to comprehensively interrogate human BM. Using both in vitro and in vivo RNA interference screens of these BMIC models, we identified SPOCK1 and TWIST2 as essential BMIC regulators. SPOCK1 in particular is a novel regulator of BMIC self-renewal, modulating tumor initiation and metastasis from the lung to the brain. A prospective cohort of primary lung cancer specimens showed that SPOCK1 was overexpressed only in patients who ultimately developed BM. Protein-protein interaction network mapping between SPOCK1 and TWIST2 identified novel pathway interactors with significant prognostic value in lung cancer patients. Of these genes, INHBA, a TGF-ß ligand found mutated in lung adenocarcinoma, showed reduced expression in BMICs with knockdown of SPOCK1. In conclusion, we have developed a useful preclinical model of BM, which has served to identify novel putative BMIC regulators, presenting potential therapeutic targets that block the metastatic process, and transform a uniformly fatal systemic disease into a locally controlled and eminently more treatable one.

Efficient generation of patient-matched malignant and normal primary cell cultures from clear cell renal cell carcinoma patients: clinically relevant models for research and personalized medicine.

BACKGROUND: Patients with clear cell renal cell carcinoma (ccRCC) have few therapeutic options, as ccRCC is unresponsive to chemotherapy and is highly resistant to radiation. Recently targeted therapies have extended progression-free survival, but responses are variable and no significant overall survival benefit has been achieved. Commercial ccRCC cell lines are often used as model systems to develop novel therapeutic approaches, but these do not accurately recapitulate primary ccRCC tumors at the genomic and transcriptional levels. Furthermore, ccRCC exhibits significant intertumor genetic heterogeneity, and the limited cell lines available fail to represent this aspect of ccRCC. Our objective was to generate accurate preclinical in vitro models of ccRCC using tumor tissues from ccRCC patients. METHODS: ccRCC primary single cell suspensions were cultured in fetal bovine serum (FBS)-containing media or defined serum-free media. Established cultures were characterized by genomic verification of mutations present in the primary tumors, expression of renal epithelial markers, and transcriptional profiling. RESULTS: The apparent efficiency of primary cell culture establishment was high in both culture conditions, but genotyping revealed that the majority of cultures contained normal, not cancer cells. ccRCC characteristically shows biallelic loss of the von Hippel Lindau (VHL) gene, leading to accumulation of hypoxia-inducible factor (HIF) and expression of HIF target genes. Purification of cells based on expression of carbonic anhydrase IX (CA9), a cell surface HIF target, followed by culture in FBS enabled establishment of ccRCC cell cultures with an efficiency of >80 %. Culture in serum-free conditions selected for growth of normal renal proximal tubule epithelial cells. Transcriptional profiling of ccRCC and matched normal cell cultures identified up- and down-regulated networks in ccRCC and comparison to The Cancer Genome Atlas confirmed the clinical validity of our cell cultures. CONCLUSIONS: The ability to establish primary cultures of ccRCC cells and matched normal kidney epithelial cells from almost every patient provides a resource for future development of novel therapies and personalized medicine for ccRCC patients.

Measuring error rates in genomic perturbation screens: gold standards for human functional genomics.

Technological advancement has opened the door to systematic genetics in mammalian cells. Genome-scale loss-of-function screens can assay fitness defects induced by partial gene knockdown, using RNA interference, or complete gene knockout, using new CRISPR techniques. These screens can reveal the basic blueprint required for cellular proliferation. Moreover, comparing healthy to cancerous tissue can uncover genes that are essential only in the tumor; these genes are targets for the development of specific anticancer therapies. Unfortunately, progress in this field has been hampered by off-target effects of perturbation reagents and poorly quantified error rates in large-scale screens. To improve the quality of information derived from these screens, and to provide a framework for understanding the capabilities and limitations of CRISPR technology, we derive gold-standard reference sets of essential and nonessential genes, and provide a Bayesian classifier of gene essentiality that outperforms current methods on both RNAi and CRISPR screens. Our results indicate that CRISPR technology is more sensitive than RNAi and that both techniques have nontrivial false discovery rates that can be mitigated by rigorous analytical methods.

Bioequivalence of Alendronate and Vitamin D3 in a Combination Tablet Versus Corresponding-Dose Individual Tablets in Healthy Taiwanese Volunteers, Determined Using a Novel Plasma Alendronate Assay.

OBJECTIVE: This study was designed to demonstrate that alendronate (ALN)/vitamin D3 combination tablets (ALN/D5600) are bioequivalent to corresponding doses of ALN and vitamin D3 as individual tablets in healthy Taiwanese volunteers. METHODS: In this open-label, randomized, 2-period, crossover study, 68 volunteers were randomized to a single ALN/D5600 combination tablet or corresponding doses of 70 mg ALN + 5600 IU vitamin D3 (2 × 2800 IU), followed by a 12-day washout period and administration of the alternate formulation. Plasma ALN levels were measured using a newly developed assay. Geometric mean ratios of ALN AUC0-last, AUC0-∞, and Cmax, and unadjusted vitamin D3 AUC0-80h and Cmax were compared and considered bioequivalent if the 90% CI was within 0.8 to 1.25. RESULTS: The geometric mean ratios were: AUC0-last, 1.084 (90% CI, 0.937-1.253); AUC0-∞, 1.081 (90% CI, 0.935-1.249); and Cmax, 1.112 (90% CI, 0.959-1.289) for ALN, and AUC0-80h 0.953 (90% CI, 0.827-1.098) and Cmax, 0.982 (90% CI, 0.854-1.130) for vitamin D3 unadjusted for endogenous levels. CONCLUSIONS: The combination tablet was considered bioequivalent to coadministration based on ALN AUC0-∞ and unadjusted vitamin D3 parameters. Slight differences for ALN AUC0-last and Cmax (upper 90% CIs outside the bounds) were not considered clinically significant. The combination tablet was well tolerated. No serious adverse experiences were reported. © 2015. The Authors. Published by Elsevier Inc. All rights reserved.

Identification of novel ATP13A2 interactors and their role in α-synuclein misfolding and toxicity.

Lysosomes are responsible for degradation and recycling of bulky cell material, including accumulated misfolded proteins and dysfunctional organelles. Increasing evidence implicates lysosomal dysfunction in several neurodegenerative disorders, including Parkinson's disease and related synucleinopathies, which are characterized by the accumulation of α-synuclein (α-syn) in Lewy bodies. Studies of lysosomal proteins linked to neurodegenerative disorders present an opportunity to uncover specific molecular mechanisms and pathways that contribute to neurodegeneration. Loss-of-function mutations in a lysosomal protein, ATP13A2 (PARK9), cause Kufor-Rakeb syndrome that is characterized by early-onset parkinsonism, pyramidal degeneration and dementia. While loss of ATP13A2 function plays a role in α-syn misfolding and toxicity, the normal function of ATP13A2 in the brain remains largely unknown. Here, we performed a screen to identify ATP13A2 interacting partners, as a first step toward elucidating its function. Utilizing a split-ubiquitin membrane yeast two-hybrid system that was developed to identify interacting partners of full-length integral membrane proteins, we identified 43 novel interactors that primarily implicate ATP13A2 in cellular processes such as endoplasmic reticulum (ER) translocation, ER-to-Golgi trafficking and vesicular transport and fusion. We showed that a subset of these interactors modified α-syn aggregation and α-syn-mediated degeneration of dopaminergic neurons in Caenorhabditis elegans, further suggesting that ATP13A2 and α-syn are functionally linked in neurodegeneration. These results implicate ATP13A2 in vesicular trafficking and provide a platform for further studies of ATP13A2 in neurodegeneration.

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.

The TSC22D, WNK, and NRBP gene families exhibit functional buffering and evolved with Metazoa for cell volume regulation.

The ability to sense and respond to osmotic fluctuations is critical for the maintenance of cellular integrity. We used gene co-essentiality analysis to identify an unappreciated relationship between TSC22D2, WNK1, and NRBP1 in regulating cell volume homeostasis. All of these genes have paralogs and are functionally buffered for osmo-sensing and cell volume control. Within seconds of hyperosmotic stress, TSC22D, WNK, and NRBP family members physically associate into biomolecular condensates, a process that is dependent on intrinsically disordered regions (IDRs). A close examination of these protein families across metazoans revealed that TSC22D genes evolved alongside a domain in NRBPs that specifically binds to TSC22D proteins, which we have termed NbrT (NRBP binding region with TSC22D), and this co-evolution is accompanied by rapid IDR length expansion in WNK-family kinases. Our study reveals that TSC22D, WNK, and NRBP genes evolved in metazoans to co-regulate rapid cell volume changes in response to osmolarity.

Survival-based CRISPR genetic screens across a panel of permissive cell lines identify common and cell-specific SARS-CoV-2 host factors.

SARS-CoV-2 depends on host cell components for infection and replication. Identification of virus-host dependencies offers an effective way to elucidate mechanisms involved in viral infection and replication. If druggable, host factor dependencies may present an attractive strategy for anti-viral therapy. In this study, we performed genome wide CRISPR knockout screens in Vero E6 cells and four human cell lines including Calu-3, UM-UC-4, HEK-293 and HuH-7 to identify genetic regulators of SARS-CoV-2 infection. Our findings identified only ACE2, the cognate SARS-CoV-2 entry receptor, as a common host dependency factor across all cell lines, while other host genes identified were largely cell line specific, including known factors TMPRSS2 and CTSL. Several of the discovered host-dependency factors converged on pathways involved in cell signalling, immune-related pathways, and chromatin modification. Notably, the chromatin modifier gene KMT2C in Calu-3 cells had the strongest impact in preventing SARS-CoV-2 infection when perturbed.