Integrated genomics and proteomics define huntingtin CAG length-dependent networks in mice

Altered activity-regulated H3K9 acetylation at TGF-beta signaling genes during egocentric memory in Huntington's disease

Molecular mechanisms underlying cognitive deficits in Huntington's disease (HD), a striatal neurodegenerative disorder, are unknown. Here, we generated ChIPseq, 4Cseq and RNAseq data on striatal tissue of HD and control mice during striatum-dependent egocentric memory process. Multi-omics analyses showed altered activity-dependent epigenetic gene reprogramming of neuronal and glial genes regulating striatal plasticity in HD mice, which correlated with memory deficit. First, our data reveal that spatial chromatin re-organization and transcriptional induction of BDNF-related markers, regulating neuronal plasticity, were reduced since memory acquisition in the striatum of HD mice. Second, our data show that epigenetic memory implicating H3K9 acetylation, which established during late phase of memory process (e.g. during consolidation/recall) and contributed to glia-mediated, TGFβ-dependent plasticity, was compromised in HD mouse striatum. Specifically, memory-dependent regulation of H3K9 acetylation was impaired at genes controlling extracellular matrix and myelination. Our study investigating the interplay between epigenetics and memory identifies H3K9 acetylation and TGFβ signaling as new targets of striatal plasticity, which might offer innovative leads to improve HD.

Postnatal Conditional Deletion of Bcl11b in Striatal Projection Neurons Mimics the Transcriptional Signature of Huntington's Disease

The dysregulation of striatal gene expression and function is linked to multiple diseases, including Huntington's disease (HD), Parkinson's disease, X-linked dystonia-parkinsonism (XDP), addiction, autism, and schizophrenia. Striatal medium spiny neurons (MSNs) make up 90% of the neurons in the striatum and are critical to motor control. The transcription factor, Bcl11b (also known as Ctip2), is required for striatal development, but the function of Bcl11b in adult MSNs in vivo has not been investigated. We conditionally deleted Bcl11b specifically in postnatal MSNs and performed a transcriptomic and behavioral analysis on these mice. Multiple enrichment analyses showed that the D9-Cre-Bcl11btm1.1Leid transcriptional profile was similar to the HD gene expression in mouse and human data sets. A Gene Ontology enrichment analysis linked D9-Cre-Bcl11btm1.1Leid to calcium, synapse organization, specifically including the dopaminergic synapse, protein dephosphorylation, and HDAC-signaling, commonly dysregulated pathways in HD. D9-Cre-Bcl11btm1.1Leid mice had decreased DARPP-32/Ppp1r1b in MSNs and behavioral deficits, demonstrating the dysregulation of a subtype of the dopamine D2 receptor expressing MSNs. Finally, in human HD isogenic MSNs, the mislocalization of BCL11B into nuclear aggregates points to a mechanism for BCL11B loss of function in HD. Our results suggest that BCL11B is important for the function and maintenance of mature MSNs and Bcl11b loss of function drives, in part, the transcriptomic and functional changes in HD.

Identification of a transcriptional signature found in multiple models of ASD and related disorders

Epigenetic regulation plays a critical role in many neurodevelopmental disorders (NDDs), including autism spectrum disorder (ASD). In particular, many such disorders are the result of mutations in genes that encode chromatin-modifying proteins. However, although these disorders share many features, it is unclear whether they also share gene expression disruptions resulting from the aberrant regulation of chromatin. We examined five chromatin modifiers that are all linked to ASD despite their different roles in regulating chromatin. Specifically, we depleted ASH1L, CHD8, CREBBP, EHMT1, and NSD1 in parallel in a highly controlled neuronal culture system. We then identified sets of shared genes, or transcriptional signatures, that are differentially expressed following loss of multiple ASD-linked chromatin modifiers. We examined the functions of genes within the transcriptional signatures and found an enrichment in many neurotransmitter transport genes and activity-dependent genes. In addition, these genes are enriched for specific chromatin features such as bivalent domains that allow for highly dynamic regulation of gene expression. The down-regulated transcriptional signature is also observed within multiple mouse models of NDDs that result in ASD, but not those only associated with intellectual disability. Finally, the down-regulated transcriptional signature can distinguish between control and idiopathic ASD patient iPSC-derived neurons as well as postmortem tissue, demonstrating that this gene set is relevant to the human disorder. This work identifies a transcriptional signature that is found within many neurodevelopmental syndromes, helping to elucidate the link between epigenetic regulation and the underlying cellular mechanisms that result in ASD.

A survey of protein interactions and posttranslational modifications that influence the polyglutamine diseases

The presence and aggregation of misfolded proteins has deleterious effects in the nervous system. Among the various diseases caused by misfolded proteins is the family of the polyglutamine (polyQ) disorders. This family comprises nine members, all stemming from the same mutation—the abnormal elongation of a polyQ repeat in nine different proteins—which causes protein misfolding and aggregation, cellular dysfunction and disease. While it is the same type of mutation that causes them, each disease is distinct: it is influenced by regions and domains that surround the polyQ repeat; by proteins with which they interact; and by posttranslational modifications they receive. Here, we overview the role of non-polyQ regions that control the pathogenicity of the expanded polyQ repeat. We begin by introducing each polyQ disease, the genes affected, and the symptoms experienced by patients. Subsequently, we provide a survey of protein-protein interactions and posttranslational modifications that regulate polyQ toxicity. We conclude by discussing shared processes and pathways that bring some of the polyQ diseases together and may serve as common therapeutic entry points for this family of incurable disorders.

Precision machine learning to understand micro-RNA regulation in neurodegenerative diseases

Micro-RNAs (miRNAs) are short (∼21 nt) non-coding RNAs that regulate gene expression through the degradation or translational repression of mRNAs. Accumulating evidence points to a role of miRNA regulation in the pathogenesis of a wide range of neurodegenerative (ND) diseases such as, for example, Alzheimer’s disease, Parkinson’s disease, amyotrophic lateral sclerosis and Huntington disease (HD). Several systems level studies aimed to explore the role of miRNA regulation in NDs, but these studies remain challenging. Part of the problem may be related to the lack of sufficiently rich or homogeneous data, such as time series or cell-type-specific data obtained in model systems or human biosamples, to account for context dependency. Part of the problem may also be related to the methodological challenges associated with the accurate system-level modeling of miRNA and mRNA data. Here, we critically review the main families of machine learning methods used to analyze expression data, highlighting the added value of using shape-analysis concepts as a solution for precisely modeling highly dimensional miRNA and mRNA data such as the ones obtained in the study of the HD process, and elaborating on the potential of these concepts and methods for modeling complex omics data.

Huntington's disease age at motor onset is modified by the tandem hexamer repeat in TCERG1

Huntington's disease is caused by an expanded CAG tract in HTT. The length of the CAG tract accounts for over half the variance in age at onset of disease, and is influenced by other genetic factors, mostly implicating the DNA maintenance machinery. We examined a single nucleotide variant, rs79727797, on chromosome 5 in the TCERG1 gene, previously reported to be associated with Huntington's disease and a quasi-tandem repeat (QTR) hexamer in exon 4 of TCERG1 with a central pure repeat. We developed a method for calling perfect and imperfect repeats from exome-sequencing data, and tested association between the QTR in TCERG1 and residual age at motor onset (after correcting for the effects of CAG length in the HTT gene) in 610 individuals with Huntington's disease via regression analysis. We found a significant association between age at onset and the sum of the repeat lengths from both alleles of the QTR (p = 2.1 × 10-9), with each added repeat hexamer reducing age at onset by one year (95% confidence interval [0.7, 1.4]). This association explained that previously observed with rs79727797. The association with age at onset in the genome-wide association study is due to a QTR hexamer in TCERG1, translated to a glutamine/alanine tract in the protein. We could not distinguish whether this was due to cis-effects of the hexamer repeat on gene expression or of the encoded glutamine/alanine tract in the protein. These results motivate further study of the mechanisms by which TCERG1 modifies onset of HD.

HTT-OMNI: A Web-based Platform for Huntingtin Interaction Exploration and Multi-omics Data Integration

