SPA-STOCSY: an automated tool for identifying annotated and non-annotated metabolites in high-throughput NMR spectra.

MOTIVATION: Nuclear magnetic resonance spectroscopy (NMR) is widely used to analyze metabolites in biological samples, but the analysis requires specific expertise, it is time-consuming, and can be inaccurate. Here, we present a powerful automate tool, SPatial clustering Algorithm-Statistical TOtal Correlation SpectroscopY (SPA-STOCSY), which overcomes challenges faced when analyzing NMR data and identifies metabolites in a sample with high accuracy. RESULTS: As a data-driven method, SPA-STOCSY estimates all parameters from the input dataset. It first investigates the covariance pattern among datapoints and then calculates the optimal threshold with which to cluster datapoints belonging to the same structural unit, i.e. the metabolite. Generated clusters are then automatically linked to a metabolite library to identify candidates. To assess SPA-STOCSY's efficiency and accuracy, we applied it to synthesized spectra and spectra acquired on Drosophila melanogaster tissue and human embryonic stem cells. In the synthesized spectra, SPA outperformed Statistical Recoupling of Variables (SRV), an existing method for clustering spectral peaks, by capturing a higher percentage of the signal regions and the close-to-zero noise regions. In the biological data, SPA-STOCSY performed comparably to the operator-based Chenomx analysis while avoiding operator bias, and it required <7 min of total computation time. Overall, SPA-STOCSY is a fast, accurate, and unbiased tool for untargeted analysis of metabolites in the NMR spectra. It may thus accelerate the use of NMR for scientific discoveries, medical diagnostics, and patient-specific decision making. AVAILABILITY AND IMPLEMENTATION: The codes of SPA-STOCSY are available at https://github.com/LiuzLab/SPA-STOCSY.

Functional screening of lysosomal storage disorder genes identifies modifiers of alpha-synuclein neurotoxicity.

Heterozygous variants in the glucocerebrosidase (GBA) gene are common and potent risk factors for Parkinson's disease (PD). GBA also causes the autosomal recessive lysosomal storage disorder (LSD), Gaucher disease, and emerging evidence from human genetics implicates many other LSD genes in PD susceptibility. We have systemically tested 86 conserved fly homologs of 37 human LSD genes for requirements in the aging adult Drosophila brain and for potential genetic interactions with neurodegeneration caused by α-synuclein (αSyn), which forms Lewy body pathology in PD. Our screen identifies 15 genetic enhancers of αSyn-induced progressive locomotor dysfunction, including knockdown of fly homologs of GBA and other LSD genes with independent support as PD susceptibility factors from human genetics (SCARB2, SMPD1, CTSD, GNPTAB, SLC17A5). For several genes, results from multiple alleles suggest dose-sensitivity and context-dependent pleiotropy in the presence or absence of αSyn. Homologs of two genes causing cholesterol storage disorders, Npc1a / NPC1 and Lip4 / LIPA, were independently confirmed as loss-of-function enhancers of αSyn-induced retinal degeneration. The enzymes encoded by several modifier genes are upregulated in αSyn transgenic flies, based on unbiased proteomics, revealing a possible, albeit ineffective, compensatory response. Overall, our results reinforce the important role of lysosomal genes in brain health and PD pathogenesis, and implicate several metabolic pathways, including cholesterol homeostasis, in αSyn-mediated neurotoxicity.

Functional variants identify sex-specific genes and pathways in Alzheimer's Disease.

The incidence of Alzheimer's Disease in females is almost double that of males. To search for sex-specific gene associations, we build a machine learning approach focused on functionally impactful coding variants. This method can detect differences between sequenced cases and controls in small cohorts. In the Alzheimer's Disease Sequencing Project with mixed sexes, this approach identified genes enriched for immune response pathways. After sex-separation, genes become specifically enriched for stress-response pathways in male and cell-cycle pathways in female. These genes improve disease risk prediction in silico and modulate Drosophila neurodegeneration in vivo. Thus, a general approach for machine learning on functionally impactful variants can uncover sex-specific candidates towards diagnostic biomarkers and therapeutic targets.

