Total loss of VHL gene function impairs neuroendocrine cancer cell fitness due to excessive HIF2α activity.

Loss-of-function germline von Hippel-Lindau (VHL) tumor suppressor mutations cause VHL disease, which predisposes individuals to kidney cancer, hemangioblastomas, and paragangliomas. The risk that a given VHL disease family will manifest some or all these tumor types is profoundly influenced by the VHL allele it carries. For example, almost all VHL disease families that develop paraganglioma have missense VHL mutations. VHL families with null VHL alleles develop kidney cancer and hemangioblastomas without a high risk of paraganglioma. The latter is surprising because the VHL gene product, pVHL, suppresses the HIF2 transcription factor and gain-of-function HIF2 mutations are also linked to paraganglioma. Paragangliomas arise from the sympathetic or parasympathetic nervous system. Given the lack of human paraganglioma cell lines, we studied the effects of inactivating VHL in neuroblastoma cell lines, which also arise from the sympathetic nervous system. We found that total loss of pVHL function profoundly impairs the fitness of neuroblastoma cell lines in a HIF2-dependent manner both ex vivo and in vivo. This fitness defect can be rescued by pVHL variants linked to paraganglioma, but not by pVHL variants associated with a low risk of paraganglioma. These findings suggest that HIF2 activity above a critical threshold prevents the development of paraganglioma.

Monosynaptic inputs to ventral tegmental area glutamate and GABA co-transmitting neurons.

A unique population of ventral tegmental area (VTA) neurons co-transmits glutamate and GABA. However, the circuit inputs to VTA VGluT2+VGaT+ neurons are unknown, limiting our understanding of their functional capabilities. By coupling monosynaptic rabies tracing with intersectional genetic targeting in male and female mice, we found that VTA VGluT2+VGaT+ neurons received diverse brain-wide inputs. The largest numbers of monosynaptic inputs to VTA VGluT2+VGaT+ neurons were from superior colliculus, lateral hypothalamus, midbrain reticular nucleus, and periaqueductal gray, whereas the densest inputs relative to brain region volume were from dorsal raphe nucleus, lateral habenula, and ventral tegmental area. Based on these and prior data, we hypothesized that lateral hypothalamus and superior colliculus inputs were glutamatergic neurons. Optical activation of glutamatergic lateral hypothalamus neurons activated VTA VGluT2+VGaT+ neurons regardless of stimulation frequency and resulted in flee-like ambulatory behavior. In contrast, optical activation of glutamatergic superior colliculus neurons activated VTA VGluT2+VGaT+ neurons for a brief period of time at high stimulation frequency and resulted in head rotation and arrested ambulatory behavior (freezing). Stimulation of glutamatergic lateral hypothalamus neurons, but not glutamatergic superior colliculus neurons, was associated with VTA VGluT2+VGaT+ footshock-induced activity. In addition, inhibition of lateral hypothalamus glutamatergic neurons disrupted VTA VGluT2+VGaT+ tailshock-induced activity. We interpret these results such that inputs to VTA VGluT2+VGaT+ neurons may integrate diverse signals related to the detection and processing of motivationally-salient outcomes.Significance Statement VTA glutamate neurons have roles in motivated behavior and unique neurotransmission capabilities. A specific VTA glutamate neuron subtype, those that co-transmit glutamate and GABA, have unique outcome signaling properties compared to other VTA cell-types. However, the circuits that regulate these neurons are unclear. We identified the whole-brain inputs to VTA glutamate and GABA co-transmitting neurons. We also identify two distinct glutamatergic inputs that activate VTA glutamate and GABA co-transmitting neurons and result in different behavioral repertoires suggestive of threat processing. Together, these results provide novel insights into the circuit and cell-type specific influences on VTA glutamate and GABA co-transmitting neuronal activity as integrators of motivationally salient outcomes.

Genome-wide CRISPR-Cas9 knockout screens identify DNMT1 as a druggable dependency in sonic hedgehog medulloblastoma.

