Germline variation contributes to false negatives in CRISPR-based experiments with varying burden across ancestries.

Reducing disparities is vital for equitable access to precision treatments in cancer. Socioenvironmental factors are a major driver of disparities, but differences in genetic variation likely also contribute. The impact of genetic ancestry on prioritization of cancer targets in drug discovery pipelines has not been systematically explored due to the absence of pre-clinical data at the appropriate scale. Here, we analyze data from 611 genome-scale CRISPR/Cas9 viability experiments in human cell line models to identify ancestry-associated genetic dependencies essential for cell survival. Surprisingly, we find that most putative associations between ancestry and dependency arise from artifacts related to germline variants. Our analysis suggests that for 1.2-2.5% of guides, germline variants in sgRNA targeting sequences reduce cutting by the CRISPR/Cas9 nuclease, disproportionately affecting cell models derived from individuals of recent African descent. We propose three approaches to mitigate this experimental bias, enabling the scientific community to address these disparities.

TYK2 as a novel therapeutic target in Alzheimer's Disease with TDP-43 inclusions.

Neuroinflammation is a pathological feature of many neurodegenerative diseases, including Alzheimer's disease (AD)1,2 and amyotrophic lateral sclerosis (ALS)3, raising the possibility of common therapeutic targets. We previously established that cytoplasmic double-stranded RNA (cdsRNA) is spatially coincident with cytoplasmic pTDP-43 inclusions in neurons of patients with C9ORF72-mediated ALS4. CdsRNA triggers a type-I interferon (IFN-I)-based innate immune response in human neural cells, resulting in their death4. Here, we report that cdsRNA is also spatially coincident with pTDP-43 cytoplasmic inclusions in brain cells of patients with AD pathology and that type-I interferon response genes are significantly upregulated in brain regions affected by AD. We updated our machine-learning pipeline DRIAD-SP (Drug Repurposing In Alzheimer's Disease with Systems Pharmacology) to incorporate cryptic exon (CE) detection as a proxy of pTDP-43 inclusions and demonstrated that the FDA-approved JAK inhibitors baricitinib and ruxolitinib that block interferon signaling show a protective signal only in cortical brain regions expressing multiple CEs. Furthermore, the JAK family member TYK2 was a top hit in a CRISPR screen of cdsRNA-mediated death in differentiated human neural cells. The selective TYK2 inhibitor deucravacitinib, an FDA-approved drug for psoriasis, rescued toxicity elicited by cdsRNA. Finally, we identified CCL2, CXCL10, and IL-6 as candidate predictive biomarkers for cdsRNA-related neurodegenerative diseases. Together, we find parallel neuroinflammatory mechanisms between TDP-43 associated-AD and ALS and nominate TYK2 as a possible disease-modifying target of these incurable neurodegenerative diseases.

CRISPR-Cas9 screens reveal common essential miRNAs in human cancer cell lines.