Tau polarizes an aging transcriptional signature to excitatory neurons and glia.

Aging is a major risk factor for Alzheimer's disease (AD), and cell-type vulnerability underlies its characteristic clinical manifestations. We have performed longitudinal, single-cell RNA-sequencing in Drosophila with pan-neuronal expression of human tau, which forms AD neurofibrillary tangle pathology. Whereas tau- and aging-induced gene expression strongly overlap (93%), they differ in the affected cell types. In contrast to the broad impact of aging, tau-triggered changes are strongly polarized to excitatory neurons and glia. Further, tau can either activate or suppress innate immune gene expression signatures in a cell-type-specific manner. Integration of cellular abundance and gene expression pinpoints nuclear factor kappa B signaling in neurons as a marker for cellular vulnerability. We also highlight the conservation of cell-type-specific transcriptional patterns between Drosophila and human postmortem brain tissue. Overall, our results create a resource for dissection of dynamic, age-dependent gene expression changes at cellular resolution in a genetically tractable model of tauopathy.

Upregulation of the ESCRT pathway and multivesicular bodies accelerates degradation of proteins associated with neurodegeneration.

Many neurodegenerative diseases, including Huntington's disease (HD) and Alzheimer's disease (AD), occur due to an accumulation of aggregation-prone proteins, which results in neuronal death. Studies in animal and cell models show that reducing the levels of these proteins mitigates disease phenotypes. We previously reported a small molecule, NCT-504, which reduces cellular levels of mutant huntingtin (mHTT) in patient fibroblasts as well as mouse striatal and cortical neurons from an HdhQ111 mutant mouse. Here, we show that NCT-504 has a broader potential, and in addition reduces levels of Tau, a protein associated with Alzheimer's disease, as well as other tauopathies. We find that in untreated cells, Tau and mHTT are degraded via autophagy. Notably, treatment with NCT-504 diverts these proteins to multivesicular bodies (MVB) and the ESCRT pathway. Specifically, NCT-504 causes a proliferation of endolysosomal organelles including MVB, and an enhanced association of mHTT and Tau with endosomes and MVB. Importantly, depletion of proteins that act late in the ESCRT pathway blocked NCT-504 dependent degradation of Tau. Moreover, NCT-504-mediated degradation of Tau occurred in cells where Atg7 is depleted, which indicates that this pathway is independent of canonical autophagy. Together, these studies reveal that upregulation of traffic through an ESCRT-dependent MVB pathway may provide a therapeutic approach for neurodegenerative diseases.

SPA-STOCSY: An Automated Tool for Identification of Annotated and Non-Annotated Metabolites in High-Throughput NMR Spectra.

Nuclear Magnetic Resonance (NMR) spectroscopy is widely used to analyze metabolites in biological samples, but the analysis can be cumbersome and inaccurate. Here, we present a powerful automated tool, SPA-STOCSY (Spatial Clustering Algorithm - Statistical Total Correlation Spectroscopy), which overcomes the challenges by identifying metabolites in each sample with high accuracy. As a data-driven method, SPA-STOCSY estimates all parameters from the input dataset, first investigating the covariance pattern and then calculating the optimal threshold with which to cluster data points belonging to the same structural unit, i.e. metabolite. The generated clusters are then automatically linked to a compound library to identify candidates. To assess SPA-STOCSY’s efficiency and accuracy, we applied it to synthesized and real NMR data obtained from Drosophila melanogaster brains and human embryonic stem cells. In the synthesized spectra, SPA outperforms Statistical Recoupling of Variables, an existing method for clustering spectral peaks, by capturing a higher percentage of the signal regions and the close-to-zero noise regions. In the real spectra, SPA-STOCSY performs comparably to operator-based Chenomx analysis but avoids operator bias and performs the analyses in less than seven minutes of total computation time. Overall, SPA-STOCSY is a fast, accurate, and unbiased tool for untargeted analysis of metabolites in the NMR spectra. As such, it might accelerate the utilization of NMR for scientific discoveries, medical diagnostics, and patient-specific decision making.

