Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

Similar striatal gene expression profiles in the striatum of the YAC128 and HdhQ150 mouse models of Huntington’s disease are not reflected in mutant Huntingtin inclusion prevalence

BMC Genomics volume 16, Article number: 1079 (2015) Cite this article

Abstract

Background

The YAC128 model of Huntington’s disease (HD) shows substantial deficits in motor, learning and memory tasks and alterations in its transcriptional profile. We examined the changes in the transcriptional profile in the YAC128 mouse model of HD at 6, 12 and 18 months and compared these with those seen in other models and human HD caudate.

Results

Differential gene expression by genotype showed that genes related to neuronal function, projection outgrowth and cell adhesion were altered in expression. A Time-course ANOVA revealed that genes downregulated with increased age in wild-type striata were likely to be downregulated in the YAC128 striata. There was a substantial overlap of concordant gene expression changes in the YAC128 striata compared with those in human HD brain. Changes in gene expression over time showed fewer striatal YAC128 RNAs altered in abundance than in the HdhQ150 striata but there was a very marked overlap in transcriptional changes at all time points. Despite the similarities in striatal expression changes at 18 months the HdhQ150 mice showed widespread mHTT and ubiquitin positive inclusion staining in the striatum whereas this was absent in the YAC128 striatum.

Conclusions

The gene expression changes in YAC128 striata show a very closely matched profile to that of HdhQ150 striata and are already significantly different between genotypes by six months of age, implying that the temporal molecular gene expression profiles of these models match very closely, despite differences in the prevalence of brain inclusion formation between the models. The YAC128 gene expression changes appear to correlate well with gene expression differences caused by ageing. A relatively small number of genes showed significant differences in expression between the striata of the two models and these could explain some of the phenotypic differences between the models.

Background

Huntington’s disease (HD) is a fatal progressive neurodegeneration with motor, cognitive and psychiatric manifestations. It is caused by an expansion of a CAG triplet repeat in exon 1 of the HTT gene, which is translated to give an expanded glutamine tract at the N-terminus of the protein, huntingtin (HTT) [1]. A series of genetic mouse models of the disease have been generated using various technologies to give transgenic and knock in models of the disease that include both truncated and full-length Htt [2–8]. These models have been tested in multiple behavioural paradigms and show deficits in tests of motor ability and in cognitive and behavioural assessments [5, 9–16].

Changes in gene expression have also been seen in the brains and other tissues of the mouse models of HD. These changes show a substantial overlap between the various mouse models tested and also overlap with gene expression changes seen in human HD brain [17]. The profiles of genes with reduced expression appear to overlap to greater extent between models and between models and human brain than those with increased expression [17, 18]. Treatments that alleviate the decreased expression can improve the phenotype in mouse models, whether directed at transcriptional mechanisms or not [19–23].

The YAC128 model of HD carries a full length human HTT gene in a GSE70656 artificial chromosome, and extensive behavioural testing has shown that it displays substantial phenotypes that correlate with human HD symptoms from a relatively early age [13, 14, 16, 24, 25] including affective symptoms [26]. Deficits on the rotarod were observed from 4 months of age and persisted throughout life, and deficits on the balance beam occurred from 8 months [27]. In the water maze, reduced ability to find the hidden platform was seen at 8–10 months with reversal learning showing deficits from 4 months [27]. Deficits in reversal learning in a set-shifting task were seen by 6 months and of extra-dimensional set-shifting at 16 months [28], although no implicit learning deficits were observed in these mice [29]. Despite these early changes in the behavioural phenotype of the mice, frank neuronal inclusions were only visible from 15 months of age, although diffuse staining with antibodies that detect inclusions was seen from 12 months of age [30].

As HD affects the caudate and putamen earliest and most profoundly and shows substantial gene expression dysregulation [18] we chose to examine mouse striatum, as the nearest equivalent, in order to investigate whether dysregulated gene expression also occurred in the YAC128 striatum. We examined global gene expression changes in the striata from wild-type (WT) and transgenic mice from the YAC128 mouse line at 6, 12 and 18 months. We observed alterations in gene expression at all time points, which overlapped with changes seen in human HD brain and in other mouse models of the disease. We observed that genes downregulated with age in normal mice tended to be even more decreased in expression in the YAC128 striata.

Results

The effects of age on gene expression

Between 6 and 12 months, 2469 mRNAs (probesets: 1272 distinct genes) are altered in abundance in the WT animals and 1266 mRNAs (800 genes) in the YAC128 animals (nominal p <0.05). Of these, 241 mRNAs (156 genes), more than expected by chance, occur in both sets (9.8 % and 19.0 % of probesets respectively, p = 0.045). Between 12 and 18 months, more mRNAs are altered: 2579 mRNA 1352 genes) in WT mice and 3019 (2211 genes) in YAC128 mice, of which 522 mRNA (398 genes) are common to both cohorts (20.2 % and 17.3 % of probesets respectively, p <10−4). The overlapping probesets and corresponding genes are given in Additional file 1: Table S1. A GO term enrichment analysis (Table 1, Additional file 2: Table S2) shows that few pathways show an over-representation of significantly differentially expressed genes: only translation (GO:0006412) in the WT 6–12 month striata and cell adhesion (GO:0007155) in the YAC128 6–12 month striata are close to significance. Translation does not appear as even nominally significant in the YAC128 data, nor cell adhesion in the WT data (Additional file 2: Table S2). Between 12 – 18 months several processes are significantly over-represented in both WT and YAC128 striata (Table 1), including neuronal processes in both cohorts. Cell adhesion (GO:0007155) is significantly over-represented in the YAC128 striata over this time period.

Table 1 Pathways altered with age in mouse striatum in YAC128 mice
Full size table

Effects of genotype on gene expression

