- Short report
- Open Access
Interleukin-6 deficiency exacerbates Huntington’s disease model phenotypes
Molecular Neurodegeneration volume 15, Article number: 29 (2020)
Huntington’s disease (HD) is an incurable neurodegenerative disorder caused by CAG trinucleotide expansions in the huntingtin gene. Markers of both systemic and CNS immune activation and inflammation have been widely noted in HD and mouse models of HD. In particular, elevation of the pro-inflammatory cytokine interleukin-6 (IL-6) is the earliest reported marker of immune activation in HD, and this elevation has been suggested to contribute to HD pathogenesis. To test the hypothesis that IL-6 deficiency would be protective against the effects of mutant huntingtin, we generated R6/2 HD model mice that lacked IL-6. Contrary to our prediction, IL-6 deficiency exacerbated HD-model associated behavioral phenotypes. Single nuclear RNA Sequencing (snRNA-seq) analysis of striatal cell types revealed that IL-6 deficiency led to the dysregulation of various genes associated with synaptic function, as well as the BDNF receptor Ntrk2. These data suggest that IL-6 deficiency exacerbates the effects of mutant huntingtin through dysregulation of genes of known relevance to HD pathobiology in striatal neurons, and further suggest that modulation of IL-6 to a level that promotes proper regulation of genes associated with synaptic function may hold promise as an HD therapeutic target.
The molecular mechanisms that link huntingtin mutation to neuronal cell death in Huntington’s disease (HD) are still not fully understood . However, due to extensive reports of systemic and CNS innate immune and inflammatory activation in human HD patients and HD mouse models [2,3,4,5,6,7,8,9], it has been proposed that innate immune activation may play a role in mediating the pathogenic effects of mutant huntingtin (mHTT). Upregulation of interleukin-6 (IL-6) is the earliest reported marker of immune activation in HD, as early as 16 years before the predicted onset of clinical symptoms . As the toxic properties of mHTT have been linked to elevation of nuclear factor-kappaB (NF-κB) activity , which is a potent inducer of IL-6 gene expression , mHTT-mediated upregulation of this pro-inflammatory cytokine in particular may lead to activation of neurotoxic innate immune signaling in HD. However, preclinical work in HD rodent models has suggested both a protective and negative role for IL-6 in HD [13, 14].
To investigate the role of IL-6 in protecting against or promoting the pathogenic effects of mHTT from early phenotypic stages, we crossed R6/2 transgenic exon 1 mHTT model mice  to mice that that were deficient in IL-6  and performed behavioral assays of HD motor phenotype progression (Schematic: Fig. 1). Consistent with previous reports that IL-6 −/− mice do not demonstrate differences in their spontaneous motor activity , we did not observe any differences in the rotarod performance or open field and rearing and climbing behaviors between IL-6 −/− and IL-6 +/+ wild type control mice (WT, those that did not carry the R6/2 exon1 mHTT transgene) (Fig. 2a-h and S1A-F). And consistent with reports that IL-6 exerts anti-obesity effects in rodents , IL6 −/− mice that carried the R6/2 transgene had slightly higher overall weights than their IL-6 +/+ R6/2 transgenic littermates (Fig. 2a). However, IL-6 −/− mice that carried the R6/2 transgene had more severe HD-associated behavioral symptoms than their IL-6 +/+ R6/2 transgenic littermates as assessed by rotarod performance (Fig. 2b and S1A), open field assay measurements (Fig. 2c-h), and rearing and climbing activity (Fig. S1B-F). The only assay in which IL-6 −/− mice that carried the R6/2 transgene did not have more severe HD-associated symptoms than their IL-6 +/+ R6/2 transgenic littermates was in the grip strength test (Fig. S1G). And although the IL-6 −/− mice that carried the R6/2 transgene had slightly higher overall body weights than their IL-6 +/+ R6/2 transgenic littermates (Fig. 2a), there was no significant correlation between mouse body weight and performance on the rotarod test (Fig. S2). Finally, since the R6/2 model mice are usually bred and tested in a mixed strain background (CBA x C57BL/6), we performed single nucleotide polymorphism (SNP) array genotyping to assess the strain characteristics of the F2 mice used in behavioral testing. We observed no significant correlation between the percentage similarity to the congenic C57BL6 strain background and motor testing performance (Fig. S3), and thus we conclude that whatever differences in strain background that exist among the tested F2 mice do not account for the differences in behavioral testing results, and that these are rather a result of the IL-6 genotype. Thus together our data reveal that constitutive IL-6 loss exacerbates several HD-like behavioral symptoms in the R6/2 exon1 mHTT model.
In order to investigate the molecular basis for the aggravation of HD model phenotypes upon IL-6 KO, we performed single nuclear RNA-sequencing (snRNA-seq, n = 3 per group; Methods) on nuclei isolated from the striatal tissue of the same mice that were used for behavioral testing, harvested at 10 weeks of age. We used the ACTIONet framework  to identify major cell types in the striatum across the replicate samples in each group. Using a curated set of marker genes, we recovered major expected cell types in the striatum, including the two major types of striatal neurons, direct and indirect pathway spiny projection neurons (dSPNs and iSPNs), as well as astroglia, microglia, oligodendrocytes, oligodendrocyte precursor cells, Chat-expressing cholinergic interneurons, Sst/Npy-expressing GABAergic interneurons, Pval/Th-expressing GABAergic interneurons, and endothelial and mural cells of the blood brain barrier. Since the subventricular zone was included in the striatal dissection, we also recovered ciliated ependymal cells, secretory ependymal cells, and migrating neuroblasts (Fig. 3a and Methods). Additionally, we recovered a distinct Foxp2/Olfm3-expressing neuron cluster representing a distinct neuron subtype that is characterized by expression of Olfm3, Foxp2, Adarb2, and Otof that likely represents the same striatal neuronal subtype recently characterized in other studies as expressing Otof and Olfm3 [20, 21]. Differential gene expression analysis of the most abundant identified cell types (dSPNs, iSPNs, astroglia, and oligodendrocytes) (Methods) revealed changes to gene expression occurring upon IL-6 KO in non-carrier control (WT) mice and in the R6/2 mice (Fig. 3b and Tables S1-S2). Pathway analysis of the differentially expressed genes revealed that among the top gene pathways altered in SPNs by IL-6 KO in both the R6/2 carrier and WT mice were various terms related to synaptic transmission, including terms related to the glutamatergic synapse and long-term potentiation (Fig. 4a). Chromatin Enrichment Analysis (ChEA)  for predicted regulators of these changes in SPNs revealed that Stat3 was among the top predicted regulators of these gene expression changes in SPNs (Fig. 4b). In support of Stat3 directly linking IL-6 deficiency to alterations in genes involved in synaptic transmission, Stat3 is not only one of the canonical transcription factors that transduces IL-6 receptor signaling  but also has been shown to have a role in the regulation synaptic plasticity . In addition, we noted that the brain-derived neurotrophic factor (BDNF) receptor Ntrk2 was downregulated in WT and R6/2 SPNs that lacked IL-6, and several synaptic protein-encoding genes such as Nrxn1, Dlg2, Cntnap2, and Gabrg3 were among the most downregulated genes in R6/2 SPNs that lacked IL-6 (Fig. 3b and Tables S1-S2). Although we did not recover enough microglial cells to determine with high confidence genes that were differentially expressed in striatal microglial cells upon IL-6 deficiency, analysis of gene pathways and predicted regulators of gene expression changes that were detected in striatal astroglia and oligodendrocytes did not show evidence of alteration of glial inflammatory or innate immune signaling pathways (Figs. S4). Although it is possible that glial populations in mouse models of HD may not fully recapitulate changes to gene expression that are observed in human glial populations in HD, our gene expression data suggest that IL-6 deficiency does not alter innate immune activation in an HD context, but rather may aggravate HD model phenotypes in part by dysregulation of genes related to synaptic transmission and neurotrophin signaling, two pathways that have been linked to mHTT pathogenesis [25, 26].
In conclusion, our data reveal that a constitutive KO of IL-6, which is normally expressed in both neurons and glia in the CNS , including in the striatum , exacerbates several HD-related behavioral phenotypes in the R6/2 exon1 mHTT model of HD. Combined with reports of a protective effect of IL-6 in models of traumatic brain injury [29, 30], Parkinson’s disease [31, 32], and in modulating Aβ deposition in Alzheimer’s disease models  it is thus possible that the clinical reports of IL-6 elevation in HD at early disease timepoints may reflect a protective rather than pathogenic alteration, and that elevation of IL-6 levels, in certain ranges and timepoints, may have therapeutic benefit in HD.
All animal experiments were conducted with the approval of the Massachusetts Institute of Technology Animal Care and Use Committee. Mice were housed with food and water provided ad libitum on a standard 12 h light/12 h dark cycle. All mice were obtained from the Jackson Laboratory (Bar Harbor, ME). R6/2 HD model mice (B6CBA-Tg (HDexon1)62Gpb/1 J, Jackson Laboratory stock #002810) were crossed with IL-6 knockout mice maintained the C57BL/6 background (B6.129S2-Il6tm1Kopf/J, Jackson Laboratory stock #002650). F1 heterozygous IL-6 −/+ x R6/2 hemizygous mice were crossed to obtain the F2 generation, which was used for experiments between 4 and 10 weeks of age. Given the two-strain mixed background of the original R6/2 breeders (CBA and C57BL/6), strain characteristics for each mouse used in behavioral testing were assayed by single nucleotide polymorphism (SNP) genotyping (Transnetyx, Cordova, TN, Full Strain Genetic Monitoring, 120 SNPs with markers on all autosomes); these data are presented in Fig. S3. Beginning at 4 weeks of age, animals were monitored weekly for body condition score and weight. Both male and female mice were used for the behavioral and biochemical analyses. Behavioral testing experiments were performed on mice at the time points outlined in Fig. 1 and listed in the Behavioral Testing section below, and mice were otherwise naïve to the testing or prior analyses.
Behavioral analyses were performed on the F2 animals as previously described . Behavioral testing was conducted and analyzed by an investigator blinded to genotype. Statistical testing was performed using GraphPad Prism 8 with the number of mice (n) per group and the statistical measures performed as reported in detail in the figure legends. Behavioral data is presented as mean ± standard error of the mean unless otherwise stipulated. Mice were identified as outliers and excluded from further analysis if they scored ±2 standard deviations from the mean on multiple behavioral tests.
Rotarod testing was performed using an accelerating rotarod (Med Associates, St. Albans, VT). Mice were trained using 3 consecutive 5 min sessions with a fixed speed of 20 rotations per minute (RPM) and gently placed back on the rotarod after each fall. Mice were given 1 min of rest between training trials. At 5 weeks of age, testing was performed on 5 consecutive days with the rotarod accelerating from 5 RPM to 40 RPM over the course of 5 min. Latency to fall was measured as the time from the beginning of a trial until the mouse fell off the rod or completed two or more passive rotations. On subsequent weeks mice were re-tested on the accelerating rotarod on a single day to assess their motor phenotype over time.
Open Field testing was performed at 7 and 10 weeks of age in 60-min sessions each week per mouse using an infrared photobeam open field chamber. The field had 16 infrared beams spaced regularly along the x, y, and z axes (#MED-OFAS-RSU, Med Associates, St. Albans, VT). Data was analyzed for distance traveled, time traveling, vertical activity and vertical time.
Rearing and Climbing analysis was performed by placing mice under an overturned black metal mesh pencil cup (Rolodex #82406) 4.375 in. in diameter and 5.5 in. in height. The latency until the mouse reared on hind legs and touched the mesh with its front paws as well as the latency until the mouse climbed on the mesh with all paws off of the ground were recorded. All subsequent rearing and climbing events were counted, and summed together to get the total vertical activity. Mice were given one 5 min trial at 8 weeks of age.
Grip strength was performed using a grip strength meter (Ugo Basile,Varese, Italy). Briefly. The mice were suspended by their tails and allowed to grab the measurement bar. They were then pulled away from the bar by the tail until they released the bar and the maximum force (g) was recorded. Each mouse was given 5 trials. Trials where the mouse failed to grasp the bar with two hands were excluded from subsequent analyses. The average maximum force of the 5 trials was used for each mouse.
EnrichR pathway analysis
Pathway and chromatin enrichment analysis was performed using the EnrichR package [35, 36] considering only protein-coding genes. Significant pathways were identified by Fisher’s exact test with adjusted p-value < 0.05.
Single nuclear (snRNA) RNA sequencing and analysis
Nuclei isolation protocol was adapted from . Briefly, striata were dissected and flash frozen in liquid nitrogen. Frozen tissue was homogenized in 700 μL of homogenization buffer with a 2 mL KIMBLE Dounce tissue grinder (MilliporeSigma, Burlington MA) using 10 strokes with loose pestle followed by 10 strokes with tight pestle. Homogenized tissue was filtered through a 40 μm cell strainer and mixed with 450 μL of working solution (50% OptiPrep density gradient medium (MilliporeSigma, Burlington MA), 5 mM CaCl2, 3 mM Mg(CH3COO)2, 10 mM Tris HCl [pH 7.8], 0.1 mM EDTA [pH 8.0], and 1 mM β-mercaptoethanol). Nuclei were pelleted at the interface of an OptiPrep density gradient containing 750 μL of 30% OptiPrep Solution on top of 300 μL of 40% OptiPrep Solution inside a Sorenson Dolphin microcentrifuge tube (MilliporeSigma, Burlington MA) by centrifugation at 10,000 x g for 5 min at 4 °C using a fixed angle rotor (FA-45-24-11-Kit). The nuclear pellet was collected at the interface and washed with 2% BSA (in 1x PBS) containing 0.12 U/μL SUPERase In RNase Inhibitor. The nuclei were pelleted by centrifugation at 300 x g for 3 min at 4 °C using a swing-bucket rotor (S-24-11-AT). Nuclei were washed three times with 2% BSA and centrifuged under the same conditions. The nuclear pellet was re-suspended in 100 μL of 2% BSA.
Droplet-based snRNA sequencing libraries were prepared using the Chromium Single Cell 3′ Reagent Kit v3 (10x Genomics, Pleasanton CA) according to the manufacturer’s protocol and sequenced on a NovaSeq 6000 at the Broad Institute Genomics Platform. FASTQ files were aligned to the pre-mRNA annotated Mus musculus reference genome version GRCm38.
The R package batchelor , was used to correct for batch effects observed across biological replicates within each experimental group from the count matrix. Batch-corrected data was used as input to the archetypal analysis for cell type identification (ACTION) algorithm  to identify a set of landmark cells or ‘archetypes’, each representing a potential underlying cell state. Using ACTION-decompositions with varying numbers of archetypes, we employed our recently developed ACTION-based network (ACTIONet) framework  to create a multi-resolution nearest neighbor graph. ACTIONet graphs were visualized using a modified version of the stochastic gradient descent-based layout method in the uniform manifold approximation and projection (UMAP) algorithm . A curated set of known cell type-specific markers was used to annotate individual cells with their expected cell type and assign a confidence score to each annotation. Network connectivity was used to correct low-confidence annotations. Multiple iterations of this process were performed to identify and prune low quality cells, such as those with ambiguous profiles resembling dissimilar cell types (generally corresponding to doublet nuclei), or cells corresponding to nodes with a low k-core number in the network (generally corresponding to high ambient RNA content or doublet nuclei).
Cell-wise gene counts were normalized and log-scaled using the R package scran using cell type and genotype as normalization factors. Differential gene expression analysis was performed using Wilcoxon rank-sum test with the R package presto. Genes were considered differentially expressed if they had an absolute log-fold change > 0.1 with FDR < 0.001 vs. the respective control in each statistical test.
GraphPad Prism 8 was used to perform the statistical analysis of the behavioural testing data experiments. Specific information on the number (n) of values used as well as the statistical tests applied to the data can be found in the figures and/or figure legends.
Availability of data and materials
The datasets generated during and analyzed in the current study are available in the NCBI GEO repository, under accession # GSE143355.
Brain-derived neurotrophic factor
Bovine serum albumin
Cyclic adenosine monophosphate
Chromatin enrichment analysis
Central nervous system
Direct pathway spiny projection neuron
First filial generation
Second filial generation
False discovery rate
Indirect pathway spiny projection neuron
Kyoto encyclopedia of genes and genomes
Neurotrophic receptor tyrosine kinase 2
Phosphate buffered saline
Restricted maximum likelihood
Revolutions per minute
single nuclear RNA sequencing
Single nucleotide polymorphism
Spiny projection neuron
Signal transducer and activator of transcription 3
Uniform manifold approximation and projection
Unique molecular identifiers
Jimenez-Sanchez M, Licitra F, Underwood BR, Rubinsztein DC. Huntington's disease: mechanisms of pathogenesis and therapeutic strategies. Cold Spring Harbor perspectives in medicine. 2017;7.
Neueder A, Bates GP. A common gene expression signature in Huntington's disease patient brain regions. BMC Med Genet. 2014;7:60.
Leblhuber F, Walli J, Jellinger K, Tilz GP, Widner B, Laccone F, Fuchs D. Activated immune system in patients with Huntington's disease. Clin Chem Lab Med. 1998;36:747–50.
Labadorf A, Hoss AG, Lagomarsino V, Latourelle JC, Hadzi TC, Bregu J, MacDonald ME, Gusella JF, Chen JF, Akbarian S, et al. RNA sequence analysis of human Huntington disease brain reveals an extensive increase in inflammatory and developmental gene expression. PLoS One. 2015;10:e0143563.
Hensman Moss DJ, Flower MD, Lo KK. Miller JR, van Ommen GB, t Hoen PA, stone TC, Guinee a, Langbehn DR, Jones L, et al: Huntington's disease blood and brain show a common gene expression pattern and share an immune signature with Alzheimer's disease. Sci Rep. 2017;7:44849.
Dalrymple A, Wild EJ, Joubert R, Sathasivam K, Bjorkqvist M, Petersen A, Jackson GS, Isaacs JD, Kristiansen M, Bates GP, et al. Proteomic profiling of plasma in Huntington's disease reveals neuroinflammatory activation and biomarker candidates. J Proteome Res. 2007;6:2833–40.
Hodges A, Strand AD, Aragaki AK, Kuhn A, Sengstag T, Hughes G, Elliston LA, Hartog C, Goldstein DR, Thu D, et al. Regional and cellular gene expression changes in human Huntington's disease brain. Hum Mol Genet. 2006;15:965–77.
Wild E, Magnusson A, Lahiri N, Krus U, Orth M, Tabrizi SJ, Bjorkqvist M. Abnormal peripheral chemokine profile in Huntington's disease. PLoS currents. 2011;3:RRN1231.
Chang KH, Wu YR, Chen YC, Chen CM. Plasma inflammatory biomarkers for Huntington's disease patients and mouse model. Brain Behav Immun. 2015;44:121–7.
Bjorkqvist M, Wild EJ, Thiele J, Silvestroni A, Andre R, Lahiri N, Raibon E, Lee RV, Benn CL, Soulet D, et al. A novel pathogenic pathway of immune activation detectable before clinical onset in Huntington's disease. J Exp Med. 2008;205:1869–77.
Khoshnan A, Ko J, Watkin EE, Paige LA, Reinhart PH, Patterson PH. Activation of the IkappaB kinase complex and nuclear factor-kappaB contributes to mutant huntingtin neurotoxicity. J Neurosci. 2004;24:7999–8008.
Libermann TA, Baltimore D. Activation of interleukin-6 gene expression through the NF-kappa B transcription factor. Mol Cell Biol. 1990;10:2327–34.
Bensadoun JC, de Almeida LP, Dreano M, Aebischer P, Deglon N. Neuroprotective effect of interleukin-6 and IL6/IL6R chimera in the quinolinic acid rat model of Huntington's syndrome. Eur J Neurosci. 2001;14:1753–61.
Bouchard J, Truong J, Bouchard K, Dunkelberger D, Desrayaud S, Moussaoui S, Tabrizi SJ, Stella N, Muchowski PJ. Cannabinoid receptor 2 signaling in peripheral immune cells modulates disease onset and severity in mouse models of Huntington's disease. J Neurosci. 2012;32:18259–68.
Mangiarini L, Sathasivam K, Seller M, Cozens B, Harper A, Hetherington C, Lawton M, Trottier Y, Lehrach H, Davies SW, Bates GP. Exon 1 of the HD gene with an expanded CAG repeat is sufficient to cause a progressive neurological phenotype in transgenic mice. Cell. 1996;87:493–506.
Kopf M, Baumann H, Freer G, Freudenberg M, Lamers M, Kishimoto T, Zinkernagel R, Bluethmann H, Kohler G. Impaired immune and acute-phase responses in interleukin-6-deficient mice. Nature. 1994;368:339–42.
Braida D, Sacerdote P, Panerai AE, Bianchi M, Aloisi AM, Iosue S, Sala M. Cognitive function in young and adult IL (interleukin)-6 deficient mice. Behav Brain Res. 2004;153:423–9.
Wallenius V, Wallenius K, Ahren B, Rudling M, Carlsten H, Dickson SL, Ohlsson C, Jansson JO. Interleukin-6-deficient mice develop mature-onset obesity. Nat Med. 2002;8:75–9.
Mohammadi S, Davila-Velderrain J, Kellis M: Multi-resolution single-cell state characterization via joint archetypal/network analysis. bioRxiv 2019.
Saunders A, Macosko EZ, Wysoker A, Goldman M, Krienen FM, de Rivera H, Bien E, Baum M, Bortolin L, Wang S, et al. Molecular diversity and specializations among the cells of the adult mouse brain. Cell. 2018;174:1015–30 e1016.
Martin A, Calvigioni D, Tzortzi O, Fuzik J, Warnberg E, Meletis K. A Spatiomolecular map of the striatum. Cell Rep. 2019;29:4320–33 e4325.
Lachmann A, Xu H, Krishnan J, Berger SI, Mazloom AR, Ma'ayan A. ChEA: transcription factor regulation inferred from integrating genome-wide ChIP-X experiments. Bioinformatics. 2010;26:2438–44.
Zhong Z, Wen Z, Darnell JE Jr. Stat3: a STAT family member activated by tyrosine phosphorylation in response to epidermal growth factor and interleukin-6. Science. 1994;264:95–8.
Nicolas CS, Peineau S, Amici M, Csaba Z, Fafouri A, Javalet C, Collett VJ, Hildebrandt L, Seaton G, Choi SL, et al. The Jak/STAT pathway is involved in synaptic plasticity. Neuron. 2012;73:374–90.
Zuccato C, Cattaneo E. Brain-derived neurotrophic factor in neurodegenerative diseases. Nat Rev Neurol. 2009;5:311–22.
Plotkin JL, Surmeier DJ. Corticostriatal synaptic adaptations in Huntington's disease. Curr Opin Neurobiol. 2015;33:53–62.
Erta M, Quintana A, Hidalgo J. Interleukin-6, a major cytokine in the central nervous system. Int J Biol Sci. 2012;8:1254–66.
Gadient RA, Otten U. Expression of interleukin-6 (IL-6) and interleukin-6 receptor (IL-6R) mRNAs in rat brain during postnatal development. Brain Res. 1994;637:10–4.
Swartz KR, Liu F, Sewell D, Schochet T, Campbell I, Sandor M, Fabry Z. Interleukin-6 promotes post-traumatic healing in the central nervous system. Brain Res. 2001;896:86–95.
Penkowa M, Giralt M, Lago N, Camats J, Carrasco J, Hernandez J, Molinero A, Campbell IL, Hidalgo J. Astrocyte-targeted expression of IL-6 protects the CNS against a focal brain injury. Exp Neurol. 2003;181:130–48.
Bolin LM, Strycharska-Orczyk I, Murray R, Langston JW, Di Monte D. Increased vulnerability of dopaminergic neurons in MPTP-lesioned interleukin-6 deficient mice. J Neurochem. 2002;83:167–75.
Akaneya Y, Takahashi M, Hatanaka H. Interleukin-1 beta enhances survival and interleukin-6 protects against MPP+ neurotoxicity in cultures of fetal rat dopaminergic neurons. Exp Neurol. 1995;136:44–52.
Chakrabarty P, Jansen-West K, Beccard A, Ceballos-Diaz C, Levites Y, Verbeeck C, Zubair AC, Dickson D, Golde TE, Das P. Massive gliosis induced by interleukin-6 suppresses Abeta deposition in vivo: evidence against inflammation as a driving force for amyloid deposition. FASEB journal : official publication of the Federation of American Societies for Experimental Biology. 2010;24:548–59.
Hachigian LJ, Carmona V, Fenster RJ, Kulicke R, Heilbut A, Sittler A. Pereira de Almeida L, Mesirov JP, Gao F, Kolaczyk ED, Heiman M: control of Huntington's disease-associated phenotypes by the striatum-enriched transcription factor Foxp2. Cell Rep. 2017;21:2688–95.
Kuleshov MV, Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, Koplev S, Jenkins SL, Jagodnik KM, Lachmann A, et al. Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016;44:W90–7.
Chen EY, Tan CM, Kou Y, Duan Q, Wang Z, Meirelles GV, Clark NR, Ma'ayan A. Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC bioinformatics. 2013;14:128.
Mathys H, Davila-Velderrain J, Peng Z, Gao F, Mohammadi S, Young JZ, Menon M, He L, Abdurrob F, Jiang X, et al. Single-cell transcriptomic analysis of Alzheimer's disease. Nature. 2019;570:332–7.
Haghverdi L, Lun ATL, Morgan MD, Marioni JC. Batch effects in single-cell RNA-sequencing data are corrected by matching mutual nearest neighbors. Nat Biotechnol. 2018;36:421–7.
Mohammadi S, Ravindra V, Gleich DF, Grama A. A geometric approach to characterize the functional identity of single cells. Nat Commun. 2018;9:1516.
Becht E, Mcinnes L, Healy J, Dutertre CA, Kwok IWH, Ng LG, Ginhoux F, Newell EW: Dimensionality reduction for visualizing single-cell data using UMAP. Nature biotechnology 2018.
RNA sequencing was performed at the Broad Institute Genomics Platform, and we thank Vanessa Lau for assisting with sample preparation for genotyping.
This work was supported by NIH/NINDS award 1 R01 NS085880 (M. Heiman), an award from the JPB Foundation (M. Heiman), an award from the Bev Hartig Huntington’s Disease Foundation (M. Heiman), a Fay/Frank Seed Award from the Brain Research Foundation (M. Heiman), the Jeptha H. and Emily V. Wade Award (M. Heiman), NIH awards R01 AG058002, U01 NS110453, UG3 NS115064 (M.K.), and fellowships from the JPB Foundation and the Hereditary Disease Foundation (M.H.W.), the NIH (5T32EB019940–05, S.S.P.), and the JPB Foundation (H.L.).
Ethics approval and consent to participate
All animal experiments were conducted with the approval of the Massachusetts Institute of Technology Animal Care and Use Committee.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Wertz, M.H., Pineda, S.S., Lee, H. et al. Interleukin-6 deficiency exacerbates Huntington’s disease model phenotypes. Mol Neurodegeneration 15, 29 (2020). https://doi.org/10.1186/s13024-020-00379-3
- Huntington’s disease