Evolutionarily conserved regulators of tau identify targets for new therapies.

Tauopathies are neurodegenerative diseases that involve the pathological accumulation of tau proteins; in this family are Alzheimer disease, corticobasal degeneration, and chronic traumatic encephalopathy, among others. Hypothesizing that reducing this accumulation could mitigate pathogenesis, we performed a cross-species genetic screen targeting 6,600 potentially druggable genes in human cells and Drosophila. We found and validated 83 hits in cells and further validated 11 hits in the mouse brain. Three of these hits (USP7, RNF130, and RNF149) converge on the C terminus of Hsc70-interacting protein (CHIP) to regulate tau levels, highlighting the role of CHIP in maintaining tau proteostasis in the brain. Knockdown of each of these three genes in adult tauopathy mice reduced tau levels and rescued the disease phenotypes. This study thus identifies several points of intervention to reduce tau levels and demonstrates that reduction of tau levels via regulation of this pathway is a viable therapeutic strategy for Alzheimer disease and other tauopathies.

An altered extracellular matrix-integrin interface contributes to Huntington's disease-associated CNS dysfunction in glial and vascular cells.

Astrocytes and brain endothelial cells are components of the neurovascular unit that comprises the blood-brain barrier (BBB) and their dysfunction contributes to pathogenesis in Huntington's disease (HD). Defining the contribution of these cells to disease can inform cell-type-specific effects and uncover new disease-modifying therapeutic targets. These cells express integrin (ITG) adhesion receptors that anchor the cells to the extracellular matrix (ECM) to maintain the integrity of the BBB. We used HD patient-derived induced pluripotent stem cell (iPSC) modeling to study the ECM-ITG interface in astrocytes and brain microvascular endothelial cells and found ECM-ITG dysregulation in human iPSC-derived cells that may contribute to the dysfunction of the BBB in HD. This disruption has functional consequences since reducing ITG expression in glia in an HD Drosophila model suppressed disease-associated CNS dysfunction. Since ITGs can be targeted therapeutically and manipulating ITG signaling prevents neurodegeneration in other diseases, defining the role of ITGs in HD may provide a novel strategy of intervention to slow CNS pathophysiology to treat HD.

Integration of transcriptome-wide association study with neuronal dysfunction assays provides functional genomics evidence for Parkinson's disease genes.

Genome-wide association studies (GWAS) have markedly advanced our understanding of the genetics of Parkinson's disease (PD), but they currently do not account for the full heritability of PD. In many cases it is difficult to unambiguously identify a specific gene within each locus because GWAS does not provide functional information on the identified candidate loci. Here we present an integrative approach that combines transcriptome-wide association study (TWAS) with high-throughput neuronal dysfunction analyses in Drosophila to discover and validate candidate PD genes. We identified 160 candidate genes whose misexpression is associated with PD risk via TWAS. Candidates were validated using orthogonal in silico methods and found to be functionally related to PD-associated pathways (i.e. endolysosome). We then mimicked these TWAS-predicted transcriptomic alterations in a Drosophila PD model and discovered that 50 candidates can modulate α-Synuclein(α-Syn)-induced neurodegeneration, allowing us to nominate new genes in previously known PD loci. We also uncovered additional novel PD candidate genes within GWAS suggestive loci (e.g. TTC19, ADORA2B, LZTS3, NRBP1, HN1L), which are also supported by clinical and functional evidence. These findings deepen our understanding of PD, and support applying our integrative approach to other complex trait disorders.

Identification of risk genes for Alzheimer's disease by gene embedding.