Analysis between genotypes shows that 2557 probesets (1821 genes) are dysregulated between YAC128 and WT animals when using data from all time points (nominal p < 0.05), of which 1151 are up-regulated and 1106 are down-regulated in the YAC128 striata. Using an FDR threshold of p > 0.05, 87 probesets are dysregulated (Table 2), 53 down- and 34 up-regulated. Analysing the time points individually, at 6 months 1287 (856 genes), at 12 months 885 (482 genes) and at 18 months 2484 probesets (1514 genes) are altered in expression (nominal p < 0.05; Additional file 3: Table S3). Between 6 and 12 months 148 probesets (118 genes) are altered in abundance at both times (11.5 % and 6.0 % respectively; p = 0.02) and between 12 and 18 months there are 223 such probesets (164 genes, 25.2 % and 9.0 % respectively; p = 0.039) (Additional file 3: Table S3).

Table 2 Probesets dysregulated between genotypes
Full size table

The TANOVA analysis (Fig. 1, Table 3) highlights a set of 146 probesets (125 genes) dysregulated between genotypes but also includes the trajectory of change in expression over time. Seventy-six genes are common to the TANOVA and LIMMA FDR p < 0.05 analyses. The TANOVA results show that if the trajectory of gene expression with increased age is towards decreasing mRNA abundance, then these genes are much more likely to be decreased in expression in the striata of YAC128 mice compared with wild-type: 39 such differentially expressed genes decreased in expression in YAC128 mice compared with 4 increased in expression. The opposite is true for genes increased in expression with age in the YAC128 mouse striatum: these are more likely to be increased in expression in the YAC128 striata: 48 increased compared with 22 decreased mRNAs. The full TANOVA plots of gene expression over time are given in Additional file 4: Figure S1.

Fig. 1
figure 1

Patterns of expression identified using TANOVA. Each different pattern is illustrated by one gene showing a significant change in expression over time and between genotypes. The gene ID and the Affymetrix probeset ID are given. Expression is given as Log-2 fold change and the time points are in months. WT = wild-type and HD = YAC128

Full size image
Table 3 Differentially expressed genes between WT and YAC128 striata identified using TANOVA
Full size table

An over-representation analysis using EASE and DAVID [31, 32] of the differentially expressed genes by TANOVA does not reveal any over-represented categories, probably due to the small numbers of genes within each pattern of change. However, an over-representation analysis of the full sets of genes identified by LIMMA analysis (Table 4) shows that the largest number of pathways are identified by the down-regulated genes and that these are mostly related to G-protein and other intracellular signalling pathways. An examination of the genes that contribute to these significant processes shows that there are substantial overlaps of genes amongst these significant categories and that the most specific pathway highlighted (the smallest) is GO:0019226, transmission of nerve impulse. Categories relating to nervous transmission and synaptic events appear as nominally significant even if they are not signficant once the FDR correction has been applied (Additional file 5: Table S4). In order to highlight the most significant functional relationships in the data we conducted a DAVID analysis [31, 32] visualised in cytoscape in Fig. 2. The most interconnected ontological term is membrane and all the other processes are related to membrane events: these include cell adhesion, neuronal projections, synaptic functions and transmission of nerve impulse.

Table 4 GO terms showing over-representation amongst genes differentially expressed in YAC128 compared with WT striata
Full size table
Fig. 2
figure 2

Pathways implicated by the gene expression differences in YAC128 mouse striata. List of differentially expressed genes were generated by t-test using a nominal p < 0.01; DAVID was used for the pathway analysis with a pathway filter q < 0.05. The pathway node size is proportional to gene membership and the edges joining nodes are weighted by gene overlap between nodes. Terms have been clustered into groups that contain 90 % genetic similarity on average. The most significant term of the cluster has been displayed. These terms have been further clustered at a level of 50 % average genetic similarity and colour-coded by cluster

Full size image

Comparison with human HD and other mouse models

A straightforward examination of the overlap between the YAC128 and HdhQ150 FDR controlled lists of differentially expressed genes shows that more genes are common to the two lists than might be expected (p <10−4; Fig. 3a). However this only examines a very small part of each genelist and the relationship of alterations in gene expression extends over a much more substantial proportion of the genes surveyed: in addition the simple analysis above does not take direction of change in gene expression into account. It is also hard to compare across species and gene expression platforms. To gain a more detailed understanding of the relationship between the striatal gene expression changes between the YAC128, human brain and other mouse model striata we examined the overlap by direction in ranked bins of the gene expression differences. The overlap between the top 1000 probesets altered in abundance between the YAC128 striata and human caudate is substantial (179/1000 probesets, p = 0.018, 80/674 orthologous genes (that are coding genes in both species)) [18]. There is also a significant overlap with human cerebellum (117/1000 probesets, p = 0.008, 47/674 orthologous genes), human BA4 cortex (124/1000 probesets, p =0.016, 55/674 orthologous genes) and BA9 cortex (71/1000 probesets, p = =0.025, 32/674 genes). The YAC128 striatal genes altered in abundance overlapped significantly with those seen in the HdhQ150 striata (272/1000 probesets, p <10−4, 132/779 genes) and also R6/1 (113/1000 probesets, p = 0.001, 103/779 genes) and R6/2 whole brain (65/1000 probesets, p = 0.012, 62/779 genes) [19, 33]. The direction of these changes is also largely concordant (Fig. 3b). Conducting the comparison at the individual time-points reveals that the substantial overlap between the expression profile of these tissues is already significant at 6 months and becomes increasingly significant over time (Additional file 6: Figure S2). A direct comparison of the YAC128 and HdhQ150 gene expression profiles at the time points for which gene expression data are available shows that concordance and overlap increase over time although they are already substantial by 6 months (Fig. 3c). To ensure that the overlaps observed were not just a result of non-specific striatal pathology, we compared the YAC128 gene expression data in this study to that from a study of MPTP effects on gene expression which also generated gene expression profiles in mouse striatal tissue [34]. We detected no significant overlap of genes differentially expressed between the two studies (19/1000 genes, NS).

Fig. 3
figure 3