BACKGROUND: Genome-wide functional screening using the CRISPR-Cas9 system is a powerful tool to uncover tumor-specific and common genetic dependencies across cancer cell lines. Current CRISPR-Cas9 knockout libraries, however, primarily target protein-coding genes. This limits functional genomics-based investigations of miRNA function. METHODS: We designed a novel CRISPR-Cas9 knockout library (lentiG-miR) of 8107 distinct sgRNAs targeting a total of 1769 human miRNAs and benchmarked its single guide RNA (sgRNA) composition, predicted on- and off-target activity, and screening performance against previous libraries. Using a total of 45 human cancer cell lines, representing 16 different tumor entities, we performed negative selection screens to identify miRNA fitness genes. Fitness miRNAs in each cell line were scored using a combination of supervised and unsupervised essentiality classifiers. Common essential miRNAs across distinct cancer cell lines were determined using the 90th percentile method. For subsequent validation, we performed knockout experiments for selected common essential miRNAs in distinct cancer cell lines and gene expression profiling. RESULTS: We found significantly lower off-target activity for protein-coding genes and a higher miRNA gene coverage for lentiG-miR as compared to previously described miRNA-targeting libraries, while preserving high on-target activity. A minor fraction of miRNAs displayed robust depletion of targeting sgRNAs, and we observed a high level of consistency between redundant sgRNAs targeting the same miRNA gene. Across 45 human cancer cell lines, only 217 (12%) of all targeted human miRNAs scored as a fitness gene in at least one model, and fitness effects for most miRNAs were confined to small subsets of cell lines. In contrast, we identified 49 common essential miRNAs with a homogenous fitness profile across the vast majority of all cell lines. Transcriptional profiling verified highly consistent gene expression changes in response to knockout of individual common essential miRNAs across a diverse set of cancer cell lines. CONCLUSIONS: Our study presents a miRNA-targeting CRISPR-Cas9 knockout library with high gene coverage and optimized on- and off-target activities. Taking advantage of the lentiG-miR library, we define a catalogue of miRNA fitness genes in human cancer cell lines, providing the foundation for further investigation of miRNAs in human cancer.

Salience signaling and stimulus scaling of ventral tegmental area glutamate neuron subtypes.

Ventral tegmental area (VTA) glutamatergic neurons participate in reward, aversion, drug-seeking, and stress. Subsets of VTA VGluT2+ neurons are capable of co-transmitting glutamate and GABA (VGluT2+VGaT+ neurons), transmitting glutamate without GABA (VGluT2+VGaT- neurons), or co-transmitting glutamate and dopamine (VGluT2+TH+ neurons), but whether these molecularly distinct subpopulations show behavior-related differences is not wholly understood. We identified that neuronal activity of each VGluT2+ subpopulation is sensitive to reward value but signaled this in different ways. The phasic maximum activity of VGluT2+VGaT+ neurons increased with sucrose concentration, whereas VGluT2+VGaT- neurons increased maximum and sustained activity with sucrose concentration, and VGluT2+TH+ neurons increased sustained but not maximum activity with sucrose concentration. Additionally, VGluT2+ subpopulations signaled consummatory preferences in different ways. VGluT2+VGaT- neurons and VGluT2+TH+ neurons showed a signaling preference for a behaviorally-preferred fat reward over sucrose, but in temporally-distinct ways. In contrast, VGluT2+VGaT+ neurons uniquely signaled a less behaviorally-preferred sucrose reward compared with fat. Further experiments suggested that VGluT2+VGaT+ consummatory reward-related activity was related to sweetness, partially modulated by hunger state, and not dependent on caloric content or behavioral preference. All VGluT2+ subtypes increased neuronal activity following aversive stimuli but VGluT2+VGaT+ neurons uniquely scaled their magnitude and sustained activity with footshock intensity. Optogenetic activation of VGluT2+VGaT+ neurons during low intensity footshock enhanced fear-related behavior without inducing place preference or aversion. We interpret these data such that VTA glutamatergic subpopulations signal different elements of rewarding and aversive experiences and highlight the unique role of VTA VGluT2+VGaT+ neurons in enhancing the salience of behavioral experiences.

Targeting TRIP13 in favorable histology Wilms tumor with nuclear export inhibitors synergizes with doxorubicin.

Wilms tumor (WT) is the most common renal malignancy of childhood. Despite improvements in the overall survival, relapse occurs in ~15% of patients with favorable histology WT (FHWT). Half of these patients will succumb to their disease. Identifying novel targeted therapies remains challenging in part due to the lack of faithful preclinical in vitro models. Here we establish twelve patient-derived WT cell lines and demonstrate that these models faithfully recapitulate WT biology using genomic and transcriptomic techniques. We then perform loss-of-function screens to identify the nuclear export gene, XPO1, as a vulnerability. We find that the FDA approved XPO1 inhibitor, KPT-330, suppresses TRIP13 expression, which is required for survival. We further identify synergy between KPT-330 and doxorubicin, a chemotherapy used in high-risk FHWT. Taken together, we identify XPO1 inhibition with KPT-330 as a potential therapeutic option to treat FHWTs and in combination with doxorubicin, leads to durable remissions in vivo.