Most disease-gene association methods do not account for gene-gene interactions, even though these play a crucial role in complex, polygenic diseases like Alzheimer's disease (AD). To discover new genes whose interactions may contribute to pathology, we introduce GeneEMBED. This approach compares the functional perturbations induced in gene interaction network neighborhoods by coding variants from disease versus healthy subjects. In two independent AD cohorts of 5,169 exomes and 969 genomes, GeneEMBED identified novel candidates. These genes were differentially expressed in post mortem AD brains and modulated neurological phenotypes in mice. Four that were differentially overexpressed and modified neurodegeneration in vivo are PLEC, UTRN, TP53, and POLD1. Notably, TP53 and POLD1 are involved in DNA break repair and inhibited by approved drugs. While these data show proof of concept in AD, GeneEMBED is a general approach that should be broadly applicable to identify genes relevant to risk mechanisms and therapy of other complex diseases.

Cross-species genetic screens identify transglutaminase 5 as a regulator of polyglutamine-expanded ataxin-1.

Many neurodegenerative disorders are caused by abnormal accumulation of misfolded proteins. In spinocerebellar ataxia type 1 (SCA1), accumulation of polyglutamine-expanded (polyQ-expanded) ataxin-1 (ATXN1) causes neuronal toxicity. Lowering total ATXN1, especially the polyQ-expanded form, alleviates disease phenotypes in mice, but the molecular mechanism by which the mutant ATXN1 is specifically modulated is not understood. Here, we identified 22 mutant ATXN1 regulators by performing a cross-species screen of 7787 and 2144 genes in human cells and Drosophila eyes, respectively. Among them, transglutaminase 5 (TG5) preferentially regulated mutant ATXN1 over the WT protein. TG enzymes catalyzed cross-linking of ATXN1 in a polyQ-length-dependent manner, thereby preferentially modulating mutant ATXN1 stability and oligomerization. Perturbing Tg in Drosophila SCA1 models modulated mutant ATXN1 toxicity. Moreover, TG5 was enriched in the nuclei of SCA1-affected neurons and colocalized with nuclear ATXN1 inclusions in brain tissue from patients with SCA1. Our work provides a molecular insight into SCA1 pathogenesis and an opportunity for allele-specific targeting for neurodegenerative disorders.

Suppression of toxicity of the mutant huntingtin protein by its interacting compound, desonide.

Identifying inhibitors of pathogenic proteins is the major strategy of targeted drug discoveries. This strategy meets challenges in targeting neurodegenerative disorders such as Huntington's disease (HD), which is mainly caused by the mutant huntingtin protein (mHTT), an "undruggable" pathogenic protein with unknown functions. We hypothesized that some of the chemical binders of mHTT may change its conformation and/or stability to suppress its downstream toxicity, functioning similarly to an "inhibitor" under a broader definition. We identified 21 potential mHTT selective binders through a small-molecule microarray­based screening. We further tested these compounds using secondary phenotypic screens for their effects on mHTT-induced toxicity and revealed four potential mHTT-binding compounds that may rescue HD-relevant phenotypes. Among them, a Food and Drug Administration­approved drug, desonide, was capable of suppressing mHTT toxicity in HD cellular and animal models by destabilizing mHTT through enhancing its polyubiquitination at the K6 site. Our study reveals the therapeutic potential of desonide for HD treatment and provides the proof of principle for a drug discovery pipeline: target-binder screens followed by phenotypic validation and mechanistic studies.

Dynamics of huntingtin protein interactions in the striatum identifies candidate modifiers of Huntington disease.

Huntington disease (HD) is a monogenic neurodegenerative disorder with one causative gene, huntingtin (HTT). Yet, HD pathobiology is multifactorial, suggesting that cellular factors influence disease progression. Here, we define HTT protein-protein interactions (PPIs) perturbed by the mutant protein with expanded polyglutamine in the mouse striatum, a brain region with selective HD vulnerability. Using metabolically labeled tissues and immunoaffinity purification-mass spectrometry, we establish that polyglutamine-dependent modulation of HTT PPI abundances and relative stability starts at an early stage of pathogenesis in a Q140 HD mouse model. We identify direct and indirect PPIs that are also genetic disease modifiers using in-cell two-hybrid and behavioral assays in HD human cell and Drosophila models, respectively. Validated, disease-relevant mHTT-dependent interactions encompass mediators of synaptic neurotransmission (SNAREs and glutamate receptors) and lysosomal acidification (V-ATPase). Our study provides a resource for understanding mHTT-dependent dysfunction in cortico-striatal cellular networks, partly through impaired synaptic communication and endosomal-lysosomal system. A record of this paper's Transparent Peer Review process is included in the supplemental information.