Correlation of direction of expression changes in YAC128 striata, other mouse models of HD and human HD caudate. a shows the overlap of genes within the FDR adjusteddifferentially expressed genes from the YAC128 compared with the HdhQ150 mouse striata. In B and C frequency  represents the fraction of the top 200 HdhQ150 expression changes that map to a particular bin of ranked data (1000 genes per bin) in the other dataset, which is then split to concordant or discordant direction of expression change. A higher frequency of concordant (green) rather than discordant (red) in the first bins indicates a similarity between the YAC128 caudate and other model or human HD gene expression signature. b shows the comparison of YAC128 striatal gene expression with human HD brain regions [18] and the comparison with other mouse models and c shows the comparison at the same time points with the HdhQ150 striatal gene expression.

Full size image

Despite the substantial overlap in gene expression profiles, there are genes whose variation in expression differs between the two models. Differences between genotype and model were assessed using ANOVA. This was done by fitting a regression model containing the main effects of genotype and model, together with their interaction term. A significant (p < 0.05) interaction term was taken as evidence of a difference in gene expression between models. There were 869 probesets with a significant interaction term (p < 0.05) (Table 5, Additional file 7: Table S5). Examination of enrichment does not reveal any significant functional pathways identified by these genes, but the most significant individual gene is Htt itself.

Table 5 Genes that show significantly different patterns of differential expression between YAC128 and HdhQ150 striata
Full size table

Comparison of inclusion prevalence between mouse models

To assess whether the gene expression changes correlated with HTT-specific inclusions in the YAC128 and HdhQ150 mouse models we examined mutant HTT and ubiquitin distribution in mouse brain using immunohistochemistry (Fig. 4). All regions of the HdhQ150 mouse brain showed widespread S830 mutant HTT positive inclusions (Fig. 4a). In contrast the YAC128 mouse brain at 18 months showed few frank inclusions (Fig. 4c), though there is nuclear filling in all regions and a few small inclusions in the hippocampus. The HTT in the Q150 mice is entirely mouse and all mutant, whereas YAC128s have human mutant HTT and endogenous mouse HTT. In order to ascertain that the S830 was detecting similar inclusions we also assessed inclusion load using anti-ubiquitin (Fig. 4b and d). The anti-ubiquitin antibody also detects the inclusions in the Q150 mouse brain, though they are less frequent than the S830 positive inclusions, showing that the inclusions shown by the HdhQ150 brain are not a direct result of the relevant mutant HTT epitope only being available to the S830 antibody in the HdhQ150 but not the YAC128 brains.

Fig. 4
figure 4

Comparative mutant HTT pathology in the HdhQ150 (panels a, b, e, f) and YAC128 (panels c, d, g, h)  mouse brain. Comparative immunohistochemistry in brains from 18 month old YAC128 or HdhQ150 mice using either the mutant HTT specific antibody S830 or the anti-ubiquitin at 100x magnification

Full size image

Discussion

The data demonstrate that the gene expression changes in the YAC128 mouse striata are similar to those in the HdhQ150 striata at all ages [35]. These results are directly comparable as samples from mice of the same age were arrayed on the same chip. The similarity increases with age which is most likely to indicate a convergent gene expression phenotype in the striata of these models as the effects of the mutation become more marked and overcome differences between the models. It also demonstrates that the trajectory of molecular changes in the striata of these two differently constructed model lines parallel each other very well, despite differences apparent in the onset of their phenotypic changes [27, 28, 36].

The similarities in RNA changes are also paralleled in the human HD brain, especially in the caudate [18]. It is unsurprising that the human caudate gene expression profile should more closely parallel that of the mouse striatum than those of the human cortex or cerebellum, as the mouse caudate is contained within the striatal tissue analysed. These data therefore indicate that the gene expression profiles of these two mouse models of HD are convergent over time, and this is reinforced by comparisons with other models and with human caudate.

One major difference between these two models of HD is in the development of HTT-positive inclusions [30, 37]. The YAC128 model develops frank nuclear inclusions in the brain relatively late, from 15 months of age onwards, whereas the HdhQ150 model has such inclusions present from 5 months of age. The reason for this difference is not understood, though it is possible that the human and mouse proteins show differences in reactivity to the S830 antibody: however, the ubiquitin immunohistochemistry also shows that intranuclear inclusions are present in the HdhQ150 but not the YAC128 mouse brain at 5 months. Fewer ubiquitin-positive inclusions are observed than S830 positive inclusions which is consistent with previous data indicating that ubiquitin positive staining is a later event than mutant-HTT positive staining of inclusions in mouse brain [38]. However, the trajectory of nuclear filling with mHTT immunoreactivity followed by inclusion formation is common to both models and thus it seems most likely that this is a true difference in mHTT aggregation in these models. The concentrations of mHTT present in the YAC128 brain are probably similar to that in the HdhQ150 model as we noted reduced expression of mHTT in the HdhQ150 model such that it probably only expresses around half the WT endogenous level in brain [13, 39]. However, the presence of normal mouse HTT might inhibit inclusion formation. If endogenous mouse HTT is recruited into the inclusions then the presence of heterozygous protein differences is known to slow aggregation of the cognate proteins in other neurodegenerative diseases where protein inclusions are present in the disease, for instance in prion-related disease [40, 41]. While we have noted there are strong similarities in the striatal gene expression changes between these models, there are also differences, although the differences do not highlight any specific functional pathways: nevertheless these differences might contribute to the difference in inclusions observed. Discordant inclusion formation and pathogenic effects have been noted previously in a number of different systems [42, 43] though the connection with similar gene expression profiles in the face of differential inclusion formation has not been made previously.