Three million images and morphological profiles of cells treated with matched chemical and genetic perturbations.

The identification of genetic and chemical perturbations with similar impacts on cell morphology can elucidate compounds' mechanisms of action or novel regulators of genetic pathways. Research on methods for identifying such similarities has lagged due to a lack of carefully designed and well-annotated image sets of cells treated with chemical and genetic perturbations. Here we create such a Resource dataset, CPJUMP1, in which each perturbed gene's product is a known target of at least two chemical compounds in the dataset. We systematically explore the directionality of correlations among perturbations that target the same protein encoded by a given gene, and we find that identifying matches between chemical and genetic perturbations is a challenging task. Our dataset and baseline analyses provide a benchmark for evaluating methods that measure perturbation similarities and impact, and more generally, learn effective representations of cellular state from microscopy images. Such advancements would accelerate the applications of image-based profiling of cellular states, such as uncovering drug mode of action or probing functional genomics.

Comprehensive mutational scanning of EGFR reveals TKI sensitivities of extracellular domain mutants.

The epidermal growth factor receptor, EGFR, is frequently activated in lung cancer and glioblastoma by genomic alterations including missense mutations. The different mutation spectra in these diseases are reflected in divergent responses to EGFR inhibition: significant patient benefit in lung cancer, but limited in glioblastoma. Here, we report a comprehensive mutational analysis of EGFR function. We perform saturation mutagenesis of EGFR and assess function of ~22,500 variants in a human EGFR-dependent lung cancer cell line. This approach reveals enrichment of erlotinib-insensitive variants of known and unknown significance in the dimerization, transmembrane, and kinase domains. Multiple EGFR extracellular domain variants, not associated with approved targeted therapies, are sensitive to afatinib and dacomitinib in vitro. Two glioblastoma patients with somatic EGFR G598V dimerization domain mutations show responses to dacomitinib treatment followed by within-pathway resistance mutation in one case. In summary, this comprehensive screen expands the landscape of functional EGFR variants and suggests broader clinical investigation of EGFR inhibition for cancers harboring extracellular domain mutations.

Translation of non-canonical open reading frames as a cancer cell survival mechanism in childhood medulloblastoma.

A hallmark of high-risk childhood medulloblastoma is the dysregulation of RNA translation. Currently, it is unknown whether medulloblastoma dysregulates the translation of putatively oncogenic non-canonical open reading frames (ORFs). To address this question, we performed ribosome profiling of 32 medulloblastoma tissues and cell lines and observed widespread non-canonical ORF translation. We then developed a stepwise approach using multiple CRISPR-Cas9 screens to elucidate non-canonical ORFs and putative microproteins implicated in medulloblastoma cell survival. We determined that multiple lncRNA-ORFs and upstream ORFs (uORFs) exhibited selective functionality independent of main coding sequences. A microprotein encoded by one of these ORFs, ASNSD1-uORF or ASDURF, was upregulated, associated with MYC-family oncogenes, and promoted medulloblastoma cell survival through engagement with the prefoldin-like chaperone complex. Our findings underscore the fundamental importance of non-canonical ORF translation in medulloblastoma and provide a rationale to include these ORFs in future studies seeking to define new cancer targets.

Functional screening of amplification outlier oncogenes in organoid models of early tumorigenesis.