Huntington's disease (HD) is a progressive neurological disorder that is caused by polyglutamine expansion of the huntingtin (HTT) protein. With the hope to uncover key modifiers of disease, a focus of the field of HD research has been on characterizing HTT-interacting proteins (HIPs) and the effect of the HTT polyglutamine expansion on the cellular omics landscape. However, while hundreds of studies have uncovered over 3000 potential HIPs to date, a means to interrogate these complementary interaction and omics datasets does not exist. The lack of a unified platform for exploring this breadth of potential HIPs and associated omics data represents a substantial barrier toward understanding the impact of HTT polyQ expansion and identifying interactions proximal to HD pathogenesis. Here, we describe the development of a web-based platform called HTT-OMNI (HTT OMics and Network Integration). This application facilitates the visualization and exploration of ∼3400 potential HTT interactors (from the HINT database) and their associated polyQ-dependent omics measurements, such as transcriptome and proteome abundances. Additionally, HTT-OMNI allows for the integration of user-generated datasets with existing HIPs and omic measurements. We first demonstrate the utility of HTT-OMNI for filtering existing HTT PPIs based on a variety of experimental metadata parameters, highlighting its capacity to select for HIPs detected in specific model organisms and tissues. Next, we leverage our application to visualize the relationships between HTT PPIs, genetic disease modifiers, and their multiomic landscape. Finally, we generate and analyze a previously unreported dataset of HTT PPIs, aimed at defining tissue-specific HTT interactions and the polyQ-dependent modulation of their relative stabilities in the cortex and striatum of HD mouse models.

Huntingtin and Its Partner Huntingtin-Associated Protein 40: Structural and Functional Considerations in Health and Disease

Since the discovery of the mutation causing Huntington’s disease (HD) in 1993, it has been debated whether an expanded polyglutamine (polyQ) stretch affects the properties of the huntingtin (HTT) protein and thus contributes to the pathological mechanisms responsible for HD. Here we review the current knowledge about the structure of HTT, alone (apo-HTT) or in a complex with Huntingtin-Associated Protein 40 (HAP40), the influence of polyQ-length variation on apo-HTT and the HTT-HAP40 complex, and the biology of HAP40. Phylogenetic analyses suggest that HAP40 performs essential functions. Highlighting the relevance of its interaction with HTT, HAP40 is one of the most abundant partners copurifying with HTT and is rapidly degraded, when HTT levels are reduced. As the levels of both proteins decrease during disease progression, HAP40 could also be a biomarker for HD. Whether declining HAP40 levels contribute to disease etiology is an open question. Structural studies have shown that the conformation of apo-HTT is less constrained but resembles that adopted in the HTT-HAP40 complex, which is exceptionally stable because of extensive interactions between HAP40 and the three domains of HTT. The complex— and to some extent apo-HTT— resists fragmentation after limited proteolysis. Unresolved regions of apo-HTT, constituting about 25% of the protein, are the main sites of post-translational modifications and likely have major regulatory functions. PolyQ elongation does not substantially alter the structure of HTT, alone or when associated with HAP40. Particularly, polyQ above the disease length threshold does not induce drastic conformational changes in full-length HTT. Therefore, models of HD pathogenesis stating that polyQ expansion drastically alters HTT properties should be reconsidered.

Modelling the Human Blood-Brain Barrier in Huntington Disease

While blood–brain barrier (BBB) dysfunction has been described in neurological disorders, including Huntington’s disease (HD), it is not known if endothelial cells themselves are functionally compromised when promoting BBB dysfunction. Furthermore, the underlying mechanisms of BBB dysfunction remain elusive given the limitations with mouse models and post mortem tissue to identify primary deficits. We established models of BBB and undertook a transcriptome and functional analysis of human induced pluripotent stem cell (iPSC)-derived brain-like microvascular endothelial cells (iBMEC) from HD patients or unaffected controls. We demonstrated that HD-iBMECs have abnormalities in barrier properties, as well as in specific BBB functions such as receptor-mediated transcytosis.

Potential disease-modifying therapies for Huntington's disease: lessons learned and future opportunities

Huntington's disease is the most frequent autosomal dominant neurodegenerative disorder; however, no disease-modifying interventions are available for patients with this disease. The molecular pathogenesis of Huntington's disease is complex, with toxicity that arises from full-length expanded huntingtin and N-terminal fragments of huntingtin, which are both prone to misfolding due to proteolysis; aberrant intron-1 splicing of the HTT gene; and somatic expansion of the CAG repeat in the HTT gene. Potential interventions for Huntington's disease include therapies targeting huntingtin DNA and RNA, clearance of huntingtin protein, DNA repair pathways, and other treatment strategies targeting inflammation and cell replacement. The early termination of trials of the antisense oligonucleotide tominersen suggest that it is time to reflect on lessons learned, where the field stands now, and the challenges and opportunities for the future.

Brain microvascular endothelial cell dysfunction in an isogenic juvenile iPSC model of Huntington's disease

Huntington's disease (HD) is an inherited neurodegenerative disease caused by expansion of cytosine-adenine-guanine (CAG) repeats in the huntingtin gene, which leads to neuronal loss and decline in cognitive and motor function. Increasing evidence suggests that blood-brain barrier (BBB) dysfunction may contribute to progression of the disease. Studies in animal models, in vitro models, and post-mortem tissue find that disease progression is associated with increased microvascular density, altered cerebral blood flow, and loss of paracellular and transcellular barrier function. Here, we report on changes in BBB phenotype due to expansion of CAG repeats using an isogenic pair of induced pluripotent stem cells (iPSCs) differentiated into brain microvascular endothelial-like cells (iBMECs). We show that CAG expansion associated with juvenile HD alters the trajectory of iBMEC differentiation, producing cells with ~ two-fold lower percentage of adherent endothelial cells. CAG expansion is associated with diminished transendothelial electrical resistance and reduced tight junction protein expression, but no significant changes in paracellular permeability. While mutant huntingtin protein (mHTT) aggregates were not observed in HD iBMECs, widespread transcriptional dysregulation was observed in iBMECs compared to iPSCs. In addition, CAG expansion in iBMECs results in distinct responses to pathological and therapeutic perturbations including angiogenic factors, oxidative stress, and osmotic stress. In a tissue-engineered BBB model, iBMECs show subtle changes in phenotype, including differences in cell turnover and immune cell adhesion. Our results further support that CAG expansion in BMECs contributes to BBB dysfunction during HD.

CK2 alpha prime and alpha-synuclein pathogenic functional interaction mediates synaptic dysregulation in huntington's disease

