- Research article
- Open Access
Elucidating a normal function of huntingtin by functional and microarray analysis of huntingtin-null mouse embryonic fibroblasts
BMC Neuroscience volume 9, Article number: 38 (2008)
The polyglutamine expansion in huntingtin (Htt) protein is a cause of Huntington's disease (HD). Htt is an essential gene as deletion of the mouse Htt gene homolog (Hdh) is embryonic lethal in mice. Therefore, in addition to elucidating the mechanisms responsible for polyQ-mediated pathology, it is also important to understand the normal function of Htt protein for both basic biology and for HD.
To systematically search for a mouse Htt function, we took advantage of the Hdh +/- and Hdh-floxed mice and generated four mouse embryonic fibroblast (MEF) cells lines which contain a single copy of the Hdh gene (Hdh-HET) and four MEF lines in which the Hdh gene was deleted (Hdh-KO). The function of Htt in calcium (Ca2+) signaling was analyzed in Ca2+ imaging experiments with generated cell lines. We found that the cytoplasmic Ca2+ spikes resulting from the activation of inositol 1,4,5-trisphosphate receptor (InsP3R) and the ensuing mitochondrial Ca2+ signals were suppressed in the Hdh-KO cells when compared to Hdh-HET cells. Furthermore, in experiments with permeabilized cells we found that the InsP3-sensitivity of Ca2+ mobilization from endoplasmic reticulum was reduced in Hdh-KO cells. These results indicated that Htt plays an important role in modulating InsP3R-mediated Ca2+ signaling. To further evaluate function of Htt, we performed genome-wide transcription profiling of generated Hdh-HET and Hdh-KO cells by microarray. Our results revealed that 106 unique transcripts were downregulated by more than two-fold with p < 0.05 and 173 unique transcripts were upregulated at least two-fold with p < 0.05 in Hdh-KO cells when compared to Hdh-HET cells. The microarray results were confirmed by quantitative real-time PCR for a number of affected transcripts. Several signaling pathways affected by Hdh gene deletion were identified from annotation of the microarray results.
Functional analysis of generated Htt-null MEF cells revealed that Htt plays a direct role in Ca2+ signaling by modulating InsP3R sensitivity to InsP3. The genome-wide transcriptional profiling of Htt-null cells yielded novel and unique information about the normal function of Htt in cells, which may contribute to our understanding and treatment of HD.
Huntington's disease (HD) is an autosomal-dominant neurodegenerative disorder which is caused by polyglutamine (polyQ) expansion in the amino-terminus of huntingtin (Htt). Htt is a soluble protein of 3,144 amino acids that has no sequence homology with other proteins. Except for the extreme amino-terminus, with its adjacent polyQ region and proline-rich segments, the entire ~350-kD protein is predicted to be composed of 36 α-helical HEAT repeats. Increasing evidence indicates that Htt functions as a molecular scaffold that is able to organize a variety of signaling complexes [1, 2]. Htt is expressed ubiquitously in humans and rodents, with the highest levels found in CNS neurons and the testes [3–5]. Intracellularly, Htt is associated with various organelles, including the nucleus, endoplasmic reticulum (ER) and Golgi complex [6–8]. This widespread subcellular localization does not facilitate the definition of its function. Hdh is evolutionary conserved – a single copy of the Htt gene is expressed in all vertebrates (from fish to humans) . The Htt gene is also present in D. melanogaster genome, but absent in the C. elegans and S. cerevisiae genomes . All vertebrate isoforms of Htt, but not Drosophila Htt, contain an amino-terminal polyQ region.
Complete knockout of the mouse Htt gene (Hdh) causes embryonic death before day 8.5 (E8.5, before gastrulation and the formation of the nervous system) [10–12]. After gastrulation, Htt becomes important for neurogenesis – mice carrying a <50% dose of wild-type Htt display profound malformations of the cortex and striatum . Another study has shown that greatly reduced Htt levels are insufficient to support normal mouse development . In addition to its function in development, Htt may play a role in the regulation of apoptosis, control of BDNF production, vesicular and mitochondrial transport, neuronal gene transcription, and synaptic transmission (reviewed in ). Despite all of these efforts and results, the exact function of Htt in cells still remains largely unknown.
In addition to answering an academic question concerning the normal function of Htt, knowledge of its function is important for understanding HD pathogenesis and for the treatment of Huntington's disease (HD). Although the HD mutation is considered to be a "gain of function" mutation, it has been suggested that the loss of normal Htt function might also contribute to the pathogenesis of HD . Approaches that are based on reducing mutant Htt expression such as RNA interference  and the use of intrabodies [16, 17] are currently considered to be promising strategies for HD treatment. It is likely that these agents will cause inactivation or impair normal function of both mutant and wild type Htt alleles. One can envision a therapy that combines such Htt-inactivating agents with drugs that restore the function of targets and pathways downstream from wild-type Htt. However, because both normal Hdh function is not known and downstream pathways have not been identified, such a combined therapy approach is not feasible at the moment.
To systematically search for Htt's normal function, we used Hdh +/-  and Hdh-floxed mice  to generate immortalized mouse embryonic fibroblasts (MEF) which contain a single functional copy of Hdh gene (Hdh-HET) or lack Hdh completely (Hdh-KO). We compared inositol 1,4,5-trisphosphate receptor (InsP3R)-mediated Ca2+ signals in these cells. We then performed a genome-wide gene transcription profiling of Hdh-HET and Hdh-KO MEF cells using microarrays to obtain novel, unique, and unbiased information about the normal function of Htt in fibroblasts, which may contribute to our understanding and treatment of HD.
Generation of Hdh-HET and Hdh-KO MEF cell lines
To generate cell lines lacking Htt expression, we employed a conditional mutagenesis strategy based on the in vitro recombination of an Hdh(flox) allele in cultured fibroblasts that are also carrying either a wild-type (+) or null (standard knock-out) Hdh allele. The Hdh(flox/+) and Hdh(flox/-) fibroblasts were obtained from embryos derived from a cross between Hdh +/- and Hdh-floxed/floxed mice (Fig 1). Primary fibroblasts were prepared and plated separately from each embryo as described in Methods. After two days in culture, the primary fibroblasts from all embryos with identical genotype (Hdh floxed/+ or Hdh floxed/-) were pooled together and transfected with a linearized SV40 plasmid. Transfected cells were then cultured for four to six weeks until immortalized Hdh floxed/+ and Hdh floxed/- mouse embryonic fibroblasts (MEFs) were obtained (Fig 1). To recombine the Hdh floxed allele, immortalized Hdh floxed/+ and Hdh floxed/- MEFs were infected with Lenti-NLS-GFP-Cre virus encoding nuclear-targeted GFP-Cre fusion protein  (Fig 1). Using the procedure described above, we generated four Hdh-HET (lines 1,2,3,5) and four Hdh-KO (lines 11, 12, 16, 27) MEF cell lines. The expression of Htt in the generated MEF lines was assessed by Western blotting of whole cell lysates using anti-Htt monoclonal antibody. Quantification of Western blotting data verified similar levels of Htt expression in all 4 Hdh-HET lines (data not shown). Consistent with the genotype of the generated cells, we detected a protein of predicted size (~350 kD) in lysates from the Hdh-HET cells, but not in lysates from Hdh-KO lines (Fig 2). The same samples were probed with monoclonal antibodies against β-actin as a loading control (Fig 2). Thus, we concluded that we successfully generated four Hdh-HET and four Hdh-KO MEF lines on similar genetic background. We reasoned that comparison of resulting MEF lines may reveal clues about normal function of Htt protein in cells.
Intracellular calcium signaling in Hdh-HET and Hdh-KO MEF cell lines
Previous studies have implicated impaired calcium signaling in the pathogenesis of HD [20–22]. The Htt directly binds to the inositol 1,4,5-trisposphate receptor (InsP3R), an intracellular Ca2+ release channel [23, 24]. The expression of mutated Htt has been shown to affect the InsP3R activity  and mitochondrial Ca2+ signals and bioenergetics [25–30]. Since Htt and mutated Htt directly targets both ER and mitochondrial sites, it is possible that Htt may have some relevance for the physical and local Ca2+ coupling between ER and mitochondria [31, 32]. To test this idea and to determine a role played by Htt in intracellular Ca2+ signaling, we performed a series of cytosolic and mitochondrial Ca2+ imaging experiments with generated MEF lines. Two Hdh-HET (HET1 and HET 5) and two Hdh-KO MEF cell lines (KO12 and KO27) (Fig 2) were selected for Ca2+ imaging studies.
In these experiments Hdh-HET and Hdh-KO MEF cell lines were challenged by ATP, an agonist of InsP3 signaling pathway in fibroblasts. Before stimulation with ATP the incubation medium was switched to a Ca2+ free buffer to prevent Ca2+ entry. The cytosolic [Ca2+]c and mitochondrial [Ca2+]m levels were monitored simultaneously as described in Methods. We found that pre-stimulation [Ca2+]c was higher in the Hdh-HET cells (HET1, 183 ± 11 nM (n = 42) and HET5 139 ± 9 nM (n = 33)) than that in the Hdh-KO cells (KO12, 102 ± 12 nM (n = 15) and KO27, 106 ± 6 nM (n = 48), p < 0.003). Stimulation with a suboptimal dose of ATP (2 μM) elicited large [Ca2+]c spikes in Hdh-HET cells, whereas only a small and slow [Ca2+]c rise was evoked in the Hdh-KO cells (Fig 3, Δ[Ca2+]c: HET1, 374 ± 24 nM; HET5, 386 ± 40 nM; KO12, 266 ± 47 nM and KO27, 284 ± 25 nM; p < 0.03). Many Hdh-KO cells did not show any [Ca2+]c rise in response to 2 μM ATP (not shown). Mobilization of the residual Ca2+ by optimal ATP (100 μM) evoked a relatively large [Ca2+]c elevation in the Hdh-KO cells (Fig 3). The ATP-induced [Ca2+]c spikes were closely followed by a [Ca2+]m elevation in the Hdh-HET cells, whereas the [Ca2+]m rise was modest in the Hdh-KO cells (Fig 3). Thus, the ATP-induced intracellular Ca2+ mobilization was suppressed and desensitized in the Hdh-KO cells.
Figure 3 has also showed that the [Ca2+]c signal was effectively propagated to the mitochondria in the Hdh-HET cells but less Ca2+ was transferred into the mitochondria in the Hdh-KO cells. To determine whether this resulted from the attenuated ER Ca2+ mobilization or from impaired ER-mitochondrial Ca2+ coupling, the [Ca2+]m rise was plotted against the [Ca2+]c increase for each cell., The relationship between [Ca2+]c elevations and the ensuing [Ca2+]m signals was similar in both the Hdh-HET and in the Hdh-KO cells (not shown). This result indicated that the local coupling between ER and mitochondria is maintained in the Hdh-KO cells.
The attenuated ATP-induced [Ca2+]c signal in the Hdh-KO cells could be due to reduced InsP3 generation, due to reduced sensitivity of InsP3R or due to depleted ER Ca2+ pool. To discriminate between these possibilities the InsP3-induced Ca2+ mobilization was quantified in suspensions of permeabilized MEF cells. The steady state [Ca2+]c was similar in both Hdh-HET and in Hdh-KO cells (Fig 4A, HET1, n = 12; HET5, n = 12; KO12, n = 13; KO27, n = 10 measurements). The Ca2+ pool size for both the ER and the ionophore-sensitive compartment was larger in the Hdh-KO than in the control cells (Fig 4A), whereas the uncoupler-sensitive mitochondrial compartment showed no difference (n = 3; not shown). Sequential application of a suboptimal and optimal InsP3 revealed lesser InsP3-sensitivity in the Hdh-KO cells than in the HET cells (Fig 4B, 4C, p < 0.01). However, the InsP3 sensitive fraction of the ER was approximately 75% in both control and Hdh-HET cells (Fig 4D). Furthermore, neither the passive Ca2+ buffering nor the mitochondrial Ca2+ uptake was altered in the cells lacking the Hdh (Fig 4E, 4F). Collectively, the data obtained in permeabilized cells suggest that the InsP3 sensitivity of the InsP3 receptor is attenuated in the Hdh-deficient cells, providing a mechanism to underlie the suppression of the InsP3-linked [Ca2+]c signaling in the Hdh-KO cells. The effect of Hdh on the InsP3 sensitivity is likely to be mediated by direct association between Htt and InsP3R [23, 24].
Microarray analysis of transcripts expressed in Hdh-KO and Hdh-HET MEFs
The results described in the previous section suggested that Htt plays a direct role in Ca2+ signaling by modulating InsP3R function. Many studies suggested that Htt also plays a major role in control of gene transcription [33, 34]. To uncover potential gene expression changes we performed genome-wide transcription profiling of Hdh-HET and Hdh-KO MEF cells. Using the procedures described in Methods, we isolated total RNA from Hdh-HET MEF lines 1, 2 (in duplicate), 3, and 5 (in duplicate) and from Hdh-KO MEF lines 11, 12, 16 (in duplicate), 27 (in duplicate). The resulting 12 samples were provided to the UT Southwestern Microarray Core Facility (MCF) for genome-wide expression profiling using Sentrix Mouse-6 Expression Bead Chips (Illumina) (see Methods for details).
All 12 arrays produced highly reproducible and consistent gene expression data. The microarray results have been deposited in NCBIs Gene Expression Omnibus (GEO), and are accessible through GEO Series accession number GSE11139 . Cluster analysis of the over-all gene expression data from the 12 samples demonstrated that all six Hdh-HET samples and all six Hdh-KO samples were clustered together, forming two clearly distinct groups (Fig 5). We also found that the duplicate samples (HET5.1 and HET5.2, HET2.1 and HET2.2, KO16.1 and KO16.2, and KO 27.1 and KO27.2) were most similar to each other when compared to other samples (Fig 5), as should be expected. Thus, we concluded that we obtained a high quality dataset of transcripts expressed in Hdh-HET and Hdh-KO MEF lines.
In the next level of analysis, we combined all results obtained in six arrays with the Hdh-HET samples (HET group) and in six arrays with the Hdh-KO samples (KO group). Using Illumina BeadStudio software, we performed a statistical analysis to identify the differentially expressed genes between the Hdh-KO and Hdh-HET groups. We found that from 45992 probes existing on the Mouse-6 BeadChips arrays, 14,065 probes were present (with detection p-value < 0.01) in at least one of these two groups. Statistical analysis (t-test) has identified 821 transcripts that were significantly different between Hdh-HET and Hdh-KO groups (p < 0.05) (Fig 6). Among these 821 targets, 455 were up-regulated in the KO group and 366 were down-regulated (Fig 6). Thus, we concluded that inactivation of Hdh expression has a very significant effect on the transcriptional profile of MEF cells. The average signal intensities for each probe in Hdh-HET and Hdh-KO groups are included in the Excel format for the 821 differentially expressed targets (Additional files 1 and 2).
Analysis and annotation of microarray data
In the annotation of microarray data, we focused on 279 unique genes whose expression differed by at least 2-fold between Hdh-HET and Hdh-KO groups. From these genes 173 genes (104 of which are annotated) were up-regulated and 106 genes (65 of which were annotated) were down-regulated in Hdh-KO cells when compared to Hdh-HET cells. In order to extract biological information from these genes, we performed annotation of the array results using the Ingenuity Pathway Analysis platform  and the GoStat software . We found that several significant GO categories can be pulled from the data: developmental process (35), nervous system development and function (19), lipid metabolism (18), glucosamine metabolic process (2), transcription regulator activity (5), cellular_component:plasma membrane (16), regulation of cellular process (29), endocytosis (2), mitochondrion (3), extracellular matrix (2), cytoskeleton (5), others (33), indicating the possible disruption of several functional pathways in the absence of Htt (Additional files 1 and 2). We have not observed significant changes in gene expression levels of most proteins involved in Ca2+ signaling pathways (Additional files 1 and 2), indicating that Ca2+ signaling changes observed in our functional experiments (Figs 2, 3, 4) are likely to be due to post-translational effects, such as for example changes in InsP3R gating properties.
Confirmation of gene expression with real-time PCR
In order to confirm our microarray results, we performed quantitative real-time PCR (qPCR) analysis for several of the candidate genes. For these experiments we choose three genes that were significantly down-regulated in Hdh-KO MEFs (Hdh, Sox-2, Tcf2) and three genes that were significantly up-regulated in Hdh-KO MEFs (Cart1, Esm1, Pitx2) (Additional files 1 and 2). We observed a good correlation between averaged microarray results and qPCR data for all six genes evaluated (Table 1). Moreover, we observed a good correlation between microarray and qPCR data for results obtained with cDNA samples from individual Hdh-HET and Hdh-KO MEF cell lines (Fig 7).
Comparison with gene expression profiling of Hdh-null ES cells
When our paper was in preparation, another group independently reported genome-wide expression profiling of Hdh-null embryonic stem (ES) cells . It is of interest to compare our findings with Hdh-KO MEF cells with results obtained for Hdh-null ES cells. Strehlow at al (2007) reported that expression of 16 known transcripts was significantly affected in Hdh-null ES cells when compared to wild type ES cells. The affected transcripts were divided into several classes of interest: protein degradation (5), extracellular matrix (4), cell division (4), and patterning/development (3) . To compare these data with our results, we attempted to locate the genes highlighted by the study by Strehlow at al (2007) among the genes which expression differs at p < 0.05 between Hdh-KO and Hdh-HET MEF cells in our experiments. We determined that two of the three genes in the patterning/development category (Otx2 and Pem) do not present in MEF cell lines, Leftb/Lefty is present in MEF cell lines, but did not show a significant difference between HET and KO MEFs. All five genes in the protein degradation category are present in MEF cell lines but did not show significant difference in our experiments. We found that Adam23 is not present in MEF cell lines. B3galt6, col4a1 and clo4a2 are present in MEF cells, but did not show significant difference. We further found that one out of four genes in cell division category (Ccdc5) is also significantly affected in our experiments. Ccdc5 is up-regulated in both Hdh-null ES cells and in Hdh-KO MEF cells. Plk1 is also strongly upregulated (8.7-fold) in Hdh-null ES cells  but did not show significant difference in Hdh-KO MEF cells in our experiments. Aurkc is not present in MEF cells, while Pttg1 is present but did not show significant difference.
Strehlow at al (2007) induced neuronal differentiation of Hdh-null ES cells by application of retinoic acid and performed microarray analysis of the in vitro differentiated neurons at 6, 8, and 10 days post-differentiation . The expression of transcripts in the number of categories was different between Hdh-null and wild type in vitro differentiated neurons . Changes in only a few genes have been consistently observed in the two studies (Table 2).
Despite the importance of understanding the normal function of Htt for both basic biology and for HD, its function remains largely unknown [9, 34]. The generation of Hdh-null MEF cell lines described in our study provides a new and unbiased approach to search for novel Htt functions.
The Htt directly binds to the InsP3R, an intracellular Ca2+ release channel [23, 24]. The expression of mutated Htt has been shown to affect the InsP3R activity  and mitochondrial Ca2+ signals and bioenergetics [25–30]. Thus, in the first series of experiments we evaluated InsP3R-mediated Ca2+ cytosolic and mitochondrial Ca2+ signals in Hdh-null MEF cells. As a result of these experiments we found that InsP3R sensitivity to stimulation by InsP3 was reduced in the absence of Hdh (Figs 3 and 4). We further found that Htt appears to be dispensable for ER-mitochondrial Ca2+ coupling (data not shown). Thus the altered InsP3R-induced cytoplasmic and mitochondrial calcium signaling in the Hdh-null MEF cells may result from the lack of Hdh by itself and does not necessarily require a secondary change in gene regulation.
Interestingly, a large number of Ca2+-related genes, such as CACNA2D3 (calcium channel, voltage-dependent, alpha2/delta subunit 3), ITPR1 (inositol 1,4,5-trisphosphate receptor 1), HOMER1 (homer homolog 1), ATP2A2 (ATPase, Ca2+ transporting, cardiac muscle, slow twitch 2), DRD2 (dopamine receptor 2), PRKCB1 (protein kinase C beta 1), PDE1B (phosphodiesterase 1B, Ca2+-calmodulin dependent), ATP2B2 (ATPase, Ca2+ transporting, plasma membrane 2), CAMK2B (calcium/calmodulin-dependent protein kinase II, beta), PLCB1 (phospholipase C, beta 1), RGS4 (regulator of G-protein signaling 4), and CAMK2A (calcium/calmodulin-dependent protein kinase II alpha), have been recently reported to be consistently and significantly downregulated in a striatal region of symptomatic human HD patients and aging HD mouse models . These results are in agreement with the "Ca2+ hypothesis of HD"  and with a direct role of Htt in intracellular Ca2+ signaling supported by our experiments.
Many studies suggested that Htt also plays a major role in control of gene transcription [33, 34]. To search for changes in gene transcription resulting from deletion of Htt gene, we performed a genome-wide comparison of transcription profiles in MEF cells expressing a single copy of Hdh (Hdh-HET cells) and in MEF cells which lack Hdh expression (Hdh-KO cells). To minimize sources of variability, the Hdh-HET and Hdh-KO MEF cells were generated in parallel experiments and on identical genetic background. From our annotation analysis, we found that a large group of affected genes play a role in embryonic development (Additional file 1). This result was not unexpected because Htt is essential for embryonic development, and complete inactivation of Htt expression in knock-out mice causes early embryonic lethality [10–12]. The functions of these genes may provide additional clues about the mechanism responsible for embryonic lethality in Hdh knockout mice, for example there are some similar phenotypic manifestations between Hdh nullizygous embryos and knockout mutants lacking fibroblast growth factor receptor1(fgfr1) . Interestingly, we found fgfrl1 message is downregulated approximately 2-fold in Hdh-KO MEF cells when compared to Hdh-HET cells (Additional files 1 and 2).
After gastrulation, Htt is important for neurogenesis – mice carrying a <50% dose of wild-type Htt display profound malformations of the cortex and striatum . The neuronal inactivation of Htt during mid- to late gestation, for example, leads to neurological abnormalities and progressive degeneration (apoptotic cells in the hippocampus, cortex and striatum, and a lack of axons) . Our analysis revealed a number of genes involved in nervous system development and function which were affected in Hdh-KO MEF cell lines (Additional file 1), and this list should also provide useful information to guide further studies of Htt's normal function in the nervous system. For example, Sox-2 expression was absent from Hdh-KO MEFs. Sox (Sry-related HMG box) genes encode transcription factors regulating crucial developmental decisions in different systems. Sox2 is expressed in, and is essential for, totipotent inner cell mass stem cells and other early multipotent cell lineages, and its ablation causes early embryonic lethality . In many different species, Sox2 is a marker of the nervous system from the beginning of its development, it maintains a stem-cell like state and actively inhibits neuronal differentiation, Sox2 deficiency causes neurodegeneration and impaired neurogenesis in the adult mouse brain. Does the absence of Sox-2 play some role in the early embryonic lethality and neurodegeneration in Hdh knock-out mice and conditional knock-out mice respectively? Further studies are required to answer these questions. Sox-11 is another sox family gene which changes dramatically in Hdh-KO MEF cell lines (increases about 10-fold, see Additional files 1 and 2). The widespread expression of sox-11 in both the central and peripheral nervous system suggests that sox-11 plays a general role in neuronal development, and its changes in Hdh-KO cells merit further investigation.
Another group of genes whose expression was significantly affected in Hdh-KO MEF cells are the genes related to lipid metabolism. It has been reported in another microarray analysis using clonal striata-derived cells, that genes involved in lipid metabolism were affected after expressing different N-terminal 548-amino-acid Htt fragments . Moreover, recent biochemical data indicated that Htt binds to caveolin and plays a direct role in cholesterol metabolism . All these data suggested that Htt plays an important role in lipid metabolism, which may be affected by HD mutation. Indeed, RXRG (Retinoic acid receptor RXR-gamma) and RBP4 (retinol binding protein 4) are consistently downregulated in a striatal region of symptomatic human HD patients and aging HD mouse models .
From our analysis we also found calcium channel voltage-dependent alpha2/delta subunit 1 (Cacna2d1) was down-regulated about 2 fold (Additional files 1 and 2), interestingly the same CACNA2D1 protein has been recently identified as novel Htt-binding partner in unbiased mass-spectroscopy screen . A closely related alpha2/delta subunit 3 (Cacna2d3) was reported on the 3rd place on the list of the genes which are consistently and significantly downregulated in a striatal region of symptomatic human HD patients and aging HD mouse models . As discussed above, these results indicate that Htt may play a role in regulation of Ca2+ channel activity and Ca2+ signaling in cells, consistent with Ca2+ hypothesis of HD [22, 43].
In conclusion, we generated four Hdh-HET and four Hdh-KO MEF cell lines and performed functional analysis of these cells by Ca2+ imaging methods and genome-wide transcription profiling of these cell lines using a microarray approach. Our results indicated that Htt plays a direct role in intracellular Ca2+ signaling by directly modulating InsP3R function in cells. The results of microarray analysis provided a novel and unique information resource for exploring normal function of Htt in cells. The microarray results have been deposited in NCBIs Gene Expression Omnibus (GEO), and are accessible through GEO Series accession number GSE11139 . The Hdh-KO cell lines will also serve as a useful tool for future follow-up experiments aimed at elucidating Htt functions in vivo.
Generation of Hdh-HET and Hdh-KO MEF cell lines
Generation and characterization of the Hdh +/- and Hdh-floxed mice have been described previously [12, 18]. E13.5 embryos obtained from a cross between the Hdh-floxed and Hdh +/- mice were first eviscerated and decapitated, and then the carcasses were finely minced using scissors. The tissue obtained from each embryo was digested with 0.25% Trypsin-EDTA at 37°C for 10 min, washed once with 10% FBS in DMEM, and the cell suspensions from each embryo were plated separately in 10% FBS-DMEM culture medium in order to obtain cultures of primary fibroblasts. Following plating of the cells, the genotype of each embryo was determined by PCR [12, 18]. After two days in culture, the primary fibroblasts from each genotype (Hdh floxed/+ and Hdh floxed/-) were pooled together. The pooled cells were then plated on six-well tissue culture plates, grown to 60–80% confluence and transfected with SV40-Large T-antigen plasmid in pcDNA3-Zeo vector (linearized with PvuI) using the Fugene-6 lipofection reagent (Roche). Transfected cells were cultured for four-six weeks until immortalized Hdh floxed/+ and Hdh floxed/- mouse embryonic fibroblasts (MEFs) were obtained. Hdh floxed/+ and Hdh floxed/- immortalized MEFs were then infected with Lenti-NLS-GFP-Cre virus  encoding a nuclear-targeted GFP-Cre fusion protein. Expression of NLS-GFP-Cre in infected cells catalyzed recombination at the Hdh(flox) loxP sites leading to excision of the promoter and exon 1 sequences of Hdh(flox) allele . Following infection with Lenti-NLS-GFP-Cre, the Hdh floxed/+ and Hdh floxed/- MEFs were plated in 10% FBS-DMEM culture medium in 96 well plates at an average density of one cell/well for clonal selection. Successful recombination with NLS-GFP-Cre converts the Hdh floxed/+ MEFs to Hdh +/- (Hdh-HET) MEFs and the Hdh floxed/- MEFs to Hdh-/-(Hdh-KO) MEFs. After four weeks in culture, several potential Hdh-HET and Hdh-KO MEF lines were selected, expanded and analyzed by Western blotting with monoclonal antibodies against Hdh (Chemicon MAB2166, 1:500) and monoclonal antibodies against actin (Chemicon MAB1501, 1:2000).
Cytosolic and mitochondrial Ca2+imaging
For measurements of mitochondrial matrix [Ca2+] ([Ca2+]m), the cells were transfected with a mitochondrial matrix targeted inverse pericam construct  by electroporation 48–72 h prior to the imaging experiment. Before use, the cells were preincubated in an extracellular medium (2% BSA/ECM) consisting of 121 mM NaCl, 5 mM NaHCO3, 10 mM Na-HEPES, 4.7 mM KCl, 1.2 mM KH2PO4, 1.2 mM MgSO4, 2 mM CaCl2, 10 mM glucose and 2% bovine serum albumin (BSA), pH 7.4. To monitor [Ca2+]c cells were loaded with 2.5 μM Fura2/AM for 20–30 min in the presence of 200 μM sulfinpyrazone and 0.003% (w/v) pluronic acid at room temperature. Before start of the measurement the buffer was replaced by a Ca2+ free 0.25%BSA/ECM ([Ca2+] <1 μM).
Coverslips were mounted on the thermo stated stage (35°C) of an Olympus IX70 inverted microscope fitted with a 40× (UApo, NA 1.35) oil immersion objective. Fluorescence images were collected using a cooled CCD camera (Pluto, Pixelvision). Excitation was rapidly switched among 340 and 380 nm for fura2 and 495 nm for pericam and a 510 nm longpass dichroic mirror and a 520 nm longpass emission filter were used.
For evaluation of cytoplasmic [Ca2+] ([Ca2+]c), Fura2 fluorescence was calculated for the total area of individual cells. [Ca2+]c was calibrated in terms of nM using in vitro dye calibration. For evaluation of [Ca2+]m the pericam-mt signal was masked. Recordings obtained from all transfected cells on the field (2–10 cells) were averaged for comparison in each experiment. Experiments were carried out with at least four different cell preparations. Significance of differences from the relevant controls was calculated by ANOVA.
Measurements of cytosolic Ca2+in permeabilized cells
The cells (2 mg protein/1.5 ml) were permeabilized in an intracellular medium (KCl 120 mM, NaCl 10 mM, KH2PO4 1 mM, Tris-HEPES 20 mM, and antipain, leupeptin and pepstatin 1 μg/ml each at pH 7.2) supplemented with 40 μg/ml digitonin and with fura2/FA (1.5–3 μM) or furaFF/FA (0.5 μM) and TMRE 2 μM in a stirred thermo stated cuvette at 35°C. All the measurements were carried out in the presence of succinate 2 mM, 2 mM MgATP and ATP regenerating system composed of 5 mM phosphocreatine, 5 U/ml creatine kinase. Fura2/FA or FuraFF/FA was added to monitor [Ca2+] in the intracellular medium that exchanges readily with the cytosolic space and so [Ca2+]fura2 was abbreviated as [Ca2+]c. Fluorescence was monitored in a fluorometer (Delta-RAM, PTI) using 340 nm, 380 nm excitation and 500 nm emission for fura2FF and 540 nm excitation and 580 nm emission for TMRE. Calibration of the fura fluorescence was carried out at the end of each measurement as described previously . Experiments were with at least five different cell preparations in multiple parallels. Significance of differences from the relevant controls was calculated by ANOVA.
Total RNA was isolated from fibroblast cultures using the TRIZOL reagent according to manufacturer's instructions (Invitrogen). Briefly, the MEF cells were grown to 60–80% confluence in T25 tissue culture flask, the culture medium was aspired and 1 ml of TRIZOL reagent was added to each flask. The cells were incubated with TRIZOL at room temperature for 5 min. The resulting lysates were collected from each flask, mixed with 0.2 ml of chloroform and centrifuged at 12,000 × g for 15 min at 4°C. The supernatants were collected, mixed with an equal volume of 70% ethanol at room temperature and immediately transferred to RNAeasy mini spin columns for RNA purification according to the manufacturer's (Qiagen) instructions. The final RNA samples were eluted from the RNeasy mini spin columns in 30 μl of DEPC-treated water. Using the procedures described above, we isolated total RNA from Hdh-HET MEF lines 1, 2 (in duplicate), 3, and 5 (in duplicate), and from Hdh-KO MEF lines 11, 12, 16 (in duplicate), 27 (in duplicate). The resulting 12 samples were submitted to the UT Southwestern Microarray Core Facility for microarray analysis. Biotinylated cRNA was prepared using the Illumina RNA Amplification Kit (Ambion, Inc., Austin, TX) according to the manufacturer's directions starting with ~200 ng total RNA. Samples were purified using the RNeasy kit (Qiagen, Valencia, CA). Hybridization to the Sentrix Mouse-6 Expression BeadChip (Illumina, Inc., San Diego, CA), washing and scanning were performed according to the Illumina BeadStation 500× manual (revision C). Two BeadChips were used, each one containing 6 arrays. For each chip, three HET and three KO samples were analyzed to minimize the effects of chip-to-chip variability. Arrays were scanned with an Illumina Bead array Reader confocal scanner and the data was analyzed using Illumina's BeadStudio software(Version 3). The raw data was background subtracted and normalized using "cubic spline" method in the software. The detection p values were computed using a dynamically constructed normal model based on the intensities of 700 negative controls. For differential analysis, the six arrays in HET group were compared with the six arrays in KO group using t-test with computing false discovery rate algorithm. The genes with p value <0.05 were considered differentially expressed and subject for further analysis.
Quantitative Real Time PCR
Quantitative real time (qRT)-PCR was performed using an Applied Biosystems Prism 7900HT sequence detection system using SYBR green chemistry. Briefly, total RNA was treated with DNase I (RNase-free, Roche Molecular Biochemicals), and reverse-transcribed with random hexamers using SuperScript II reverse transcriptase (Invitrogen) to generate cDNA as previously described . Primers were designed using Primer Express Software (PerkinElmer Life Sciences) and validated by analysis of template titration and dissociation curves. Each qRT-PCR contained (final volume of 10 μl) 25 ng of reverse-transcribed RNA, each primer at 150 nM, and 5 μl of 2× SYBR Green PCR Master Mix (Applied Biosystems), and each sample was analyzed in triplicate. Results were evaluated by the comparative CT method (User Bulletin No. 2, PerkinElmer Life Sciences) using cyclophilin as the invariant control gene.
MacDonald ME: Huntingtin: alive and well and working in middle management. Sci STKE. 2003, 2003: pe48.
Li SH, Li XJ: Huntingtin-protein interactions and the pathogenesis of Huntington's disease. Trends Genet. 2004, 20: 146-154. 10.1016/j.tig.2004.01.008.
Trottier Y, Lutz Y, Stevanin G, Imbert G, Devys D, Cancel G, Saudou F, Weber C, David G, Tora L, et al: Polyglutamine expansion as a pathological epitope in Huntington's disease and four dominant cerebellar ataxias. Nature. 1995, 378: 403-406. 10.1038/378403a0.
Ferrante RJ, Gutekunst CA, Persichetti F, McNeil SM, Kowall NW, Gusella JF, MacDonald ME, Beal MF, Hersch SM: Heterogeneous topographic and cellular distribution of huntingtin expression in the normal human neostriatum. J Neurosci. 1997, 17: 3052-3063.
Fusco FR, Chen Q, Lamoreaux WJ, Figueredo-Cardenas G, Jiao Y, Coffman JA, Surmeier DJ, Honig MG, Carlock LR, Reiner A: Cellular localization of huntingtin in striatal and cortical neurons in rats: lack of correlation with neuronal vulnerability in Huntington's disease. J Neurosci. 1999, 19: 1189-1202.
DiFiglia M, Sapp E, Chase K, Schwarz C, Meloni A, Young C, Martin E, Vonsattel JP, Carraway R, Reeves SA, et al: Huntingtin is a cytoplasmic protein associated with vesicles in human and rat brain neurons. Neuron. 1995, 14: 1075-1081. 10.1016/0896-6273(95)90346-1.
Velier J, Kim M, Schwarz C, Kim TW, Sapp E, Chase K, Aronin N, DiFiglia M: Wild-type and mutant huntingtins function in vesicle trafficking in the secretory and endocytic pathways. Exp Neurol. 1998, 152: 34-40. 10.1006/exnr.1998.6832.
Kegel KB, Meloni AR, Yi Y, Kim YJ, Doyle E, Cuiffo BG, Sapp E, Wang Y, Qin ZH, Chen JD, Nevins JR, Aronin N, DiFiglia M: Huntingtin is present in the nucleus, interacts with the transcriptional corepressor C-terminal binding protein, and represses transcription. J Biol Chem. 2002, 277: 7466-7476. 10.1074/jbc.M103946200.
Cattaneo E, Zuccato C, Tartari M: Normal huntingtin function: an alternative approach to Huntington's disease. Nat Rev Neurosci. 2005, 6: 919-930. 10.1038/nrn1806.
Nasir J, Floresco SB, O'Kusky JR, Diewert VM, Richman JM, Zeisler J, Borowski A, Marth JD, Phillips AG, Hayden MR: Targeted disruption of the Huntington's disease gene results in embryonic lethality and behavioral and morphological changes in heterozygotes. Cell. 1995, 81: 811-823. 10.1016/0092-8674(95)90542-1.
Duyao MP, Auerbach AB, Ryan A, Persichetti F, Barnes GT, McNeil SM, Ge P, Vonsattel JP, Gusella JF, Joyner AL, et al: Inactivation of the mouse Huntington's disease gene homolog Hdh. Science. 1995, 269: 407-410. 10.1126/science.7618107.
Zeitlin S, Liu JP, Chapman DL, Papaioannou VE, Efstratiadis A: Increased apoptosis and early embryonic lethality in mice nullizygous for the Huntington's disease gene homologue. Nat Genet. 1995, 11: 155-163. 10.1038/ng1095-155.
White JK, Auerbach W, Duyao MP, Vonsattel JP, Gusella JF, Joyner AL, MacDonald ME: Huntingtin is required for neurogenesis and is not impaired by the Huntington's disease CAG expansion. Nat Genet. 1997, 17: 404-410. 10.1038/ng1297-404.
Auerbach W, Hurlbert MS, Hilditch-Maguire P, Wadghiri YZ, Wheeler VC, Cohen SI, Joyner AL, MacDonald ME, Turnbull DH: The HD mutation causes progressive lethal neurological disease in mice expressing reduced levels of huntingtin. Hum Mol Genet. 2001, 10: 2515-2523. 10.1093/hmg/10.22.2515.
Harper SQ, Staber PD, He X, Eliason SL, Martins IH, Mao Q, Yang L, Kotin RM, Paulson HL, Davidson BL: RNA interference improves motor and neuropathological abnormalities in a Huntington's disease mouse model. Proc Natl Acad Sci U S A. 2005, 102: 5820-5825. 10.1073/pnas.0501507102.
Lecerf JM, Shirley TL, Zhu Q, Kazantsev A, Amersdorfer P, Housman DE, Messer A, Huston JS: Human single-chain Fv intrabodies counteract in situ huntingtin aggregation in cellular models of Huntington's disease. Proc Natl Acad Sci U S A. 2001, 98: 4764-4769. 10.1073/pnas.071058398.
Colby DW, Chu Y, Cassady JP, Duennwald M, Zazulak H, Webster JM, Messer A, Lindquist S, Ingram VM, Wittrup KD: Potent inhibition of huntingtin aggregation and cytotoxicity by a disulfide bond-free single-domain intracellular antibody. Proc Natl Acad Sci U S A. 2004, 101: 17616-17621. 10.1073/pnas.0408134101.
Dragatsis I, Levine MS, Zeitlin S: Inactivation of Hdh in the brain and testis results in progressive neurodegeneration and sterility in mice. Nat Genet. 2000, 26: 300-306. 10.1038/81593.
Ho A, Morishita W, Atasoy D, Liu X, Tabuchi K, Hammer RE, Malenka RC, Sudhof TC: Genetic analysis of Mint/X11 proteins: essential presynaptic functions of a neuronal adaptor protein family. J Neurosci. 2006, 26: 13089-13101. 10.1523/JNEUROSCI.2855-06.2006.
Mattson MP: Calcium and neurodegeneration. Aging Cell. 2007, 6: 337-350. 10.1111/j.1474-9726.2007.00275.x.
Bossy-Wetzel E, Schwarzenbacher R, Lipton SA: Molecular pathways to neurodegeneration. Nat Med. 2004, 10 (Suppl): S2-9. 10.1038/nm1067.
Bezprozvanny I, Hayden MR: Deranged neuronal calcium signaling and Huntington disease. Biochem Biophys Res Commun. 2004, 322: 1310-1317. 10.1016/j.bbrc.2004.08.035.
Tang TS, Tu H, Chan EY, Maximov A, Wang Z, Wellington CL, Hayden MR, Bezprozvanny I: Huntingtin and huntingtin-associated protein 1 influence neuronal calcium signaling mediated by inositol-(1,4,5) triphosphate receptor type 1. Neuron. 2003, 39: 227-239. 10.1016/S0896-6273(03)00366-0.
Kaltenbach LS, Romero E, Becklin RR, Chettier R, Bell R, Phansalkar A, Strand A, Torcassi C, Savage J, Hurlburt A, Cha GH, Ukani L, Chepanoske CL, Zhen Y, Sahasrabudhe S, Olson J, Kurschner C, Ellerby LM, Peltier JM, Botas J, Hughes RE: Huntingtin interacting proteins are genetic modifiers of neurodegeneration. PLoS Genet. 2007, 3: e82-10.1371/journal.pgen.0030082.
Panov AV, Gutekunst CA, Leavitt BR, Hayden MR, Burke JR, Strittmatter WJ, Greenamyre JT: Early mitochondrial calcium defects in Huntington's disease are a direct effect of polyglutamines. Nat Neurosci. 2002, 5: 731-736.
Choo YS, Johnson GV, MacDonald M, Detloff PJ, Lesort M: Mutant huntingtin directly increases susceptibility of mitochondria to the calcium-induced permeability transition and cytochrome c release. Hum Mol Genet. 2004, 13: 1407-1420. 10.1093/hmg/ddh162.
Oliveira JM, Jekabsons MB, Chen S, Lin A, Rego AC, Goncalves J, Ellerby LM, Nicholls DG: Mitochondrial dysfunction in Huntington's disease: the bioenergetics of isolated and in situ mitochondria from transgenic mice. J Neurochem. 2007, 101: 241-249. 10.1111/j.1471-4159.2006.04361.x.
Brustovetsky N, LaFrance R, Purl KJ, Brustovetsky T, Keene CD, Low WC, Dubinsky JM: Age-dependent changes in the calcium sensitivity of striatal mitochondria in mouse models of Huntington's Disease. J Neurochem. 2005, 93: 1361-1370. 10.1111/j.1471-4159.2005.03036.x.
Seong IS, Ivanova E, Lee JM, Choo YS, Fossale E, Anderson M, Gusella JF, Laramie JM, Myers RH, Lesort M, MacDonald ME: HD CAG repeat implicates a dominant property of huntingtin in mitochondrial energy metabolism. Hum Mol Genet. 2005, 14: 2871-2880. 10.1093/hmg/ddi319.
Oliveira JM, Chen S, Almeida S, Riley R, Goncalves J, Oliveira CR, Hayden MR, Nicholls DG, Ellerby LM, Rego AC: Mitochondrial-dependent Ca2+ handling in Huntington's disease striatal cells: effect of histone deacetylase inhibitors. J Neurosci. 2006, 26: 11174-11186. 10.1523/JNEUROSCI.3004-06.2006.
Rizzuto R, Duchen MR, Pozzan T: Flirting in little space: the ER/mitochondria Ca2+ liaison. Sci STKE. 2004, 2004: re1.
Csordas G, Renken C, Varnai P, Walter L, Weaver D, Buttle KF, Balla T, Mannella CA, Hajnoczky G: Structural and functional features and significance of the physical linkage between ER and mitochondria. J Cell Biol. 2006, 174: 915-921. 10.1083/jcb.200604016.
Zuccato C, Tartari M, Crotti A, Goffredo D, Valenza M, Conti L, Cataudella T, Leavitt BR, Hayden MR, Timmusk T, Rigamonti D, Cattaneo E: Huntingtin interacts with REST/NRSF to modulate the transcription of NRSE-controlled neuronal genes. Nat Genet. 2003, 35: 76-83. 10.1038/ng1219.
Truant R, Atwal RS, Burtnik A: Nucleocytoplasmic trafficking and transcription effects of huntingtin in Huntington's disease. Prog Neurobiol. 2007
record GEO: . [http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE11139]
Strehlow AN, Li JZ, Myers RM: Wild-type huntingtin participates in protein trafficking between the Golgi and the extracellular space. Hum Mol Genet. 2007, 16: 391-409. 10.1093/hmg/ddl467.
Kuhn A, Goldstein DR, Hodges A, Strand AD, Sengstag T, Kooperberg C, Becanovic K, Pouladi MA, Sathasivam K, Cha JH, Hannan AJ, Hayden MR, Leavitt BR, Dunnett SB, Ferrante RJ, Albin R, Shelbourne P, Delorenzi M, Augood SJ, Faull RL, Olson JM, Bates GP, Jones L, Luthi-Carter R: 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-1861. 10.1093/hmg/ddm133.
Avilion AA, Nicolis SK, Pevny LH, Perez L, Vivian N, Lovell-Badge R: Multipotent cell lineages in early mouse development depend on SOX2 function. Genes Dev. 2003, 17: 126-140. 10.1101/gad.224503.
Sipione S, Rigamonti D, Valenza M, Zuccato C, Conti L, Pritchard J, Kooperberg C, Olson JM, Cattaneo E: Early transcriptional profiles in huntingtin-inducible striatal cells by microarray analyses. Hum Mol Genet. 2002, 11: 1953-1965. 10.1093/hmg/11.17.1953.
Trushina E, Singh RD, Dyer RB, Cao S, Shah VH, Parton RG, Pagano RE, McMurray CT: Mutant huntingtin inhibits clathrin-independent endocytosis and causes accumulation of cholesterol in vitro and in vivo. Hum Mol Genet. 2006, 15: 3578-3591. 10.1093/hmg/ddl434.
Tang TS, Slow EJ, Lupu V, Stavrovskaya IG, Sugimori M, Llinas R, Kristal BS, Hayden MR, Bezprozvanny I: Disturbed Ca2+ signaling and apoptosis of medium spiny neurons in Huntington's disease. Proc Natl Acad Sci U S A. 2005, 102: 2602-2607. 10.1073/pnas.0409402102.
Nagai T, Sawano A, Park ES, Miyawaki A: Circularly permuted green fluorescent proteins engineered to sense Ca2+. Proc Natl Acad Sci U S A. 2001, 98: 3197-3202. 10.1073/pnas.051636098.
Csordas G, Hajnoczky G: Sorting of calcium signals at the junctions of endoplasmic reticulum and mitochondria. Cell Calcium. 2001, 29: 249-262. 10.1054/ceca.2000.0191.
Kurrasch DM, Huang J, Wilkie TM, Repa JJ: Quantitative real-time polymerase chain reaction measurement of regulators of G-protein signaling mRNA levels in mouse tissues. Methods Enzymol. 2004, 389: 3-15.
We are grateful to the personal of UTSW Microarray Core for assistance with these experiments and to Janet Young for administrative assistance. We are thankful to Thomas Südhof and Katsuhiko Tabuchi for a generous gift of Lenti-NLS-GFP-Cre virus. This study was supported by the Hereditary Disease Foundation and NINDS R01 NS38082 and R01 NS056224 (IB), NINDS R01 NS043466 (SZ), the NIH GM59419 (GH), and the Ara Parseghian Medical Research Foundation (JR).
HZ generated the Hdh-KO and Hdh-HET cell lines, carried out microarray experiments and prepared results for publication, SD and GH performed Ca2+ imaging studies and prepared results for publication, QZL assisted with microarray experiments and analysis and prepared results for publication, ID and SZ generated critical reagents and help with interpretation of obtained results, JJR assisted with qPCR experiments and analysis, IB conceived the study, and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Hdh -KO MEF cells. The genes are classified into significant GO categories (classes) as explained in the text. The official gene symbol and protein name are shown. Also shown is average intensity of microarray signal for Hdh-HET and Hdh-KO samples and the calculated Diff score and p-values. Negative Diff scores correspond to genes underexpressed in Hdh-KO cells, positive Diff Scores correspond to genes overexpressed in Hdh-KO cells. The fold change is calculated by dividing average Hdh-HET and Hdh-KO signals. The transcripts which are not detected in Hdh-HET samples are shown as "up", the samples not detected in Hdh-KO samples are shown as "down". The Genebank accession numbers are also shown. (DOC 226 KB)
Additional file 2: . The microarray data for 821 differentially expressed probes are presented in the Excell file. The data listed for each probe are as follows: an Illumina beads target ID, an average signal from 6 chips with Hdh-HET samples, the probability of detection of a probe in Hdh-HET samples, an average signal from 6 chips with Hdh-KO samples, the probability of detection of a probe in Hdh-KO samples, the Diff Score obtained by comparing Hdh-HET and Hdh-HET samples, the P value calculated from the Diff score, an official gene symbol and definition (when available), the sequence of the probe, and the genebank accession number (when available). The results are sorted by the Diff score with negative Diff scores corresponding to the probes downregulated in Hdh-KO cells (366 probes) and positive Diff scores corresponding to the probes upregulated in Hdh-KO cells (455 probes). (XLS 287 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Zhang, H., Das, S., Li, QZ. et al. Elucidating a normal function of huntingtin by functional and microarray analysis of huntingtin-null mouse embryonic fibroblasts. BMC Neurosci 9, 38 (2008). https://doi.org/10.1186/1471-2202-9-38
- Mouse Embryonic Fibroblast
- Mouse Embryonic Fibroblast Cell
- Early Embryonic Lethality
- Mouse Embryonic Fibroblast Cell Line
- Wild Type Embryonic Stem Cell