Human Molecular Genetics Advance Access originally published online on December 8, 2004
Human Molecular Genetics 2005 14(3):373-384; doi:10.1093/hmg/ddi033
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Human Molecular Genetics, Vol. 14, No. 3 © Oxford University Press 2005; all rights reserved
The cerebellar transcriptome during postnatal development of the Ts1Cje mouse, a segmental trisomy model for Down syndrome
1Unité Mixte de Recherche 7637, Centre National de la Recherche Scientifique, 2Equipe de Statistique Appliquée, Ecole Supérieure de Physique et de Chimie Industrielles, 10 rue Vauquelin, 75005 Paris, France, 3Institut National de la Santé Et de la Recherche Médicale Unité 549, Institut Paul Broca, 2ter, rue d'Alésia, 75014 Paris, France, 4Department of Genetic Medicine and Development, University of Geneva Medical School, 1 rue Michel Servet, 1211 Geneva, Switzerland and 5Department of Pediatrics, UCSF, San Francisco, CA 94143-0748, USA
* To whom correspondence should be addressed. Email: marie-claude.potier{at}espci.fr
Received September 16, 2004; Revised November 15, 2004; Accepted November 26, 2004
| ABSTRACT |
|---|
|
|
|---|
The central nervous system of persons with Down syndrome presents cytoarchitectural abnormalities that likely result from gene-dosage effects affecting the expression of key developmental genes. To test this hypothesis, we have investigated the transcriptome of the cerebellum of the Ts1Cje mouse model of Down syndrome during postnatal development using microarrays and quantitative PCR (qPCR). Genes present in three copies were consistently overexpressed, with a mean ratio relative to euploid of 1.52 as determined by qPCR. Out of 63 three-copy genes tested, only five, nine and seven genes had ratios >2 or <1.2 at postnatal days 0 (P0), P15 and P30, respectively. This gene-dosage effect was associated with a dysregulation of the expression of some two-copy genes. Out of 8258 genes examined, the Ts1Cje/euploid ratios differed significantly from 1.0 for 406 (80 and 154 with ratios above 1.5 and below 0.7, respectively), 333 (11 above 1.5 and 55 below 0.7) and 246 genes (59 above 1.5 and 69 below 0.7) at P0, P15 and P30, respectively. Among the two-copy genes differentially expressed in the trisomic cerebellum, six homeobox genes, two belonging to the Notch pathway, were severely repressed. Overall, at P0, transcripts involved in cell differentiation and development were over-represented among the dysregulated genes, suggesting that cell differentiation and migration might be more altered than cell proliferation. Finally, global gene profiling revealed that transcription in Ts1Cje mice is more affected by the developmental changes than by the trisomic state, and that there is no apparent detectable delay in the postnatal development of the cerebellum of Ts1Cje mice.
| INTRODUCTION |
|---|
|
|
|---|
Down syndrome results from trisomy of human chromosome 21 and is the most frequent genetic cause of mental retardation, occurring in
1 in 800 newborns. Most cases are full trisomies, and about 300 genes are triplicated (1
Because of the prominence of mental retardation, brain has been the subject of particular interest in research on Down syndrome, but a systematic study of gene expression in brain during development, particularly at postnatal stages, is only possible in animal models. Mouse models for Down syndrome are trisomic for single gene, for several genes or for a large segment corresponding to the distal part of chromosome 16, orthologous to a large portion of human chromosome 21 (8
). Two models of segmental trisomy 16 have been generated, the Ts65Dn mice and the Ts1Cje mice (9
11
). These models have revealed previously unknown phenotypes that may be relevant to Down syndrome, such as a substantial loss of glutamatergic granule cells in the internal layer of the cerebellum in adult (12
,13
).
We studied the transcriptome during postnatal development of the cerebellum for two reasons. First, there are massive developmental changes during the first 20 days after birth. During this period, granule cells, which represent
40% of the total number of neurons of the mature cerebellum, proliferate, migrate and differentiate from the external to the internal layer, and Purkinje cells develop their dense dendritic trees making connections with other cells. Secondly, gene profiling during postnatal development of the cerebellum in mice is well documented (14
,15
), Furthermore, previously investigated mutants of the cerebellum could serve as a reference to our study (14
).
The Ts1Cje Down syndrome model was used in the present study after backcrossing the mice onto a pure genetic background to reduce variability in gene expression. Ts1Cje mice carry a segmental duplication of the syntenic region orthologous to human chromosome 21 from Sod1 to Znf295, including about 95 genes (8
,10
,11
). Differential gene expression was studied at postnatal days 0 (P0), P15 and P30 in the cerebellum using Affymetrix U74Av2 microarrays. Expression of genes was also measured in parallel by quantitative PCR (qPCR).
Two groups of genes have been characterized: one that contains the triplicated genes and shows a consistent gene-dosage effect (approximately 1.5 overexpression ratio in Ts1Cje mice when compared with euploid mice) during postnatal development, and the second that contains euploid genes that are significantly differentially expressed (over or under) and are mainly involved in cell differentiation and development.
| RESULTS |
|---|
|
|
|---|
Experimental data and analysis: validation
Using Affymetrix murine genome U74Av2 microarrays, we measured the expression levels of 12 488 genes (half known genes and half ESTs) in the cerebellum of Ts1Cje and euploid mice each at P0, P15 and P30 during postnatal development. The mean expression value for each gene was calculated from duplicate experiments performed on the six groups of three animals each. According to the Affymetrix algorithm, 8287 genes were found to be expressed at least at one time point (see Materials and Methods for details on the filtering step and Supplementary Material, Table S1 for a complete list of the 8287 genes). Two statistical analyses were applied to the data: variability analysis of microarray data (VARAN), a freely available web server performing a signal intensity-based analysis of the log 2 expression ratios variability deduced from DNA microarray raw data (16
i and ßi parameters, respectively, in Table 1).
|
To validate our experimental data and analysis, we compared genes differentially expressed between P0 and P15 or between P15 and P30 (
i parameters) in either euploid mice or Ts1Cje mice with genes differentially expressed along normal postnatal development of the cerebellum as reported by Diaz et al. (14
i parameters). Many genes differentially expressed during postnatal development of the cerebellum, such as GABA-
6 and GABA-
, myelin oligodendrocyte glycoprotein and parvalbumin (Supplemental Material, Table S2), belonged to the cell differentiation class B of Diaz et al. (14
Gene expression in the Ts1Cje cerebellum: triplicated versus euploid genes from mouse chromosome 16
Of the 8287 expressed genes, 145 were from mouse chromosome 16, and among them 29 were triplicated in Ts1Cje mice. At P0, P15 and P30, mean expression ratios of triplicated genes between Ts1Cje and controls were 1.66, 1.32 and 1.32 at P0, P15 and P30, respectively, whereas with euploid genes the ratios were 1.08, 1.12 and 1.02 at P0, P15 and P30, respectively (Fig. 1). After ANOVA, 12, six and nine triplicated genes were found to be statistically overexpressed (
=1%) in Ts1Cje mice when compared with euploid mice at P0, P15 and P30, respectively (Table 2). Among these genes, four (Dscr3, HMGP14, Donson, C21orf4) were overexpressed at all three time points.
|
|
qPCR on 78 genes from chromosome 16 (15 euploid and 63 triplicated) showed overexpression of the triplicated genes with mean ratios of 1.56, 1.43 and 1.55 at P0, P15 and P30, respectively, whereas the euploid genes from chromosome 16 had ratios very close to 1.0 (1.04, 1.00 and 1.01 at P0, P15 and P30, respectively) (Fig. 2; Supplementary Material, Table S3). The 18 control genes mapping to chromosomes other than 16 and described previously (17
|
qPCR and microarray measurements were highly correlated (P<0.01) for 45 transcripts tested that were present in both microarray and qPCR studies. A single and expected exception was the Sod1 transcript. The wild-type Sod1 gene is composed of five exons, but in Ts1Cje mice a neomycin cassette is inserted into exon 3, thereby generating a shorter transcript with only exons 1, 2 and 5 (18
Gene expression in the Ts1Cje cerebellum: euploid genes from all chromosomes (statistical analysis and validation by qPCR)
In addition to the triplicated genes from chromosome 16 that were found to be consistently overexpressed in Ts1Cje mice when compared with euploid mice, euploid genes from chromosomes other than 16 were also found to be differentially expressed at P0, P15 and P30. A first statistical analysis of the data using VARAN selected a total of 95 genes with ratios between Ts1Cje and euploid ranging from 0.1 to 0.8 and from 1.2 to 70.7, respectively. Five genes were from the triplicated segment of chromosome 16, and 90 were euploid genes. Table 3 displays the number of genes differentially expressed deduced from the four independent analysis at each time point (see Data Analysis in Materials and Methods). Overlaps between the four analyses gave the final number of differentially expressed genes: 13, 57 and 25 at P0, P15 and P30, respectively (Table 3). These genes are listed in Table 4. No differentially expressed genes were in common for all three developmental stages studied. However, three genes were in common between P0 and P15 and three other genes between P15 and P30 (bold in Table 4).
|
|
A second statistical analysis, ANOVA, was applied to the data. Table 5 gives the number of genes that are differentially expressed (over or under) and the test of the nullity of each parameter for several values of the risk of type I error (1, 0.1 and 0.01%). For a 1% risk, the number of genes differentially expressed is much higher between the stages of development (24.6 and 27.5% between P0 and P15 and between P15 and P30, respectively) than between the trisomic versus euploid state (5.1, 4.1 and 3.1% at P0, P15 and P30, respectively). A complete list of genes is included in Supplementary Material, Table S4. All genes differentially expressed from VARAN were also detected by ANOVA and corresponded to the short list in Table 4, indicating that analysis with VARAN is more stringent than with ANOVA. At P0, 86 and 154 genes had ratios above 1.5 and below 0.7, respectively. At P15, more genes had ratios below 0.7 than above 1.5 (154 and 86 genes, respectively). At P30, 63 and 69 genes had ratios above 1.5 and below 0.7, respectively. Thirty-one genes were in common between all three stages (Supplementary Material, Table S4). Figure 3 represents gene expression ratios between Ts1Cje and euploid mice with respect to the value of their statistic (Fisher), and shows that variations are more significant (higher Fisher) at P0 and P15 than at P30. Differentially expressed genes are either up- or down-regulated, and nearly all triplicated genes belong to the up-regulated genes.
|
|
Classification of the differentially expressed genes from ANOVA according to gene ontology showed a significant enrichment only at P0 for GO categories GO : 0007275 (P=0.00079) and GO : 0030154 (P=0.00051), corresponding to genes involved in cell development and differentiation, respectively. Among the 419 differentially expressed genes at P0 (Table 5, Ts1Cje/euploid at P0 with
=1%), there were 297 genes with annotated ontology of which six and 17 were involved in cell differentiation and development, respectively. To confirm the reliability of microarray data, seven differentially expressed genes from the short list of 95 genes selected by VARAN and ANOVA were analyzed by qPCR at P0, P15 or P30 (Table 6). From a total of eight experiments conducted across three developmental stages, six variations were confirmed. Differential expression of Glutamine synthase was found to be different when measured by qPCR compared with microarray, thus suggesting a lack of specificity of the microarray probe set.
|
Principal component analysis
For gene profiling across all conditions, trisomic and euploid at P0, P15 and P30, we performed PCA. With the 8287 genes detected on the Affymetrix U74Av2 microarray, the projection onto the subspace spanned by the two first principal components clearly grouped the animals according to their stage of development, thus demonstrating that the effect of development on gene expression is more important than the trisomic effect (Fig. 4A). PCA was then performed on subsets of genes of increasing size, the genes being ranked according to their discriminatory power for trisomy using Student statistic on data from either euploid mice or trisomic mice whatever their stage of development (see Materials and Methods). For the first 100 genes, PCA grouped the animals according to the trisomic versus euploid state (Fig. 4D). However on larger subsets (3000 genes), animals were grouped according to their developmental stage (Fig. 4B). With 3000 genes, within each stage of development, Ts1Cje animals could be distinguished from the euploids. The use of 2000 genes resulted in an intermediate situation where the most likely separation was made according to the trisomic versus euploid state, although one could also still get a separation according to the stage of development (Fig. 4C).
|
Comparison of microarray data from the literature
Our data set was compared to those of two microarray studies published recently: that of Saran et al. (6
|
Using the same statistical analysis of microarray data (VARAN), three genes/ESTs from the euploid genomic region were found to be differentially expressed in the cerebellum of both Ts1Cje mice at P30 (our study) and Ts65Dn mice at 34 months (6
| DISCUSSION |
|---|
|
|
|---|
Gene-dosage effect of the triplicated genes from mouse chromosome 16 (MMU16) during postnatal development of the cerebellum in Ts1Cje mice
Study of the cerebellar transcriptome during postnatal development in the Ts1Cje model of Down syndrome clearly shows an overexpression of the three copy-genes, with mean ratios over postnatal development of 1.43 and 1.52 as determined by microarray and qPCR experiments, respectively. These results are in agreement with other studies published recently (6
Development has a greater effect than trisomy on global gene expression
Global gene expression profiling during postnatal development of the cerebellum enables mice to be classified according to their stage of development rather than to their trisomic versus eusomic status. This indicates that development has more impact on the expression profile than does gene-dosage (Fig. 4A). Furthermore, gene expression profiles of the Ts1Cje cerebellum at P15 and P30 were not closer to the euploid animals at P0 and P15, respectively, thus indicating that there was no developmental delay in Ts1Cje mice (Fig. 4A). Figure 3 showed that there were more genes differentially expressed with significant statistical values at P0 and P15 than at P30. In addition, qPCR on triplicated genes indicated that there was a highly significant correlation (P<0.01) between the P15 and P30 time points, but not between P0 and P15 or P30, suggesting that the developmental changes occurred between P0 and P15, and then stabilized.
Gene expression dysregulation of euploid genes
Statistical analysis revealed a group of dysregulated genes with ratios ranging from 0.1 to 0.8 and from 1.2 to 70.7 (short list in Table 4). The analysis was validated by qPCR on a group of seven genes, and the differential expression of six were confirmed. Comparative analysis of microarray data [our study on Ts1Cje and study of Saran et al. (6
) on the cerebellum of adult Ts65Dn mice] revealed euploid genes that are severely differentially expressed (up or down) with large statistical values (Fig. 5A and C). In the transcriptome study on whole brain from Ts1Cje at P0, most differentially expressed genes were from the triplicated genomic segment, with ratios below 2 (7
) (Fig. 5B). Because whole brain is a mixture of many different cell types, including cerebellar cells, cell-specific regulation is diluted and below the level of detection. Reducing cellular complexity by the use of cerebellum increases the chance of seeing genes that are differentially expressed in a particular cell type. Using even smaller brain structures may reveal differentially expressed genes that are specific for one region of the brain. The higher number of dysregulated genes found in our study (5.1, 4.1 and 3.1% at P0, P15 and P30, respectively, Table 5) and in the study of Saran et al. (6
) in adult Ts65Dn cerebellum (5.35%) as compared to 0.95% in the study of Amano et al. (7
) on whole brain of Ts1Cje at P0 could relate to the complexity of the tissue used. In a similar way, pooling the same tissue (cerebellum) at different time points (P0, P15 and P30) decreased the number of differentially expressed genes to 0.61% (Fig. 5A). From VARAN analysis, two euploid genes were selected as differentially expressed both from our study (Ts1Cje at P30) and from the data of Saran et al. (6
) (Ts65Dn at 34 months). One, corresponding to hemoglobin ß1 chain, was repressed, and the other, corresponding to erythroid differentiation regulator, was overexpressed. Supplementary Material, Table S4 shows that several hemoglobin
- and ß-subunits were highly repressed at P30 and P15, thus suggesting a decrease in the amount of meninges that surround the cerebellum. Additional experiments done on animals bred on the same genetic background and at the same time points will be necessary to draw a final conclusion on particular genes.
Pathways modified during postnatal development of the cerebellum in Ts1Cje mice: cell differentiation and development
We used gene ontology to categorize the differentially expressed genes from our microarray experiments and found that at P0 there was an enrichment of genes involved in cell differentiation and development. This suggests that cell differentiation and migration of cells might be more affected by trisomy than cell proliferation. At P0, of the 17 genes known to be involved in development, six are homeobox genes, all of which were repressed. In addition, at P15 and P30, respectively, 12 and one homeobox genes were found to be differentially expressed. Homeoproteins are involved in the early patterning of the nervous system and possibly at later stages of neuronal differentiation (23
). Among them, Homeo box A5 was found to be severely down-regulated (to 0.12) at P0 in our study, as well as in two pools of cerebellum from adult Ts65Dn mice (0.47 and 0.54) (6
). In the study of Amano et al. (7
) on whole brain of Ts1Cje mice, Homeo box A5 was not detected. Homeo box A5 is mostly expressed during the neonatal period, but there are studies relating to its expression in the adult, particularly in Purkinje cells, where it activates the transcription of Purkinje cell protein 2 (24
).
At P15, Dlx1 was found to be overexpressed in Ts1Cje mice when compared with euploid mice (Table 6: 80-fold and 70.8-fold by microarray and qPCR, respectively). This very high ratio probably indicates that Dlx1 is expressed only in trisomic cerebellum. Dlx1 is not known to be involved in postnatal development of the cerebellum. However, its role in the lineage of a particular type of GABAergic interneurons in human neocortex has been demonstrated recently (25
). In addition, an increase of Dlx1 and Dlx2 negatively regulates the Notch signalling to specify a later subset of neuronal progenitors and promote their terminal differentiation (26
). Dlx2 was also found to be increased at P15 (7-fold; Supplementary Material, Table S4). It is therefore likely that increase of Dlx1 and Dlx2 may be associated with cell differentiation. We cannot exclude significant changes in the expression of transcription factors important for development of the cerebellum, such as Math1, Mash1 and Shh, that have been described in microarray studies (14
), as they were not present on the Affymetrix U74Av2 microarray. More data on these particular genes will be useful for understanding cerebellum development, as Math1 has recently been shown to be involved in the control of cerebellar granule cell differentiation by regulating multiple components of the Notch signalling pathway (27
). Study of gene expression of Notch receptors and ligands and transcription factor targets will confirm whether the Notch pathway is altered in Down syndrome.
Gene profiling from the study of Diaz et al. (14
) on cerebellum mutants (weaver and lurcher) were compared to gene profiling in Ts1Cje mice using the same statistical analysis (ANOVA with the parameterization described in Materials and Methods). In weaver mice, granule cells are severely reduced in the internal layer of the adult cerebellum, whereas in lurcher mice, Purkinje cells are reduced (28
,29
). Comparison of Ts1Cje, weaver and lurcher revealed more similarity between Ts1cje and weaver mutants than between Ts1Cje and lurcher mutants (unpublished data). This is consistent with morphological data (13
) suggesting that granule cells, rather than Purkinje cells, are altered in Ts1Cje mice.
The list of genes differentially expressed also includes transthyretin, a transporter of thyroid hormone and vitamin A, which has been shown to be a carrier of amyloid beta peptide in the cerebrospinal fluid. Transthyretin prevents formation of amyloid fibrils and Apo-E-induced accumulation of amyloid beta (30
,31
). In Ts1Cje mice, transthyretin was repressed on microarrays and by qPCR at P0 (Table 6). Repression of transthyretin, if present in Down syndrome, might contribute to the amyloid beta deposition process in brain in addition to APP overexpression.
Another gene expression study in Down syndrome has shown that REST (neuron-specific silencer factor) and its targets (SCG10, L1, synapsin) are down-regulated in neuronal precursor cells from Down syndrome fetal cortex (32
). From our study, expression of SCG10, synapsin and L1, but not of REST, was detected during postnatal development of the cerebellum. Only L1 was found to be differentially expressed, with a ratio of 0.52 between Ts1Cje and euploid at P0, although with low statistical value (Fisher=4.6 under the 1% threshold). Our data do not, therefore, provide evidence for or against the hypothesis of down-regulation of REST pathway in neuronal precursor cells proposed by Bahn et al. (32
).
In conclusion, in the Ts1Cje mouse model of Down syndrome, the presence of three copies of
85 genes having orthologues on human chromosome 21 induces a consistent change in expression of these triplicated genes throughout postnatal development. In addition to this general gene-dosage effect, where the expected 1.5 overexpression ratio is largely but not totally respected, a number of euploid genes are differentially expressed in a time-dependant manner, with an over-representation of transcripts involved in cell differentiation and development. Finally, gene profiling reveals that development has more impact on the transcriptome than the trisomic/eusomic status, and that there is no detectable delay in the postnatal development of the cerebellum of the segmental trisomy Ts1Cje mice.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Mice
Ts1Cje mice carry a segmental duplication of the MMU16 region from Sod1 to Znf295. Ts1Cje mice were bred on a C57BL/6 background (11 backcrosses). For gene expression analysis, three pairs of euploid (C57BL/6) and trisomic Ts1Cje male sibs from two litters were used for each developmental stage P0, P15 and P30.
Genotyping
Because the trisomic segment of MMU16 has a truncated Sod1 allele with the neomycin resistance cassette, the mice were first typed on genomic DNA from the tail by PCR using Neo primers (10
).
In addition, after individual RNA extraction (discussed subsequently), all animals were tested for the presence of both Sod1 wild-type and truncated alleles by PCR on cDNAs using Sod1 primers: upCAATGTGACTGCTGGAAAGG and lowATCCCAATCACTCCACAGGC. For each mouse, 500 ng of total RNA from cerebellum were reverse-transcribed overnight at 37°C in the presence of 200 U of Superscript II (InVitrogen), 1x first-strand buffer, 100 µM DTT, 5 µM random hexamers oligonucleotides (Pasteur Institute) and 500 µM dNTPs. PCRs were performed using 2.5 U of HotstartTaq® DNA polymerase (Qiagen) with 1x PCR buffer containing 1.5 mM MgCl2, 0.3 µM of each primer and 200 µM dNTPs for 40 cycles of 94°C for 30 s, 58°C for 30 s, 72°C for 40 s. Because one of the P30 euploid mice showed a pattern with an unexpected number of bands, we decided not to use it, and the euploid P30 pool was therefore composed of only two animals instead of three.
RNA extraction and microarray hybridizations
Total RNA was extracted from frozen individual cerebellum and treated with DNase using RNeasy Midi kit (Qiagen). The quality of each RNA sample was then checked using the Agilent 2100 Bioanalyzer (Agilent Technologies, Inc., Palo Alto, CA, USA).
For the three postnatal stages, euploid and trisomic samples were prepared by pooling an equal amount of each individual RNA. An aliquot 20 µg of total RNA from each pool was converted to cDNA and then to biotinylated cRNA, and were hybridized to Affymetrix Murine Genome U74A version 2 microarrays (Affymetrix, Inc., Santa Clara, CA, USA) on the Curie Institute Affymetrix microarray platform (Paris, France), according to the Affymetrix procedures. Each pool was hybridized in duplicate on independent U74Av2 microarrays.
Data analysis: filtering, normalization, VARAN and ANOVA and ontology
After hybridization, signal intensities were calculated using the Affymetrix GeneChip® software MAS 5.0. The software generates a detection P-value and a detection call to decide whether a gene is expressed, and then attributes to each gene a status: present, absent or marginal. For further analysis, we used only the 8287 genes that were defined as Present or Marginal in Ts1Cje or euploid mice at least at one time point. Genes called Absent in all conditions were not used for further analysis.
We have deposited all the raw data (accession number GSE1611) in the NCBI public database Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/).
Filtered data (8287 genes) were submitted to VARAN (http://www.bionet.espci.fr,) for normalization and differential expression analysis (16
). The microarray measurements were grouped into six pairs of replicates: euploids at P0, P15, and P30, and Ts1Cje mice at P0, P15, and P30. The expression levels were corrected in order to constrain the mean ratio of each pair to be equal to one using a Lowess fit on the MA plot of the log-expression levels (16
). In the numerical analysis, the logarithm of the expression levels was used. For each experiment replicates on euploid mice (i.e. at P0, P15 and P30), a reference MA plot was established. Assuming that the variance of the gene expression level was a function of their mean, and that the expression levels were Gaussian, the 1% threshold values for the difference as a function of the mean were estimated.
For each stage, the MA plots of the four possible combinations of euploid and Ts1Cje mice were established, thus providing four analysis at each time point, and the genes whose difference of expression between euploid and trisomic exceeded the 1% threshold values were selected. For the ANOVA, we applied a regression model for gene expression levels. Table 1 indicates the parameters {
i}, which model the effects of development in euploid mice, and {ßi}, which model the effects of trisomy at each development stage. The number n of parameters equals six, and the number N of microarray measurements equals 12.
Assuming a Gaussian distribution of gene expression levels with common variance s2, the sum of squared residuals of the model was
2 distributed: SSR/
2
2(Nn). If trisomy had no effect at a given stage, we assumed the nullity of the corresponding parameter (null hypothesis H0). If H0 is true, (SSR0SSR)(Nn)/SSR
Fisher(1, Nn), where SSR0 corresponded to the sum of squared residuals of the associated submodel. Hence the test rejects H0 if the Fisher statistic (SSR0SSR) (Nn)/SSR>f
(1, Nn), where
denoted the chosen significance level and f
(1, Nn) denoted the corresponding quantile of the Fisher distribution. The hypergeometric distribution was used to calculate the probability of observing individual GO terms as enriched in sets of differentially expressed genes relative to the reference set of expressed genes on the Affymetrix U74Av2 microarray.
Principal component analysis
PCA was first performed on the 8287 expressed genes and then on subsets of genes of increasing size, the genes being according to their discriminatory power for trisomy. The discriminatory power of gene i was measured by the absolute value of its Student statistic ti:
![]() |
and
the mean expression of gene i for the trisomic and euploid mice, respectively, and si the intraclass variance estimate.
Real-time quantitative PCR
From each euploid and trisomic pool, 250 ng and 6 µg of total RNA were individually converted into cDNAs overnight at 37°C in the presence of 200 U of Superscript II (In vitrogen), 1x first-strand buffer, 100 µM DTT, 5 µM random hexamers oligonucleotides (Pasteur Institute) and 500 µM dNTPs.
The cDNAs were then diluted 1 : 5 and 1 : 14 for the 250 ng and 6 µg RNA samples, respectively, for qPCR according to the following protocols: qPCR on the 78 MMU16 genes and 18 control genes mapping to chromosomes other than 16 were performed as Taqman assays in an ABI 7900 Sequence Detection System (Applied Biosystems). Assay sequences are available in Lyle et al. (17
). Each expression value corresponds to the mean of five replicates.
qPCR on the seven genes outside MMU16 that were found to be differentially expressed by microarray analysis (ALDR, Astrotactin, Dlx1, Glutamine synthase, Homeo box A5, transthyretin, Uncx4.1) were performed in a Lightcycler system (Roche Molecular Biochemicals), in the presence of 0.5 µM of each specific primer (designed by Oligo4 software) and 1x QuantitectTM SYBR® Green PCR master mix (Qiagen) containing 2.5 mM MgCl2, HotstartTaq® Polymerase, dNTP mix and the fluorescent dye SYBR Green I. Each reaction was performed in triplicate. The list of primers is given in Supplementary Material, Table S5. For each sample tested, three control genes (GAPD, tubulin and Ppox) that showed no difference between Ts1Cje and euploids were used to normalize raw data, and the mean relative ratios Ts/Eu were calculated according to geNorm (33
) as described in Lyle et al. (17
).
| ACKNOWLEDGEMENTS |
|---|
We wish to thank David Gentien from Curie Institute Paris for assistance in microarray experiments on the Affymetrix platform, Michel Volovitch and Alain Aurias for helpful discussions, John Ngai and Jean Yee Hwa Yang from UCSF and Nidhi Saran and Roger Reeves from Johns Hopkins Medical School for providing raw data from their published work. This work was supported by EEC grant 00816 to M.C.P. and S.E.A., Fondation Jérôme Lejeune grants to M.C.P. and P.M.S. and to S.E.A., grant from the Ministère de la Recherche to M.C.P., NIH grant HD-31498 to C.J.E., the Swiss National Science Foundation, NCCR Frontiers in Genetics, the ChildCare Foundation to S.E.A. and the Université Paris 5. G.G. had a fellowship from the Fondation Jérôme Lejeune France and RXM a C.J. Martin fellowship from NHMRC Australia.
| SUPPLEMENTARY MATERIAL |
|---|
|
|
|---|
Supplementary Material is available on HMG Online.
| REFERENCES |
|---|
|
|
|---|
- Gardiner, K., Fortna, A., Bechtel, L. and Davisson, M.T. (2003) Mouse models of Down syndrome: how useful can they be? Comparison of the gene content of human chromosome 21 with orthologous mouse genomic regions. Gene, 318, 137147.[CrossRef][Web of Science][Medline]
- Epstein, C.J., Korenberg, J.R., Anneren, G., Antonarakis, S.E., Ayme, S., Courchesne, E., Epstein, L.B., Fowler, A., Groner, Y., Huret, J.L. et al. (1991) Protocols to establish genotypephenotype correlations in Down syndrome. Am. J. Hum. Genet., 49, 207235.[Web of Science][Medline]
-
Chrast, R., Scott, H.S., Papasavvas, M.P., Rossier, C., Antonarakis, E.S., Barras, C., Davisson, M.T., Schmidt, C., Estivill, X., Dierssen, M. et al. (2000) The mouse brain transcriptome by SAGE: differences in gene expression between P30 brains of the partial trisomy 16 mouse model of Down syndrome (Ts65Dn) and normals. Genome Res., 10, 20062021.
[Abstract/Free Full Text] -
FitzPatrick, D.R., Ramsay, J., McGill, N.I., Shade, M., Carothers, A.D. and Hastie, N.D. (2002) Transcriptome analysis of human autosomal trisomy. Hum. Mol. Genet., 11, 32493256.
[Abstract/Free Full Text] - Mao, R., Zielke, C.L., Zielke, H.R. and Pevsner, J. (2003) Global up-regulation of chromosome 21 gene expression in the developing Down syndrome brain. Genomics, 81, 457467.[CrossRef][Web of Science][Medline]
-
Saran, N.G., Pletcher, M.T., Natale, J.E., Cheng, Y. and Reeves, R.H. (2003) Global disruption of the cerebellar transcriptome in a Down syndrome mouse model. Hum. Mol. Genet., 12, 20132019.
[Abstract/Free Full Text] -
Amano, K., Sago, H., Uchikawa, C., Suzuki, T., Kotliarova, S.E., Nukina, N., Epstein, C.J. and Yamakawa, K. (2004) Dosage-dependent over-expression of genes in the trisomic region of Ts1Cje mouse model for Down syndrome. Hum. Mol. Genet., 13, 13331340.
[Abstract/Free Full Text] - Antonarakis, S.E., Lyle, R., Dermitzakis, E.T., Reymond, A. and Deutsch, S. (2004) Chromosome 21 and Down syndrome: from genomics to pathophysiology. Nat. Rev. Genet., 5, 114.
- Reeves, R.H., Irving, N.G., Moran, T.H., Wohn, A., Kitt, C., Sisodia, S.S., Schmidt, C., Bronson, R.T. and Davisson, M.T. (1995) A mouse model for Down syndrome exhibits learning and behaviour deficits. Nat. Genet., 11, 177184.[CrossRef][Web of Science][Medline]
-
Sago, H., Carlson, E.J., Smith, D.J., Kilbridge, J., Rubin, E.M., Mobley, W.C., Epstein, C.J. and Huang, T.T. (1998) Ts1Cje, a partial trisomy 16 mouse model for Down syndrome, exhibits learning and behavioral abnormalities. Proc. Natl Acad. Sci. USA, 95, 62566261.
[Abstract/Free Full Text] - Sago, H., Carlson, E.J., Smith, D.J., Rubin, E.M., Crnic, L.S., Huang, T.T. and Epstein, C.J. (2000) Genetic dissection of region associated with behavioral abnormalities in mouse models for Down syndrome. Pediatr. Res., 48, 606613.[Web of Science][Medline]
-
Baxter, L.L., Moran, T.H., Richtsmeier, J.T., Troncoso, J. and Reeves, R.H. (2000) Discovery and genetic localization of Down syndrome cerebellar phenotypes using the Ts65Dn mouse. Hum. Mol. Genet., 9, 195202.
[Abstract/Free Full Text] - Olson, L.E., Roper, R.J., Baxter, L.L., Carlson, E.J., Epstein, C.J. and Reeves, R.H. (2004) Down syndrome mouse models Ts65Dn, Ts1Cje, and Ms1Cje/Ts65Dn exhibit variable severity of cerebellar phenotypes. Dev. Dyn., 230, 581589.[CrossRef][Web of Science][Medline]
- Diaz, E., Ge, Y., Yang, Y.H., Loh, K.C., Serafini, T.A., Okazaki, Y., Hayashizaki, Y., Speed, T.P., Ngai, J. and Scheiffele, P. (2002) Molecular analysis of gene expression in the developing pontocerebellar projection system. Neuron, 36, 417434.[CrossRef][Web of Science][Medline]
- Kagami, Y. and Furuichi, T. (2001) Investigation of differentially expressed genes during the development of mouse cerebellum. Gene Expr. Patterns, 1, 3959.[CrossRef][Medline]
-
Golfier, G., Dang, M.T., Dauphinot, L., Graison, E., Rossier, J. and Potier, MC. (2004) VARAN: a web server for variability analysis of DNA microarray experiments. Bioinformatics, 20, 16411643.
[Abstract/Free Full Text] -
Lyle, R., Gehrig, C., Neergaard-Henrichsen, C., Deutsch, S. and Antonarakis, S.E. (2004) Gene expression from the aneuploid chromosome in a trisomy mouse model of down syndrome. Genome Res., 14, 12681274.
[Abstract/Free Full Text] - Huang, T.T., Yasunami, M., Carlson, E.J., Gillespie, A.M., Reaume, A.G., Hoffman, E.K., Chan, P.H., Scott, R.W. and Epstein, C.J. (1997) Superoxide-mediated cytotoxicity in superoxide dismutase-deficient fetal fibroblasts. Arch. Biochem. Biophys., 344, 424432.[CrossRef][Web of Science][Medline]
-
Kahlem, P., Sultan, M., Herwig, R., Steinfath, M., Balzereit, D., Eppens, B., Saran, N.G., Pletcher, M.T., South, S.T., Stetten, G. et al. (2004) Transcript level alterations reflect gene dosage effects across multiple tissues in a mouse model of down syndrome. Genome Res., 14, 12581267.
[Abstract/Free Full Text] -
Pollack, J.R., Sorlie, T., Perou, C.M., Rees, C.A., Jeffrey, S.S., Lonning, P.E., Tibshirani, R., Botstein, D., Borresen-Dal, A.L. and Brown, P.O. (2002) Microarray analysis reveals a major direct role of DNA copy number alteration in the transcriptional program of human breast tumors. Proc. Natl Acad. Sci. USA, 99, 1296312968.
[Abstract/Free Full Text] - Zhou, X., Cole, S.W., Hu, S. and Wong, D.T. (2004) Detection of DNA copy number abnormality by microarray expression analysis. Hum. Genet., 114, 464467.[CrossRef][Web of Science][Medline]
- Janoueix-Lerosey, I., Novikov, E., Monteiro, M., Gruel, N., Schleiermacher, G., Loriod, B., Nguyen, C. and Delattre, O. (2004) Gene expression profiling of 1p35-36 genes in neuroblastoma. Oncogene, 23, 59125922.[CrossRef][Web of Science][Medline]
- Chatelin, L., Volovitch, M., Joliot, A.H., Perez, F. and Prochiantz, A. (1996) Transcription factor hoxa-5 is taken up by cells in culture and conveyed to their nuclei. Mech. Dev., 55, 111117.[CrossRef][Web of Science][Medline]
- Sanlioglu, S., Zhang, X., Baader, S.L. and Oberdick, J. (1998) Regulation of a Purkinje cell-specific promoter by homeodomain proteins: repression by engrailed-2 vs. synergistic activation by Hoxa5 and Hoxb7. J. Neurobiol., 36, 559571.[CrossRef][Web of Science][Medline]
- Letinic, K., Zoncu, R. and Rakic, P. (2002) Origin of GABAergic neurons in the human neocortex. Nature, 417, 645649.[CrossRef][Medline]
-
Yun, K., Fischman, S., Johnson, J., Hrabe de Angelis, M., Weinmaster, G. and Rubenstein, J.L. (2002) Modulation of the notch signaling by Mash1 and Dlx1/2 regulates sequential specification and differentiation of progenitor cell types in the subcortical telencephalon. Development, 129, 50295040.
[Abstract/Free Full Text] -
Gazit, R., Krizhanovsky, V. and Ben-Arie, N. (2004) Math1 controls cerebellar granule cell differentiation by regulating multiple components of the Notch signaling pathway. Development, 131, 903913.
[Abstract/Free Full Text] - Zuo, J., De Jager, P.L., Takahashi, K.A., Jiang, W., Linden, D.J. and Heintz, N. (1997) Neurodegeneration in lurcher mice caused by mutation in delta2 glutamate receptor gene. Nature, 388, 769773.[CrossRef][Medline]
- Patil, N., Cox, D.R., Bhat, D., Faham, M., Myers, R.M. and Peterson, A.S. (1995) A potassium channel mutation in weaver mice implicates membrane excitability in granule cell differentiation. Nat. Genet., 11, 126129.[CrossRef][Web of Science][Medline]
-
Schwarzman, A.L., Gregori, L., Vitek, M.P., Lyubski, S., Strittmatter, J., Enghilde, J.J., Bhasin, R., Silverman, J., Weisgraber, K.H., Coyle, P.K. et al. (1994) Transthyretin sequesters amyloid ß protein and prevents amyloid formation. Proc. Natl. Acad. Sci. USA, 91, 83688372.
[Abstract/Free Full Text] -
Li, M.D., Kane, J.K., Matta, S.G., Blaner, W.S. and Sharp, B.M. (2000) Nicotine enhances the biosynthesis and secretion of transthyretin from the choroid plexus in rats: implications for beta-amyloid formation. J. Neurosci., 20, 13181323.
[Abstract/Free Full Text] - Bahn, S., Mimmack, M., Ryan, M., Caldwell, M.A., Jauniaux, E., Starkey, M., Svendsen, C.N. and Emson, P. (2002) Neuronal target genes of the neuron-restrictive silencer factor in neurospheres derived from fetuses with Down's syndrome: a gene expression study. Lancet, 359, 310315.[CrossRef][Web of Science][Medline]
-
Vandesompele, J., De Preter, K., Pattyn, F., Poppe, B., Van Roy, N., De Paepe, A. and Speleman, F. (2002) Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol., 3, RESEARCH0034.
This article has been cited by other articles:
![]() |
D. K. Slonim, K. Koide, K. L. Johnson, U. Tantravahi, J. M. Cowan, Z. Jarrah, and D. W. Bianchi Functional genomic analysis of amniotic fluid cell-free mRNA suggests that oxidative stress is significant in Down syndrome fetuses PNAS, June 9, 2009; 106(23): 9425 - 9429. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Fernandez-Martinez, E. M. Vela, M. Tora-Ponsioen, O. H. Ocana, M. A. Nieto, and J. Galceran Attenuation of Notch signalling by the Down-syndrome-associated kinase DYRK1A J. Cell Sci., May 15, 2009; 122(10): 1574 - 1583. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Esposito, J. Imitola, J. Lu, D. De Filippis, C. Scuderi, V. S. Ganesh, R. Folkerth, J. Hecht, S. Shin, T. Iuvone, et al. Genomic and functional profiling of human Down syndrome neural progenitors implicates S100B and aquaporin 4 in cell injury Hum. Mol. Genet., February 1, 2008; 17(3): 440 - 457. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. Besson, V. Brault, A. Duchon, D. Togbe, J.-C. Bizot, V. F.J. Quesniaux, B. Ryffel, and Y. Herault Modeling the monosomy for the telomeric part of human chromosome 21 reveals haploinsufficient genes modulating the inflammatory and airway responses Hum. Mol. Genet., September 1, 2007; 16(17): 2040 - 2052. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Chang, P. Ostling, M. Akerfelt, D. Trouillet, M. Rallu, Y. Gitton, R. El Fatimy, V. Fardeau, S. Le Crom, M. Morange, et al. Role of heat-shock factor 2 in cerebral cortex formation and as a regulatorof p35 expression. Genes & Dev., April 1, 2006; 20(7): 836 - 847. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||