Somatic copy number gains are pervasive across cancer types, yet their roles in oncogenesis are insufficiently evaluated. This inadequacy is partly due to copy gains spanning large chromosomal regions, obscuring causal loci. Here, we employed organoid modeling to evaluate candidate oncogenic loci identified via integrative computational analysis of extreme copy gains overlapping with extreme expression dysregulation in The Cancer Genome Atlas. Subsets of "outlier" candidates were contextually screened as tissue-specific cDNA lentiviral libraries within cognate esophagus, oral cavity, colon, stomach, pancreas, and lung organoids bearing initial oncogenic mutations. Iterative analysis nominated the kinase DYRK2 at 12q15 as an amplified head and neck squamous carcinoma oncogene in p53-/- oral mucosal organoids. Similarly, FGF3, amplified at 11q13 in 41% of esophageal squamous carcinomas, promoted p53-/- esophageal organoid growth reversible by small molecule and soluble receptor antagonism of FGFRs. Our studies establish organoid-based contextual screening of candidate genomic drivers, enabling functional evaluation during early tumorigenesis.

Vps18 contributes to phagosome membrane integrity in Mycobacterium tuberculosis-infected macrophages.

Mycobacterium tuberculosis (Mtb) has evolved to be exquisitely adapted to survive within host macrophages. The capacity to damage the phagosomal membrane has emerged as central to Mtb virulence. While Mtb factors driving membrane damage have been described, host factors that repair that damage to contain the pathogen remain largely unknown. We used a genome-wide CRISPR screen to identify novel host factors required to repair Mtb-damaged phagosomal membranes. Vacuolar protein sorting-associated protein 18 (Vps18), a member of the HOPS and CORVET trafficking complexes, was among the top hits. Vps18 colocalized with Mtb in macrophages beginning shortly after infection, and Vps18-knockout macrophages demonstrated increased damage of Mtb-containing phagosomes without impaired autophagy. Mtb grew more robustly in Vps18-knockout cells, and the first-line anti-tuberculosis antibiotic pyrazinamide was less effective. Our results identify Vps18 as required for phagosomal membrane integrity in Mtb-infected cells and suggest that modulating phagosome integrity may hold promise for improving the efficacy of antibiotic treatment for TB.

Bed nucleus of the stria terminalis GABA neurons are necessary for changes in foraging behaviour following an innate threat.

Foraging is a universal behaviour that has co-evolved with predation pressure. We investigated the role of the bed nucleus of the stria terminalis (BNST) GABA neurons in robotic and live predator threat processing and their consequences in post-threat encounter foraging. Both robotic and live predator interactions increased BNST GABA neuron activity. Mice were trained to procure food in a laboratory-based foraging apparatus in which food pellets were placed at incrementally greater distances from a nest zone. After mice learned to forage, they were exposed to a robotic or live predator threat, while BNST GABA neurons were chemogenetically inhibited. Post-robotic threat encounter, mice spent more time in the nest zone, but other foraging parameters were unchanged compared with pre-encounter behaviour. Inhibition of BNST GABA neurons had no effect on foraging behaviour post-robotic threat encounter. Following live predator exposure, control mice spent significantly more time in the nest zone, increased their latency to successfully forage, and significantly altered their overall foraging performance. Inhibition of BNST GABA neurons during live predator exposure prevented changes in foraging behaviour from developing after a live predator threat. BNST GABA neuron inhibition did not alter foraging behaviour during robotic or live predator threats. We conclude that these results demonstrate that while both robotic and live predator encounters effectively intrude on foraging behaviour, the perceived risk and behavioural consequences of the threat are distinguishable. Additionally, BNST GABA neurons may play a role in the integration of prior innate predator threat experience that results in hypervigilance during post-encounter foraging behaviour.

Combining RNAscope and immunohistochemistry to visualize inflammatory gene products in neurons and microglia.

