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

To test the hypotheses that mutant huntingtin protein length and wild-type huntingtin dosage have important effects on disease-related transcriptional dysfunction, we compared the changes in mRNA in seven genetic mouse models of Huntington's disease (HD) and postmortem human HD caudate. Transgenic models expressing short N-terminal fragments of mutant huntingtin (R6/1 and R6/2 mice) exhibited the most rapid effects on gene expression, consistent with previous studies. Although changes in the brains of knock-in and full-length transgenic models of HD took longer to appear, 15- and 22-month CHL2(Q150/Q150), 18-month Hdh(Q92/Q92) and 2-year-old YAC128 animals also exhibited significant HD-like mRNA signatures. Whereas it was expected that the expression of full-length huntingtin transprotein might result in unique gene expression changes compared with those caused by the expression of an N-terminal huntingtin fragment, no discernable differences between full-length and fragment models were detected. In addition, very high correlations between the signatures of mice expressing normal levels of wild-type huntingtin and mice in which the wild-type protein is absent suggest a limited effect of the wild-type protein to change basal gene expression or to influence the qualitative disease-related effect of mutant huntingtin. The combined analysis of mouse and human HD transcriptomes provides important temporal and mechanistic insights into the process by which mutant huntingtin kills striatal neurons. In addition, the discovery that several available lines of HD mice faithfully recapitulate the gene expression signature of the human disorder provides a novel aspect of validation with respect to their use in preclinical therapeutic trials.

Congenital hypothyroidism (cretinism) in neuroD2-deficient mice

Mice lacking neuroD2, a basic helix-loop-helix transcription factor involved in brain development, show growth retardation and other abnormalities consistent with hypothalamic-pituitary-thyroid (HPT) axis dysfunction. neuroD2 is expressed in the paraventricular hypothalamic nuclei, the anterior lobe of pituitary, and the thyroid gland. In neuroD2-deficient mice, thyrotropin-releasing hormone, thyroid-stimulating hormone, and thyroid hormone are decreased in these three regions, respectively. neuroD2-null mice typically die 2 to 3 weeks after birth, but those treated with replacement doses of thyroxine survived more than 8 weeks. These data indicate that neuroD2 is expressed throughout the HPT axis and that all levels of the axis are functionally affected by its absence in mice.

Suberoylanilide hydroxamic acid is effective in preclinical studies of medulloblastoma

PURPOSE

Suberoylanilide hydroxamic acid (SAHA) has been studied in adult solid and hematologic malignancies. However, little information has been reported on the effects of SAHA on central nervous system (CNS) tumors including medulloblastoma, the most common malignant brain tumor in children. We investigated SAHA in preclinical medulloblastoma models to determine its anti-cancer efficacy as well as its ability to affect intracranial lesions when administered systemically.

EXPERIMENTAL DESIGN AND RESULTS

Tissue culture studies were performed treating primary human fibroblasts, established medulloblastoma cell lines, and primary human medulloblastoma tumors with SAHA. At 10 microM concentration, SAHA had little effect on normal fibroblasts but caused >90% apoptosis in cultured medulloblastoma cells. Primary medulloblastomas from patients were sensitive to SAHA compared to vehicle alone in ex vivo studies. In athymic mice with medulloblastoma xenograft tumors, oral SAHA resulted in apoptosis of tumor tissue and significantly slowed tumor growth. In the ND2:Smo transgenic mouse medulloblastoma model, SAHA treatment caused significant apoptosis in these cerebellar tumors.

CONCLUSIONS

SAHA effectively induces cell death in established medulloblastoma cell lines, human patient primary tumor cultures, medulloblastoma xenografts and intracranial spontaneous medulloblastomas. Fibroblasts in culture and mice treated with SAHA did not reveal prohibitive toxicity profiles. These findings support the advancement of SAHA to pediatric clinical trials.

Assessment of the relationship between pre-chip and post-chip quality measures for Affymetrix GeneChip expression data

BACKGROUND

Gene expression microarray experiments are expensive to conduct and guidelines for acceptable quality control at intermediate steps before and after the samples are hybridised to chips are vague. We conducted an experiment hybridising RNA from human brain to 117 U133A Affymetrix GeneChips and used these data to explore the relationship between 4 pre-chip variables and 22 post-chip outcomes and quality control measures.

RESULTS

We found that the pre-chip variables were significantly correlated with each other but that this correlation was strongest between measures of RNA quality and cRNA yield. Post-mortem interval was negatively correlated with these variables. Four principal components, reflecting array outliers, array adjustment, hybridisation noise and RNA integrity, explain about 75% of the total post-chip measure variability. Two significant canonical correlations existed between the pre-chip and post-chip variables, derived from MAS 5.0, dChip and the Bioconductor packages affy and affyPLM. The strongest (CANCOR 0.838, p < 0.0001) correlated RNA integrity and yield with post chip quality control (QC) measures indexing 3'/5' RNA ratios, bias or scaling of the chip and scaling of the variability of the signal across the chip. Post-mortem interval was relatively unimportant. We also found that the RNA integrity number (RIN) could be moderately well predicted by post-chip measures B_ACTIN35, GAPDH35 and SF.

CONCLUSION

We have found that the post-chip variables having the strongest association with quantities measurable before hybridisation are those reflecting RNA integrity. Other aspects of quality, such as noise measures (reflecting the execution of the assay) or measures reflecting data quality (outlier status and array adjustment variables) are not well predicted by the variables we were able to determine ahead of time. There could be other variables measurable pre-hybridisation which may be better associated with expression data quality measures. Uncovering such connections could create savings on costly microarray experiments by eliminating poor samples before hybridisation.

A statistical method for predicting splice variants between two groups of samples using GeneChip expression array data

BACKGROUND

Alternative splicing of pre-messenger RNA results in RNA variants with combinations of selected exons. It is one of the essential biological functions and regulatory components in higher eukaryotic cells. Some of these variants are detectable with the Affymetrix GeneChip that uses multiple oligonucleotide probes (i.e. probe set), since the target sequences for the multiple probes are adjacent within each gene. Hybridization intensity from a probe correlates with abundance of the corresponding transcript. Although the multiple-probe feature in the current GeneChip was designed to assess expression values of individual genes, it also measures transcriptional abundance for a sub-region of a gene sequence. This additional capacity motivated us to develop a method to predict alternative splicing, taking advance of extensive repositories of GeneChip gene expression array data.

