Development of quantitative high-throughput screening assays to identify, validate, and optimize small-molecule stabilizers of misfolded ß-glucocerebrosidase with therapeutic potential for Gaucher disease and Parkinson's disease.

Glucocerebrosidase (GCase) is implicated in both a rare, monogenic disorder (Gaucher disease, GD) and a common, multifactorial condition (Parkinson's disease); hence, it is an urgent therapeutic target. To identify correctors of severe protein misfolding and trafficking obstruction manifested by the pathogenic L444P-variant of GCase, we developed a suite of quantitative, high-throughput, cell-based assays. First, we labeled GCase with a small pro-luminescent HiBiT peptide reporter tag, enabling quantitation of protein stabilization in cells while faithfully maintaining target biology. TALEN-based gene editing allowed for stable integration of a single HiBiT-GBA1 transgene into an intragenic safe-harbor locus in GBA1-knockout H4 (neuroglioma) cells. This GD cell model was amenable to lead discovery via titration-based quantitative high-throughput screening and lead optimization via structure-activity relationships. A primary screen of 10,779 compounds from the NCATS bioactive collections identified 140 stabilizers of HiBiT-GCase-L444P, including both pharmacological chaperones (ambroxol and non-inhibitory chaperone NCGC326) and proteostasis regulators (panobinostat, trans-ISRIB, and pladienolide B). Two complementary high-content imaging-based assays were deployed to triage hits: the fluorescence-quenched substrate LysoFix-GBA captured functional lysosomal GCase activity, while an immunofluorescence assay featuring antibody hGCase-1/23 provided direct visualization of GCase lysosomal translocation. NCGC326 was active in both secondary assays and completely reversed pathological glucosylsphingosine accumulation. Finally, we tested the concept of combination therapy, by demonstrating synergistic actions of NCGC326 with proteostasis regulators in enhancing GCase-L444P levels. Looking forward, these physiologically-relevant assays can facilitate the identification, pharmacological validation, and medicinal chemistry optimization of new chemical matter targeting GCase, ultimately leading to a viable therapeutic for two protein-misfolding diseases.

Comparative chemical genomics in Babesia species identifies the alkaline phosphatase PhoD as a determinant of antiparasitic resistance.

Babesiosis is an emerging zoonosis and widely distributed veterinary infection caused by 100+ species of Babesia parasites. The diversity of Babesia parasites and the lack of specific drugs necessitate the discovery of broadly effective antibabesials. Here, we describe a comparative chemogenomics (CCG) pipeline for the identification of conserved targets. CCG relies on parallel in vitro evolution of resistance in independent populations of Babesia spp. (B. bovis and B. divergens). We identified a potent antibabesial, MMV019266, from the Malaria Box, and selected for resistance in two species of Babesia. After sequencing of multiple independently derived lines in the two species, we identified mutations in a membrane-bound metallodependent phosphatase (phoD). In both species, the mutations were found in the phoD-like phosphatase domain. Using reverse genetics, we validated that mutations in bdphoD confer resistance to MMV019266 in B. divergens. We have also demonstrated that BdPhoD localizes to the endomembrane system and partially with the apicoplast. Finally, conditional knockdown and constitutive overexpression of BdPhoD alter the sensitivity to MMV019266 in the parasite. Overexpression of BdPhoD results in increased sensitivity to the compound, while knockdown increases resistance, suggesting BdPhoD is a pro-susceptibility factor. Together, we have generated a robust pipeline for identification of resistance loci and identified BdPhoD as a resistance mechanism in Babesia species.

Multiplex single-cell chemical genomics reveals the kinase dependence of the response to targeted therapy.

Chemical genetic screens are a powerful tool for exploring how cancer cells' response to drugs is shaped by their mutations, yet they lack a molecular view of the contribution of individual genes to the response to exposure. Here, we present sci-Plex-Gene-by-Environment (sci-Plex-GxE), a platform for combined single-cell genetic and chemical screening at scale. We highlight the advantages of large-scale, unbiased screening by defining the contribution of each of 522 human kinases to the response of glioblastoma to different drugs designed to abrogate signaling from the receptor tyrosine kinase pathway. In total, we probed 14,121 gene-by-environment combinations across 1,052,205 single-cell transcriptomes. We identify an expression signature characteristic of compensatory adaptive signaling regulated in a MEK/MAPK-dependent manner. Further analyses aimed at preventing adaptation revealed promising combination therapies, including dual MEK and CDC7/CDK9 or nuclear factor κB (NF-κB) inhibitors, as potent means of preventing transcriptional adaptation of glioblastoma to targeted therapy.

Discovery and Characterization of Fluopipamine, a Putative Cellulose Synthase 1 Antagonist within Arabidopsis.

Herbicide-resistant weeds are increasingly a problem in crop fields when exposed to similar chemistry over time. To avoid future yield losses, identifying herbicidal chemistry needs to be accelerated. We screened 50,000 small molecules using a liquid-handling robot and light microscopy focusing on pre-emergent herbicides in the family of cellulose biosynthesis inhibitors. Through phenotypic, chemical, genetic, and in silico methods we uncovered 6-{[4-(2-fluorophenyl)-1-piperazinyl]methyl}-N-(2-methoxy-5-methylphenyl)-1,3,5-triazine-2,4-diamine (fluopipamine). Symptomologies support fluopipamine as a putative antagonist of cellulose synthase enzyme 1 (CESA1) from Arabidopsis (Arabidopsis thaliana). Ectopic lignification, inhibition of etiolation, phenotypes including loss of anisotropic cellular expansion, swollen roots, and live cell imaging link fluopipamine to cellulose biosynthesis inhibition. Radiolabeled glucose incorporation of cellulose decreased in short-duration experiments when seedlings were incubated in fluopipamine. To elucidate the mechanism, ethylmethanesulfonate mutagenized M2 seedlings were screened for fluopipamine resistance. Two loci of genetic resistance were linked to CESA1. In silico docking of fluopipamine, quinoxyphen, and flupoxam against various CESA1 mutations suggests that an alternative binding site at the interface between CESA proteins is necessary to preserve cellulose polymerization in compound presence. These data uncovered potential fundamental mechanisms of cellulose biosynthesis in plants along with feasible leads for herbicidal uses.

Randomly barcoded transposon mutant libraries for gut commensals II: Applying libraries for functional genetics.

The critical role of the intestinal microbiota in human health and disease is well recognized. Nevertheless, there are still large gaps in our understanding of the functions and mechanisms encoded in the genomes of most members of the gut microbiota. Genome-scale libraries of transposon mutants are a powerful tool to help us address this gap. Recent advances in barcoded transposon mutagenesis have dramatically lowered the cost of mutant fitness determination in hundreds of in vitro and in vivo experimental conditions. In an accompanying review, we discuss recent advances and caveats for the construction of pooled and arrayed barcoded transposon mutant libraries in human gut commensals. In this review, we discuss how these libraries can be used across a wide range of applications, the technical aspects involved, and expectations for such screens.

Colorectal Cancer Is Borrowing Blueprints from Intestinal Ontogenesis.

Colorectal tumors are heterogenous cellular systems harboring small populations of self-renewing and highly tumorigenic cancer stem cells (CSCs). Understanding the mechanisms fundamental to the emergence of CSCs and colorectal tumor initiation is crucial for developing effective therapeutic strategies. Two recent studies have highlighted the importance of developmental gene expression programs as potential therapeutic targets to suppress pro-oncogenic stem cell populations in the colonic epithelium. Specifically, a subset of aberrant stem cells was identified in preneoplastic intestinal lesions sharing significant transcriptional similarities with fetal gut development. In such aberrant stem cells, Sox9 was shown as a cornerstone for altered cell plasticity, the maintenance of premalignant stemness, and subsequent colorectal tumor initiation. Independently, chemical genomics was used to identify FDA-approved drugs capable of suppressing neoplastic self-renewal based on the ontogenetic root of a target tumor and transcriptional programs embedded in pluripotency. Here, we discuss the joint conclusions from these two approaches, underscoring the importance of developmental networks in CSCs as a novel paradigm for identifying therapeutics targeting colorectal cancer stemness.

Chemical genomics reveals targetable programs of human cancers rooted in pluripotency.

Overlapping principles of embryonic and tumor biology have been described, with recent multi-omics campaigns uncovering shared molecular profiles between human pluripotent stem cells (hPSCs) and adult tumors. Here, using a chemical genomic approach, we provide biological evidence that early germ layer fate decisions of hPSCs reveal targets of human cancers. Single-cell deconstruction of hPSCs-defined subsets that share transcriptional patterns with transformed adult tissues. Chemical screening using a unique germ layer specification assay for hPSCs identified drugs that enriched for compounds that selectively suppressed the growth of patient-derived tumors corresponding exclusively to their germ layer origin. Transcriptional response of hPSCs to germ layer inducing drugs could be used to identify targets capable of regulating hPSC specification as well as inhibiting adult tumors. Our study demonstrates properties of adult tumors converge with hPSCs drug induced differentiation in a germ layer specific manner, thereby expanding our understanding of cancer stemness and pluripotency.

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.

Comparative chemical genomic profiling across plant-based hydrolysate toxins reveals widespread antagonism in fitness contributions.

The budding yeast Saccharomyces cerevisiae has been used extensively in fermentative industrial processes, including biofuel production from sustainable plant-based hydrolysates. Myriad toxins and stressors found in hydrolysates inhibit microbial metabolism and product formation. Overcoming these stresses requires mitigation strategies that include strain engineering. To identify shared and divergent mechanisms of toxicity and to implicate gene targets for genetic engineering, we used a chemical genomic approach to study fitness effects across a library of S. cerevisiae deletion mutants cultured anaerobically in dozens of individual compounds found in different types of hydrolysates. Relationships in chemical genomic profiles identified classes of toxins that provoked similar cellular responses, spanning inhibitor relationships that were not expected from chemical classification. Our results also revealed widespread antagonistic effects across inhibitors, such that the same gene deletions were beneficial for surviving some toxins but detrimental for others. This work presents a rich dataset relating gene function to chemical compounds, which both expands our understanding of plant-based hydrolysates and provides a useful resource to identify engineering targets.

Genomic Approaches to Antifungal Drug Target Identification and Validation.

The last several decades have witnessed a surge in drug-resistant fungal infections that pose a serious threat to human health. While there is a limited arsenal of drugs that can be used to treat systemic infections, scientific advances have provided renewed optimism for the discovery of novel antifungals. The development of chemical-genomic assays using Saccharomyces cerevisiae has provided powerful methods to identify the mechanism of action of molecules in a living cell. Advances in molecular biology techniques have enabled complementary assays to be developed in fungal pathogens, including Candida albicans and Cryptococcus neoformans. These approaches enable the identification of target genes for drug candidates, as well as genes involved in buffering drug target pathways. Here, we examine yeast chemical-genomic assays and highlight how such resources can be utilized to predict the mechanisms of action of compounds, to study virulence attributes of diverse fungal pathogens, and to bolster the antifungal pipeline.

Global analysis of Saccharomyces cerevisiae growth in mucin.

Metagenomic profiling of the human gut microbiome has discovered DNA from dietary yeasts like Saccharomyces cerevisiae. However, it is unknown if the S. cerevisiae detected by common metagenomic methods are from dead dietary sources, or from live S. cerevisiae colonizing the gut similar to their close relative Candida albicans. While S. cerevisiae can adapt to minimal oxygen and acidic environments, it has not been explored whether this yeast can metabolize mucin, the large, gel-forming, highly glycosylated proteins representing a major source of carbon in the gut mucosa. We reveal that S. cerevisiae can utilize mucin as their main carbon source, as well as perform both a transcriptome analysis and a chemogenomic screen to identify biological pathways required for this yeast to grow optimally in mucin. In total, 739 genes demonstrate significant differential expression in mucin culture, and deletion of 21 genes impact growth in mucin. Both screens suggest that mitochondrial function is required for proper growth in mucin, and through secondary assays we determine that mucin exposure induces mitogenesis and cellular respiration. We further show that deletion of an uncharacterized ORF, YCR095W-A, led to dysfunction in mitochondrial morphology and oxygen consumption in mucin. Finally, we demonstrate that Yps7, an aspartyl protease and homolog to mucin-degrading proteins in C. albicans, is important for growth on mucin. Collectively, our work serves as the initial step toward establishing how this common dietary fungus can survive in the mucus environment of the human gut.

Discovery of Novel eEF2K Inhibitors Using HTS Fingerprint Generated from Predicted Profiling of Compound-Protein Interactions.

Background: Eukaryotic elongation factor 2 kinase (eEF2K) regulates the elongation stage of protein synthesis by phosphorylating eEF2, a process related to various diseases including cancer and cardiovascular and neurodegenerative diseases. In this study, we describe the identification of novel eEF2K inhibitors using high-throughput screening fingerprints (HTSFP) generated from predicted profiling of compound-protein interactions (CPIs). Methods: We utilized computationally generated HTSFPs referred to as chemical genomics-based fingerprint (CGBFP). Generally, HTSFPs are generated from multiple biochemical or cell-based assay data. On the other hand, CGBFPs are generated from computational prediction of CPIs using the Chemical Genomics-Based Virtual Screening (CGBVS) method. Therefore, CGBFPs do not have missing information mainly caused by the absence of assay data. Results: Chemogenomics-Based Similarity Profiling (CGBSP) of the screening library (2.6 million compounds) yielded 27 compounds which were evaluated for in vitro eEF2K inhibitory activity. Three compounds with interesting results were identified. Compounds 2 (IC50 = 11.05 µM) and 4 (IC50 = 43.54 µM) are thieno[2,3-b]pyridine derivatives that have the same scaffolds with a known eEF2K inhibitor, while compound 13 (IC50 = 70.13 µM) was a new thiophene-2-amine-type eEF2K inhibitor. Conclusions: CGBSP supplied an efficient strategy in the identification of novel eEF2K inhibitors and provided useful scaffolds for optimization.

Advances in fungal chemical genomics for the discovery of new antifungal agents.

Invasive fungal infections have escalated from a rare curiosity to a major cause of human mortality around the globe. This is in part due to a scarcity in the number of antifungal drugs available to combat mycotic disease, making the discovery of novel bioactive compounds and determining their mode of action of utmost importance. The development and application of chemical genomic assays using the model yeast Saccharomyces cerevisiae has provided powerful methods to identify the mechanism of action of diverse molecules in a living cell. Furthermore, complementary assays are continually being developed in fungal pathogens, most notably Candida albicans and Cryptococcus neoformans, to elucidate compound mechanism of action directly in the pathogen of interest. Collectively, the suite of chemical genetic assays that have been developed in multiple fungal species enables the identification of candidate drug target genes, as well as genes involved in buffering drug target pathways, and genes involved in general cellular responses to small molecules. In this review, we examine current yeast chemical genomic assays and highlight how such resources provide powerful tools that can be utilized to bolster the antifungal pipeline.

A chemical genomics-aggrephagy integrated method studying functional analysis of autophagy inducers.

Macroautophagy/autophagy plays a critical role in the pathogenesis of various human diseases including neurodegenerative disorders such as Parkinson disease (PD) and Huntington disease (HD). Chemical autophagy inducers are expected to serve as disease-modifying agents by eliminating cytotoxic/damaged proteins. Although many autophagy inducers have been identified, their precise molecular mechanisms are not fully understood because of the complicated crosstalk among signaling pathways. To address this issue, we performed several chemical genomic analyses enabling us to comprehend the dominancy among the autophagy-associated pathways followed by an aggresome-clearance assay. In a first step, more than 400 target-established small molecules were assessed for their ability to activate autophagic flux in neuronal PC12D cells, and we identified 39 compounds as autophagy inducers. We then profiled the autophagy inducers by testing their effect on the induction of autophagy by 200 well-established signal transduction modulators. Our principal component analysis (PCA) and clustering analysis using a dataset of "autophagy profiles" revealed that two Food and Drug Administration (FDA)-approved drugs, memantine and clemastine, activate endoplasmic reticulum (ER) stress responses, which could lead to autophagy induction. We also confirmed that SMK-17, a recently identified autophagy inducer, induced autophagy via the PRKC/PKC-TFEB pathway, as had been predicted from PCA. Finally, we showed that almost all of the autophagy inducers tested in this present work significantly enhanced the clearance of the protein aggregates observed in cellular models of PD and HD. These results, with the combined approach, suggested that autophagy-activating small molecules may improve proteinopathies by eliminating nonfunctional protein aggregates.Abbreviations: ADK: adenosine kinase; AMPK: AMP-activated protein kinase; ATF4: activating transcription factor 4; BECN1: beclin-1; DDIT3/CHOP: DNA damage inducible transcript 3; EIF2AK3/PERK: eukaryotic translation initiation factor 2 alpha kinase 3; EIF2S1/eIF2α: eukaryotic translation initiation factor 2 subunit alpha; ER: endoplasmic reticulum; ERN1/IRE1α: endoplasmic reticulum to nucleus signaling 1; FDA: Food and Drug Administration; GSH: glutathione; HD: Huntington disease; HSPA5/GRP78: heat shock protein family A (Hsp70) member 5; HTT: huntingtin; JAK: Janus kinase, MAP1LC3B/LC3: microtubule associated protein 1 light chain 3 beta; MAP2K/MEK: mitogen-activated protein kinase kinase; MAP3K8/Tpl2: mitogen-activated protein kinase kinase kinase 8; MAPK: mitogen-activated protein kinase; MPP+: 1-methyl-4-phenylpyridinium; MTOR: mechanistic target of rapamycin kinase; MTORC: MTOR complex; NAC: N-acetylcysteine; NGF: nerve growth factor 2; NMDA: N-methyl-D-aspartate; PCA: principal component analysis; PD: Parkinson disease; PDA: pancreatic ductal adenocarcinoma; PIK3C3: phosphatidylinositol 3-kinase catalytic subunit type 3; PMA: phorbol 12-myristate 13-acetate; PRKC/PKC: protein kinase C; ROCK: Rho-associated coiled-coil protein kinase; RR: ribonucleotide reductase; SIGMAR1: sigma non-opioid intracellular receptor 1; SQSTM1/p62: sequestosome 1; STK11/LKB1: serine/threonine kinase 11; TFEB: Transcription factor EB; TGFB/TGF-ß: Transforming growth factor beta; ULK1: unc-51 like autophagy activating kinase 1; XBP1: X-box binding protein 1.

In Vivo Chemical Genomics with Random Cyclized Peptides.

We have developed and applied a novel strategy that can best be described as in vivo chemical genomics, a concept where populations of any transformable organism may be screened for consequences of novel RNAs or peptides. We created a library of ~800,000 random DNA sequences biased only by third-position nucleotide substitutions that suppress the frequency of termination codons. The sequences may be shuttled to any plant, microbial, or animal expression vector with recombination cloning. We then generated large populations of Arabidopsis thaliana plants, each expressing a randomized DNA sequence, presumably giving rise to synthetic RNA species and/or the peptides they encode. These novel molecules are produced within the context of the cell and have been shown to affect plant biology with a relatively high frequency, as evidenced by diverse phenotypes. This chapter provides the protocols necessary to construct the libraries and isolate plants expressing randomized DNA sequences.

Identification of ABA Receptor Agonists Using a Multiplexed High-Throughput Chemical Screening.

Small molecules that can activate abscisic acid (ABA) receptors represent valuable probes to study ABA perception and signaling. Additionally, these compounds have the potential to be used in the field to counteract the negative effect of drought stress on plant productivity. The PYR/PYL ABA receptors, in their ligand-bound conformation, inactivate protein phosphatases 2C (PP2Cs), triggering physiological responses that are essential for plant adaptation to environmental stresses, including drought. Based on this ligand-induced PP2C inactivation mechanism, we have developed an in vitro assay for the identification of ABA-receptor agonists by high-throughput screening of chemical libraries. The assay allows simultaneous use of different ABA receptors, increasing the chances to find new agonists and eliminates the need for parallel screening. In this chapter, we describe detailed procedures for the identification of ABA agonists using this multiplexed assay in a medium- (96-well plates) or a high-throughput (384-well plates) setup.

A Luciferase Reporter Assay to Identify Chemical Activators of ABA Signaling.

Plant stress tolerance relies on intricate signaling networks that are not fully understood. Several plant hormones are involved in the adaptation to different environmental conditions. Abscisic acid (ABA) has an essential role in stress tolerance, especially in the adaptation to drought. During the last years, chemical genomics has gained attention as an alternative approach to decipher complex traits. Additionally, chemical-based strategies have been very useful to untangle genetic redundancy, which is hard to address by other approaches such as classical genetics. Here, we describe the use of an ABA-inducible luciferase (LUC) reporter line for the high-throughput identification of chemical activators of the ABA signaling pathway. In this assay, seven-day-old pMAPKKK18-LUC+ seedlings are grown on 96-well plates and treated with test compounds. Next, the activity of the LUC reporter is quantified semiautomatically by image analysis. Candidate compounds able to activate the reporter are thus identified and subjected to a secondary screen by analyzing their effect on ABA-related phenotypes (e.g., inhibition of seed germination). This assay is fast, high-throughput, nondestructive, semiquantitative and can be applied to any other luciferase reporter lines, making it ideal for forward chemical genetic screenings.

Tau Pathology Drives Dementia Risk-Associated Gene Networks toward Chronic Inflammatory States and Immunosuppression.

To understand how neural-immune-associated genes and pathways contribute to neurodegenerative disease pathophysiology, we performed a systematic functional genomic analysis in purified microglia and bulk tissue from mouse and human AD, FTD, and PSP. We uncover a complex temporal trajectory of microglial-immune pathways involving the type 1 interferon response associated with tau pathology in the early stages, followed by later signatures of partial immune suppression and, subsequently, the type 2 interferon response. We find that genetic risk for dementias shows disease-specific patterns of pathway enrichment. We identify drivers of two gene co-expression modules conserved from mouse to human, representing competing arms of microglial-immune activation (NAct) and suppression (NSupp) in neurodegeneration. We validate our findings by using chemogenetics, experimental perturbation data, and single-cell sequencing in post-mortem brains. Our results refine the understanding of stage- and disease-specific microglial responses, implicate microglial viral defense pathways in dementia pathophysiology, and highlight therapeutic windows.

6,11-Dioxobenzo[f]pyrido[1,2-a]indoles Kill Mycobacterium tuberculosis by Targeting Iron-Sulfur Protein Rv0338c (IspQ), A Putative Redox Sensor.

Screening of a diversity-oriented compound library led to the identification of two 6,11-dioxobenzo[f]pyrido[1,2-a]indoles (DBPI) that displayed low micromolar bactericidal activity against the Erdman strain of Mycobacterium tuberculosis in vitro. The activity of these hit compounds was limited to tubercle bacilli, including the nonreplicating form, and to Mycobacterium marinum. On hit expansion and investigation of the structure activity relationship, selected modifications to the dioxo moiety of the DBPI scaffold were either neutral or led to reduction or abolition of antimycobacterial activity. To find the target, DBPI-resistant mutants of M. tuberculosis Erdman were raised and characterized first microbiologically and then by whole genome sequencing. Four different mutations, all affecting highly conserved residues, were uncovered in the essential gene rv0338c (ispQ) that encodes a membrane-bound protein, named IspQ, with 2Fe-2S and 4Fe-4S centers and putative iron-sulfur-binding reductase activity. With the help of a structural model, two of the mutations were localized close to the 2Fe-2S domain in IspQ and another in transmembrane segment 3. The mutant genes were recessive to the wild type in complementation experiments and further confirmation of the hit-target relationship was obtained using a conditional knockdown mutant of rv0338c in M. tuberculosis H37Rv. More mechanistic insight was obtained from transcriptome analysis, following exposure of M. tuberculosis to two different DBPI; this revealed strong upregulation of the redox-sensitive SigK regulon and genes induced by oxidative and thiol-stress. The findings of this investigation pharmacologically validate a novel target in tubercle bacilli and open a new vista for tuberculosis drug discovery.

Seasonal changes in NRF2 antioxidant pathway regulates winter depression-like behavior.

Seasonal changes in the environment lead to depression-like behaviors in humans and animals. The underlying mechanisms, however, are unknown. We observed decreased sociability and increased anxiety-like behavior in medaka fish exposed to winter-like conditions. Whole brain metabolomic analysis revealed seasonal changes in 68 metabolites, including neurotransmitters and antioxidants associated with depression. Transcriptome analysis identified 3,306 differentially expressed transcripts, including inflammatory markers, melanopsins, and circadian clock genes. Further analyses revealed seasonal changes in multiple signaling pathways implicated in depression, including the nuclear factor erythroid-derived 2-like 2 (NRF2) antioxidant pathway. A broad-spectrum chemical screen revealed that celastrol (a traditional Chinese medicine) uniquely reversed winter behavior. NRF2 is a celastrol target expressed in the habenula (HB), known to play a critical role in the pathophysiology of depression. Another NRF2 chemical activator phenocopied these effects, and an NRF2 mutant showed decreased sociability. Our study provides important insights into winter depression and offers potential therapeutic targets involving NRF2.