A challenge for central nervous system (CNS) tissue analysis in neuroscience research has been the difficulty to codetect and colocalize gene and protein expression in the same tissue. Given the importance of identifying gene expression relative to proteins of interest, for example, cell-type specific markers, we aimed to develop a protocol to optimize their codetection. RNAscope fluorescent in situ hybridization (FISH) combined with immunohistochemistry (IHC) in fixed (CNS) tissue sections allows for reliable quantification of gene transcripts of interest within IHC-labeled cells. This paper describes a new method for simultaneous visualization of FISH and IHC in thicker (14-µm), fixed tissue samples, using spinal cord sections. This method's effectiveness is shown by the cell-type-specific quantification of two genes, namely the proinflammatory cytokine interleukin-1beta (IL-1b) and the inflammasome NLR family pyrin domain containing 3 (NLRP3). These genes are challenging to measure accurately using immunohistochemistry (IHC) due to the nonspecificity of available antibodies and the hard-to-distinguish, dot-like visualizations of the labeled proteins within the tissue. These measurements were carried out in spinal cord sections after unilateral chronic constriction injury of the sciatic nerve to induce neuroinflammation in the spinal cord. RNAscope is used to label transcripts of genes of interest and IHC is used to label cell-type specific antigens (IBA1 for microglia, NeuN for neurons). This combination allowed for labeled RNA transcripts to be quantified within cell-type specific boundaries using confocal microscopy and standard image analysis methods. This method makes it easy to answer empirical questions that are intractable with standard IHC or in situ hybridization alone. The method, which has been optimized for spinal cord tissue and to minimize tissue preparation time and costs, is described in detail from tissue collection to image analysis. Further, the relative expression changes in inflammatory genes NLRP3 and IL-1b in spinal cord microglia vs. neurons of somatotopically relevant laminae are described for the first time.

Hyper-Dependence on NHEJ Enables Synergy between DNA-PK Inhibitors and Low-Dose Doxorubicin in Leiomyosarcoma.

PURPOSE: Leiomyosarcoma (LMS) is an aggressive sarcoma for which standard chemotherapies achieve response rates under 30%. There are no effective targeted therapies against LMS. Most LMS are characterized by chromosomal instability (CIN), resulting in part from TP53 and RB1 co-inactivation and DNA damage repair defects. We sought to identify therapeutic targets that could exacerbate intrinsic CIN and DNA damage in LMS, inducing lethal genotoxicity. EXPERIMENTAL DESIGN: We performed clinical targeted sequencing in 287 LMS and genome-wide loss-of-function screens in 3 patient-derived LMS cell lines, to identify LMS-specific dependencies. We validated candidate targets by biochemical and cell-response assays in vitro and in seven mouse models. RESULTS: Clinical targeted sequencing revealed a high burden of somatic copy-number alterations (median fraction of the genome altered =0.62) and demonstrated homologous recombination deficiency signatures in 35% of LMS. Genome-wide short hairpin RNA screens demonstrated PRKDC (DNA-PKcs) and RPA2 essentiality, consistent with compensatory nonhomologous end joining (NHEJ) hyper-dependence. DNA-PK inhibitor combinations with unconventionally low-dose doxorubicin had synergistic activity in LMS in vitro models. Combination therapy with peposertib and low-dose doxorubicin (standard or liposomal formulations) inhibited growth of 5 of 7 LMS mouse models without toxicity. CONCLUSIONS: Combinations of DNA-PK inhibitors with unconventionally low, sensitizing, doxorubicin dosing showed synergistic effects in LMS in vitro and in vivo models, without discernable toxicity. These findings underscore the relevance of DNA damage repair alterations in LMS pathogenesis and identify dependence on NHEJ as a clinically actionable vulnerability in LMS.

EZH2 Cooperates with BRD4-NUT to Drive NUT Carcinoma Growth by Silencing Key Tumor Suppressor Genes.