Huntington's disease (HD) is a neurodegenerative disorder caused by a CAG trinucleotide repeat expansion in the HTT gene for which no therapies are available. HTT mutation causes protein misfolding and aggregation, preferentially affecting medium spiny neurons (MSNs) of the basal ganglia. Transcriptional perturbations in synaptic genes and neuroinflammation are key processes that precede MSN dysfunction and motor symptom onset. Understanding the interplay between these processes is crucial to develop effective therapeutic strategies to treat HD. We investigated the role of protein kinase CK2α', a kinase upregulated in MSNs in HD and previously associated with Parkinson's disease (PD), in the regulation of neuroinflammation and synaptic function in HD. We used the heterozygous knock-in zQ175 HD mouse model and compared that to zQ175 mice lacking one allele of CK2α' (zQ175:CK2α'(±)). CK2α' haploinsufficiency in zQ175 mice resulted in decreased levels of pro-inflammatory cytokines, HTT aggregation, astrogliosis and transcriptional alterations of synaptic genes related to glutamatergic signaling. zQ175:CK2α'(±) mice also presented increased frequency of striatal miniature excitatory postsynaptic currents (mEPSCs), an indicator of synaptic activity, and improved motor coordination compared to zQ175 mice. Neuropathological and phenotypic changes mediated by CK2α' were connected to alpha-synuclein (α-syn) dysregulation and correlated with differences in α-syn serine 129 phosphorylation (pS129-α-syn), a post-translational modification involved in α-synucleinopathy and shown to be regulated by CK2 in PD. pS129-α-syn was increased in the nuclei of MSNs in zQ175 mice and in the striatum of patients with HD, and it decreased in zQ175:CK2α'(±) mice. Collectively, our data established a novel connection between CK2α', neuroinflammation and synaptic gene dysregulation with synucleinopathy in HD and suggested common molecular mechanisms of neurodegeneration between HD and PD. Our results also support CK2α' inhibition as a potential therapeutic strategy to modulate neuronal function and neuroprotection in HD.

Analysis of Huntington's Disease Modifiers Using the Hyperbolic Mapping of the Protein Interaction Network

Huntington's disease (HD) is caused by the production of a mutant huntingtin (HTT) with an abnormally long poly-glutamine (polyQ) tract, forming aggregates and inclusions in neurons. Previous work by us and others has shown that an increase or decrease in polyQ-triggered aggregates can be passive simply due to the interaction of proteins with the aggregates. To search for proteins with active (functional) effects, which might be more effective in finding therapies and mechanisms of HD, we selected among the proteins that interact with HTT a total of 49 pairs of proteins that, while being paralogous to each other (and thus expected to have similar passive interaction with HTT), are located in different regions of the protein interaction network (suggesting participation in different pathways or complexes). Three of these 49 pairs contained members with opposite effects on HD, according to the literature. The negative members of the three pairs, MID1, IKBKG, and IKBKB, interact with PPP2CA and TUBB, which are known negative factors in HD, as well as with HSP90AA1 and RPS3. The positive members of the three pairs interact with HSPA9. Our results provide potential HD modifiers of functional relevance and reveal the dynamic aspect of paralog evolution within the interaction network.

A Glimpse of Molecular Biomarkers in Huntington's Disease

Huntington's disease (HD) is a devastating neurodegenerative disorder that is caused by an abnormal expansion of CAG repeats in the Huntingtin (HTT) gene. Although the main symptomatology is explained by alterations at the level of the central nervous system, predominantly affecting the basal ganglia, a peripheral component of the disease is being increasingly acknowledged. Therefore, the manifestation of the disease is complex and variable among CAG expansion carriers, introducing uncertainty in the appearance of specific signs, age of onset and severity of disease. The monogenic nature of the disorder allows a precise diagnosis, but the use of biomarkers with prognostic value is still needed to achieve clinical management of the patients in an individual manner. In addition, we need tools to evaluate the patient's response to potential therapeutic approaches. In this review, we provide a succinct summary of the most interesting molecular biomarkers that have been assessed in patients, mostly obtained from body fluids such as cerebrospinal fluid, peripheral blood and saliva.

Modifier pathways in polyglutamine (PolyQ) diseases: from genetic screens to drug targets

Polyglutamine (PolyQ) diseases include a group of inherited neurodegenerative disorders caused by unstable expansions of CAG trinucleotide repeats in the coding region of specific genes. Such genetic alterations produce abnormal proteins containing an unusually long PolyQ tract that renders them more prone to aggregate and cause toxicity. Although research in the field in the last years has contributed significantly to the knowledge of the biological mechanisms implicated in these diseases, effective treatments are still lacking. In this review, we revisit work performed in models of PolyQ diseases, namely the yeast Saccharomyces cerevisiae, the nematode worm Caenorhabditis elegans and the fruit fly Drosophila melanogaster, and provide a critical overview of the high-throughput unbiased genetic screens that have been performed using these systems to identify novel genetic modifiers of PolyQ diseases. These approaches have revealed a wide variety of cellular processes that modulate the toxicity and aggregation of mutant PolyQ proteins, reflecting the complexity of these disorders and demonstrating how challenging the development of therapeutic strategies can be. In addition to the unbiased large-scale genetic screenings in non-vertebrate models, complementary studies in mammalian systems, closer to humans, have contributed with novel genetic modifiers of PolyQ diseases, revealing neuronal function and inflammation as key disease modulators. A pathway enrichment analysis, using the human orthologues of genetic modifiers of PolyQ diseases clustered modifier genes into major themes translatable to the human disease context, such as protein folding and transport as well as transcription regulation. Innovative genetic strategies of genetic manipulation, together with significant advances in genomics and bioinformatics, are taking modifier genetic studies to more realistic disease contexts. The characterization of PolyQ disease modifier pathways is of extreme relevance to reveal novel therapeutic possibilities to delay disease onset and progression in patients.

Group dynamics goes awry: PolyQ-expanded huntingtin gains unwanted partners

In this issue of Cell Systems, Greco et al. define high-confidence polyglutamine-dependent huntingtin interactors using AP-MS and complementary approaches and categorize them based on their interaction abundance and stability. The study reveals that a toxic gain of polyQ-dependent Htt interacting partners is a robust feature of HD pathogenesis.

Association of DNA Methylation Patterns in 7 Novel Genes With Ischemic Stroke in the Northern Chinese Population

Background: Ischemic stroke is a highly complex disorder. This study aims to identify novel methylation changes in ischemic stroke.

Methods: We carried out an epigenome-wide study of ischemic stroke using an Infinium HumanMethylation 850K array (cases:controls = 4:4). 10 CpG sites in 8 candidate genes from gene ontology analytics top-ranked pathway were selected to validate 850K BeadChip results (cases:controls = 20:20). We further qualified the methylation level of promoter regions in 8 candidate genes (cases:controls = 188:188). Besides, we performed subgroup analysis, dose-response relationship and diagnostic prediction polygenic model of candidate genes.

Results: In the discovery stage, we found 462 functional DNA methylation positions to be associated with ischemic stroke. Gene ontology analysis highlighted the “calcium-dependent cell-cell adhesion via plasma membrane cell adhesion molecules” item, including 8 candidate genes (CDH2/PCDHB10/PCDHB11/PCDHB14/PCDHB16/PCDHB3/PCDHB6/PCDHB9). In the replication stage, we identified 5 differentially methylated loci in 20 paired samples and 7 differentially methylated genes (CDH2/PCDHB10/PCDHB11/PCDHB14/PCDHB16/PCDHB3/PCDHB9) in 188 paired samples. Subgroup analysis showed that the methylation level of above 7 genes remained significantly different in the male subgroup, large-artery atherosclerosis subgroup and right hemisphere subgroup. The methylation level of each gene was grouped into quartiles, and Q4 groups of the 7 genes were associated with higher risk of ischemic stroke than Q1 groups (p < 0.05). Besides, the polygenic model showed high diagnostic specificity (0.8723), sensitivity (0.883), and accuracy (0.8777).

Conclusion: Our results demonstrate that DNA methylation plays a crucial part in ischemic stroke. The methylation of these 7 genes may be potential diagnostic biomarker for ischemic stroke.

Back-to-BACs in Huntington's disease modeling

In this issue of Neuron, Gu et al. (2022) describe a new BAC mouse model that faithfully recapitulates various aspects of Huntington's disease pathobiology and reveals important insights into the relative toxicities of mHTT-derived products.

Uninterrupted CAG repeat drives striatum-selective transcriptionopathy and nuclear pathogenesis in human Huntingtin BAC mice

In Huntington's disease (HD), the uninterrupted CAG repeat length, but not the polyglutamine length, predicts disease onset. However, the underlying pathobiology remains unclear. Here, we developed bacterial artificial chromosome (BAC) transgenic mice expressing human mutant huntingtin (mHTT) with uninterrupted, and somatically unstable, CAG repeats that exhibit progressive disease-related phenotypes. Unlike prior mHTT transgenic models with stable, CAA-interrupted, polyglutamine-encoding repeats, BAC-CAG mice show robust striatum-selective nuclear inclusions and transcriptional dysregulation resembling those in murine huntingtin knockin models and HD patients. Importantly, the striatal transcriptionopathy in HD models is significantly correlated with their uninterrupted CAG repeat length but not polyglutamine length. Finally, among the pathogenic entities originating from mHTT genomic transgenes and only present or enriched in the uninterrupted CAG repeat model, somatic CAG repeat instability and nuclear mHTT aggregation are best correlated with early-onset striatum-selective molecular pathogenesis and locomotor and sleep deficits, while repeat RNA-associated pathologies and repeat-associated non-AUG (RAN) translation may play less selective or late pathogenic roles, respectively.

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.

Development of a ligand for in vivo imaging of mutant huntingtin in Huntington's disease

Huntington's disease (HD) is a dominantly inherited neurodegenerative disorder caused by a CAG trinucleotide expansion in the huntingtin (HTT) gene that encodes the pathologic mutant HTT (mHTT) protein with an expanded polyglutamine (polyQ) tract. Whereas several therapeutic programs targeting mHTT expression have advanced to clinical evaluation, methods to visualize mHTT protein species in the living brain are lacking. Here, we demonstrate the development and characterization of a positron emission tomography (PET) imaging radioligand with high affinity and selectivity for mHTT aggregates. This small molecule radiolabeled with ¹¹C ([¹¹C]CHDI-180R) allowed noninvasive monitoring of mHTT pathology in the brain and could track region- and time-dependent suppression of mHTT in response to therapeutic interventions targeting mHTT expression in a rodent model. We further showed that in these animals, therapeutic agents that lowered mHTT in the striatum had a functional restorative effect that could be measured by preservation of striatal imaging markers, enabling a translational path to assess the functional effect of mHTT lowering.

Deletion of SUMO1 attenuates behavioral and anatomical deficits by regulating autophagic activities in Huntington disease

The CAG expansion of huntingtin (mHTT) associated with Huntington disease (HD) is a ubiquitously expressed gene, yet it prominently damages the striatum and cortex, followed by widespread peripheral defects as the disease progresses. However, the underlying mechanisms of neuronal vulnerability are unclear. Previous studies have shown that SUMO1 (small ubiquitin-like modifier-1) modification of mHtt promotes cellular toxicity, but the in vivo role and functions of SUMO1 in HD pathogenesis are unclear. Here, we report that SUMO1 deletion in Q175DN HD-het knockin mice (HD mice) prevented age-dependent HD-like motor and neurological impairments and suppressed the striatal atrophy and inflammatory response. SUMO1 deletion caused a drastic reduction in soluble mHtt levels and nuclear and extracellular mHtt inclusions while increasing cytoplasmic mHtt inclusions in the striatum of HD mice. SUMO1 deletion promoted autophagic activity, characterized by augmented interactions between mHtt inclusions and a lysosomal marker (LAMP1), increased LC3B- and LAMP1 interaction, and decreased interaction of sequestosome-1 (p62) and LAMP1 in DARPP-32-positive medium spiny neurons in HD mice. Depletion of SUMO1 in an HD cell model also diminished the mHtt levels and enhanced autophagy flux. In addition, the SUMOylation inhibitor ginkgolic acid strongly enhanced autophagy and diminished mHTT levels in human HD fibroblasts. These results indicate that SUMO is a critical therapeutic target in HD and that blocking SUMO may ameliorate HD pathogenesis by regulating autophagy activities.

An Augmented High-Dimensional Graphical Lasso Method to Incorporate Prior Biological Knowledge for Global Network Learning

Biological networks are often inferred through Gaussian graphical models (GGMs) using gene or protein expression data only. GGMs identify conditional dependence by estimating a precision matrix between genes or proteins. However, conventional GGM approaches often ignore prior knowledge about protein-protein interactions (PPI). Recently, several groups have extended GGM to weighted graphical Lasso (wGlasso) and network-based gene set analysis (Netgsa) and have demonstrated the advantages of incorporating PPI information. However, these methods are either computationally intractable for large-scale data, or disregard weights in the PPI networks. To address these shortcomings, we extended the Netgsa approach and developed an augmented high-dimensional graphical Lasso (AhGlasso) method to incorporate edge weights in known PPI with omics data for global network learning. This new method outperforms weighted graphical Lasso-based algorithms with respect to computational time in simulated large-scale data settings while achieving better or comparable prediction accuracy of node connections. The total runtime of AhGlasso is approximately five times faster than weighted Glasso methods when the graph size ranges from 1,000 to 3,000 with a fixed sample size (n = 300). The runtime difference between AhGlasso and weighted Glasso increases when the graph size increases. Using proteomic data from a study on chronic obstructive pulmonary disease, we demonstrate that AhGlasso improves protein network inference compared to the Netgsa approach by incorporating PPI information.

Analysis of LINE1 Retrotransposons in Huntington’s Disease

Transposable elements (TEs) are mobile genetic elements that made up about half the human genome. Among them, the autonomous non-LTR retrotransposon long interspersed nuclear element-1 (L1) is the only currently active TE in mammals and covers about 17% of the mammalian genome. L1s exert their function as structural elements in the genome, as transcribed RNAs to influence chromatin structure and as retrotransposed elements to shape genomic variation in somatic cells. L1s activity has been shown altered in several diseases of the nervous system. Huntington disease (HD) is a dominantly inherited neurodegenerative disorder caused by an expansion of a CAG repeat in the HTT gene which leads to a gradual loss of neurons most prominently in the striatum and, to a lesser extent, in cortical brain regions. The length of the expanded CAG tract is related to age at disease onset, with longer repeats leading to earlier onset. Here we carried out bioinformatic analysis of public RNA-seq data of a panel of HD mouse models showing that a decrease of L1 RNA expression recapitulates two hallmarks of the disease: it correlates to CAG repeat length and it occurs in the striatum, the site of neurodegeneration. Results were then experimentally validated in HttQ111 knock-in mice. The expression of L1-encoded proteins was independent from L1 RNA levels and differentially regulated in time and tissues. The pattern of expression L1 RNAs in human HD post-mortem brains showed similarity to mouse models of the disease. This work suggests the need for further study of L1s in HD and adds support to the current hypothesis that dysregulation of TEs may be involved in neurodegenerative diseases.

Identification of cyclin D1 as a major modulator of 3-nitropropionic acid-induced striatal neurodegeneration

Mitochondria dysfunction occurs in the aging brain as well as in several neurodegenerative disorders and predisposes neuronal cells to enhanced sensitivity to neurotoxins. 3-nitropropionic acid (3-NP) is a naturally occurring plant and fungal neurotoxin that causes neurodegeneration predominantly in the striatum by irreversibly inhibiting the tricarboxylic acid respiratory chain enzyme, succinate dehydrogenase (SDH), the main constituent of the mitochondria respiratory chain complex II. Significantly, although 3-NP-induced inhibition of SDH occurs in all brain regions, neurodegeneration occurs primarily and almost exclusively in the striatum for reasons still not understood. In rodents, 3-NP-induced striatal neurodegeneration depends on the strain background suggesting that genetic differences among genotypes modulate toxicant variability and mechanisms that underlie 3-NP-induced neuronal cell death. Using the large BXD family of recombinant inbred (RI) strains we demonstrate that variants in Ccnd1 - the gene encoding cyclin D1 - of the DBA/2 J parent underlie the resistance to 3-NP-induced striatal neurodegeneration. In contrast, the Ccnd1 variant inherited from the widely used C57BL/6 J parental strain confers sensitivity. Given that cellular stress triggers induction of cyclin D1 expression followed by cell-cycle re-entry and consequent neuronal cell death, we sought to determine if the C57BL/6 J and DBA/2 J Ccnd1 variants are differentially modulated in response to 3-NP. We confirm that 3-NP induces cyclin D1 expression in striatal neuronal cells of C57BL/6 J, but this response is blunted in the DBA/2 J. We further show that striatal-specific alternative processing of a highly conserved 3'UTR negative regulatory region of Ccnd1 co-segregates with the C57BL/6 J parental Ccnd1 allele in BXD strains and that its differential processing accounts for sensitivity or resistance to 3-NP. Our results indicate that naturally occurring Ccnd1 variants may play a role in the variability observed in neurodegenerative disorders involving mitochondria complex II dysfunction and point to cyclin D1 as a possible therapeutic target.

Nucleolar stress controls mutant Huntington toxicity and monitors Huntington's disease progression

Transcriptional and cellular-stress surveillance deficits are hallmarks of Huntington's disease (HD), a fatal autosomal-dominant neurodegenerative disorder caused by a pathological expansion of CAG repeats in the Huntingtin (HTT) gene. The nucleolus, a dynamic nuclear biomolecular condensate and the site of ribosomal RNA (rRNA) transcription, is implicated in the cellular stress response and in protein quality control. While the exact pathomechanisms of HD are still unclear, the impact of nucleolar dysfunction on HD pathophysiology in vivo remains elusive. Here we identified aberrant maturation of rRNA and decreased translational rate in association with human mutant Huntingtin (mHTT) expression. The protein nucleophosmin 1 (NPM1), important for nucleolar integrity and rRNA maturation, loses its prominent nucleolar localization. Genetic disruption of nucleolar integrity in vulnerable striatal neurons of the R6/2 HD mouse model decreases the distribution of mHTT in a disperse state in the nucleus, exacerbating motor deficits. We confirmed NPM1 delocalization in the gradually progressing zQ175 knock-in HD mouse model: in the striatum at a presymptomatic stage and in the skeletal muscle at an early symptomatic stage. In Huntington's patient skeletal muscle biopsies, we found a selective redistribution of NPM1, similar to that in the zQ175 model. Taken together, our study demonstrates that nucleolar integrity regulates the formation of mHTT inclusions in vivo, and identifies NPM1 as a novel, readily detectable peripheral histopathological marker of HD progression.

Modulating FKBP5/FKBP51 and autophagy lowers HTT (huntingtin) levels

Current disease-modifying therapies for Huntington disease (HD) focus on lowering mutant HTT (huntingtin; mHTT) levels, and the immunosuppressant drug rapamycin is an intriguing therapeutic for aging and neurological disorders. Rapamycin interacts with FKBP1A/FKBP12 and FKBP5/FKBP51, inhibiting the MTORC1 complex and increasing cellular clearance mechanisms. Whether the levels of FKBP (FK506 binding protein) family members are altered in HD models and if these proteins are potential therapeutic targets for HD have not been investigated. Here, we found levels of FKBP5 are significantly reduced in HD R6/2 and zQ175 mouse models and human HD isogenic neural stem cells and medium spiny neurons derived from induced pluripotent stem cells. Moreover, FKBP5 interacts and colocalizes with HTT in the striatum and cortex of zQ175 mice and controls. Importantly, when we decreased FKBP5 levels or activity by genetic or pharmacological approaches, we observed reduced levels of mHTT in our isogenic human HD stem cell model. Decreasing FKBP5 levels by siRNA or pharmacological inhibition increased LC3-II levels and macroautophagic/autophagic flux, suggesting autophagic cellular clearance mechanisms are responsible for mHTT lowering. Unlike rapamycin, the effect of pharmacological inhibition with SAFit2, an inhibitor of FKBP5, is MTOR independent. Further, in vivo treatment for 2 weeks with SAFit2, results in reduced HTT levels in both HD R6/2 and zQ175 mouse models. Our studies establish FKBP5 as a protein involved in the pathogenesis of HD and identify FKBP5 as a potential therapeutic target for HD.Abbreviations : ACTB/β-actin: actin beta; AD: Alzheimer disease; BafA1: bafilomycin A1; BCA: bicinchoninic acid; BBB: blood brain barrier; BSA: bovine serum albumin; CoIP: co-immunoprecipitation; DMSO: dimethyl sulfoxide; DTT: dithiothreitol; FKBPs: FK506 binding proteins; HD: Huntington disease; HTT: huntingtin; iPSC: induced pluripotent stem cells; MAP1LC3/LC3:microtubule associated protein 1 light chain 3; MAPT/tau: microtubule associated protein tau; MES: 2-ethanesulfonic acid; MOPS: 3-(N-morphorlino)propanesulfonic acid); MSN: medium spiny neurons; mHTT: mutant huntingtin; MTOR: mechanistic target of rapamycin kinase; NSC: neural stem cells; ON: overnight; PD: Parkinson disease; PPIase: peptidyl-prolyl cis/trans-isomerases; polyQ: polyglutamine; PPP1R1B/DARPP-32: protein phosphatase 1 regulatory inhibitor subunit 1B; PTSD: post-traumatic stress disorder; RT: room temperature; SQSTM1/p62: sequestosome 1; SDS-PAGE: sodium dodecyl sulfate-polyacrylamide gel electrophoresis; TBST:Tris-buffered saline, 0.1% Tween 20; TUBA: tubulin; ULK1: unc-51 like autophagy activating kinase 1; VCL: vinculin; WT: littermate controls.