The significant differences between the gene expression profiles of the striata from the YAC128 and HdhQ150 lines do not highlight any obvious functional differences between the two models that could account for the observed differences in inclusion prevalence. Grk4 is close to the Htt locus in mice but it is not contained within the YAC construct used to generate these mice [13]: thus this cannot explain the increased expression of this gene. Other than this the significantly altered genes are not in the same chromosomal locations, so no chromosomally specific effect can be inferred that is due to the direct action of the transgene. The gene most significantly altered in expression is Htt itself and this is expected in light of the substantial down-regulation of mHTT in the HdhQ150 homozygous knock in striata [39].

Most of the genes that are differentially expressed between the models show increased expression in YAC128 striata and it is possible that this is due to the increased expression of HTT itself over the endogenous HTT expression. Mutant HTT can rescue the lethal effects of knocking out WT HTT [4, 44] so the mutant protein may well be exerting functional effects in relation to its normal as well as its pathological function. The expression of the YAC transgene in the YAC128 animals is close to the level of expression of the endogenous mouse gene [45]. The differential changes could potentially be related to the normal function of HTT through the effects of increased overall huntingtin expression and thus enhancement of the normal function. The gain in weight of the YAC128 animals over time compared with most HD mouse models has been suggested to be the result of the third copy of HTT and higher huntingtin expression [46]. There is, however, no overlap with the genes found to be altered in Htt null cell lines [47], although this could be the result of the very different biological systems studied and the different chips used.

Htr2a is expressed more highly in the YAC128 than the HdhQ150 striata. HTR2A in humans is the major serotonin receptor in the brain and is a target of the SSRI citalopram, which downregulates its expression [48–50]. The potentially depressive-like symptoms noted in the YAC128 animals [26] might therefore relate to alterations in the expression of this receptor. The deficits in the forced swim test in the YAC128 mice are thought to be a surrogate test for a syndrome related to depression in mice and this alteration in gene expression might underlie this observation. The SSRI fluoxetine did not improve performance in the forced swim task in YAC128 mice, but this test is confounded in these mice by the motor deficit and it is not clear that fluoxetine has the same effects on Htr2a as citalopram. In addition, the improvement in symptoms seen in mice treated with SSRIs is thought to be mediated by BDNF levels: it is possible that these are so compromised in the YAC128 mice that the drugs are incapable of improving them although levels of Bdnf RNA measured in the striata are unchanged in the cohort of animals that we used. However, given the prevalence of depressive symptoms in HD patients [26, 51–53] this may shed light on their mechanism.

The differences in gene expression with age in these animals is not the same as for the HdhQ150 animals. There are fewer changes and they do not highlight as many pathways, or the same pathways, as in the HdhQ150 mouse striata, especially between 6 and 12 months [35]. Between 12 and 18 months more pathways are highlighted some of which relate to perception of smell and neurological processes and notably genes for the perception of smell were enriched between 6 and 12 months in the HdhQ150 cohorts. It is not clear why these differences exist, and it must be remembered that these are not true longitudinal data (the same measurements in the same animals) but rather pseudo-longitudinal: the animals are maintained in the same way and are on the same genetic background but are different animals at the different time points. There may be subtle differences in the housing over time that have contributed to these differences in striatal gene expression profiles with age.

Examining the effects of genotype with age using TANOVA shows that genes that are down-regulated in the WT animal striata with age are also down-regulated in the Q150 knock in homozygote striata. This may indicate that the processes underlying the down-regulation of gene expression in HD model striata are related to those in normal ageing. The genes highlight a number of differentially regulated pathways including those involved in G-protein mediated signal transduction and homeostasis, although there is no direct overlap with the enriched pathways seen in the Q150 striata, despite the strong concordant overlap in the expression of individual genes [39].

Conclusions

The related pathways of cell adhesion, neuronal projections, synaptic functions and transmission of nerve impulse seen in the DAVID analysis for enriched pathways indicates that the maintainance and regulation of connections between neurons is central to the molecular pathogenesis in YAC128 striatum.

The similarities in striatal gene expression differences between this HD model and human caudate and with other mouse models indicate that similar molecular processes are probably occurring. These appear to occur despite differences in inclusion formation. This indicates that the presence of HTT-positive inclusions is not necessary for these changes to occur. This separation of inclusion formation from the molecular consequences of the Htt mutation is worthy of further investigation as many studies in in vitro systems use inclusion formation as a measurable end-point to assess the effects of potential therapeutics.

Methods

Samples

Heterozygous YAC128 mice on a C57BL6/J background [13] were bred in house and genotypes ascertained using tail tip DNA (Laragen Inc., Los Angeles). Mice were culled by cervical dislocation at the same point in the light phase of the diurnal cycle. CAG repeat lengths in the YAC128 animals were 121 and no variation from this was detected. WT and hemizygous YAC128 animals of both sexes were used in the experiments (57 % male). The animals were housed as sex matched littermate groups and had access to food and water ad libitum. All experiments were carried out in accordance with the United Kingdom Animals (Scientific Procedures) Act of 1986, and subject to local ethical review (Project licence PPL30/1968 and PPL30/2305). The behavioural data relating to the complete cohort of mice are given in Brooks et al. [27].

Gene expression

From this experimental group, 15 hemizygous YAC128 (7 female and 8 male) and 14 WT (Hdh+/+) mice (6 female and 9 male) were used for gene expression studies. Brains from age matched mice from each genotype were harvested at 6, 12 and 18 months and micro-dissected into striatum, motor cortex, cerebellum, prefrontal cortex and hippocampus. The dissected brain samples were snap frozen in liquid nitrogen and stored at −80 °C.

Total RNA was extracted from micro-dissected striata for gene expression analysis as previously described [19, 35]. RNA quality was determined using an Agilent RNA 6000 Nano Kit and Agilent 2100 Bioanalyser (Agilent Technologies, Santa Clara, USA). Samples with RIN (RNA integrity number) values greater than 7.5 were selected for subsequent analysis. For each RNA sample, cDNA was generated from 100 ng total RNA using an Ambion® WT expression kit (Applied Biosystems Carlsbad, California, USA), followed by fragmentation, labelling and hybridisation to a Mouse GeneChip Gene 1.0 ST Array. An Affymetrix WT Terminal Labelling and Hybridisation kit was used according to the manufacturer’s protocol. Gene Chips were processed using a Fluidics station 450 and a GeneChip scanner 3000 7G (Affymetrix UK Ltd, High Wycombe UK).