The developmental roots of neurodegeneration.

Neurodegenerative disorders can alter neural circuitry long before symptoms appear, but the path from early changes to later pathologies is obscure. In this issue of Neuron, Capizzi et al. (2021) show how early axonal growth defects in Huntington's disease create vulnerability to later degeneration.

Downregulation of glial genes involved in synaptic function mitigates Huntington's disease pathogenesis.

Most research on neurodegenerative diseases has focused on neurons, yet glia help form and maintain the synapses whose loss is so prominent in these conditions. To investigate the contributions of glia to Huntington's disease (HD), we profiled the gene expression alterations of Drosophila expressing human mutant Huntingtin (mHTT) in either glia or neurons and compared these changes to what is observed in HD human and HD mice striata. A large portion of conserved genes are concordantly dysregulated across the three species; we tested these genes in a high-throughput behavioral assay and found that downregulation of genes involved in synapse assembly mitigated pathogenesis and behavioral deficits. To our surprise, reducing dNRXN3 function in glia was sufficient to improve the phenotype of flies expressing mHTT in neurons, suggesting that mHTT's toxic effects in glia ramify throughout the brain. This supports a model in which dampening synaptic function is protective because it attenuates the excitotoxicity that characterizes HD.

When a neuron dies, through injury or disease, the body loses all communication that passes through it. The brain compensates by rerouting the flow of information through other neurons in the network. Eventually, if the loss of neurons becomes too great, compensation becomes impossible. This process happens in Alzheimer's, Parkinson's, and Huntington's disease. In the case of Huntington's disease, the cause is mutation to a single gene known as huntingtin. The mutation is present in every cell in the body but causes particular damage to parts of the brain involved in mood, thinking and movement. Neurons and other cells respond to mutations in the huntingtin gene by turning the activities of other genes up or down, but it is not clear whether all of these changes contribute to the damage seen in Huntington's disease. In fact, it is possible that some of the changes are a result of the brain trying to protect itself. So far, most research on this subject has focused on neurons because the huntingtin gene plays a role in maintaining healthy neuronal connections. But, given that all cells carry the mutated gene, it is likely that other cells are also involved. The glia are a diverse group of cells that support the brain, providing care and sustenance to neurons. These cells have a known role in maintaining the connections between neurons and may also have play a role in either causing or correcting the damage seen in Huntington's disease. The aim of Onur et al. was to find out which genes are affected by having a mutant huntingtin gene in neurons or glia, and whether severity of Huntington's disease improved or worsened when the activity of these genes changed. First, Onur et al. identified genes affected by mutant huntingtin by comparing healthy human brains to the brains of people with Huntington's disease. Repeating the same comparison in mice and fruit flies identified genes affected in the same way across all three species, revealing that, in Huntington's disease, the brain dials down glial cell genes involved in maintaining neuronal connections. To find out how these changes in gene activity affect disease severity and progression, Onur et al. manipulated the activity of each of the genes they had identified in fruit flies that carried mutant versions of huntingtin either in neurons, in glial cells or in both cell types. They then filmed the flies to see the effects of the manipulation on movement behaviors, which are affected by Huntington's disease. This revealed that purposely lowering the activity of the glial genes involved in maintaining connections between neurons improved the symptoms of the disease, but only in flies who had mutant huntingtin in their glial cells. This indicates that the drop in activity of these genes observed in Huntington's disease is the brain trying to protect itself. This work suggests that it is important to include glial cells in studies of neurological disorders. It also highlights the fact that changes in gene expression as a result of a disease are not always bad. Many alterations are compensatory, and try to either make up for or protect cells affected by the disease. Therefore, it may be important to consider whether drugs designed to treat a condition by changing levels of gene activity might undo some of the body's natural protection. Working out which changes drive disease and which changes are protective will be essential for designing effective treatments.