Karyopherin abnormalities in neurodegenerative proteinopathies

Neurodegenerative proteinopathies are characterized by progressive cell loss that is preceded by the mislocalization and aberrant accumulation of proteins prone to aggregation. Despite their different physiological functions, disease-related proteins like tau, α-synuclein, TAR DNA binding protein-43, fused in sarcoma and mutant huntingtin, all share low complexity regions that can mediate their liquid-liquid phase transitions. The proteins' phase transitions can range from native monomers to soluble oligomers, liquid droplets and further to irreversible, often-mislocalized aggregates that characterize the stages and severity of neurodegenerative diseases. Recent advances into the underlying pathogenic mechanisms have associated mislocalization and aberrant accumulation of disease-related proteins with defective nucleocytoplasmic transport and its mediators called karyopherins. These studies identify karyopherin abnormalities in amyotrophic lateral sclerosis, frontotemporal dementia, Alzheimer's disease, and synucleinopathies including Parkinson's disease and dementia with Lewy bodies, that range from altered expression levels to the subcellular mislocalization and aggregation of karyopherin α and β proteins. The reported findings reveal that in addition to their classical function in nuclear import and export, karyopherins can also act as chaperones by shielding aggregation-prone proteins against misfolding, accumulation and irreversible phase-transition into insoluble aggregates. Karyopherin abnormalities can, therefore, be both the cause and consequence of protein mislocalization and aggregate formation in degenerative proteinopathies. The resulting vicious feedback cycle of karyopherin pathology and proteinopathy identifies karyopherin abnormalities as a common denominator of onset and progression of neurodegenerative disease. Pharmacological targeting of karyopherins, already in clinical trials as therapeutic intervention targeting cancers such as glioblastoma and viral infections like COVID-19, may therefore represent a promising new avenue for disease-modifying treatments in neurodegenerative proteinopathies.