NUT carcinoma is an aggressive carcinoma driven by the BRD4-NUT fusion oncoprotein, which activates chromatin to promote expression of progrowth genes. BET bromodomain inhibitors (BETi) are a promising treatment for NUT carcinoma that can impede BRD4-NUT's ability to activate genes, but the efficacy of BETi as monotherapy is limited. Here, we demonstrated that enhancer of zeste homolog 2 (EZH2), which silences genes through establishment of repressive chromatin, is a dependency in NUT carcinoma. Inhibition of EZH2 with the clinical compound tazemetostat potently blocked growth of NUT carcinoma cells. Epigenetic and transcriptomic analysis revealed that tazemetostat reversed the EZH2-specific H3K27me3 silencing mark and restored expression of multiple tumor suppressor genes while having no effect on key oncogenic BRD4-NUT-regulated genes. Indeed, H3K27me3 and H3K27ac domains were found to be mutually exclusive in NUT carcinoma cells. CDKN2A was identified as the only gene among all tazemetostat-derepressed genes to confer resistance to tazemetostat in a CRISPR-Cas9 screen. Combined inhibition of EZH2 and BET synergized to downregulate cell proliferation genes, resulting in more pronounced growth arrest and differentiation than either inhibitor alone. In preclinical models, combined tazemetostat and BETi synergistically blocked tumor growth and prolonged survival of NUT carcinoma-xenografted mice, with complete remission without relapse in one cohort. Identification of EZH2 as a dependency in NUT carcinoma substantiates the reliance of NUT carcinoma tumor cells on epigenetic dysregulation of functionally opposite, yet highly complementary, chromatin regulatory pathways to maintain NUT carcinoma growth. SIGNIFICANCE: Repression of tumor suppressor genes, including CDKN2A, by EZH2 provides a mechanistic rationale for combining EZH2 and BET inhibitors for the clinical treatment of NUT carcinoma. See related commentary by Kazansky and Kentsis, p. 3827.

EZH2 synergizes with BRD4-NUT to drive NUT carcinoma growth through silencing of key tumor suppressor genes.

NUT carcinoma (NC) is an aggressive carcinoma driven by the BRD4-NUT fusion oncoprotein, which activates chromatin to promote expression of pro-growth genes. BET bromodomain inhibitors (BETi) impede BRD4-NUT's ability to activate genes and are thus a promising treatment but limited as monotherapy. The role of gene repression in NC is unknown. Here, we demonstrate that EZH2, which silences genes through establishment of repressive chromatin, is a dependency in NC. Inhibition of EZH2 with the clinical compound tazemetostat (taz) potently blocked growth of NC cells. Epigenetic and transcriptomic analysis revealed that taz reversed the EZH2-specific H3K27me3 silencing mark, and restored expression of multiple tumor suppressor genes while having no effect on key oncogenic BRD4- NUT-regulated genes. CDKN2A was identified as the only gene amongst all taz-derepressed genes to confer resistance to taz in a CRISPR-Cas9 screen. Combined EZH2 inhibition and BET inhibition synergized to downregulate cell proliferation genes resulting in more pronounced growth arrest and differentiation than either inhibitor alone. In pre-clinical models, combined taz and BETi synergistically blocked growth and prolonged survival of NC-xenografted mice, with all mice cured in one cohort. STATEMENT OF SIGNIFICANCE: Identification of EZH2 as a dependency in NC substantiates the reliance of NC tumor cells on epigenetic dysregulation of functionally opposite, yet highly complementary chromatin regulatory pathways to maintain NC growth. In particular, repression of CDKN2A expression by EZH2 provides a mechanistic rationale for combining EZH2i with BETi for the clinical treatment of NC.

Partial gene suppression improves identification of cancer vulnerabilities when CRISPR-Cas9 knockout is pan-lethal.