Sonic hedgehog subgroup of medulloblastoma (SHH-MB) is characterized by aberrant activation of the SHH signaling pathway. An inhibition of the positive SHH regulator Smoothened (SMO) has demonstrated promising clinical efficacy. Yet, primary and acquired resistance to SMO inhibitors limit their efficacy. An understanding of underlying molecular mechanisms of resistance to therapy is warranted to bridge this unmet need. Here, we make use of genome-wide CRISPR-Cas9 knockout screens in murine SMB21 and human DAOY cells, in order to unravel genetic dependencies and drug-related genetic interactors that could serve as alternative therapeutic targets for SHH-MB. Our screens reinforce SMB21 cells as a faithful model system for SHH-MB, as opposed to DAOY cells, and identify members of the epigenetic machinery including DNA methyltransferase 1 (DNMT1) as druggable targets in SHH-dependent tumors. We show that Dnmt1 plays a crucial role in normal murine cerebellar development and is required for SHH-MB growth in vivo. Additionally, DNMT1 pharmacological inhibition alone and in combination with SMO inhibition effectively inhibits tumor growth in murine and human SHH-MB cell models and prolongs survival of SHH-MB mouse models by inhibiting SHH signaling output downstream of SMO. In conclusion, our data highlight the potential of inhibiting epigenetic regulators as a novel therapeutic avenue in SMO-inhibitor sensitive as well as resistant SHH-MBs.

Targeting Cholesterol Biosynthesis with Statins Synergizes with AKT Inhibitors in Triple-Negative Breast Cancer.

Triple-negative breast cancer (TNBC) is responsible for a disproportionate number of breast cancer patient deaths due to extensive molecular heterogeneity, high recurrence rates and lack of targeted therapies. Dysregulation of the phosphoinositide 3-kinase (PI3K)/AKT pathway occurs in approximately 50% of TNBC patients. Here, we performed a genome-wide CRISPR/Cas9 screen with PI3Kα and AKT inhibitors to find targetable synthetic lethalities in TNBC. Cholesterol homeostasis was identified as a collateral vulnerability with AKT inhibition. Disruption of cholesterol homeostasis with pitavastatin synergized with AKT inhibition to induce TNBC cytotoxicity in vitro, in mouse TNBC xenografts and in patient-derived, estrogen receptor (ER)-negative breast cancer organoids. Neither ER-positive breast cancer cell lines nor ER-positive organoids were sensitive to combined AKT inhibitor and pitavastatin. Mechanistically, TNBC cells showed impaired sterol regulatory element-binding protein 2 (SREBP-2) activation in response to single agent or combination treatment with AKT inhibitor and pitavastatin, which was rescued by inhibition of the cholesterol trafficking protein Niemann-Pick C1 (NPC1). NPC1 loss caused lysosomal cholesterol accumulation, decreased endoplasmic reticulum cholesterol levels, and promoted SREBP-2 activation. Taken together, these data identify a TNBC-specific vulnerability to the combination of AKT inhibitors and pitavastatin mediated by dysregulated cholesterol trafficking. These findings support combining AKT inhibitors with pitavastatin as a therapeutic modality in TNBC. .

Generation of a biliary tract cancer cell line atlas reveals molecular subtypes and therapeutic targets.

Biliary tract cancers (BTCs) are a group of deadly malignancies encompassing intrahepatic and extrahepatic cholangiocarcinoma, gallbladder carcinoma, and ampullary carcinoma. Here, we present the integrative analysis of 63 BTC cell lines via multi-omics clustering and genome- scale CRISPR screens, providing a platform to illuminate BTC biology and inform therapeutic development. We identify dependencies broadly enriched in BTC compared to other cancers as well as dependencies selective to the anatomic subtypes. Notably, cholangiocarcinoma cell lines are stratified into distinct lineage subtypes based on biliary or dual biliary/hepatocyte marker signatures, associated with dependency on specific lineage survival factors. Transcriptional analysis of patient specimens demonstrates the prognostic significance of these lineage subtypes. Additionally, we delineate strategies to enhance targeted therapies or to overcome resistance in cell lines with key driver gene mutations. Furthermore, clustering based on dependencies and proteomics data elucidates unexpected functional relationships, including a BTC subgroup with partial squamous differentiation. Thus, this cell line atlas reveals potential therapeutic targets in molecularly defined BTCs, unveils biologically distinct disease subtypes, and offers a vital resource for BTC research.

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.

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.

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.

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.

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.

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.

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.

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.