Modifiers of CAG/CTG Repeat Instability: Insights from Mammalian Models

At fifteen different genomic locations, the expansion of a CAG/CTG repeat causes a neurodegenerative or neuromuscular disease, the most common being Huntington's disease and myotonic dystrophy type 1. These disorders are characterized by germline and somatic instability of the causative CAG/CTG repeat mutations. Repeat lengthening, or expansion, in the germline leads to an earlier age of onset or more severe symptoms in the next generation. In somatic cells, repeat expansion is thought to precipitate the rate of disease. The mechanisms underlying repeat instability are not well understood. Here we review the mammalian model systems that have been used to study CAG/CTG repeat instability, and the modifiers identified in these systems. Mouse models have demonstrated prominent roles for proteins in the mismatch repair pathway as critical drivers of CAG/CTG instability, which is also suggested by recent genome-wide association studies in humans. We draw attention to a network of connections between modifiers identified across several systems that might indicate pathway crosstalk in the context of repeat instability, and which could provide hypotheses for further validation or discovery. Overall, the data indicate that repeat dynamics might be modulated by altering the levels of DNA metabolic proteins, their regulation, their interaction with chromatin, or by direct perturbation of the repeat tract. Applying novel methodologies and technologies to this exciting area of research will be needed to gain deeper mechanistic insight that can be harnessed for therapies aimed at preventing repeat expansion or promoting repeat contraction.

Huntington's disease mouse models: unraveling the pathology caused by CAG repeat expansion

Huntington's disease (HD) is a neurodegenerative disease that results in motor and cognitive dysfunction, leading to early death. HD is caused by an expansion of CAG repeats in the huntingtin gene (HTT). Here, we review the mouse models of HD. They have been used extensively to better understand the molecular and cellular basis of disease pathogenesis as well as to provide non-human subjects to test the efficacy of potential therapeutics. The first and best-studied in vivo rodent model of HD is the R6/2 mouse, in which a transgene containing the promoter and exon 1 fragment of human HTT with 150 CAG repeats was inserted into the mouse genome. R6/2 mice express rapid, robust behavioral pathologies and display a number of degenerative abnormalities in neuronal populations most vulnerable in HD. The first conditional full-length mutant huntingtin (mHTT) mouse model of HD was the bacterial artificial chromosome (BAC) transgenic mouse model of HD (BACHD), which expresses human full-length mHTT with a mixture of 97 CAG-CAA repeats under the control of endogenous HTT regulatory machinery. It has been useful in identifying the role of mHTT in specific neuronal populations in degenerative processes. In the knock-in (KI) model of HD, the expanded human CAG repeats and human exon 1 are inserted into the mouse Htt locus, so a chimera of the full-length mouse protein with the N-terminal human portion is expressed. Many of aspects of the pathology and behavioral deficits in the KI model better mimic disease characteristics found in HD patients than other models. Accordingly, some have proposed that these mice may be preferable models of the disease over others. Indeed, as our understanding of HD advances, so will the design of animal models to test and develop HD therapies.

Netboost: Boosting-Supported Network Analysis Improves High-Dimensional Omics Prediction in Acute Myeloid Leukemia and Huntington's Disease

State-of-the art selection methods fail to identify weak but cumulative effects of features found in many high-dimensional omics datasets. Nevertheless, these features play an important role in certain diseases. We present Netboost, a three-step dimension reduction technique. First, a boosting-based filter is combined with the topological overlap measure to identify the essential edges of the network. Second, sparse hierarchical clustering is applied on the selected edges to identify modules and finally module information is aggregated by the first principal components. We demonstrate the application of the newly developed Netboost in combination with CoxBoost for survival prediction of DNA methylation and gene expression data from 180 acute myeloid leukemia (AML) patients and show, based on cross-validated prediction error curve estimates, its prediction superiority over variable selection on the full dataset as well as over an alternative clustering approach. The identified signature related to chromatin modifying enzymes was replicated in an independent dataset, the phase II AMLSG 12-09 study. In a second application we combine Netboost with Random Forest classification and improve the disease classification error in RNA-sequencing data of Huntington's disease mice. Netboost is a freely available Bioconductor R package for dimension reduction and hypothesis generation in high-dimensional omics applications.