BACKGROUND: Hundreds of functional genomic screens have been performed across a diverse set of cancer contexts, as part of efforts such as the Cancer Dependency Map, to identify gene dependencies-genes whose loss of function reduces cell viability or fitness. Recently, large-scale screening efforts have shifted from RNAi to CRISPR-Cas9, due to superior efficacy and specificity. However, many effective oncology drugs only partially inhibit their protein targets, leading us to question whether partial suppression of genes using RNAi could reveal cancer vulnerabilities that are missed by complete knockout using CRISPR-Cas9. Here, we compare CRISPR-Cas9 and RNAi dependency profiles of genes across approximately 400 matched cancer cell lines. RESULTS: We find that CRISPR screens accurately identify more gene dependencies per cell line, but the majority of each cell line's dependencies are part of a set of 1867 genes that are shared dependencies across the entire collection (pan-lethals). While RNAi knockdown of about 30% of these genes is also pan-lethal, approximately 50% have selective dependency patterns across cell lines, suggesting they could still be cancer vulnerabilities. The accuracy of the unique RNAi selectivity is supported by associations to multi-omics profiles, drug sensitivity, and other expected co-dependencies. CONCLUSIONS: Incorporating RNAi data for genes that are pan-lethal knockouts facilitates the discovery of a wider range of gene targets than could be detected using the CRISPR dataset alone. This can aid in the interpretation of contrasting results obtained from CRISPR and RNAi screens and reinforce the importance of partial gene suppression methods in building a cancer dependency map.

Ventral tegmental area glutamate neurons establish a mu-opioid receptor gated circuit to mesolimbic dopamine neurons and regulate opioid-seeking behavior.

A two-neuron model of ventral tegmental area (VTA) opioid function classically involves VTA GABA neuron regulation of VTA dopamine neurons via a mu-opioid receptor dependent inhibitory circuit. However, this model predates the discovery of a third major type of neuron in the VTA: glutamatergic neurons. We found that about one-quarter of VTA neurons expressing the mu-opioid receptor are glutamate neurons without molecular markers of GABA co-release. Glutamate-Mu opioid receptor neurons are largely distributed in the anterior VTA. The majority of remaining VTA mu-opioid receptor neurons are GABAergic neurons that are mostly within the posterior VTA and do not express molecular markers of glutamate co-release. Optogenetic stimulation of VTA glutamate neurons resulted in excitatory currents recorded from VTA dopamine neurons that were reduced by presynaptic activation of the mu-opioid receptor ex vivo, establishing a local mu-opioid receptor dependent excitatory circuit from VTA glutamate neurons to VTA dopamine neurons. This VTA glutamate to VTA dopamine pathway regulated dopamine release to the nucleus accumbens through mu-opioid receptor activity in vivo. Behaviorally, VTA glutamate calcium-related neuronal activity increased following oral oxycodone consumption during self-administration and response-contingent oxycodone-associated cues during abstinent reinstatement of drug-seeking behavior. Further, chemogenetic inhibition of VTA glutamate neurons reduced abstinent oral oxycodone-seeking behavior in male but not female mice. These results establish 1) a three-neuron model of VTA opioid function involving a mu-opioid receptor gated VTA glutamate neuron pathway to VTA dopamine neurons that controls dopamine release within the nucleus accumbens, and 2) that VTA glutamate neurons participate in opioid-seeking behavior.

Genome-wide pooled CRISPR screening in neurospheres.

Spheroid culture systems have allowed in vitro propagation of cells unable to grow in canonical cell culturing conditions, and may capture cellular contexts that model tumor growth better than current model systems. The insights gleaned from genome-wide clustered regularly interspaced short palindromic repeat (CRISPR) screening of thousands of cancer cell lines grown in conventional culture conditions illustrate the value of such CRISPR pooled screens. It is clear that similar genome-wide CRISPR screens of three-dimensional spheroid cultures will be important for future biological discovery. Here, we present a protocol for genome-wide CRISPR screening of three-dimensional neurospheres. While many in-depth protocols and discussions have been published for more typical cell lines, few detailed protocols are currently available in the literature for genome-wide screening in spheroidal cell lines. For those who want to screen such cell lines, and particularly neurospheres, we provide a step-by-step description of assay development tests to be performed before screening, as well as for the screen itself. We highlight considerations of variables that make these screens distinct from, or similar to, typical nonspheroid cell lines throughout. Finally, we illustrate typical outcomes of neurosphere genome-wide screens, and how neurosphere screens typically produce slightly more heterogeneous signal distributions than more canonical cancer cell lines. Completion of this entire protocol will take 8-12 weeks from the initial assay development tests to deconvolution of the sequencing data.