RESULTS

We developed a two-step approach to predict alternative splicing from GeneChip data. First, we clustered the probes from a probe set into pseudo-exons based on similarity of probe intensities and physical adjacency. A pseudo-exon is defined as a sequence in the gene within which multiple probes have comparable probe intensity values. Second, for each pseudo-exon, we assessed the statistical significance of the difference in probe intensity between two groups of samples. Differentially expressed pseudo-exons are predicted to be alternatively spliced. We applied our method to empirical data generated from GeneChip Hu6800 arrays, which include 7129 probe sets and twenty probes per probe set. The dataset consists of sixty-nine medulloblastoma (27 metastatic and 42 non-metastatic) samples and four cerebellum samples as normal controls. We predicted that 577 genes would be alternatively spliced when we compared normal cerebellum samples to medulloblastomas, and predicted that thirteen genes would be alternatively spliced when we compared metastatic medulloblastomas to non-metastatic ones. We checked the consistency of some of our findings with information in UCSC Human Genome Browser.

CONCLUSION

The two-step approach described in this paper is capable of predicting some alternative splicing from multiple oligonucleotide-based gene expression array data with GeneChip technology. Our method employs the extensive repositories of gene expression array data available and generates alternative splicing hypotheses, which can be further validated by experimental studies.

Regulation of Thalamocortical Patterning and Synaptic Maturation by NeuroD2

During cortical development, both activity-dependent and genetically determined mechanisms are required to establish proper neuronal connectivity. While activity-dependent transcription may link the two processes, specific transcription factors that mediate such a process have not been identified. We identified the basic helix-loop-helix (bHLH) transcription factor Neurogenic Differentiation 2 (NeuroD2) in a screen for calcium-regulated transcription factors and report that it is required for the proper development of thalamocortical connections. In neuroD2 null mice, thalamocortical axon terminals fail to segregate in the somatosensory cortex, and the postsynaptic barrel organization is disrupted. Additionally, synaptic transmission is defective at thalamocortical synapses in neuroD2 null mice. Total excitatory synaptic currents are reduced in layer IV in the knockouts, and the relative contribution of AMPA and NMDA receptor-mediated currents to evoked responses is decreased. These observations indicate that NeuroD2 plays a critical role in regulating synaptic maturation and the patterning of thalamocortical connections.

Regional and cellular gene expression changes in human Huntington's disease brain