A Systematic Review of Transcriptional Dysregulation in Huntington's Disease Studied by RNA Sequencing

Huntington's disease (HD) is a chronic neurodegenerative disorder caused by an expansion of polyglutamine repeats in exon 1 of the Huntingtin gene. Transcriptional dysregulation accompanied by epigenetic alterations is an early and central disease mechanism in HD yet, the exact mechanisms and regulators, and their associated gene expression programs remain incompletely understood. This systematic review investigates genome-wide transcriptional studies that were conducted using RNA sequencing (RNA-seq) technology in HD patients and models. The review protocol was registered at the Open Science Framework (OSF). The biomedical literature and gene expression databases, PubMed and NCBI BioProject, Array Express, European Nucleotide Archive (ENA), European Genome-Phenome Archive (EGA), respectively, were searched using the defined terms specified in the protocol following the PRISMA guidelines. We conducted a complete literature and database search to retrieve all RNA-seq-based gene expression studies in HD published until August 2020, retrieving 288 articles and 237 datasets from PubMed and the databases, respectively. A total of 27 studies meeting the eligibility criteria were included in this review. Collectively, comparative analysis of the datasets revealed frequent genes that are consistently dysregulated in HD. In postmortem brains from HD patients, DNAJB1, HSPA1B and HSPB1 genes were commonly upregulated across all brain regions and cell types except for medium spiny neurons (MSNs) at symptomatic disease stage, and HSPH1 and SAT1 genes were altered in expression in all symptomatic brain datasets, indicating early and sustained changes in the expression of genes related to heat shock response as well as response to misfolded proteins. Specifically in indirect pathway medium spiny neurons (iMSNs), mitochondria related genes were among the top uniquely dysregulated genes. Interestingly, blood from HD patients showed commonly differentially expressed genes with a number of brain regions and cells, with the highest number of overlapping genes with MSNs and BA9 region at symptomatic stage. We also found the differential expression and predicted altered activity of a set of transcription factors and epigenetic regulators, including BCL6, EGR1, FOSL2 and CREBBP, HDAC1, KDM4C, respectively, which may underlie the observed transcriptional changes in HD. Altogether, our work provides a complete overview of the transcriptional studies in HD, and by data synthesis, reveals a number of common and unique gene expression and regulatory changes across different cell and tissue types in HD. These changes could elucidate new insights into molecular mechanisms of differential vulnerability in HD. Systematic Review Registration: https://osf.io/pm3wq.

PaIRKAT: A pathway integrated regression-based kernel association test with applications to metabolomics and COPD phenotypes

High-throughput data such as metabolomics, genomics, transcriptomics, and proteomics have become familiar data types within the "-omics" family. For this work, we focus on subsets that interact with one another and represent these "pathways" as graphs. Observed pathways often have disjoint components, i.e., nodes or sets of nodes (metabolites, etc.) not connected to any other within the pathway, which notably lessens testing power. In this paper we propose the Pathway Integrated Regression-based Kernel Association Test (PaIRKAT), a new kernel machine regression method for incorporating known pathway information into the semi-parametric kernel regression framework. This work extends previous kernel machine approaches. This paper also contributes an application of a graph kernel regularization method for overcoming disconnected pathways. By incorporating a regularized or "smoothed" graph into a score test, PaIRKAT can provide more powerful tests for associations between biological pathways and phenotypes of interest and will be helpful in identifying novel pathways for targeted clinical research. We evaluate this method through several simulation studies and an application to real metabolomics data from the COPDGene study. Our simulation studies illustrate the robustness of this method to incorrect and incomplete pathway knowledge, and the real data analysis shows meaningful improvements of testing power in pathways. PaIRKAT was developed for application to metabolomic pathway data, but the techniques are easily generalizable to other data sources with a graph-like structure.

Isoform-Specific Reduction of the Basic Helix-Loop-Helix Transcription Factor TCF4 Levels in Huntington's Disease

Huntington's disease (HD) is an inherited neurodegenerative disorder with onset of characteristic motor symptoms at midlife, preceded by subtle cognitive and behavioral disturbances. Transcriptional dysregulation emerges early in the disease course and is considered central to HD pathogenesis. Using wild-type (wt) and HD knock-in mouse striatal cell lines we observed a HD genotype-dependent reduction in the protein levels of transcription factor 4 (TCF4), a member of the basic helix-loop-helix (bHLH) family with critical roles in brain development and function. We characterized mouse Tcf4 gene structure and expression of alternative mRNAs and protein isoforms in cell-based models of HD, and in four different brain regions of male transgenic HD mice (R6/1) from young to mature adulthood. The largest decrease in the levels of TCF4 at mRNA and specific protein isoforms were detected in the R6/1 mouse hippocampus. Translating this finding to human disease, we found reduced expression of long TCF4 isoforms in the postmortem hippocampal CA1 area and in the cerebral cortex of HD patients. Additionally, TCF4 protein isoforms showed differential synergism with the proneural transcription factor ASCL1 in activating reporter gene transcription in hippocampal and cortical cultured neurons. Induction of neuronal activity increased these synergistic effects in hippocampal but not in cortical neurons, suggesting brain region-dependent differences in TCF4 functions. Collectively, this study demonstrates isoform-specific changes in TCF4 expression in HD that could contribute to the progressive impairment of transcriptional regulation and neuronal function in this disease.

Label propagation-based semi-supervised feature selection on decoding clinical phenotypes with RNA-seq data

BACKGROUND

Clinically, behavior, cognitive, and mental functions are affected during the neurodegenerative disease progression. To date, the molecular pathogenesis of these complex disease is still unclear. With the rapid development of sequencing technologies, it is possible to delicately decode the molecular mechanisms corresponding to different clinical phenotypes at the genome-wide transcriptomic level using computational methods. Our previous studies have shown that it is difficult to distinguish disease genes from non-disease genes. Therefore, to precisely explore the molecular pathogenesis under complex clinical phenotypes, it is better to identify biomarkers corresponding to different disease stages or clinical phenotypes. So, in this study, we designed a label propagation-based semi-supervised feature selection approach (LPFS) to prioritize disease-associated genes corresponding to different disease stages or clinical phenotypes.

METHODS

In this study, we pioneering put label propagation clustering and feature selection into one framework and proposed label propagation-based semi-supervised feature selection approach. LPFS prioritizes disease genes related to different disease stages or phenotypes through the alternative iteration of label propagation clustering based on sample network and feature selection with gene expression profiles. Then the GO and KEGG pathway enrichment analysis were carried as well as the gene functional analysis to explore molecular mechanisms of specific disease phenotypes, thus to decode the changes in individual behavioral and mental characteristics during neurodegenerative disease progression.

RESULTS

Large amounts of experiments were conducted to verify the performance of LPFS with Huntington's gene expression data. Experimental results shown that LPFS performs better in comparison with the-state-of-art methods. GO and KEGG enrichment analysis of key gene sets shown that TGF-beta signaling pathway, cytokine-cytokine receptor interaction, immune response, and inflammatory response were gradually affected during the Huntington's disease progression. In addition, we found that the expression of SLC4A11, ZFP474, AMBP, TOP2A, PBK, CCDC33, APSL, DLGAP5, and Al662270 changed seriously by the development of the disease.

CONCLUSIONS

In this study, we designed a label propagation-based semi-supervised feature selection model to precisely selected key genes of different disease phenotypes. We conducted experiments using the model with Huntington's disease mice gene expression data to decode the mechanisms of it. We found many cell types, including astrocyte, microglia, and GABAergic neuron, could be involved in the pathological process.

Gene expression profiles complement the analysis of genomic modifiers of the clinical onset of Huntington disease