Statistical methods

Gene expression analysis

An analysis of GeneChip expression data was undertaken using R/Bioconductor. Expression values were computed using robust multichip average (RMA) (affy package [54]), with testing for differential gene expression by age or genotype performed using moderated t-tests in LIMMA [55]. Changes in gene response over time were identified using TANOVA [56]. Genes with a false discovery rate (FDR) [57] corrected p < 0.05 were extracted and the data for these genes classified in three patterns representing an up, no change or down difference in expression over the time course (this was done separately for the WT and YAC128 animals tagging each gene with its highest correlation to theoretical expression profiles). The gene expression data are available through GEO accession number GSE70656.

Determining biological themes

The resultant gene lists from the differential gene expression, time course ANOVA and behaviour/expression correlation were analysed for over-representation of genes in pathways against GO Biological Process gene sets using the Bioconductor GOstats package with the conditional hypergeometic test (which only uses those terms that were not already significant when testing a higher order (parent) term). Changes in expression of genes in GO gene sets were assessed using Gene Set Analysis [58] against the whole dataset. The biological themes were further analysed using the DAVID database [31, 32] (with the appropriate background gene list selected for the make of microarray at the website). An input list was chosen for all probesets significantly differentially expressed at a nominal p-value < 0.01. The output from this and visualised on Cytoscape (version 2.8.3) [59, 60], via the Enrichment Map plugin.

Comparison with other gene expression data

Comparison with differentially expressed genes from human HD brain [18] and other HD models [19, 33] was calculated using hypergeometric tests on the top 1000 ranked genes in gene lists for differently expressed genes between WT and YAC128. To enable comparisons between different array platforms where the species was identical, probesets in gene lists were first converted to unique EntrezGene IDs and the overlap calculated using these. Where overlaps were made between data for different species, data was first converted to Entrez Gene IDs and then to Homologene IDs which were used to calculate the overlap between lists. In addition, a graphical representation of the overlap, along with information about the relative direction of changes was generated using the method of Kuhn [17]. Genes differentially expressed between YAC128 and HdhQ150 striata were identified by a significant interaction term in the ANOVA.

Histology

Tissue was processed and immunohistochemistry carried out as previously [61]. The ubiquitin antibody was Mouse anti-ubiquitin from Invitrogen used at a dilution of 1/1000 in TBS which stains huntingtin positive inclusions [62]. S830, which recognises mutant HTT was a kind gift from Gillian Bates [63].