Huntington's disease (HD) pathology is well understood at a histological level but a comprehensive molecular analysis of the effect of the disease in the human brain has not previously been available. To elucidate the molecular phenotype of HD on a genome-wide scale, we compared mRNA profiles from 44 human HD brains with those from 36 unaffected controls using microarray analysis. Four brain regions were analyzed: caudate nucleus, cerebellum, prefrontal association cortex [Brodmann's area 9 (BA9)] and motor cortex [Brodmann's area 4 (BA4)]. The greatest number and magnitude of differentially expressed mRNAs were detected in the caudate nucleus, followed by motor cortex, then cerebellum. Thus, the molecular phenotype of HD generally parallels established neuropathology. Surprisingly, no mRNA changes were detected in prefrontal association cortex, thereby revealing subtleties of pathology not previously disclosed by histological methods. To establish that the observed changes were not simply the result of cell loss, we examined mRNA levels in laser-capture microdissected neurons from Grade 1 HD caudate compared to control. These analyses confirmed changes in expression seen in tissue homogenates; we thus conclude that mRNA changes are not attributable to cell loss alone. These data from bona fide HD brains comprise an important reference for hypotheses related to HD and other neurodegenerative diseases.

Dysfunction of the cholesterol biosynthetic pathway in Huntington's disease

The expansion of a polyglutamine tract in the ubiquitously expressed huntingtin protein causes Huntington's disease (HD), a dominantly inherited neurodegenerative disease. We show that the activity of the cholesterol biosynthetic pathway is altered in HD. In particular, the transcription of key genes of the cholesterol biosynthetic pathway is severely affected in vivo in brain tissue from HD mice and in human postmortem striatal and cortical tissue; this molecular dysfunction is biologically relevant because cholesterol biosynthesis is reduced in cultured human HD cells, and total cholesterol mass is significantly decreased in the CNS of HD mice and in brain-derived ST14A cells in which the expression of mutant huntingtin has been turned on. The transcription of the genes of the cholesterol biosynthetic pathway is regulated via the activity of sterol regulatory element-binding proteins (SREBPs), and we found an approximately 50% reduction in the amount of the active nuclear form of SREBP in HD cells and mouse brain tissue. As a consequence, mutant huntingtin reduces the transactivation of an SRE-luciferase construct even under conditions of SREBP overexpression or in the presence of an exogenous N-terminal active form of SREBP. Finally, the addition of exogenous cholesterol to striatal neurons expressing mutant huntingtin prevents their death in a dose-dependent manner. We conclude that the cholesterol biosynthetic pathway is impaired in HD cells, mice, and human subjects, and that the search for HD therapies should also consider cholesterol levels as both a potential target and disease biomarker.

Cisplatin-based chemotherapy followed by focal, reduced-dose irradiation for pediatric primary central nervous system germinomas

The objective of this study was to evaluate retrospectively one institution's experience treating pediatric central nervous system (CNS) pure germinomas with platinum-based chemotherapy followed by focal, reduced-dose irradiation. Eight patients were identified with localized, pure CNS germinomas from 1993 to 2004 at the authors' institution. The median age at diagnosis was 13 years (range 7-19). The median follow-up was 40 months (range 8-141). The tumor location was suprasellar in four, the pineal region in three, and the third ventricle in one. Irradiation was started a median of 20 weeks (range 17-22) from diagnosis and consisted of conformal fields to the primary site as determined by the initial diagnostic MR plus a 1.5- to 2-cm margin. Six of the eight patients received a dose of 3,060 cGy; two patients received 3,600 cGy. The 5-year actuarial event free survival was 71% (56-86%, 95% CI). Two patients suffered marginal (at field edge) failures and both were salvaged using reinduction platinum-based chemotherapy followed by cranial spinal irradiation and a boost to the primary tumor. The 5-year actuarial overall survival was 100%. There were no spinal failures. These data suggest that a reduction in both volume and dose (30.6-36 Gy) retains the excellent survival rates for patients with localized, pure germinomas of the CNS. A higher rate of ventricular relapse rate is observed, although salvage of those patients is feasible.

Mutant huntingtin alters MAPK signaling pathways in PC12 and striatal cells: ERK1/2 protects against mutant huntingtin-associated toxicity

Huntington's disease (HD) is a devastating neurodegenerative disorder caused by an expanded polyglutamine (polyQ) tract within the huntingtin protein (Htt). Identifying the pathways that are altered in response to the mutant protein is crucial for understanding the cellular processes impacted by the disease as well as for the rational development of effective pharmacological interventions. Here, expression profiling of a cellular HD model identifies genes that implicate altered mitogen-activated protein kinase (MAPK) signaling. Targeted biochemical studies and pharmacological modulation of these MAPK pathways suggest that mutant Htt affects signaling at upstream points such that both ERK and JNK are activated. Modulation of the ERK pathway suggests that this pathway is associated with cell survival, whereas inhibition of JNK was found to effectively suppress pathogenesis. These studies suggest that pharmacological intervention in MAPK pathways, particularly at the level of ERK activation, may be an appropriate approach to HD therapy.

The dosage of the neuroD2 transcription factor regulates amygdala development and emotional learning

The amygdala is centrally involved in formation of emotional memory and response to fear or risk. We have demonstrated that the lateral and basolateral amygdala nuclei fail to form in neuroD2 null mice and neuroD2 heterozygotes have fewer neurons in this region. NeuroD2 heterozygous mice show profound deficits in emotional learning as assessed by fear conditioning. Unconditioned fear was also diminished in neuroD2 heterozygotes compared to wild-type controls. Several key molecular regulators of emotional learning were diminished in the brains of neuroD2 heterozygotes including Ulip1, alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptor, and GABA(A) receptor. Thus, neuroD2 is essential for amygdala development and genes involved in amygdala function are altered in neuroD2-deficient mice.

Contribution of nuclear and extranuclear polyQ to neurological phenotypes in mouse models of Huntington's disease

In postmortem Huntington's disease brains, mutant htt is present in both nuclear and cytoplasmic compartments. To dissect the impact of nuclear and extranuclear mutant htt on the initiation and progression of disease, we generated a series of transgenic mouse lines in which nuclear localization or nuclear export signal sequences have been placed N-terminal to the htt exon 1 protein carrying 144 glutamines. Our data indicate that the exon 1 mutant protein is present in the nucleus as part of an oligomeric or aggregation complex. Increasing the concentration of the mutant transprotein in the nucleus is sufficient for and dramatically accelerates the onset and progression of behavioral phenotypes. Furthermore, nuclear exon 1 mutant protein is sufficient to induce cytoplasmic neurodegeneration and transcriptional dysregulation. However, our data suggest that cytoplasmic mutant exon 1 htt, if present, contributes to disease progression.

Gene expression in Huntington's disease skeletal muscle: a potential biomarker

Huntington's disease (HD) is an incurable and fatal neurodegenerative disorder. Improvements in the objective measurement of HD will lead to more efficient clinical trials and earlier therapeutic intervention. We hypothesized that abnormalities seen in the R6/2 mouse, a greatly accelerated HD model, might highlight subtle phenotypes in other mouse models and human HD. In this paper, we identify common gene expression changes in skeletal muscle from R6/2 mice, Hdh(CAG(150)) homozygous knock-in mice and HD patients. This HD-triggered gene expression phenotype is consistent with the beginnings of a transition from fast-twitch to slow-twitch muscle fiber types. Metabolic adaptations similar to those induced by diabetes or fasting are also present but neither metabolic disorder can explain the full phenotype of HD muscle. The HD-induced gene expression changes reflect disease progression. This raises the possibility that muscle gene expression may be used as an objective biomarker to complement clinical HD-rating systems. Furthermore, an understanding of the molecular basis of muscle dysfunction in HD should provide insight into mechanisms involved in neuronal abnormalities and neurodegeneration.

Step structures on III-V phosphide (001) surfaces: how do steps and Sb affect CuPt ordering of GaInP2?

The observation of III-V phosphide (001)-(2 x 2) surfaces makes it possible to solve a long standing mystery of step structures. First-principles calculations show that a bulklike type-B step on a hydrogenated 2 x 2 surface is more stable than a rebonded one by 1.1 eV/unit step. In contrast, this energy difference for a H-free beta(2 x 4) surface is only 0.5 eV/unit step. The large difference explains why the CuPt ordering of GaInP is stronger in metal-organic chemical vapor deposition than in molecular beam epitaxy. However, a minute amount of Sb will preferentially attach to the 2 x 2 surface steps and induce additional step structures that cause ordering disruption.

A class of models for analyzing GeneChip gene expression analysis array data

Background

Various analytical methods exist that first quantify gene expression and then analyze differentially expressed genes from Affymetrix GeneChip^® gene expression analysis array data. These methods differ in the choice of probe measure (quantification of probe hybridization), summarizing multiple probe intensities into a gene expression value, and analysis of differential gene expression. Research papers that describe these methods focus on performance, and how their approaches differ from others. To better understand the common features and differences between various methods, and to evaluate their impact on the results of gene expression analysis, we describe a class of models, referred to as generalized probe models (GPMs), which encompass various currently available methods.

Results

Using an empirical dataset, we compared different formulations of GPMs, and GPMs with three other commonly used methods, i.e. MAS 5.0, dChip, and RMA. The comparison shows that, on a genome-wide scale , different methods yield similar results if the same probe measures are chosen.

Conclusion

In this paper we present a general framework, i.e. GPMs, which encompasses various methods. GPMs permit the use of a wide range of probe measures and facilitate appropriate comparison between commonly used methods. We demonstrate that the dissimilar results stem primarily from different choice of probe measures, rather than other factors.

Significance testing for small microarray experiments

Which significance test is carried out when the number of repeats is small in microarray experiments can dramatically influence the results. When in two sample comparisons both conditions have fewer than, say, five repeats traditional test statistics require extreme results, before a gene is considered statistically significant differentially expressed after a multiple comparisons correction. In the literature many approaches to circumvent this problem have been proposed. Some of these proposals use (empirical) Bayes arguments to moderate the variance estimates for individual genes. Other proposals try to stabilize these variance estimate by combining groups of genes or similar experiments. In this paper we compare several of these approaches, both on data sets where both experimental conditions are the same, and thus few statistically significant differentially expressed genes should be identified, and on experiments where both conditions do differ. This allows us to identify which approaches are most powerful without identifying many false positives. We conclude that after balancing the numbers of false positives and true positives an empirical Bayes approach and an approach which combines experiments perform best. Standard t-tests are inferior and offer almost no power when the sample size is small.

p38 MAP kinase: a convergence point in cancer therapy

Recent studies show that activation of p38 mitogen-activated protein kinase (MAPK) results in cancer cell apoptosis initiated by retinoids, cisplatin and other chemotherapeutic agents. The observation that divergent therapies act through a common signal transduction pathway raises the possibility of developing new anti-cancer agents that lack the side-effects caused by events upstream of p38 MAPK. Here, we review p38-MAPK-mediated tumor cell apoptosis and implications for cancer therapeutics.

The SmoA1 Mouse Model Reveals That Notch Signaling Is Critical for the Growth and Survival of Sonic Hedgehog-Induced Medulloblastomas

To develop a genetically faithful model of medulloblastoma with increased tumor incidence compared with the current best model we activated the Sonic Hedgehog (Shh) pathway by transgenically expressing a constitutively active form of Smoothened in mouse cerebellar granule neuron precursors (ND2:SmoA1 mice). This resulted in early cerebellar granule cell hyper-proliferation and a 48% incidence of medulloblastoma formation. Gene expression studies showed an increase in the known Shh targets Gli1 and Nmyc that correlated with increasing hyperplasia and tumor formation. Notch2 and the Notch target gene, HES5, were also significantly elevated in Smoothened-induced tumors showing that Shh pathway activation is sufficient to induce Notch pathway signaling. In human medulloblastomas reverse transcription-PCR for Shh and Notch targets revealed activation of both of these pathways in most tumors when compared with normal cerebellum. Notch pathway inhibition with soluble Delta ligand or gamma secretase inhibitors resulted in a marked reduction of viable cell numbers in medulloblastoma cell lines and primary tumor cultures. Treatment of mice with D283 medulloblastoma xenografts with a gamma secretase inhibitor resulted in decreased proliferation and increased apoptosis, confirming that Notch signaling contributes to human medulloblastoma proliferation and survival. Medulloblastomas in ND2:SmoA1 mice and humans have concomitant increase in Shh and Notch pathway activities, both of which contribute to tumor survival.

The vulnerability of values to attack: inoculation of values and value-relevant attitudes

Based on the values-as-truisms hypothesis and inoculation theory, two experiments tested whether providing cognitive defenses for the value of equality induces resistance against a message attacking this value. Experiment 1 found that participants who generated cognitive support in an active-supportive or an active-refutational defense were less persuaded by a subsequent message attacking equality than were participants who engaged in no prior defense. Experiment 2 examined the effects of an active-refutational defense and a passive-refutational defense, which simply asked participants to read reasons supporting or opposing equality. Results indicated additive effects of the active and passive defenses, such that participants were most resistant to the anti-equality message when they were given both defenses. Mediational analysis across both experiments revealed that the defenses increased counterargumentation of the anti-equality message, which led to increased post-attack importance of equality and predicted more favorable equality-relevant attitudes and values.

Regulation of neuroD2 expression in mouse brain

The basic helix-loop-helix (bHLH) transcription factor, neuroD2, induces neuronal differentiation and promotes neuronal survival. Reduced levels of neuroD2 were previously shown to cause motor deficits, ataxia, and seizure propensity. Because neuroD2 levels may be critical for brain function, we studied the regulation of neuroD2 gene in cell culture and transgenic mouse models. In transgenic mice, a 10-kb fragment of the neuroD2 promoter fully recapitulated the endogenous neuroD2 staining pattern. A 1-kb fragment of the neuroD2 promoter drove reporter gene expression in most, but not all neuroD2-positive neuronal populations. Mutation of two critical E-boxes, E4 and E5 (E4 and E5 situated 149 and 305 bp upstream of the transcriptional start site) eliminated gene expression. NeuroD2 expression was diminished in mice lacking neurogenin1 demonstrating that neurogenin1 regulates neuroD2 during murine brain development. These studies demonstrate that neuroD2 expression is highly dependent on bHLH-responsive E-boxes in the proximal promoter region, that additional distal regulatory elements are important for neuroD2 expression in a subset of cortical neurons, and that neurogenin1 regulates neuroD2 expression during mouse brain development.

Genetic heterogeneity of stably transfected cell lines revealed by expression profiling with oligonucleotide microarrays

Large-scale gene expression measurements with oligonucleotide microarrays have contributed tremendously to biological research. However, to distinguish between relevant expression changes and falsely identified positives, the source and magnitude of errors must be understood. Here, we report a source of biological variability in microarray experiments with stably transfected cell lines. Mouse embryonic fibroblast (MEF/3T3) and rat schwannoma (RT4) cell lines were generated to provide regulatable schwannomin expression. The expression levels of 29 samples from five different mouse embryonic fibroblast clonal cell lines and 18 samples from 3 RT4 cell lines were monitored with oligonucleotide microarrays. Using hierarchical clustering, we determined that the changes in gene expression induced by schwannomin overexpression were subtle when compared with those detected as a consequence of clonal selection during generation of the cell lines. The hierarchical clustering implies that significant alterations of gene expression were introduced during the transfection and selection processes. A total of 28 genes were identified by Kruskal-Wallis rank test that showed significant variation between clonal lines. Most of them were related to cytoskeletal function and signaling pathways. Based on these analyses, we recommend that replications of experiments with several selected cell lines are necessary to assess biological effects of induced gene expression.

BMP-2 mediates retinoid-induced apoptosis in medulloblastoma cells through a paracrine effect

The mechanisms of retinoid activity in tumors remain largely unknown. Here we establish that retinoids cause extensive apoptosis of medulloblastoma cells. In a xenograft model, retinoids largely abrogated tumor growth. Using receptor-specific retinoid agonists, we defined a subset of mRNAs that were induced by all active retinoids in retinoid-sensitive cell lines. We also identified bone morphogenetic protein-2 (BMP-2) as a candidate mediator of retinoid activity. BMP-2 protein induced medulloblastoma cell apoptosis, whereas the BMP-2 antagonist noggin blocked both retinoid and BMP-2-induced apoptosis. BMP-2 also induced p38 mitogen-activated protein kinase (MAPK), which is necessary for BMP-2- and retinoid-induced apoptosis. Retinoid-resistant medulloblastoma cells underwent apoptosis when treated with BMP-2 or when cultured with retinoid-sensitive medulloblastoma cells. Retinoid-induced expression of BMP-2 is thus necessary and sufficient for apoptosis of retinoid-responsive cells, and expression of BMP-2 by retinoid-sensitive cells is sufficient to induce apoptosis in surrounding retinoid-resistant cells.

A genome screen of systemic lupus erythematosus using affected-relative-pair linkage analysis with covariates demonstrates genetic heterogeneity

Systemic lupus erythematosus (SLE) appears to be the consequence of complex genetics and of only partly understood environmental contributions. Previous work by ourselves and by others has established genetic effects on 1q, 2q, 4p, 6p, and 16p using SLE as the phenotype. However, individual SLE affecteds are extraordinarily different from one another by clinical and laboratory measures. This variation may have a genetic basis; if so, it is advantageous to incorporate measures of between-family clinical variability as covariates in a genetic linkage analysis of affected relative pairs (ARPs) to allow for locus heterogeneity. This approach was applied to genome scan marker data from 160 pedigrees multiplex for SLE and containing 202 ARPs. Because the number of potential covariates was large, we used both ad hoc methods and formal principal components analysis to construct four composite covariates using the SLE classification criteria plus age of onset, ethnicity, and sex. Linkage analysis without covariates has detected evidence for linkage at 1q22-24, 2q37, 4p16, 12p12-11, and 17p13. Linkage analysis with these covariates uncovered linkage at 13p11, 17q11-25, and 20q12 and greatly improved evidence for linkage at 1q22-24, 2q37, 12p12-11, and 17p13. Follow-up analysis identified the original variables contributing to locus heterogeneity in each of these locations. In conclusion, allowing for locus heterogeneity through the incorporation of covariates in linkage analysis is a useful way to dissect the genetic contributions to SLE and uncover new genetic effects.

Estimating the statistical significance of gene expression changes observed with oligonucleotide arrays

We present a simple method to assign approximate P-values to gene expression changes detected with Affymetrix oligonucleotide arrays and software. The method pools data for groups of genes and a small number of like-to-like comparisons in order to estimate the significance of changes observed for single genes in comparisons of experimental interest. Statistical significance levels are based on the observed variability in the fractional majority of probe pairs that indicate increasing or decreasing differential expression in comparisons of technical replicates. From this reference distribution or error model, we compute the expected frequency for fractional majorities in comparisons for N > or = 2. These computed distributions are the source of P-value estimates for changes seen in the experimental comparisons. The method is intended to complement the Affymetrix software and to rationalize gene selection for experimental designs involving limited replication.

Evaluating test statistics to select interesting genes in microarray experiments

A randomization procedure to evaluate the significance level and the false-discovery rate in complex microarray experiments is proposed. A related graph can be used to compare different test statistics that can be used to analyze the same experiment. This graph is closely related to receiver operator characteristic (ROC) curves. The proposed method is applied to a subset of the data from a cell-line experiment related to Huntington's disease. A small simulation study compares the effectiveness of the proposed procedure with the significance analysis of microarrays (SAM) procedure.

Medulloblastoma growth inhibition by hedgehog pathway blockade

Constitutive Hedgehog (Hh) pathway activity is associated with initiation of neoplasia, but its role in the continued growth of established tumors is unclear. Here, we investigate the therapeutic efficacy of the Hh pathway antagonist cyclopamine in preclinical models of medulloblastoma, the most common malignant brain tumor in children. Cyclopamine treatment of murine medulloblastoma cells blocked proliferation in vitro and induced changes in gene expression consistent with initiation of neuronal differentiation and loss of neuronal stem cell-like character. This compound also caused regression of murine tumor allografts in vivo and induced rapid death of cells from freshly resected human medulloblastomas, but not from other brain tumors, thus establishing a specific role for Hh pathway activity in medulloblastoma growth.

Polyglutamine and transcription: gene expression changes shared by DRPLA and Huntington's disease mouse models reveal context-independent effects

Recent evidence indicates that transcriptional abnormalities may play an important role in the pathophysiology of polyglutamine diseases. In the present study, we have explored the extent to which polyglutamine-related changes in gene expression may be independent of protein context by comparing mouse models of dentatorubral-pallidoluysian atrophy (DRPLA) and Huntington's disease (HD). Microarray gene expression profiling was conducted in mice of the same background strain in which the same promoter was employed to direct the expression of full-length atrophin-1 or partial huntingtin transproteins (At-65Q or N171-82Q mice). A large number of overlapping gene expression changes were observed in the cerebella of At-65Q and N171-82Q mice. Six of the gene expression changes common to both huntingtin and atrophin-1 transgenic mice were also observed in the cerebella of mouse models expressing full-length mutant ataxin-7 or the androgen receptor. These results demonstrate that some of the gene expression effects of expanded polyglutamine proteins occur independently of protein context.

Increased huntingtin protein length reduces the number of polyglutamine-induced gene expression changes in mouse models of Huntington's disease

Both transcriptional dysregulation and proteolysis of mutant huntingtin (htt) are postulated to be important components of Huntington's disease (HD) pathogenesis. In previous studies, we demonstrated that transgenic mice that express short mutant htt fragments containing 171 or fewer N-terminal residues (R6/2 and N171-82Q mice) recapitulate many of the mRNA changes observed in human HD brain. To examine whether htt protein length influences the ability of its expanded polyglutamine domain to alter gene expression, we conducted mRNA profiling analyses of mice that express an extended N-terminal fragment (HD46, HD100; 964 amino acids) or full-length (YAC72; 3144 amino acids) mutant htt transprotein. Oligonucleotide microarray analyses of HD46 and YAC72 mice identified fewer differentially expressed mRNAs than were seen in transgenic mice expressing short N-terminal mutant htt fragments. Histologic analyses also detected limited changes in these mice (small decreases in adenosine A2a receptor mRNA and dopamine D2 receptor binding in HD100 animals; small increases in dopamine D1 receptor binding in HD46 and HD100 mice). Neither HD46 nor YAC72 mice exhibited altered mRNA levels similar to those observed previously in R6/2 mice, N171-82Q mice or human HD patients. These findings suggest that htt protein length influences the ability of an expanded polyglutamine domain to alter gene expression. Furthermore, our findings suggest that short N-terminal fragments of mutant htt might be responsible for the gene expression alterations observed in human HD brain.

Altered transcriptional regulation in cells expressing the expanded polyglutamine androgen receptor

Kennedy's disease is a degenerative disease of motor neurons in which the causative mutation is expansion of a CAG/polyglutamine tract near the 5' end of the androgen receptor gene. The mutant protein misfolds, aggregates, and interacts abnormally with other proteins, leading to a novel, toxic gain of function and an alteration of normal function. We used a cell culture model to explore the mechanisms underlying the alterations in androgen receptor function conferred by the mutation. Here we show that cells expressing the wild-type androgen receptor with 24 CAG repeats respond to ligand by showing trophic effects including prolonged survival in low serum, whereas cells expressing the mutant receptor with 65 CAG repeats do not show a robust trophic response. This partial loss of function correlates with decreased levels of the mutant protein due to its preferential degradation by the ubiquitin-proteasome pathway. Expression analysis using oligonucleotide arrays confirms that the mutant receptor has undergone a partial loss of function, and fails to regulate a subset of genes whose expression is normally affected by ligand activation of the wild-type receptor. The mutant receptor has also undergone several functionally important post-translational modifications in the absence of ligand that the wild-type receptor undergoes in the presence of ligand, including acetylation and phosphorylation. These modifications correlate with a ligand-independent gain of function exhibited by the mutant receptor in expression analysis. Our findings suggest that polyglutamine expansion alters androgen receptor function by promoting its degradation and by modifying its activity as a transcription factor.

A regression-based method to identify differentially expressed genes in microarray time course studies and its application in an inducible Huntington's disease transgenic model

Time-course studies with microarray technologies provide enormous potential for exploring underlying mechanisms of biological phenomena in many areas of biomedical research, but the large amount of gene expression data generated by such studies also presents great challenges to data analysis. Here we introduce a regression-based statistical modeling approach that identifies differentially expressed genes in microarray time-course studies. To illustrate this method, we applied it to data generated from an inducible Huntington's disease transgenic model. The regression method accounts for the induction process, incorporates relevant experimental information, and includes parameters that specifically address the research interest: the temporal differences in gene expression profiles between the mutant and control mice over the time course, in addition to heterogeneities that commonly exist in microarray data. Least-squares and estimating equation techniques were used to estimate parameters and variances, and inferences were made based on efficient and robust Z-statistics under a set of well-defined assumptions. A permutation test was also used to estimate the number of false-positives, providing an alternative measurement of statistical significance useful for investigators to make decisions on follow-up studies.

Dysregulation of gene expression in the R6/2 model of polyglutamine disease: parallel changes in muscle and brain

Previous analyses of gene expression in a mouse model of Huntington's disease (R6/2) indicated that an N-terminal fragment of mutant huntingtin causes downregulation of striatal signaling genes and particularly those normally induced by cAMP and retinoic acid. The present study expands the regional and temporal scope of this previous work by assessing whether similar changes occur in other brain regions affected in Huntington's disease and other polyglutamine diseases and by discerning whether gene expression changes precede the appearance of disease signs. Oligonucleotide microarrays were employed to survey the expression of approximately 11,000 mRNAs in the cerebral cortex, cerebellum and striatum of symptomatic R6/2 mice. The number and nature of gene expression changes were similar among these three regions, influenced as expected by regional differences in baseline gene expression. Time-course studies revealed that mRNA changes could only reliably be detected after 4 weeks of age, coincident with development of early pathologic and behavioral changes in these animals. In addition, we discovered that skeletal muscle is also a target of polyglutamine-related perturbations in gene expression, showing changes in mRNAs that are dysregulated in brain and also muscle-specific mRNAs. The complete dataset is available at www.neumetrix.info.

Early transcriptional profiles in huntingtin-inducible striatal cells by microarray analyses

Gene expression studies conducted with mouse models of Huntington's disease (HD) have revealed profound modifications in gene transcription. However, the complexity of in vivo tissue hampers definition of very early transcriptional modifications and does not allow discrimination between cell-autonomous changes and those resulting from intercellular activity processes. To identify early, cell-autonomous transcriptional changes, we compared gene expression profiles of clonal striata-derived cells expressing different N-terminal 548-amino-acid huntingtin fragments (with 26, 67, 105 or 118 glutamines) under the control of a doxycycline-regulated promoter. In these cells, mutant huntingtin did not form aggregates or cause cell death; therefore, the gene expression profiles report transcriptional changes reflecting early pathogenic events. We found that genes involved in cell signaling, transcription, lipid metabolism and vesicle trafficking were affected, in some cases, within 12 hours of mutant protein induction. Interestingly, this study revealed differential expression of a number of genes involved in cholesterol and fatty acid metabolism, suggesting that these metabolic pathways may play a role in HD pathogenesis.

Strain-driven mesoscopic reconstruction of the As/Ge(111) surface

Periodic arrays of large hexagonal tiles (up to 170 A in size) are observed on As/Ge(111) surfaces. First-principles total energy calculations combined with scanning tunneling microscopy reveal a (5-7-5)-ringed structure for the trenches that separate the tiles. We find that trenches form via an exothermic process. The calculated equilibrium trench spacing of approximately 104 A agrees with experiment. Comparison between first-principles calculations and continuum elasticity theory suggests that the observed mesoscopic reconstruction is driven entirely by long-range surface strain relaxation.

Locating the genes underlying a simulated complex disease by discriminant analysis

The purpose of this study was to propose and evaluate a novel multivariate approach for genetic mapping of complex binary human diseases. This approach uses the application of either of two methods of standard (stepwise) discriminant analysis to detect linkage based on the differential marker identity-by-descent distributions among the three affection groups of sib pairs (concordantly affected, discordant, and concordantly unaffected). One of the advantages of this approach is that it allows for simultaneously testing all markers, as well as other genetic and environmental factors, in a single multivariate setting. We have explored its properties and behaviors via an application to the simulated data in Genetic Analysis Workshop 12.

Comparison of three methods for obtaining principal components from family data in genetic analysis of complex disease

Three multivariate techniques used to derive principal components (PCs) from family data were compared for their ability to model family data and power to detect linkage. Using the simulated data from Genetic Analysis Workshop 12, the five quantitative traits were first adjusted for age, sex, and environmental factors 1 and 2. Then, standard PCs, PCs obtained from between-family covariance, and PCs obtained from within-family genetic covariance were derived and subjected to multivariate sib pair linkage analysis. The standard PCs obtained from the overall correlation matrix allowed identification of key features of the true genetic model more readily than did the other methods. For detection of linkage, standard PCs and PCs obtained from the between-family genetic covariance performed similarly in terms of both power and type 1 error, and both methods performed better than the PCs obtained from within-family genetic covariance.

Prediction of central nervous system embryonal tumour outcome based on gene expression

Embryonal tumours of the central nervous system (CNS) represent a heterogeneous group of tumours about which little is known biologically, and whose diagnosis, on the basis of morphologic appearance alone, is controversial. Medulloblastomas, for example, are the most common malignant brain tumour of childhood, but their pathogenesis is unknown, their relationship to other embryonal CNS tumours is debated, and patients' response to therapy is difficult to predict. We approached these problems by developing a classification system based on DNA microarray gene expression data derived from 99 patient samples. Here we demonstrate that medulloblastomas are molecularly distinct from other brain tumours including primitive neuroectodermal tumours (PNETs), atypical teratoid/rhabdoid tumours (AT/RTs) and malignant gliomas. Previously unrecognized evidence supporting the derivation of medulloblastomas from cerebellar granule cells through activation of the Sonic Hedgehog (SHH) pathway was also revealed. We show further that the clinical outcome of children with medulloblastomas is highly predictable on the basis of the gene expression profiles of their tumours at diagnosis.

Linkage analysis of human systemic lupus erythematosus-related traits: a principal component approach

OBJECTIVE

To identify chromosomal regions containing genes involved in the susceptibility to human systemic lupus erythematosus (SLE)-related traits.

METHODS

In the context of a genome scan, we analyzed 101 SLE-affected sibpairs with respect to dermatologic, renal, immunologic, hematologic, neurologic, cardiopulmonary, and arthritic characteristics. Phenotypes were redefined in terms of principal components, which are synthetic variables composed of linear combinations of the original traits. Using 9 principal components obtained from these 7 traits plus age at SLE onset and race, we analyzed genome scan data with the multivariate version of the new Haseman-Elston regression model.

RESULTS

The largest linkage for an individual trait was on chromosome 2 at 228 cM (immunologic; P = 0.00048). The most significant linkage to an individual principal component was on chromosome 4 at 208 cM (P = 0.00007). The largest multivariate linkage was on chromosome 7 at 69 cM (P = 0.0001). Of the individual organ systems, dermatologic involvement had the largest effect (P = 0.0083) at this peak at 7p13 on chromosome 7. Further analyses revealed that malar rash, a subtype of dermatologic involvement, was linked significantly (P = 0.00458) to this location.

CONCLUSION

These results provide evidence of the presence and locations of genes that are involved in the genetic susceptibility to SLE-related traits in humans.

Altered transcription in yeast expressing expanded polyglutamine

Expanded polyglutamine tracts are responsible for at least eight fatal neurodegenerative diseases. In mouse models, proteins with expanded polyglutamine cause transcriptional dysregulation before onset of symptoms, suggesting that this dysregulation may be an early event in polyglutamine pathogenesis. Transcriptional dysregulation and cellular toxicity may be due to interaction between expanded polyglutamine and the histone acetyltransferase CREB-binding protein. To determine whether polyglutamine-mediated transcriptional dysregulation occurs in yeast, we expressed polyglutamine tracts in Saccharomyces cerevisiae. Gene expression profiles were determined for strains expressing either a cytoplasmic or nuclear protein with 23 or 75 glutamines, and these profiles were compared to existing profiles of mutant yeast strains. Transcriptional induction of genes encoding chaperones and heat-shock factors was caused by expression of expanded polyglutamine in either the nucleus or cytoplasm. Transcriptional repression was most prominent in yeast expressing nuclear expanded polyglutamine and was similar to profiles of yeast strains deleted for components of the histone acetyltransferase complex Spt/Ada/Gcn5 acetyltransferase (SAGA). The promoter from one affected gene (PHO84) was repressed by expanded polyglutamine in a reporter gene assay, and this effect was mitigated by the histone deacetylase inhibitor, Trichostatin A. Consistent with an effect on SAGA, nuclear expanded polyglutamine enhanced the toxicity of a deletion in the SAGA component SPT3. Thus, an early component of polyglutamine toxicity, transcriptional dysregulation, is conserved in yeast and is pharmacologically antagonized by a histone deacetylase inhibitor. These results suggest a therapeutic approach for treatment of polyglutamine diseases and provide the potential for yeast-based screens for agents that reverse polyglutamine toxicity.

Steps on As-terminated Ge(001) revisited: theory versus experiment

We examine the double-layer B-type steps on As-terminated vicinal Ge(001) surfaces. The currently accepted structure is a chemically inert bulklike structure without any gap state, and with all the chemical bonds of the Ge and As atoms being satisfied. However, we show that the need for optimizing the p(3) pyramidal angles of the threefold coordinated As atoms drives unusual atomic rearrangement. This leads to a more stable reconstruction involving odd-membered (5-7-5) rings at the step edge. Comparison between theoretical and experimental scanning tunneling microscopy images yields excellent agreement.

The amyloid precursor protein locus and very-late-onset Alzheimer disease

Although mutations in the amyloid-beta precursor protein (APP) gene are known to confer high risk of Alzheimer disease (AD) to a small percentage of families in which it has early onset, convincing evidence of a major role for the APP locus in late-onset AD has not been forthcoming. In this report, we have used a covariate-based affected-sib-pair linkage method to analyze the chromosome 21 clinical and genetic data obtained on affected sibships by the National Institute of Mental Health Alzheimer Disease Genetics Initiative. The baseline model (without covariates) gave a LOD score of 0.02, which increases to 1.43 when covariates representing the additive effects of E2 and E4 are added. Larger increases in LOD scores were found when age at last examination/death (LOD score 5.54; P=.000002) or age at onset plus disease duration (LOD score 5.63; P=.000006) were included in the linkage model. We conclude that the APP locus may predispose to AD in the very elderly.

Ontogeny of catabolic and morphological properties of skeletal muscle of the red-winged blackbird (Agelaius phoeniceus)

Thermogenic capabilities of red-winged blackbirds improve markedly during their 10-12-day nestling period, especially between day 5 and day 8. The time course of improvements may be determined by the maturation of skeletal muscles involved in shivering thermogenesis, particularly the pectoralis muscles. To test this hypothesis, morphological and biochemical changes in pectoral and leg muscles were measured in young and adult blackbirds. Both muscles grew disproportionately relative to body mass. The pectoralis consisted entirely of fast-twitch fibers, predominantly fast oxidative glycolytic. In contrast, the gastrocnemius muscle consisted of a mixture of slow and fast fibers (predominantly fast glycolytic). Although fiber composition was constant, both cross-sectional area and density of fibers increased with age in both muscles. Catabolic capacities of the pectoralis increased significantly (approximately 7-8-fold) throughout the nestling period, most abruptly after day 3 (citrate synthase, CS) or day 4 (3-hydroxacyl-CoA-dehydrogenase, HOAD). Myofibrillar ATPase activities in the pectoralis were initially low, but increased after day 5. Further increases in CS and myofibrillar ATPase activities occurred in the pectoralis after fledging. CS and HOAD activities in the leg were much lower, but myofibrillar ATPase activities were remarkably similar in the two muscles, differing only in adults. These results are consistent with the hypothesis that the development of endothermy is dependent on the morphological and biochemical maturation of skeletal muscles important in thermogenesis.

No excess of early onset cancer in family members of Wilms tumor patients

BACKGROUND

Wilms tumor is one of the few pediatric cancers with well-defined familial and genetic components. The authors assessed the risk of early-onset cancers in first- and second-degree relatives of patients enrolled by the National Wilms Tumor Study Group.

METHODS

Using a stratified sampling scheme that targeted 530 families of patients who were believed a priori to have a genetic contribution to their disease, the authors conducted interviews regarding cancer occurrence in 4258 family members from 296 families of patients with Wilms tumor. Reports of malignant neoplasms that occurred before 55 years of age were confirmed by review of medical records wherever possible. A period of risk was defined for each family member based on calendar time and his or her relationship to the proband.

RESULTS

Ninety-nine cancers were observed, whereas 126.8 were expected by applying standard cancer rates for age and calendar period to the 120,885 person-years at risk. The standardized incidence ratio (SIR) was O-E = 0.78 with 95% confidence interval (CI) of (0.64, 0.95). In subgroup analyses, the highest relative risks were observed for parents of the index case (O/E = 21/13.0 = 1.6, 95% CI = 1.0, 2.5) and for leukemia (O/E = 9/4.9 = 1.9, 95% CI= 0.85,3.5).

CONCLUSIONS

The results of this study may provide reassurance to families of children who have had Wilms tumor. Potential sources of bias included the low (56%) rate of participation of targeted families. In general, the biases might have led to the underreporting of some cancers, especially in more distant relatives. The possibility of a slight excess of cancer in parents of Wilms tumor patients could not be excluded.

An efficient and robust statistical modeling approach to discover differentially expressed genes using genomic expression profiles

We have developed a statistical regression modeling approach to discover genes that are differentially expressed between two predefined sample groups in DNA microarray experiments. Our model is based on well-defined assumptions, uses rigorous and well-characterized statistical measures, and accounts for the heterogeneity and genomic complexity of the data. In contrast to cluster analysis, which attempts to define groups of genes and/or samples that share common overall expression profiles, our modeling approach uses known sample group membership to focus on expression profiles of individual genes in a sensitive and robust manner. Further, this approach can be used to test statistical hypotheses about gene expression. To demonstrate this methodology, we compared the expression profiles of 11 acute myeloid leukemia (AML) and 27 acute lymphoblastic leukemia (ALL) samples from a previous study (Golub et al. 1999) and found 141 genes differentially expressed between AML and ALL with a 1% significance at the genomic level. Using this modeling approach to compare different sample groups within the AML samples, we identified a group of genes whose expression profiles correlated with that of thrombopoietin and found that genes whose expression associated with AML treatment outcome lie in recurrent chromosomal locations. Our results are compared with those obtained using t-tests or Wilcoxon rank sum statistics.

NeuroD2 is necessary for development and survival of central nervous system neurons

NeuroD2 is sufficient to induce cell cycle arrest and neurogenic differentiation in nonneuronal cells. To determine whether this bHLH transcription factor was necessary for normal brain development, we used homologous recombination to replace the neuroD2 coding region with a beta-galactosidase reporter gene. The neuroD2 gene expressed the reporter in a subset of neurons in the central nervous system, including in neurons of the neocortex and hippocampus and cerebellum. NeuroD2(-/-) mice showed normal development until about day P14, when they began exhibiting ataxia and failure to thrive. Brain areas that expressed neuroD2 were smaller than normal and showed higher rates of apoptosis. Cerebella of neuroD2-null mice expressed reduced levels of genes encoding proteins that support cerebellar granule cell survival, including brain-derived neurotrophic factor (BDNF). Decreased levels of BDNF and higher rates of apoptosis in cerebellar granule cells of neuroD2(-/-) mice indicate that neuroD2 is necessary for the survival of specific populations of central nervous system neurons in addition to its known effects on cell cycle regulation and neuronal differentiation.