Huntington disease (HD) is a neurodegenerative disorder that is caused by a CAG repeat expansion in HTT. The length of this repeat, however, only explains a proportion of the variability in age of onset in patients. Genome-wide association studies have identified modifiers that contribute toward a proportion of the observed variance. By incorporating tissue-specific transcriptomic information with these results, additional modifiers can be identified. We performed a transcriptome-wide association study assessing heritable differences in genetically determined expression in diverse tissues, with genome-wide data from over 4000 patients. Functional validation of prioritized genes was undertaken in isogenic HD stem cells and patient brains. Enrichment analyses were performed with biologically relevant gene sets to identify the core pathways. HD-associated gene coexpression modules were assessed for associations with neurological phenotypes in an independent cohort and to guide drug repurposing analyses. Transcriptomic analyses identified genes that were associated with age of HD onset and displayed colocalization with gene expression signals in brain tissue (FAN1, GPR161, PMS2, SUMF2), with supporting evidence from functional experiments. This included genes involved in DNA repair, as well as novel-candidate modifier genes that have been associated with other neurological conditions. Further, cortical coexpression modules were also associated with cognitive decline and HD-related traits in a longitudinal cohort. In summary, the combination of population-scale gene expression information with HD patient genomic data identified novel modifier genes for the disorder. Further, these analyses expanded the pathways potentially involved in modifying HD onset and prioritized candidate therapeutics for future study.

Integrative enrichment analysis of gene expression based on an artificial neuron

BACKGROUND

Huntington's disease is a kind of chronic progressive neurodegenerative disease with complex pathogenic mechanisms. To data, the pathogenesis of Huntington's disease is still not fully understood, and there has been no effective treatment. The rapid development of high-throughput sequencing technologies makes it possible to explore the molecular mechanisms at the transcriptome level. Our previous studies on Huntington's disease have shown that it is difficult to distinguish disease-associated genes from non-disease genes. Meanwhile, recent progress in bio-medicine shows that the molecular origin of chronic complex diseases may not exist in the diseased tissue, and differentially expressed genes between different tissues may be helpful to reveal the molecular origin of chronic diseases. Therefore, developing integrative analysis computational methods for the multi-tissues gene expression data, exploring the relationship between differentially expressed genes in different tissues and the disease, can greatly accelerate the molecular discovery process.

METHODS

For analysis of the intra- and inter- tissues' differentially expressed genes, we designed an integrative enrichment analysis method based on an artificial neuron (IEAAN). Firstly, we calculated the differential expression scores of genes which are seen as features of the corresponding gene, using fold-change approach with intra- and inter- tissues' gene expression data. Then, we weighted sum all the differential expression scores through a sigmoid function to get differential expression enrichment score. Finally, we ranked the genes according to the enrichment score. Top ranking genes are supposed to be the potential disease-associated genes.

RESULTS

In this study, we conducted large amounts of experiments to analyze the differentially expressed genes of intra- and inter- tissues. Experimental results showed that genes differentially expressed between different tissues are more likely to be Huntington's disease-associated genes. Five disease-associated genes were selected out in this study, two of which have been reported to be implicated in Huntington's disease.

CONCLUSIONS

We proposed a novel integrative enrichment analysis method based on artificial neuron (IEAAN), which displays better prediction precision of disease-associated genes in comparison with the state-of-the-art statistical-based methods. Our comprehensive evaluation suggests that genes differentially expressed between striatum and liver tissues of health individuals are more likely to be Huntington's disease-associated genes.

A Multi-omic Huntington's Disease Transgenic Sheep-Model Database for Investigating Disease Pathogenesis

Background:

The pathological mechanism of cellular dysfunction and death in Huntington’s disease (HD) is not well defined. Our transgenic HD sheep model (OVT73) was generated to investigate these mechanisms and for therapeutic testing. One particular cohort of animals has undergone focused investigation resulting in a large interrelated multi-omic dataset, with statistically significant changes observed comparing OVT73 and control ‘omic’ profiles and reported in literature.

Objective:

Here we make this dataset publicly available for the advancement of HD pathogenic mechanism discovery.

Methods:

To enable investigation in a user-friendly format, we integrated seven multi-omic datasets from a cohort of 5-year-old OVT73 (n = 6) and control (n = 6) sheep into a single database utilising the programming language R. It includes high-throughput transcriptomic, metabolomic and proteomic data from blood, brain, and other tissues.

Results:

We present the ‘multi-omic’ HD sheep database as a queryable web-based platform that can be used by the wider HD research community (https://hdsheep.cer.auckland.ac.nz/). The database is supported with a suite of simple automated statistical analysis functions for rapid exploratory analyses. We present examples of its use that validates the integrity relative to results previously reported. The data may also be downloaded for user determined analysis.

Conclusion:

We propose the use of this online database as a hypothesis generator and method to confirm/refute findings made from patient samples and alternate model systems, to expand our understanding of HD pathogenesis. Importantly, additional tissue samples are available for further investigation of this cohort.

HAP40 protein levels are huntingtin-dependent and decrease in Huntington disease

The huntingtin-associated protein 40 (HAP40) is an abundant interactor of huntingtin (HTT). In complexes of these proteins, HAP40 tightly binds to HTT in a cleft formed by two larger domains rich in HEAT repeats, and a smaller bridge domain connecting the two. We show that HAP40 steady-state protein levels are directly dependent on HTT (both normal and mutant HTT) and that HAP40 is strongly stabilized by the interaction with HTT resulting in an at least 5-fold increase in HAP40's half-life when bound to HTT. Cellular HAP40 protein levels were reduced in primary fibroblasts and lymphoblasts of Huntington Disease (HD) patients and in brain tissue of a full-length HTT mouse model of HD, concomitant with decreased soluble HTT levels in these cell types. This data and our previous demonstration of coevolution between HTT and HAP40 and evolutionary conservation of their interaction suggest that HAP40 is an obligate interaction partner of HTT. Our observation of reduced HAP40 levels in HD invites further studies, whether HAP40 loss-of-function contributes to the pathophysiology of HD.

HDinHD: A Rich Data Portal for Huntington's Disease Research

HDinHD (Huntington's Disease in High Definition; HDinHD.org) is an open online portal for the HD research community that presents a synthesized view of HD-related scientific data. Here, we present a broad overview of HDinHD and highlight the newly launched HDinHD Explorer tool that enables researchers to discover and explore a wide range of diverse yet interconnected HD-related data. We demonstrate the utility of HDinHD Explorer through data mining of a single collection of newly released in vivo therapeutic intervention study reports alongside previously published reports.

Shared and oppositely regulated transcriptomic signatures in Huntington's disease and brain ischemia confirm known and unveil novel potential neuroprotective genes

Huntington's disease and subcortical vascular dementia display similar dementing features, shaped by different degrees of striatal atrophy, deep white matter degeneration and tau pathology. To investigate the hypothesis that Huntington's disease transcriptomic hallmarks may provide a window into potential protective genes upregulated during brain acute and subacute ischemia, we compared RNA sequencing signatures in the most affected brain areas of 2 widely used experimental mouse models: Huntington's disease, (R6/2, striatum and cortex and Q175, hippocampus) and brain ischemia-subcortical vascular dementia (BCCAS, striatum, cortex and hippocampus). We identified a cluster of 55 shared genes significantly differentially regulated in both models and we screened these in 2 different mouse models of Alzheimer's disease, and 96 early-onset familial and apparently sporadic small vessel ischemic disease patients. Our data support the prevalent role of transcriptional regulation upon genetic coding variability of known neuroprotective genes (Egr2, Fos, Ptgs2, Itga5, Cdkn1a, Gsn, Npas4, Btg2, Cebpb) and provide a list of potential additional ones likely implicated in different dementing disorders and worth further investigation.

Single-Nucleus RNA-Seq Reveals Dysregulation of Striatal Cell Identity Due to Huntington's Disease Mutations

Huntington's disease (HD) is a dominantly inherited neurodegenerative disorder caused by a trinucleotide expansion in exon 1 of the huntingtin (HTT) gene. Cell death in HD occurs primarily in striatal medium spiny neurons (MSNs), but the involvement of specific MSN subtypes and of other striatal cell types remains poorly understood. To gain insight into cell type-specific disease processes, we studied the nuclear transcriptomes of 4524 cells from the striatum of a genetically precise knock-in mouse model of the HD mutation, HttQ175/+, and from wild-type controls. We used 14- to 15-month-old male mice, a time point at which multiple behavioral, neuroanatomical, and neurophysiological changes are present but at which there is no known cell death. Thousands of differentially expressed genes (DEGs) were distributed across most striatal cell types, including transcriptional changes in glial populations that are not apparent from RNA-seq of bulk tissue. Reconstruction of cell type-specific transcriptional networks revealed a striking pattern of bidirectional dysregulation for many cell type-specific genes. Typically, these genes were repressed in their primary cell type, yet de-repressed in other striatal cell types. Integration with existing epigenomic and transcriptomic data suggest that partial loss-of-function of the polycomb repressive complex 2 (PRC2) may underlie many of these transcriptional changes, leading to deficits in the maintenance of cell identity across virtually all cell types in the adult striatum.SIGNIFICANCE STATEMENT Huntington's disease (HD) is a dominantly inherited neurodegenerative disorder characterized by specific loss of medium spiny neurons (MSNs) in the striatum, accompanied by more subtle changes in many other cell types. It is thought that changes in transcriptional regulation are an important underlying mechanism, but cell type-specific gene expression changes are not well understood, particularly at time points relevant to the onset of disease-related symptoms. Single-nucleus (sn)RNA-seq in a genetically precise mouse model enabled us to identify novel patterns of transcriptional dysregulation because of HD mutations, including bidirectional dysregulation of many cell type identity genes that may be driven by partial loss-of-function of the polycomb repressive complex (PRC). Identifying these regulators of transcriptional dysregulation in HD can be leveraged to design novel disease-modifying therapeutics.

ADAM10 hyperactivation acts on piccolo to deplete synaptic vesicle stores in Huntington's disease

Synaptic dysfunction and cognitive decline in Huntington's disease (HD) involve hyperactive A disintegrin and metalloproteinase domain-containing protein 10 (ADAM10). To identify the molecular mechanisms through which ADAM10 is associated with synaptic dysfunction in HD, we performed an immunoaffinity purification-mass spectrometry (IP-MS) study of endogenous ADAM10 in the brains of wild-type and HD mice. We found that proteins implicated in synapse organization, synaptic plasticity, and vesicle and organelles trafficking interact with ADAM10, suggesting that it may act as hub protein at the excitatory synapse. Importantly, the ADAM10 interactome is enriched in presynaptic proteins and ADAM10 co-immunoprecipitates with piccolo (PCLO), a key player in the recycling and maintenance of synaptic vesicles. In contrast, reduced ADAM10/PCLO immunoprecipitation occurs in the HD brain, with decreased density of synaptic vesicles in the reserve and docked pools at the HD presynaptic terminal. Conditional heterozygous deletion of ADAM10 in the forebrain of HD mice reduces active ADAM10 to wild-type level and normalizes ADAM10/PCLO complex formation and synaptic vesicle density and distribution. The results indicate that presynaptic ADAM10 and PCLO are a relevant component of HD pathogenesis.

Transcriptional Dysregulation in Huntington's Disease: The Role in Pathogenesis and Potency for Pharmacological Targeting

Huntington's disease (HD) is an inherited neurodegenerative disorder caused by a mutation in the gene that encodes a critical cell regulatory protein, huntingtin (Htt). The expansion of cytosine-adenine-guanine (CAG) trinucleotide repeats causes improper folding of functional proteins and is an initial trigger of pathological changes in the brain. Recent research has indicated that the functional dysregulation of many transcription factors underlies the neurodegenerative processes that accompany HD. These disturbances are caused not only by the loss of wild-type Htt (WT Htt) function but also by the occurrence of abnormalities that result from the action of mutant Htt (mHtt). In this review, we aim to describe the role of transcription factors that are currently thought to be strongly associated with HD pathogenesis, namely, RE1-silencing transcription factor, also known as neuron-restrictive silencer factor (REST/NRSF), forkhead box proteins (FOXPs), peroxisome proliferator-activated receptor gamma coactivator-1a (PGC1α), heat shock transcription factor 1 (HSF1), and nuclear factor κ light-chain-enhancer of activated B cells (NF- κB). We also take into account the role of these factors in the phenotype of HD as well as potential pharmacological interventions targeting the analyzed proteins. Furthermore, we considered whether molecular manipulation resulting in changes in transcription factor function may have clinical potency for treating HD.

A meta-analysis of transcriptomic profiles of Huntington's disease patients

Description of robust transcriptomic alterations in Huntington’s disease is essential to identify targets for biochemical studies and drug development. We analysed publicly available transcriptome data from the brain and blood of 220 HD patients and 241 healthy controls and identified 737 and 661 genes with robustly altered mRNA levels in the brain and blood of HD patients, respectively. In the brain, a subnetwork of 320 genes strongly correlated with HD and was enriched in transport-related genes. Bioinformatical analysis of this subnetwork highlighted CDC42, PAK1, YWHAH, NFY, DLX1, HMGN3, and PRMT3. Moreover, we found that CREB1 can regulate 78.0% of genes whose mRNA levels correlated with HD in the blood of patients. Alterations in protein transport, metabolism, transcriptional regulation, and CDC42-mediated functions are likely central features of HD. Further our data substantiate the role of transcriptional regulators that have not been reported in the context of HD (e.g. DLX1, HMGN3 and PRMT3) and strongly suggest dysregulation of NFY and its target genes across tissues. A large proportion of the identified genes such as CDC42 were also altered in Parkinson’s (PD) and Alzheimer’s disease (AD). The observed dysregulation of CDC42 and YWHAH in samples from HD, AD and PD patients indicates that those genes and their upstream regulators may be interesting therapeutic targets.

Taming the Huntington's Disease Proteome: What Have We Learned?

Mass spectrometry (MS) is a physical technique used to identify specific chemicals and molecules by precise analysis of their mass and charge; this technology has been adapted for biological sciences applications. Investigators have used MS to identify differential expressions of proteins in Huntington’s disease (HD), to discover Huntingtin (HTT) interacting proteins and to analyze HTT proteoforms. Using systems biology and computational approaches, data from MS screens have been leveraged to find differentially expressed pathways. This review summarizes the data from most of the MS studies done in the HD field in the last 20 years and compares it to the protein data reported before the use of MS technology. The MS results validate early findings in the field such as differential expression of PDE10a and DARPP-32 and identify new changes. We offer a perspective on the MS approach in HD, particularly for identification of disease pathways, the challenges in interpreting data across different studies, and its application to protein studies moving forward.

Integrated analysis on transcriptome and behaviors defines HTT repeat-dependent network modules in Huntington's disease

Huntington's disease (HD) is caused by a CAG repeat expansion in the huntingtin (HTT) gene. Knock-in mice carrying a CAG repeat-expanded Htt will develop HD phenotypes. Previous studies suggested dysregulated molecular networks in a CAG length genotype- and the age-dependent manner in brain tissues from knock-in mice carrying expanded Htt CAG repeats. Furthermore, a large-scale phenome analysis defined a behavioral signature for HD genotype in knock-in mice carrying expanded Htt CAG repeats. However, an integrated analysis correlating phenotype features with genotypes (CAG repeat expansions) was not conducted previously. In this study, we revealed the landscape of the behavioral features and gene expression correlations based on 445 mRNA samples and 445 microRNA samples, together with behavioral features (396 PhenoCube behaviors and 111 NeuroCube behaviors) in Htt CAG-knock-in mice. We identified 37 behavioral features that were significantly associated with CAG repeat length including the number of steps and hind limb stand duration. The behavioral features were associated with several gene coexpression groups involved in neuronal dysfunctions, which were also supported by the single-cell RNA sequencing data in the striatum and the spatial gene expression in the brain. We also identified 15 chemicals with significant responses for genes with enriched behavioral features, most of them are agonist or antagonist for dopamine receptors and serotonin receptors used for neurology/psychiatry. Our study provides further evidence that abnormal neuronal signal transduction in the striatum plays an important role in causing HD-related phenotypic behaviors and provided rich information for the further pharmacotherapeutic intervention possibility for HD.

Epigenetic regulation in Huntington's disease

Huntington's disease (HD) is a devastating and fatal monogenic neurodegenerative disorder characterized by progressive loss of selective neurons in the brain and is caused by an abnormal expansion of CAG trinucleotide repeats in a coding exon of the huntingtin (HTT) gene. Progressive gene expression changes that begin at premanifest stages are a prominent feature of HD and are thought to contribute to disease progression. Increasing evidence suggests the critical involvement of epigenetic mechanisms in abnormal transcription in HD. Genome-wide alterations of a number of epigenetic modifications, including DNA methylation and multiple histone modifications, are associated with HD, suggesting that mutant HTT causes complex epigenetic abnormalities and chromatin structural changes, which may represent an underlying pathogenic mechanism. The causal relationship of specific epigenetic changes to early transcriptional alterations and to disease pathogenesis require further investigation. In this article, we review recent studies on epigenetic regulation in HD with a focus on DNA and histone modifications. We also discuss the contribution of epigenetic modifications to HD pathogenesis as well as potential mechanisms linking mutant HTT and epigenetic alterations. Finally, we discuss the therapeutic potential of epigenetic-based treatments.