References

  1. The Huntington's disease collaborative research group. A novel gene containing a trinucleotide repeat that is expanded and unstable on Huntington's disease chromosomes. Cell 1993, 72:971–983.

  2. Crook ZR, Housman D. Huntington’s disease: can mice lead the way to treatment? Neuron. 2011;69:423–35.

    Article  PubMed  CAS  Google Scholar 

  3. Gray M, Shirasaki DI, Cepeda C, Andre VM, Wilburn B, Lu XH, et al. Full-length human mutant huntingtin with a stable polyglutamine repeat can elicit progressive and selective neuropathogenesis in BACHD mice. J Neurosci. 2008;28:6182–95.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  4. Lin CH, Tallaksen-Greene S, Chien WM, Cearley JA, Jackson WS, Crouse AB, et al. Neurological abnormalities in a knock-in mouse model of Huntington's disease. Hum Mol Genet. 2001;10:137–44.

    Article  PubMed  CAS  Google Scholar 

  5. Mangiarini L, Sathasivam K, Seller M, Cozens B, Harper A, Hetherington C, et al. 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.

    Article  PubMed  CAS  Google Scholar 

  6. Menalled LB, Sison JD, Wu Y, Olivieri M, Li XJ, Li H, et al. Early motor dysfunction and striosomal distribution of huntingtin microaggregates in Huntington's disease knock-in mice. J Neurosci. 2002;22:8266–76.

    PubMed  CAS  Google Scholar 

  7. Wheeler VC, Auerbach W, White JK, Srinidhi J, Auerbach A, Ryan A, et al. Length-dependent gametic CAG repeat instability in the Huntington's disease knock-in mouse. Hum Mol Genet. 1999;8:115–22.

    Article  PubMed  CAS  Google Scholar 

  8. Bowles KR, Brooks SP, Dunnett SB, Jones L. Gene Expression and Behaviour in Mouse Models of Huntington's Disease. Brain Res Bull. 2012;88:276–84.

    Article  PubMed  CAS  Google Scholar 

  9. Trueman RC, Brooks SP, Jones L, Dunnett SB. The operant serial implicit learning task reveals early onset motor learning deficits in the Hdh knock-in mouse model of Huntington's disease. Eur J Neurosci. 2007;25:551–8.

    Article  PubMed  Google Scholar 

  10. Brooks SP, Betteridge H, Trueman RC, Jones L, Dunnett SB. Selective extra-dimensional set shifting deficit in a knock-in mouse model of Huntington's disease. Brain Res Bull. 2006;69:452–7.

    Article  PubMed  CAS  Google Scholar 

  11. Carter RJ, Lione LA, Humby T, Mangiarini L, Mahal A, Bates GP, et al. Characterization of progressive motor deficits in mice transgenic for the human Huntington's disease mutation. J Neurosci. 1999;19:3248–57.

    PubMed  CAS  Google Scholar 

  12. Lione LA, Carter RJ, Hunt MJ, Bates GP, Morton AJ, Dunnett SB. Selective discrimination learning impairments in mice expressing the human Huntington's disease mutation. J Neurosci. 1999;19:10428–37.

    PubMed  CAS  Google Scholar 

  13. Slow EJ, van Raamsdonk J, Rogers D, Coleman SH, Graham RK, Deng Y, et al. Selective striatal neuronal loss in a YAC128 mouse model of Huntington disease. Hum Mol Genet. 2003;12:1555–67.

    Article  PubMed  CAS  Google Scholar 

  14. Van Raamsdonk JM, Pearson J, Slow EJ, Hossain SM, Leavitt BR, Hayden MR. Cognitive dysfunction precedes neuropathology and motor abnormalities in the YAC128 mouse model of Huntington's disease. J Neurosci. 2005;25:4169–80.

    Article  PubMed  Google Scholar 

  15. Brooks S, Higgs G, Jones L, Dunnett SB. Longitudinal analysis of the behavioural phenotype in Hdh(Q92) Huntington's disease knock-in mice. Brain Res Bull. 2012;88:148–55.

    Article  PubMed  CAS  Google Scholar 

  16. Brooks S, Higgs G, Jones L, Dunnett SB. Longitudinal analysis of the behavioural phenotype in Hdh((CAG)150) Huntington's disease knock-in mice. Brain Res Bull. 2012;88:182–8.

    Article  PubMed  Google Scholar 

  17. Kuhn A, Goldstein DR, Hodges A, Strand AD, Sengstag T, Kooperberg C, et al. Mutant huntingtin's effects on striatal gene expression in mice recapitulate changes observed in human Huntington's disease brain and do not differ with mutant huntingtin length or wild-type huntingtin dosage. Hum Mol Genet. 2007;16:1845–61.

    Article  PubMed  CAS  Google Scholar 

  18. Hodges A, Strand AD, Aragaki AK, Kuhn A, Sengstag T, Hughes G, et al. Regional and cellular gene expression changes in human Huntington's disease brain. Hum Mol Genet. 2006;15:965–77.

    Article  PubMed  CAS  Google Scholar 

  19. Morton AJ, Hunt MJ, Hodges AK, Lewis PD, Redfern AJ, Dunnett SB, et al. A combination drug therapy improves cognition and reverses gene expression changes in a mouse model of Huntington's disease. European J Neuroscience. 2005;21:855–70.

    Article  Google Scholar 

  20. Hockly E, Richon VM, Woodman B, Smith DL, Zhou X, Rosa E, et al. Suberoylanilide hydroxamic acid, a histone deacetylase inhibitor, ameliorates motor deficits in a mouse model of Huntington's disease. Proc Natl Acad Sci U S A. 2003;100:2041–6.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  21. Ferrante RJ, Kubilus JK, Lee J, Ryu H, Beesen A, Zucker B, et al. Histone deacetylase inhibition by sodium butyrate chemotherapy ameliorates the neurodegenerative phenotype in Huntington's disease mice. J Neurosci. 2003;23:9418–27.

    PubMed  CAS  Google Scholar 

  22. Gardian G, Browne SE, Choi DK, Klivenyi P, Gregorio J, Kubilus JK, et al. Neuroprotective effects of phenylbutyrate in the N171-82Q transgenic mouse model of Huntington's disease. J Biol Chem. 2005;280:556–63.

    Article  PubMed  CAS  Google Scholar 

  23. Thomas EA, Coppola G, Desplats PA, Tang B, Soragni E, Burnett R, et al. The HDAC inhibitor 4b ameliorates the disease phenotype and transcriptional abnormalities in Huntington's disease transgenic mice. Proc Natl Acad Sci U S A. 2008;105:15564–9.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  24. Brooks SP, Janghra N, Higgs GV, Bayram-Weston Z, Heuer A, Jones L, et al. Selective cognitive impairment in the YAC128 Huntington's disease mouse. Brain Res Bull. 2012;88:121–9.

    Article  PubMed  Google Scholar 

  25. Brooks SP, Jones L, Dunnett SB. Longitudinal analyses of operant performance on the serial implicit learning task (SILT) in the YAC128 Huntington's disease mouse line. Brain Res Bull. 2012;88:130–6.

    Article  PubMed  Google Scholar 

  26. Pouladi MA, Graham RK, Karasinska JM, Xie Y, Santos RD, Petersen A, et al. Prevention of depressive behaviour in the YAC128 mouse model of Huntington disease by mutation at residue 586 of huntingtin. Brain. 2009;132:919–32.

    Article  PubMed  Google Scholar 

  27. Brooks S, Higgs G, Janghra N, Jones L, Dunnett SB: Longitudinal analysis of the behavioural phenotype in YAC128 (C57BL/6 J) Huntington’s disease transgenic mice. Brain Res Bull in press.

  28. Brooks SP, Janghra N, Higgs GV, Bayram-Weston Z, Heuer A, Jones L, Dunnett SB: Selective cognitive impairment in the YAC128 Huntington's disease mouse. Brain Res Bull in press.

  29. Brooks SP, Jones L, Dunnett SB: Longitudinal analyses of operant performance on the serial implicit learning task (SILT) in the YAC128 Huntington's disease mouse line. Brain Res Bull in press.

  30. Bayram-Weston Z, Jones L, Dunnett SB, Brooks SP: Light and electron microscopic characterization of the evolution of cellular pathology in YAC128 Huntington's disease transgenic mice. Brain Res Bull in press.

  31. da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.

    Article  CAS  Google Scholar 

  32. da Huang W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37:1–13.

    Article  PubMed Central  Google Scholar 

  33. Hodges A, Hughes G, Brooks S, Elliston L, Holmans P, Dunnett SB, et al. Brain gene expression correlates with changes in behavior in the R6/1 mouse model of Huntington's disease. Genes Brain Behav. 2008;7:288–99.

    Article  PubMed  CAS  Google Scholar 

  34. Strand AD, Baquet ZC, Aragaki AK, Holmans P, Yang L, Cleren C, et al. Expression profiling of Huntington's disease models suggests that brain-derived neurotrophic factor depletion plays a major role in striatal degeneration. J Neurosci. 2007;27:11758–68.

    Article  PubMed  CAS  Google Scholar 

  35. Giles P, Elliston L, Higgs GV, Brooks SP, Dunnett SB, Jones L. Gene expression changes over time and correlated with behaviour in the HdhQ150 mouse model of HD. Brain Res Bull. 2012;88:199–209.

    Article  PubMed  CAS  Google Scholar 

  36. Brooks S, Higgs G, Jones L, Dunnett SB: Longitudinal analysis of the behavioural phenotype in Hdh((CAG)150) Huntington's disease knock-in mice. Brain Res Bull in press.

  37. Bayram-Weston Z, Torres EM, Jones L, Dunnett SB, Brooks SP: Light and electron microscopic characterization of the evolution of cellular pathology in the Hdh((CAG)150) Huntington's disease knock-in mouse. Brain Res Bull in press.

  38. Gong B, Kielar C, Morton AJ. Temporal separation of aggregation and ubiquitination during early inclusion formation in transgenic mice carrying the Huntington's disease mutation. PLoS One. 2012;7, e41450.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  39. Giles P, Elliston L, Higgs GV, Brooks SP, Dunnett SB, Jones L: Gene expression changes over time and correlated with behaviour in the HdhQ150 mouse model of HD. Brain Res Bull in press.

  40. Collinge J, Whitfield J, McKintosh E, Beck J, Mead S, Thomas DJ, et al. Kuru in the 21st century--an acquired human prion disease with very long incubation periods. Lancet. 2006;367:2068–74.

    Article  PubMed  Google Scholar 

  41. Mead S, Whitfield J, Poulter M, Shah P, Uphill J, Beck J, et al. Genetic susceptibility, evolution and the kuru epidemic. Philos Trans R Soc Lond B Biol Sci. 2008;363:3741–6.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Miller J, Arrasate M, Shaby BA, Mitra S, Masliah E, Finkbeiner S. Quantitative relationships between huntingtin levels, polyglutamine length, inclusion body formation, and neuronal death provide novel insight into huntington's disease molecular pathogenesis. J Neurosci. 2010;30:10541–50.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  43. Arrasate M, Mitra S, Schweitzer ES, Segal MR, Finkbeiner S. Inclusion body formation reduces levels of mutant huntingtin and the risk of neuronal death. Nature. 2004;431:805–10.

    Article  PubMed  CAS  Google Scholar 

  44. White JK, Auerbach W, Duyao MP, Vonsattel JP, Gusella JF, Joyner AL, et al. Huntingtin is required for neurogenesis and is not impaired by the Huntington's disease CAG expansion. Nat Genet. 1997;17:404–10.

    Article  PubMed  CAS  Google Scholar 

  45. Graham RK, Slow EJ, Deng Y, Bissada N, Lu G, Pearson J, et al. Levels of mutant huntingtin influence the phenotypic severity of Huntington disease in YAC128 mouse models. Neurobiol Dis. 2006;21:444–55.

    Article  PubMed  CAS  Google Scholar 

  46. Van Raamsdonk JM, Gibson WT, Pearson J, Murphy Z, Lu G, Leavitt BR, et al. Body weight is modulated by levels of full-length huntingtin. Hum Mol Genet. 2006;15:1513–23.

    Article  PubMed  Google Scholar 

  47. Jacobsen JC, Gregory GC, Woda JM, Thompson MN, Coser KR, Murthy V, et al. HD CAG-correlated gene expression changes support a simple dominant gain of function. Hum Mol Genet. 2011;20:2846–60.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  48. McMahon FJ, Buervenich S, Charney D, Lipsky R, Rush AJ, Wilson AF, et al. Variation in the gene encoding the serotonin 2A receptor is associated with outcome of antidepressant treatment. Am J Hum Genet. 2006;78:804–14.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  49. Peremans K, Audenaert K, Hoybergs Y, Otte A, Goethals I, Gielen I, et al. The effect of citalopram hydrobromide on 5-HT2A receptors in the impulsive-aggressive dog, as measured with 123I-5-I-R91150 SPECT. Eur J Nucl Med Mol Imaging. 2005;32:708–16.

    Article  PubMed  CAS  Google Scholar 

  50. Papolos DF, Yu YM, Rosenbaum E, Lachman HM. Modulation of learned helplessness by 5-hydroxytryptamine2A receptor antisense oligodeoxynucleotides. Psychiatry Res. 1996;63:197–203.

    Article  PubMed  CAS  Google Scholar 

  51. Duff K, Paulsen JS, Beglinger LJ, Langbehn DR, Stout JC. Psychiatric symptoms in Huntington's disease before diagnosis: the predict-HD study. Biol Psychiatry. 2007;62:1341–6.

    Article  PubMed  Google Scholar 

  52. Vaccarino AL, Sills T, Anderson KE, Bachoud-Levi AC, Borowsky B, Craufurd D, et al. Assessment of depression, anxiety and apathy in prodromal and early huntington disease. PLoS Curr. 2011;3, RRN1242.

    PubMed  PubMed Central  Google Scholar 

  53. Gargiulo M, Lejeune S, Tanguy ML, Lahlou-Laforet K, Faudet A, Cohen D, et al. Long-term outcome of presymptomatic testing in Huntington disease. Eur J Hum Genet. 2009;17:165–71.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Gautier L, Cope L, Bolstad BM, Irizarry RA. affy--analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004;20:307–15.

    Article  PubMed  CAS  Google Scholar 

  55. Smyth GK. Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Statistical Applications in Genetics and Molecular Biology. 2004;Vol. 3(Iss. 1):Article 3.

    Google Scholar 

  56. Zhou B, Xu W, Herndon D, Tompkins R, Davis R, Xiao W, et al. Analysis of factorial time-course microarrays with application to a clinical study of burn injury. Proc Natl Acad Sci U S A. 2010;107:9923–8.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  57. Benjamini Y, Hochberg Y. Controlling the false discovery rate:A practical and powerful appraoch to multiple testing. J R Statist Soc B. 1995;57:289–300.

    Google Scholar 

  58. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005;102:15545–50.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  59. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  60. Merico D, Isserlin R, Stueker O, Emili A, Bader GD. Enrichment map: a network-based method for gene-set enrichment visualization and interpretation. PLoS One. 2010;5, e13984.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Bayram-Weston Z, Jones L, Dunnett SB, Brooks SP. Light and electron microscopic characterization of the evolution of cellular pathology in YAC128 Huntington's disease transgenic mice. Brain Res Bull. 2012;88:137–47.

    Article  PubMed  Google Scholar 

  62. Holbert S, Dedeoglu A, Humbert S, Saudou F, Ferrante RJ, Neri C. Cdc42-interacting protein 4 binds to huntingtin: neuropathologic and biological evidence for a role in Huntington's disease. Proc Natl Acad Sci U S A. 2003;100:2712–7.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  63. Landles C, Sathasivam K, Weiss A, Woodman B, Moffitt H, Finkbeiner S, et al. Proteolysis of mutant huntingtin produces an exon 1 fragment that accumulates as an aggregated protein in neuronal nuclei in Huntington disease. J Biol Chem. 2010;285:8808–23.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