Dual targeting of brain region-specific kinases potentiates neurological rescue in Spinocerebellar ataxia type 1.

A critical question in neurodegeneration is why the accumulation of disease-driving proteins causes selective neuronal loss despite their brain-wide expression. In Spinocerebellar ataxia type 1 (SCA1), accumulation of polyglutamine-expanded Ataxin-1 (ATXN1) causes selective degeneration of cerebellar and brainstem neurons. Previous studies revealed that inhibiting Msk1 reduces phosphorylation of ATXN1 at S776 as well as its levels leading to improved cerebellar function. However, there are no regulators that modulate ATXN1 in the brainstem-the brain region whose pathology is most closely linked to premature death. To identify new regulators of ATXN1, we performed genetic screens and identified a transcription factor-kinase axis (ZBTB7B-RSK3) that regulates ATXN1 levels. Unlike MSK1, RSK3 is highly expressed in the human and mouse brainstems where it regulates Atxn1 by phosphorylating S776. Reducing Rsk3 rescues brainstem-associated pathologies and deficits, and lowering Rsk3 and Msk1 together improves cerebellar and brainstem function in an SCA1 mouse model. Our results demonstrate that selective vulnerability of brain regions in SCA1 is governed by region-specific regulators of ATXN1, and targeting multiple regulators could rescue multiple degenerating brain areas.

Spinocerebellar Ataxia Type 1 protein Ataxin-1 is signaled to DNA damage by ataxia-telangiectasia mutated kinase.

Spinocerebellar Ataxia Type 1 (SCA1) is an autosomal dominant neurodegenerative disorder caused by a polyglutamine expansion in the ataxin-1 protein. Recent genetic correlational studies have implicated DNA damage repair pathways in modifying the age at onset of disease symptoms in SCA1 and Huntington's Disease, another polyglutamine expansion disease. We demonstrate that both endogenous and transfected ataxin-1 localizes to sites of DNA damage, which is impaired by polyglutamine expansion. This response is dependent on ataxia-telangiectasia mutated (ATM) kinase activity. Further, we characterize an ATM phosphorylation motif within ataxin-1 at serine 188. We show reduction of the Drosophila ATM homolog levels in a ATXN1[82Q] Drosophila model through shRNA or genetic cross ameliorates motor symptoms. These findings offer a possible explanation as to why DNA repair was implicated in SCA1 pathogenesis by past studies. The similarities between the ataxin-1 and the huntingtin responses to DNA damage provide further support for a shared pathogenic mechanism for polyglutamine expansion diseases.

Harnessing the paradoxical phenotypes of APOE ɛ2 and APOE ɛ4 to identify genetic modifiers in Alzheimer's disease.

The strongest genetic risk factor for idiopathic late-onset Alzheimer's disease (LOAD) is apolipoprotein E (APOE) ɛ4, while the APOE ɛ2 allele is protective. However, there are paradoxical APOE ɛ4 carriers who remain disease-free and APOE ɛ2 carriers with LOAD. We compared exomes of healthy APOE ɛ4 carriers and APOE ɛ2 Alzheimer's disease (AD) patients, prioritizing coding variants based on their predicted functional impact, and identified 216 genes with differential mutational load between these two populations. These candidate genes were significantly dysregulated in LOAD brains, and many modulated tau- or ß42-induced neurodegeneration in Drosophila. Variants in these genes were associated with AD risk, even in APOE ɛ3 homozygotes, showing robust predictive power for risk stratification. Network analyses revealed involvement of candidate genes in brain cell type-specific pathways including synaptic biology, dendritic spine pruning and inflammation. These potential modifiers of LOAD may constitute novel biomarkers, provide potential therapeutic intervention avenues, and support applying this approach as larger whole exome sequencing cohorts become available.