Download references

Acknowledgements

This work was funded by CHDI and supported by an MRC studentship to ZB-W. We thank Central Biotechnology Services at the School of Medicine, Cardiff University, for technical support.

Author information

Author notes

    Authors and Affiliations

    1. MRC Centre for Neuropsychiatric Genetics and Genomics, School of Medicine, Cardiff University, Cardiff, CF24 4HQ, UK

      Timothy C. Stone, Peter Giles, Linda Elliston, Peter A. Holmans & Lesley Jones

    2. Brain Research Group, School of Bioscience, Cardiff University, Cardiff, CF10 4AX, UK

      Zubeyde Bayram-Weston, Nari Janghra, Gemma V. Higgs, Stephen B. Dunnett & Simon P. Brooks

    Authors
    1. Zubeyde Bayram-Weston
      View author publications

      You can also search for this author in PubMed Google Scholar

    2. Timothy C. Stone
      View author publications

      You can also search for this author in PubMed Google Scholar

    3. Peter Giles
      View author publications

      You can also search for this author in PubMed Google Scholar

    4. Linda Elliston
      View author publications

      You can also search for this author in PubMed Google Scholar

    5. Nari Janghra
      View author publications

      You can also search for this author in PubMed Google Scholar

    6. Gemma V. Higgs
      View author publications

      You can also search for this author in PubMed Google Scholar

    7. Peter A. Holmans
      View author publications

      You can also search for this author in PubMed Google Scholar

    8. Stephen B. Dunnett
      View author publications

      You can also search for this author in PubMed Google Scholar

    9. Simon P. Brooks
      View author publications

      You can also search for this author in PubMed Google Scholar

    10. Lesley Jones
      View author publications

      You can also search for this author in PubMed Google Scholar

    Corresponding author

    Correspondence to Lesley Jones.

    Additional information

    Competing interests

    The authors have no competing interests.

    Authors’ contributions

    LJ and SBD designed the study, SPB, GVH and NJ conducted the animal work and tissue dissections, ZB-W carried out the pathology, LE processed the samples, TCS, PG, PAH and LJ conducted the data analysis. LJ wrote and all authors edited the MS. All authors read and approved the final manuscript.

    Zubeyde Bayram-Weston, Timothy C. Stone and Peter Giles contributed equally to this work.

    Additional files

    Additional file 1: Table S1.

    Genes altered in expression with age common to both WT and YAC128 caudate. (XLS 75 kb)

    Additional file 2: Table S2.

    Over-representation analysis of the gene expression differences with age in YAC128 and WT striata (XLS 116 kb)

    Additional file 3: Table S3.

    Genes differentially expressed between WT and YAC128 striata. (XLS 393 kb)

    Additional file 4: Figure S1.

    TANOVA patterns of genes differentially expressed between WT and YAC128 caudate. Plots are sorted by pattern of changed expression and each of the 9 patterns is indicated at the end of the relevant section. The gene ID and the Affymetrix probeset ID are given. Expression is given as Log-2 fold change and the time points are in months. WT = wild-type and HD = YAC128. (PDF 65 kb)

    Additional file 5: Table S4.

    Over-representation analysis of the 2557 genes nominally significantly differentially expressed between YAC128 and WT striata (XLS 61 kb)

    Additional file 6: Figure S2.

    Analysis of the overlap of changes in gene expression in the YAC128 striatum at A. 6 months of age B. 12 months of age and C. 18 months of age compared with other HD model mouse striata and human brain. Frequency represents the fraction of the top 200 HdhQ150 expression changes that map to a particular bin of ranked data (1000 genes per bin) in the other dataset, which is then split to identify concordant or discordant direction of expression change. A higher frequency of concordant (green) rather than discordant (red) in the first bins indicates a similarity between the YAC128 caudate and other model or human HD gene expression signature. (PDF 59 kb)

    Additional file 7: Table S5.

    The 869 probesets that are differentially expressed between YAC128 and HdhQ150 striata. (XLS 137 kb)

    Rights and permissions

    Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

    Reprints and permissions

    About this article

    Check for updates. Verify currency and authenticity via CrossMark

    Cite this article

    Bayram-Weston, Z., Stone, T.C., Giles, P. et al. Similar striatal gene expression profiles in the striatum of the YAC128 and HdhQ150 mouse models of Huntington’s disease are not reflected in mutant Huntingtin inclusion prevalence. BMC Genomics 16, 1079 (2015). https://doi.org/10.1186/s12864-015-2251-4

    Download citation

    • Received:

    • Accepted:

    • Published:

    • DOI: https://doi.org/10.1186/s12864-015-2251-4

    Keywords

    BMC Genomics

    ISSN: 1471-2164

    Contact us