Skip to main content

Constitutional trisomy 8 mosaicism as a model for epigenetic studies of aneuploidy

Abstract

Background

To investigate epigenetic patterns associated with aneuploidy we used constitutional trisomy 8 mosaicism (CT8M) as a model, enabling analyses of single cell clones, harboring either trisomy or disomy 8, from the same patient; this circumvents any bias introduced by using cells from unrelated, healthy individuals as controls. We profiled gene and miRNA expression as well as genome-wide and promoter specific DNA methylation and hydroxymethylation patterns in trisomic and disomic fibroblasts, using microarrays and methylated DNA immunoprecipitation.

Results

Trisomy 8-positive fibroblasts displayed a characteristic expression and methylation phenotype distinct from disomic fibroblasts, with the majority (65%) of chromosome 8 genes in the trisomic cells being overexpressed. However, 69% of all deregulated genes and non-coding RNAs were not located on this chromosome. Pathway analysis of the deregulated genes revealed that cancer, genetic disorder, and hematopoiesis were top ranked. The trisomy 8-positive cells displayed depletion of 5-hydroxymethylcytosine and global hypomethylation of gene-poor regions on chromosome 8, thus partly mimicking the inactivated X chromosome in females.

Conclusions

Trisomy 8 affects genes situated also on other chromosomes which, in cooperation with the observed chromosome 8 gene dosage effect, has an impact on the clinical features of CT8M, as demonstrated by the pathway analysis revealing key features that might explain the increased incidence of hematologic malignancies in CT8M patients. Furthermore, we hypothesize that the general depletion of hydroxymethylation and global hypomethylation of chromosome 8 may be unrelated to gene expression regulation, instead being associated with a general mechanism of chromatin processing and compartmentalization of additional chromosomes.

Background

Constitutional trisomy 8 mosaicism (CT8M) is a relatively rare chromosomal disorder with an estimated frequency of approximately 1/25,000 to 1/50,000 [1]. However, since the phenotypes of individuals with CT8M vary quite extensively, ranging from severe malformations with impaired cognitive functioning to rather discrete dysmorphic changes [2], the true prevalence may well be higher. Characteristic clinical features of CT8M include elongated facial features, abnormally shaped ears, strabismus, camptodactyly, clinodactyly, deep plantar and palmar skin furrows, vertebral/hip anomalies, cardiovascular and urogenital malformations and mild to moderate mental retardation [3]. As regards acquired trisomy 8, it is one of the most common abnormalities in malignant myeloid disorders, such as acute myeloid leukemia, myelodysplastic syndromes and myeloproliferative neoplasms [4, 5]. Interestingly, based on several case reports on myeloid malignancies in patients with CT8M [see Additional file 1: Table S1], CT8M seems to be associated with an increased risk of these disorders. It has even been suggested, albeit in a small patient series, that the ‘acquired’ +8 in myeloid diseases in some instances could represent an unrecognized CT8M [6]. This, however, remains to be confirmed or refuted in larger studies.

There is ample evidence that CT8M arises post-zygotically through a mitotic error, with no preferential parental origin of the gained chromosome 8 (chr8) [7–11]. The pathogenetically important consequence of CT8M is most likely, as for numerical abnormalities in general, a dosage effect, with overexpression of genes on the gained chromosome [12]. However, it is unknown whether all chr8 genes are upregulated in CT8M. In a previous study [10], a single gene (GSR at 8p) was investigated in two patients and shown to have an increased activity. To date, no global gene expression analyses of CT8M cases have been reported, but considering that myeloid malignancies with acquired trisomy 8 overexpress most, but not all, genes on this chromosome [13–15], it is plausible that the same holds true also for CT8M. Furthermore, it is also unclear whether CT8M is associated with epigenetic changes involving chr8, akin to the methylation changes on chr21 in Down syndrome (DS) [16, 17]. We have previously shown that most tri- and tetrasomic chromosomes in high hyperdiploid (51 to 67 chromosomes) acute lymphoblastic leukemia (ALL), one of the most common subtypes of childhood ALL [18], are globally less methylated than their disomic counterparts and suggested that this could be a general phenomenon associated with gains of chromosomes, irrespective of whether they are constitutional or acquired [19]. However, this has, until the present study, not been investigated.

In order to further our understanding of gene expression patterns and epigenetic mechanisms associated with gain of chr8, we used CT8M as a model, allowing analyses of single cell-derived fibroblast clones, harboring either trisomy or disomy 8, from the same patient. This model enables comparisons between identical cell types (for example, fibroblasts), differing only in respect to chr8 content, from a single individual and, hence, circumvents the bias introduced when comparing the epigenetic findings with unrelated, healthy individuals. We used this model to profile global gene and microRNA (miRNA) expression as well as genome-wide and gene promoter specific DNA methylation and hydroxymethylation patterns in trisomic and disomic fibroblasts from a patient with CT8M, using microarrays and methylated DNA immunoprecipitation (MeDIP).

Results

Unsupervised HCA of gene and miRNA expression

A standard unsupervised hierarchical cluster analysis (HCA), after using a variance filter with a cut-off value of 50% of the highest standard deviating genes to limit the number of probes and a Pearson correlation test [20] for average linkage clustering, was used to investigate similarities and differences in global gene and miRNA expression patterns among disomy 8, trisomy 8, and references. Unsupervised HCA of only gene expression and of both gene and miRNA expression data clustered the trisomy 8 cultures together but did not differentiate between disomy 8 and reference cultures [see Additional file 2: Figure S1], whereas unsupervised HCA of only miRNA expression data was less accurate in clustering the various culture groups [see Additional file 2: Figure S1]. To test for a possible bias introduced by in vitro culturing, we performed long-term culturing (16 weeks) of disomy 8 and trisomy 8 cells and repeated the analysis. Unsupervised HCA clustered the long-term cultures together with their short-term counterparts, indicating a minimal impact of in vitro effects on gene and miRNA expression [see Additional file 2: Figure S1]. It should be stressed that we cannot exclude the possibility of an effect of culturing as such during the first eight weeks compared with uncultured cells. However, since analyses could not be performed without expansion of single cell clones, this could not be addressed in the present study.

PCA and supervised HCA of gene and miRNA expression

Principal component analysis (PCA) identified 1,650 unique probes (5% of all 33,297 probes analyzed in our expression dataset) that displayed significant expression differences among the trisomy 8, disomy 8, and reference groups. Supervised HCA of these probes correctly clustered the three groups, with the trisomy 8 cultures being placed in the same branch as the disomy 8 cultures (Figure 1A). When performing PCA only of chr8 probes (n = 1,087), 87 significant genes (12% of the 728 known HG19 protein-coding genes on this chromosome) were identified. Supervised HCA of these genes again accurately clustered the three groups, this time with the trisomy 8 cultures being placed in a single branch (Figure 1B).

Figure 1
figure 1

PCA and supervised HCA of gene and miRNA expression accurately cluster the trisomy 8, disomy 8, and reference cultures. (A) The PCA identified 1,650 probes that displayed significant expression differences among the trisomy 8, disomy 8, and reference groups (left). Supervised HCA of these probes accurately clustered the subgroups, as illustrated in the heat map to the right. (B) A similar pattern was observed when analyzing only chr8 data, with the PCA identifying 87 genes (left) and HCA accurately clustering the subgroups (right), placing the trisomy 8 cultures in a single branch. (C) Median-centered gene expression of chr8 demonstrated that the majority of the genes on this chromosome in the trisomy 8 cultures were upregulated compared with the disomy 8 cultures (left), whereas few differences were observed when comparing the disomy 8 cultures with the reference cultures (right). (D) The PCA identified 108 miRNA probes that displayed significant expression differences between the trisomy 8, disomy 8, and reference groups (left). Supervised HCA of these probes grouped the cultures correctly (right). chr8, chromosome 8; HCA, hierarchical cluster analysis; miRNA, microRNA; PCA, principal component analysis.

To illustrate the expression levels of genes on chr8, the average median-centered gene expression in the trisomy 8, disomy 8, and reference cultures was plotted against the location of each gene on this chromosome, revealing a global overexpression of chr8 genes in the trisomy 8 cultures compared with both the disomy 8 and reference cultures (Figure 1C and Additional file 3: Figure S2); in contrast, the expression levels in the disomy 8 and the reference cultures were very similar (Figure 1C). In total, among the 728 protein-coding genes on chr8, 476 (65%) were overexpressed in the trisomy 8 cells compared with the disomy 8 cells, as ascertained by measuring the median expression values in the two culture groups.

PCA and supervised HCA of global miRNA expression, based on 108 significant probes, also correctly clustered the three groups, placing the trisomy 8 cultures in the same branch as the disomy 8 cultures (Figure 1D). Analysis of only chr8 miRNA expression was uninformative because of too few differentially expressed probes (data not shown).

Genes and miRNAs possibly contributing to the CT8M phenotype and their clinical and biological functions

A multiple testing approach identified 25 protein-coding genes, one small nucleolar RNA, and one miRNA that were significantly overexpressed and four genes and one miRNA that were significantly underexpressed in the trisomy 8 cultures as compared with both the disomy 8 and reference cultures. The over- or underexpression of seven of the genes identified was validated by real-time quantitative PCR (qPCR) [see Additional file 4: Figure S3]. As seen in Table 1, 10 (31%) of these 32 genes/miRNAs are located on chr8. In addition, a t-test of the global gene expression patterns in the trisomy 8 versus the disomy 8 cultures revealed 502 differentially expressed genes. The top clinical conditions and biological features associated with these genes, as identified using an Ingenuity pathway analysis (Ingenuity Systems, Redwood City, CA, USA), are provided in Additional file 5: Table S2.

Table 1 Genes and miRNAs significantly associated with the presence of trisomy 8

Hypermethylated and hyperhydroxymethylated promoters/CpG islands in relation to gene expression

An edge-preserving smoother analysis identified 150 to 200 genes that were highly enriched for 5-methylcytosine (5 mC) and/or 5-hydroxymethylcytosine (5 hmC) in the disomy 8, trisomy 8, and reference groups, respectively. As seen in Additional file 6: Table S3, approximately 80% of the hypermethylated genes were underexpressed, whereas only a minority of the hyperhydroxymethylated genes was under- or overexpressed, as ascertained by a t-test. When analyzing only genes on chr8, a similar pattern emerged – 17 (89%) of 19 hypermethylated genes were underexpressed, whereas only 8 (32%) of 25 hyperhydroxymethylated genes were underexpressed.

The patterns of hypermethylated and hyperhydroxymethylated promoters/CpG islands

The levels of promoter-specific methylation on the autosomes in the trisomy 8 group did not differ significantly from the corresponding levels in the disomy 8 and reference groups combined [see Additional file 7: Figure S4], nor did the levels of 5 hmC, with the notable exception of chr8 on which the 5 hmC levels were significantly lower (P <0.05; t-test) in the trisomy 8 group (Figure 2 and Additional file 8: Figure S5). As regards the X chromosome, its promoters/CpG islands were generally (irrespective of gender and culture group) significantly (P <0.05; t-test) less hydroxymethylated than all ‘autosomal’ promoters/CpG islands [see Additional file 9: Figure S6]. The female reference displayed a higher level of promoter-specific methylation on the X chromosome (Figure 3A) as well as a lower level of promoter-specific hydroxymethylation (Figure 3B) compared with the XY cultures.

Figure 2
figure 2

Promoter-specific hydroxymethylation levels on chromosome 8 are significantly lower in the trisomy 8 compared with the disomy 8/reference cultures. The levels of average promoter-specific hydroxymethylation on chr8 in the trisomy 8 and in the disomy 8 and reference cultures combined were log2-converted, median-centered, and plotted against genomic positions. Significantly lower levels (P < 0.05; t-test) were observed in the trisomy 8 cultures. chr8, chromosome 8.

Figure 3
figure 3

Promoter-specific methylation and hydroxymethylation on the X chromosome in the female culture is higher and lower, respectively, compared with the male cultures. A rank plot of median-centered promoter-specific methylation and hydroxymethylation log2 levels on the X chromosome reveals that the female reference culture generally displays higher levels (P <0.05; t-test) of methylation compared with (A) the average levels in the male cultures and (B) lower levels (P <0.05; t-test) of hydroxymethylation compared with the male cultures.

Genome-wide methylation patterns

PCA identified 781 unique clones, comprising 2.4% of the 32,433 bacterial artificial chromosomes (BACs) in the array, that displayed significant (P <0.05) differences in methylation between the trisomy 8, disomy 8, and reference groups. Supervised HCA of these clones clustered the three groups, with the trisomy 8 cultures being placed in the same branch as the disomy 8 cultures (Figure 4A). PCA of only BACs on chr8 revealed 81 significant clones (6.2% of the 1,312 BACs covering chr8), and supervised HCA of these BACs again accurately clustered the three groups, this time with the trisomy 8 cultures constituting a single branch (Figure 4B).

Figure 4
figure 4

Accurate clustering of the trisomy 8, disomy 8, and reference cultures by PCA and supervised HCA of genome-wide methylation patterns and significant hypomethylation of gene-poor regions on chromosome 8 only in the trisomy 8 cultures. (A) The PCA identified 781 BACs that displayed significant differences between the trisomy 8, disomy 8, and reference cultures (left). Supervised HCA of these BACs accurately clustered the subgroups, as illustrated in the heat map (right). (B) A similar pattern was observed when analyzing only chr8 data, with the PCA identifying 81 BACs (left) and HCA accurately clustering the subgroups (right). (C) Mean log2 ratio of the three culture subgroups plotted against the number of genes on the BACs for chr8 revealed a significant hypomethylation of gene-poor regions in trisomy 8 cultures. The black curve shows the proportions that the BAC clones with different number of genes constitute of all BAC clones on chr 8; for example, the BACs with no genes comprise 45% of all the BAC clones on this chromosome. (D) An independent and weighted ANOVA and a Tukey HSD test verified the global hypomethylation of gene-poor regions on chr8 in trisomy 8 cultures but not in the disomy 8 or reference cultures. (E) Mean log2 ratio of the three culture subgroups plotted against the number of genes on the BACs for chromosome 7 (used as an example of another autosome) revealed no significant differences in genome-wide methylation. (F) Mean log2 ratio of female and male cultures plotted against the number of genes on the BACs for chromosome X revealed significant hypomethylation of gene-poor regions only in the female culture. ANOVA, analysis of variance; BAC, bacterial artificial chromosome; chr8, chromosome 8; HCA, hierarchial cluster analysis; HSD, honestly significant difference; PCS, principal component analysis.

An independent and weighted ANOVA and a Tukey honestly significant difference (HSD) test identified a significant (P <0.01) global hypomethylation of gene-poor regions on chr8 in trisomy 8 cultures but not in the disomy 8 or reference cultures (Figure 4C and D). This pattern was not observed on any of the other autosomes (Figure 4E). However, a similar global hypomethylation of gene-poor regions was seen on the X chromosome in the female reference compared with all XY cultures, irrespective of sample type (Figure 4F).

Discussion

Aneuploidy is a well-studied phenomenon in human cancer as well as in various model organisms, for example Arabidopsis thaliana, Drosophila melanogaster, Mus musculus, and Saccharomyces cerevisiae, that in recent years has received increased attention as regards its effects on cellular physiology, transcriptional consequences, and genomic instability [21–24]. Such studies have revealed that aneuploidy is associated with a general proliferative disadvantage, species-conserved stress related gene expression patterns regardless of gained chromosome as well as increased and defect mitotic recombination and DNA repair, resulting in genomic instability. In humans, constitutional aneuploidy, for example, loss or gain of sex chromosomes or trisomy of chromosomes 13, 18, or 21, is a major cause of miscarriage and developmental disturbances. Prior to the present study, gene expression patterns in patients with constitutional trisomies have exclusively been compared with expression signatures derived from disomic control cells from unrelated healthy individuals [25–31]. The differences previously observed might, hence, not only represent the presence of a specific trisomy but also interindividual variations and comparisons between non-identical cell types, with the latter introducing severe interpretation difficulties. For this reason, we used single cell cloning of fibroblasts from an individual with CT8M to ascertain, in an unbiased manner, the epigenetic changes associated with gain of chr8. Because several human chromosomal disorders occasionally, or even frequently, occur as postzygotic events, resulting in two or more chromosomally different cell lines developed from a single zygote, the present approach is applicable also for epigenetic studies of other constitutional mosaicisms, with well-known examples including, apart from CT8M, triplications of 8p, 9p, 12p, and 18p, trisomies for chromosomes 9, 13, 14, 16, 18, 20, 21, and 22, and Turner and Klinefelter syndromes [32].

By single cell cloning with subsequent gene expression and methylation analyses, we here show that trisomy 8-positive fibroblasts from a patient with CT8M display a characteristic expression and methylation phenotype that is clearly distinct from both disomy 8 cells and normal reference fibroblasts. In fact, PCA and HCA of global as well as chr8-specific gene/miRNA expression and methylation profiles grouped all the trisomy 8 cultures in the same cluster branch, with no discernible effect by long-term culturing [see Additional file 2: Figure S1]. When investigating which genes were the major ‘expression contributors to the observed clustering, we could, not surprisingly, demonstrate that the majority of the genes on chr8 in the trisomy 8 cultures were overexpressed compared with the chr8 genes in the disomy 8/reference groups (Figure 1C and Additional file 3: Figure S2 and Additional file 4: Figure S3). Hence, it seems safe to assume that the clinical features as well as the elevated risk for malignant myeloid disorders associated with CT8M [3, 4] [see Additional file 1: Table S1] is, at least partly, due to a direct chr8 gene dosage effect. However, far from all deregulated genes and non-coding RNAs in the trisomy 8-positive cells were located on chr8 (Table 1), nor were all genes on this chromosome up-regulated (Figure 1C). Thus, all genes on chr8 are not dosage sensitive and trisomy 8 directly or indirectly affects genes situated on other chromosomes as well. That not every gene on a constitutional trisomic chromosome displays increased expression has previously been clearly demonstrated in DS; in fact, only a minority of the chr21 genes are overexpressed in DS [25–28].

Pathway analysis of the deregulated genes in the trisomy 8 cells revealed that ‘cancer’ and ‘genetic disorder’ were the top two biological functions in the diseases and disorders category and that ‘hematological system development and function’ and ‘hematopoiesis’ were top ranked in the physiological system development and function category [see Additional file 5: Table S2]. The latter may explain the increased incidence of myeloid neoplasms in CT8M patients [4]. Furthermore, by using a set of stringent inclusion criteria, several candidate genes (Table 1) for some of the clinical manifestations were identified [see Additional file 10: Table S4]. For example, the presence of congenital heart disease could possibly be ascribed to the deregulation of the DES, LDB3, POPDC2 and SLC4A3 genes [see Additional file 11: Figure S7], all of which have previously been associated with cardiac function and pathology, such as cardiomyopathy. In addition, eight of the candidate genes (CASP4, COX6C, FNBP1, LRP12, RAB2A, RPS6KA1, TNFRSF10A and TNKS; Table 1) have been associated with various solid tumors. These genes also displayed a joint functional association [see Additional file 11: Figure S7]. The overexpression of MIR151, which was observed in the trisomy 8 cultures only, is noteworthy because this non-coding RNA has been suggested to play a role in metastases [see Additional file 10: Table S4]. Furthermore, the deregulation of AGPAT6, CPXM2, CRYAB, EDNRA, ENTPD4, RAB2A and RAB3GAP2 as well as their functional network [see Additional file 11: Figure S7] may contribute to the mild to moderate mental retardation seen in CT8M patients since these genes have all previously been associated with neurological development and function. In addition, ELP3 and SNORD115-32, involved in small nucleolar RNA processing and transcript elongation, were upregulated in the trisomy 8 cells. They could also play a role in the cognitive impairment, considering that ELP3 regulates the maturation of projection neurons and that SNORD115-32 has been implicated in splicing defects of the serotonin receptor 2C in Prader-Willi syndrome [see Additional file 10: Table S4].

The present analyses of global methylation and hydroxymethylation patterns in trisomy 8-positive cells may provide insights into the epigenetic consequences of aneuploidy. As a proof-of-principle, we initially investigated methylation and hydroxymethylation of the X chromosome in XX and XY cultures (Figure 3). As expected, unlike its active counterpart, the inactivated X chromosome (Xi) in the XX cells exhibited hypomethylation of gene-poor regions but hypermethylation of a substantial proportion of promoters (Figures 3A and 4F); the latter is considered a sex chromosome dosage compensation mechanism [33–35]. Furthermore, and in agreement with a previous study [36], we observed a general depletion of 5 hmC on the X chromosome compared with all autosomes in both males and females [see Additional file 9: Figure S6]. Interestingly, similar methylation and hydroxymethylation profiles of chr8 were observed in the trisomy 8 cells, which – compared with the disomic cells – displayed depletion of 5 hmC (Figure 2) and a global hypomethylation of gene-poor regions (Figure 4C). However, unlike the Xi, no elevated levels of promoter/CpG island 5 mC content on chr8 were detected [see Additional file 7: Figure S4]. It is interesting to note that the methylation status of chr8 observed in the present study is reminiscent of the hypomethylation patterns of acquired chromosomal gains reported in neoplasia, that is, hypomethylation of gene-poor regions of the tri-/tetrasomic chromosomes in high hyperdiploid childhood ALL and of trisomies 7 and 14 in colon cancer [19, 37]. Thus, we suggest that an association between aneuploidy and globally lowered levels of 5 mC on the gained chromosomes is valid for both constitutional and acquired numerical changes. It should be emphasized, however, that the hypomethylation seems to be confined to genomic regions with no or few genes and this phenomenon may, hence, be unrelated to gene expression changes. Instead, it is possible that the hypomethylation of gene-poor regions may be involved in chromatin compartmentalization. In fact, it has been demonstrated that absence of methylation influences chromatin organization by increasing the nuclear clustering of pericentric heterochromatin [38]. It is hence feasible that the hypomethylation of gene-poor regions results in altered nuclear positioning of trisomic/tetrasomic chromosomes. However, such an effect on the nuclear architecture remains to be addressed experimentally.

The biological impact of the decreased hydroxymethylation levels of the gene promoters/CpG islands on chr8 in the trisomy 8 cells and on the X chromosome is currently unclear. In fact, it is debated whether the main cellular role of 5 hmC is to be a demethylation intermediate or if it constitutes a bona fide epigenetic regulator, attracting chromatin and transcriptional modifiers by itself. To date, no specific 5 hmC-binding protein has been found, arguing against the latter mechanism [39]. Whether a demethylation intermediate or not, it is clear that 5 hmC is situated primarily in the gene-bodies of actively transcribed genes [40], playing an important role in cellular lineage commitment as well as in tissue-specific gene expression [41, 42]. In support of this, it has recently been demonstrated [43] that the most important predictor of 5 hmC content in a cell is the tissue type from which it is derived and that the levels of hydroxymethylation tend to decrease rapidly over time when cells are cultured in vitro. This possible bias was circumvented in our study by the use of identical cell types as well as culture times. Interestingly, global loss of 5 hmC has been associated with a variety of human malignancies [44, 45] with an inverse relationship between 5 hmC and cell proliferation [46]. Hence, the hypohydroxymethylation on chr8 in the trisomy 8-positive cells may contribute to the propensity for myeloid malignancies in CT8M patients.

Conclusions

In this study we used a mosaic genetic syndrome, namely constitutional trisomy 8 mosaicism, to investigate the epigenetic mechanisms associated with aneuploidy using the patient’s own diploid cells as control. This new approach circumvents the interpretation difficulties introduced when comparing trisomic cells with diploid cells from an unrelated control. We profiled gene and miRNA expression as well as DNA methylation and hydroxymethylation patterns. The trisomic fibroblasts displayed a characteristic expression and methylation patterns, with the majority of chromosome 8 genes in the trisomic cells being overexpressed. However, not all deregulated genes were located on this chromosome. Thus, trisomy 8 affects genes situated also on other chromosomes, which in cooperation with the gene dosage effect might influence the clinical features and elevated cancer risk seen in this syndrome. Interestingly, trisomy 8-positive cells also displayed a general depletion of hydroxymethylation and a global hypomethylation of gene-poor regions on chromosome 8, thus partly mimicking the patterns seen in the inactivated X-chromosome. This novel finding could be associated with a general mechanism of chromatin processing of additional chromosomes.

Methods

Patient history and cytogenetic analyses

In 2002, a cytogenetic analysis of peripheral blood from a 10-year-old boy with a history of developmental delay, cryptorchidism and pectus excavatum revealed trisomy 8 in three of 25 analyzed metaphases. To confirm the diagnosis of CT8M, a skin biopsy was analyzed in 2003, revealing the karyotype 47,XY,+8[5]/46,XY[20].

Single-cell cloning of fibroblasts from the patient with CT8M

The primary culture for cytogenetic analysis of fibroblasts was established using standard methods in our department, including initial mechanical disaggregation and subsequent culture in RPM1 1640 (Invitrogen, Carlsbad, CA, USA). After the chromosome analysis, the fibroblasts were viably frozen in liquid nitrogen. For the present study, the primary culture was thawed and cultured in standard conditions with RPMI 1640 medium (Invitrogen) containing 10% FBS (Sigma-Aldrich, St. Louis, MO, USA) in a 50 ml (25 cm2) cell culture flask (BD Falcon, Stockholm, Sweden). The culture was then, at the time of subconfluence, passaged 1:3 by treatment with 0.5% trypsin. One subculture was refrozen in liquid nitrogen, one was analyzed by interphase fluorescence in situ hybridization (FISH) using CEP-8 (Abbott Molecular, Des Plaines, IL, USA) to evaluate the proportions of cells with either trisomy or disomy 8, and one was used for further passages. After a short culture time of the latter, the cells were trypsinized, counted in a Bürker chamber (Digital Bio, Seoul, Korea), and diluted so that approximately one cell could be plated per well in 96-well microtiter plates (BD Falcon). Viable confluent cell cultures were used for further passages at a ratio of 1:1 to 24-well microtiter plates (BD Falcon). After reaching confluence, all cell cultures were transferred 1:1 to 6-well microtiter plates (BD Falcon). The cell cultures were then resuspended, after which a fraction from each was used for interphase FISH analysis with dual (Cy3 and Cy5) CEP-8 probes (Abbott Molecular). A minimum of 285 (95%) of 300 normal or abnormal nuclei was used as a cut-off to designate a culture as being either 46,XY or 47,XY,+8. Three cultures with disomy 8 and three with trisomy 8 were then passaged to 250 ml (75 cm2) cell culture flasks (BD Falcon). Finally, these six cell cultures were harvested simultaneously, and RNA and DNA were extracted. The total time in culture was eight weeks. Additional file 12: Figure S8 shows a schematic overview of the culture process. The regional ethics board at Lund University approved the study and informed consent was obtained from the patient’s guardians in accordance with the Declaration of Helsinki.

Reference fibroblasts

Primary dermal fibroblasts from a male neonate and a female adult were obtained from the American Type Culture Collection (Manassas, VA, USA). The cells were maintained and passaged in RPMI 1640 medium (Invitrogen) containing 10% FBS (Sigma-Aldrich) in 250 ml cell culture flasks (BD Falcon). The cultures were cytogenetically characterized using conventional G-banding analysis, revealing normal male and female karyotypes, respectively.

Global gene and miRNA expression analysis

The mirVana kit (Invitrogen) was used to isolate RNA, including miRNA, from the fibroblast cultures from the CT8M patient and the references. The RNA concentration was measured with a NanoDrop 1000 spectrophotometer (Saveen & Werner AB, Malmö, Sweden), and the quality of the RNA was assessed with a 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). Total RNA was then hybridized to the Affymetrix Human Gene 1.0 ST and the GeneChip miRNA arrays (Affymetrix, Santa Clara, CA, USA) according to the manufacturer’s protocols. Quality control and initial analyses were performed using the Expression Console and the miRNA QC Tool (Affymetrix).

Bioinformatic analyses of gene and miRNA expression

HCA, PCA, multiclass significance analysis of microarrays (SAM), and t-tests were used to compare expression profiles between fibroblasts with either disomy 8 (from the CT8M patient and from the references) or trisomy 8 (from the CT8M patient) and to identify differentially expressed genes and miRNAs in the trisomic cells. Filtering, HCA, and PCA of log2 transformed expression data were performed using the Qlucore Omics Explorer v 2.0 (Qlucore, Lund, Sweden). For PCA, based on the Pearson correlation matrix, the data were normalized by setting the mean to 0 and σ to 1 and by filtering the variance on the basis of σ/σmax. The Benjamini-Hochberg method was applied for error correction by q-value calculation. A P-value of <0.05 was used as a cut-off. HCA of the genes identified by PCA was performed on the basis of Euclidean metric distance (samples) and Pearson correlation average linkage test (genes).

Differentially expressed genes and miRNAs were identified using the correlation matrix-based PCA in combination with SAM and t-tests (MultiExperiment Viewer v 4.2) [20]. A multiple testing approach was applied, with the criteria for inclusion being set to a P-value of <0.05 for the t-tests and a q-value of ≤10% in the SAM analysis, when comparing trisomy 8 with disomy 8 (CT8M only) and when comparing trisomy 8 with disomy 8 (CT8M and references grouped together). Only genes and miRNAs significantly differentially expressed in all the above-mentioned statistical analyses were considered to be associated with CT8M. Functional association networks of these genes were identified by the GeneMANIA analysis (http://www.genemania.org/) [47].

qPCR validation

The qPCR analyses were performed using standard protocols for predesigned TaqMan probes (Applied Biosystems, Carlsbad, CA, USA) to validate over- or underexpression of genes significantly associated with the presence of trisomy 8 in the global gene expression analyses (Table 1). The expression of transcripts of the following seven genes was ascertained: AGPAT6 (Hs00410940_m1), CASP4 (Hs01031951_m1), CRYAB (Hs00157107_m1), CPXM2 (Hs00406866_m1), EDNRA (Hs03988672_m1), ENTPD4 (Hs01555223_m1), LRP12 (Hs00940649_m). Fold expression was calculated relative to endogenous control normalization of TBP (Hs.590872_m1), encoding a TATA box binding protein. The rationale for using TBP as control was its uniform expression in all cultures in the global gene expression analyses. All qPCR analyses were performed in triplicate and assayed on a 7500 real-time PCR system (Applied Biosystems) and ΔΔCT values were calculated using the comparative CT method [48].

MeDIP

From each sample (fibroblasts from the CT8M patient and the references), 12 μg of DNA was diluted in standard Tris-ethylenediaminetetraacetic acid ( TE) buffer to a concentration of 75 ng/μl. The DNA was sonicated for 90 seconds with a Covaris S220 (Covaris, Woburn, MA, USA) set to a 5% duty cycle, intensity value 3, and 200 cycles per burst to generate fragments ranging from 200 to 1,000 bp, as verified by gel electrophoresis. The fragmented DNA was then diluted to a concentration of 13.4 ng/μl and denatured for 10 minutes at 95°C, after which two input fractions (approximately 250 ng each), serving as references in subsequent microarray hybridizations (see below), were collected and stored at −80°C. The remaining DNA aliquots were immunoprecipitated overnight with monoclonal antibodies against either 5 mC or 5 hmC using magnetic bead-based MeDIP and hMeDIP kits (Diagenode, Liège, Belgium) to purify and isolate the precipitated DNA. Both the input fraction and the immunoprecipitated DNAs were then amplified using a WGA2 kit (Sigma-Aldrich), generating sufficient amounts (approximately 6 μg) for the various microarray hybridizations detailed below.

Genome-wide methylation profiling

Of the methylated and input fraction DNAs, 1.5 μg was differentially labeled with Cy3 and Cy5, respectively, and hybridized to BAC array slides containing 32,433 tiling clones covering at least 98% of the human genome. The slides were produced at the SCIBLU DNA microarray resource center at Lund University, Sweden. Labeling of DNA, slide preparation, hybridization and analysis were performed as described previously [49].

Classification as relative enrichment of methylation was based on the log2 ratio for each of the BAC clones smoothed with the factor of 0.33, with a positive ratio indicating a higher degree of methylation. PCA and HCA were performed as described above. The three groups (disomy 8, trisomy 8, and references) were analyzed in relation to the gene density on each BAC clone by ascertaining, using a custom-made script, the number of annotated RefSeq genes on each consecutive clone. All BACs on each chromosome were grouped in relation to the number of genes they harbor, which most often were 0 to 3 genes per BAC (total range 0 to 20). The average log2 ratio for each group was then plotted against their gene content. For every chromosome, each gene count group, comprising at least 0.1% of the total number of BAC clones on the chromosome, was analyzed using an independent and weighted ANOVA to ascertain whether there were statistically significant differences in 5mC levels between disomy 8, trisomy 8, and references. The Tukey HSD test was applied to investigate between which types of culture the significant differences were present in. The analyses were performed using the VassarStats website (http://vassarstats.net/).

Gene-specific methylation and hydroxymethylation profiling

A total of 4 μg of the immunoprecipitated and input fraction DNAs were analyzed using the Human DNA Methylation 2.1 M Deluxe Promoter Arrays (Roche NimbleGen, Madison, WI, USA) by the NimbleGen Service Group, who performed the quality control, labeling, hybridization, and scanning. This array platform covers all 27,867 University of California Santa Cruz (UCSC)-annotated CpG islands and promoter regions for all known RefSeq genes as well as 730 annotated miRNA promoters.

Signal intensity data were extracted from the scanned images, with each feature on the array having a corresponding log2 ratio, that is, the ratio of the input signals for the precipitated and the input fraction DNAs. PCA and HCA were performed as above. In addition, the log2 ratio was scaled in order to center the ratio data around zero by subtracting the bi-weight mean for the log2 ratio values for all features on the array from each log2 ratio value. A modified algorithm for capturing microarray enrichment (ACME) [50], where a fixed-length window (750 bp) is placed around each consecutive probe, and the one-sided Kolmogorov-Smirnov test were then applied on the scaled log2 ratio data to determine whether the probes were drawn from a significantly more positive distribution (peak) of intensity log2 ratios than the other probes in the array. The resulting score for each probe is the -log10 P- value from the windowed Kolmogorov-Smirnov test around that probe.

For every group (disomy 8, trisomy 8, and reference), all genes with significant CpG island/promoter methylation scores were ascertained. Then, an edge-preserving smoother analysis [51], using the Pott’s filter with the penalty parameter set to 2 and the least allowed aberration size limited to two clones, was performed to identify genes that were highly enriched for 5 mC or 5 hmC in the three different culture groups. The expression of each gene was then investigated, after which possible associations between hypermethylated and hyperhydroxymethylated promoters/CpG islands and gene expression levels were compared using a t-test. Underexpression and overexpression were defined as fold changes ≤ −0.5 and ≥0.5, respectively; expression levels between these cut-offs were considered ‘intermediate’.

All microarray data are minimum information about a microarray experiment (MIAME) compliant and available to download at the Gene Expression Omnibus Archive (http://www.ncbi.nlm.nih.gov/geo/) under the accession number GSE40321.

Abbreviations

ALL:

Acute lymphoblastic leukemia

bp:

Base pair

CT8M:

Constitutional trisomy 8 mosaicism

Chr8:

Chromosome 8

DS:

Down syndrome

FBS:

Fetal bovine serum

FISH:

Fluorescence in situ hybridization

HCA:

Hierarchical cluster analysis

HSD:

Honestly significant difference

MeDIP:

Methylated DNA immunoprecipitation

miRNA:

MicroRNA

PCA:

Principal component analysis

PCR:

Polymerase chain reaction

qPCR:

Quantitative real-time PCR

SAM:

Significance analysis of microarrays

Xi:

Inactivated X chromosome.

References

  1. Gorlin RJ, Cohen MM, Hennekam RC: Syndromes of the head and neck. 2001, New York: Oxford University Press

    Google Scholar 

  2. Hale NE, Keane JF: Piecing together a picture of trisomy 8 mosaicism syndrome. J Am Osteopath Assoc. 2010, 110: 21-23.

    PubMed  Google Scholar 

  3. Riccardi VM: Trisomy 8: an international study of 70 patients. Birth Defects Orig Artic Ser. 1977, 13: 171-184.

    CAS  PubMed  Google Scholar 

  4. Paulsson K, Johansson B: Trisomy 8 as the sole chromosomal aberration in acute myeloid leukemia and myelodysplastic syndromes. Pathol Biol (Paris). 2007, 55: 37-48. 10.1016/j.patbio.2006.04.007.

    Article  CAS  Google Scholar 

  5. Mitelman F, Johansson B, Mertens F: Mitelman Database of Chromosome Aberrations and Gene Fusions in Cancer. 2012,http://cgap.nci.nih.gov/Chromosomes/Mitelman.

    Google Scholar 

  6. Maserati E, Aprili F, Vinante F, Locatelli F, Amendola G, Zatterale A, Milone G, Minelli A, Bernardi F, Lo Curto F, Pasquali F: Trisomy 8 in myelodysplasia and acute leukemia is constitutional in 15-20% of cases. Genes Chromosomes Cancer. 2002, 33: 93-97. 10.1002/gcc.1214.

    Article  CAS  PubMed  Google Scholar 

  7. DeBrasi D, Genardi M, D'Agostino A, Calvieri F, Tozzi C, Varrone S, Neri G: Double autosomal/gonosomal mosaic aneuploidy: study of nondisjunction in two cases with trisomy of chromosome 8. Hum Genet. 1995, 95: 519-525.

    Article  CAS  PubMed  Google Scholar 

  8. Robinson WP, Binkert F, Bernasconi F, Lorda-Sanchez I, Werder EA, Schinzel AA: Molecular studies of chromosomal mosaicism: relative frequency of chromosome gain or loss and possible role of cell selection. Am J Hum Genet. 1995, 56: 444-451.

    PubMed Central  CAS  PubMed  Google Scholar 

  9. James RS, Jacobs PA: Molecular studies of the aetiology of trisomy 8 in spontaneous abortions and the liveborn population. Hum Genet. 1996, 97: 283-286. 10.1007/BF02185754.

    Article  CAS  PubMed  Google Scholar 

  10. Seghezzi L, Maserati E, Minelli A, Dellavecchia C, Addis P, Locatelli F, Angioni A, Balloni P, Miano C, Cavalli P, Danesino C, Pasquali F: Constitutional trisomy 8 as first mutation in multistep carcinogenesis: clinical, cytogenetic, and molecular data on three cases. Genes Chromosomes Cancer. 1996, 17: 94-101. 10.1002/(SICI)1098-2264(199610)17:2<94::AID-GCC4>3.0.CO;2-W.

    Article  CAS  PubMed  Google Scholar 

  11. Karadima G, Bugge M, Nicolaidis P, Vassilopoulos D, Avramopoulos D, Grigoriadou M, Albrecht B, Passarge E, Annerén G, Blennow E, Clausen N, Galla-Voumvouraki A, Tsezou A, Kitsiou-Tzeli S, Hahnemann JM, Hertz JM, Houge G, Kuklík M, Macek M, Lacombe D, Miller K, Moncla A, López Pajares I, Patsalis PC, Petersen MB: Origin of nondisjunction in trisomy 8 and trisomy 8 mosaicism. Eur J Hum Genet. 1998, 6: 432-438. 10.1038/sj.ejhg.5200212.

    Article  CAS  PubMed  Google Scholar 

  12. Veitia RA, Bottani S, Birchler JA: Cellular reactions to gene dosage imbalance: genomic, transcriptomic and proteomic effects. Trends Genet. 2008, 24: 390-397. 10.1016/j.tig.2008.05.005.

    Article  CAS  PubMed  Google Scholar 

  13. Virtaneva K, Wright FA, Tanner SM, Yuan B, Lemon WJ, Caligiuri MA, Bloomfield CD, de La Chapelle A, Krahe R: Expression profiling reveals fundamental biological differences in acute myeloid leukemia with isolated trisomy 8 and normal cytogenetics. Proc Natl Acad Sci U S A. 2001, 98: 1124-1129. 10.1073/pnas.98.3.1124.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Vey N, Mozziconacci MJ, Groulet-Martinec A, Debono S, Finetti P, Carbuccia N, Beillard E, Devilard E, Arnoulet C, Coso D, Sainty D, Xerri L, Stoppa AM, Lafage-Pochitaloff M, Nguyen C, Houlgatte R, Blaise D, Maraninchi D, Birg F, Birnbaum D, Bertucci F: Identification of new classes among acute myelogenous leukaemias with normal karyotype using gene expression profiling. Oncogene. 2004, 23: 9381-9391. 10.1038/sj.onc.1207910.

    Article  CAS  PubMed  Google Scholar 

  15. Schoch C, Kohlmann A, Dugas M, Kern W, Schnittger S, Haferlach T: Impact of trisomy 8 on expression of genes located on chromosome 8 in different AML subgroups. Genes Chromosomes Cancer. 2006, 45: 1164-1168. 10.1002/gcc.20380.

    Article  CAS  PubMed  Google Scholar 

  16. Chango A, Abdennebi-Najar L, Tessier F, Ferré S, Do S, Guéant JL, Nicolas JP, Willequet F: Quantitative methylation-sensitive arbitrarily primed PCR method to determine differential genomic DNA methylation in Down Syndrome. Biochem Biophys Res Commun. 2006, 349: 492-496. 10.1016/j.bbrc.2006.08.038.

    Article  CAS  PubMed  Google Scholar 

  17. Kerkel K, Schupf N, Hatta K, Pang D, Salas M, Kratz A, Minden M, Murty V, Zigman WB, Mayeux RP, Jenkins EC, Torkamani A, Schork NJ, Silverman W, Croy BA, Tycko B: Altered DNA methylation in leukocytes with trisomy 21. PLoS Genet. 2010, 6: e1001212-10.1371/journal.pgen.1001212.

    Article  PubMed Central  PubMed  Google Scholar 

  18. Paulsson K, Johansson B: High hyperdiploid childhood acute lymphoblastic leukemia. Genes Chromosomes Cancer. 2009, 48: 637-660. 10.1002/gcc.20671.

    Article  CAS  PubMed  Google Scholar 

  19. Davidsson J, Lilljebjorn H, Andersson A, Veerla S, Heldrup J, Behrendtz M, Fioretos T, Johansson B: The DNA methylome of pediatric acute lymphoblastic leukemia. Hum Mol Genet. 2009, 18: 4054-4065. 10.1093/hmg/ddp354.

    Article  CAS  PubMed  Google Scholar 

  20. Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, Braisted J, Klapa M, Currier T, Thiagarajan M, Sturn A, Snuffin M, Rezantsev A, Popov D, Ryltsov A, Kostukovich E, Borisovsky I, Liu Z, Vinsavich A, Trush V, Quackenbush J: TM4: a free, open-source system for microarray data management and analysis. Biotechniques. 2003, 34: 374-378.

    CAS  PubMed  Google Scholar 

  21. Torres EM, Sokolsky T, Tucker CM, Chan LY, Boselli M, Dunham MJ, Amon A: Effects of aneuploidy on cellular physiology and cell division in haploid yeast. Science. 2007, 317: 916-924. 10.1126/science.1142210.

    Article  CAS  PubMed  Google Scholar 

  22. Sheltzer JM, Torres EM, Dunham MJ, Amon A: Transcriptional consequences of aneuploidy. Proc Natl Acad Sci U S A. 2012, 109: 12644-12649. 10.1073/pnas.1209227109.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Sheltzer JM, Blank HM, Pfau SJ, Tange Y, George BM, Humpton TJ, Brito IL, Hiraoka Y, Niwa O, Amon A: Aneuploidy drives genomic instability in yeast. Science. 2011, 333: 1026-1030. 10.1126/science.1206412.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Birchler JA: Aneuploidy in plants and flies: the origin of studies of genomic imbalance. Semin Cell Dev Biol. 2013, 24: 315-319. 10.1016/j.semcdb.2013.02.004.

    Article  CAS  PubMed  Google Scholar 

  25. Mao R, Wang X, Spitznagel EL, Frelin LP, Ting JC, Ding H, Kim JW, Ruczinski I, Downey TJ, Pevsner J: Primary and secondary transcriptional effects in the developing human Down syndrome brain and heart. Genome Biol. 2005, 6: R107-10.1186/gb-2005-6-13-r107.

    Article  PubMed Central  PubMed  Google Scholar 

  26. Li CM, Guo M, Salas M, Schupf N, Silverman W, Zigman WB, Husain S, Warburton D, Thaker H, Tycko B: Cell type-specific over-expression of chromosome 21 genes in fibroblasts and fetal hearts with trisomy 21. BMC Med Genet. 2006, 7: 24.

    Article  PubMed Central  PubMed  Google Scholar 

  27. Prandini P, Deutsch S, Lyle R, Gagnebin M, Delucinge Vivier C, Delorenzi M, Gehrig C, Descombes P, Sherman S, Dagna Bricarelli F, Baldo C, Novelli A, Dallapiccola B, Antonarakis SE: Natural gene-expression variation in Down syndrome modulates the outcome of gene-dosage imbalance. Am J Hum Genet. 2007, 8: 252-263.

    Article  Google Scholar 

  28. Rozovski U, Jonish-Grossman A, Bar-Shira A, Ochshorn Y, Goldstein M, Yaron Y: Genome-wide expression analysis of cultured trophoblast with trisomy 21 karyotype. Hum Reprod. 2007, 22: 2538-2545. 10.1093/humrep/dem214.

    Article  CAS  PubMed  Google Scholar 

  29. FitzPatrick DR, Ramsay J, McGill NI, Shade M, Carothers AD, Hastie ND: Transcriptome analysis of human autosomal trisomy. Hum Mol Genet. 2002, 11: 3249-3256. 10.1093/hmg/11.26.3249.

    Article  CAS  PubMed  Google Scholar 

  30. Altug-Teber Ö, Bonin M, Walter M, Mau-Holzmann UA, Dufke A, Stappert H, Tekesin I, Heilbronner H, Nieselt K, Riess O: Specific transcriptional changes in human fetuses with autosomal trisomies. Cytogenet Genome Res. 2007, 119: 171-184. 10.1159/000112058.

    Article  CAS  PubMed  Google Scholar 

  31. Biancotti JC, Narwani K, Buehler N, Mandefro B, Golan-Lev T, Yanuka O, Clark A, Hill D, Benvenisty N, Lavon N: Human embryonic stem cells as models for aneuploid chromosomal syndromes. Stem Cells. 2010, 28: 1530-1540. 10.1002/stem.483.

    Article  CAS  PubMed  Google Scholar 

  32. Schinzel A: Catalogue of unbalanced chromosome aberrations In. 2001, New York: Walter de Gruyter

    Google Scholar 

  33. Hellman A, Chess A: Gene body-specific methylation on the active X chromosome. Science. 2007, 315: 1141-1143. 10.1126/science.1136352.

    Article  CAS  PubMed  Google Scholar 

  34. Yasukochi Y, Maruyama O, Mahajan MC, Padden C, Euskirchen GM, Schulz V, Hirakawa H, Kuhara S, Pan XH, Newburger PE, Snyder M, Weissman SM: X chromosome-wide analyses of genomic DNA methylation states and gene expression in male and female neutrophils. Proc Natl Acad Sci U S A. 2010, 107: 3704-3709. 10.1073/pnas.0914812107.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Sharp AJ, Stathaki E, Migliavacca E, Brahmachary M, Montgomery SB, Dupre Y, Antonarakis SE: DNA methylation profiles of human active and inactive X chromosomes. Genome Res. 2011, 21: 1592-1600. 10.1101/gr.112680.110.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Szulwach KE, Li X, Li Y, Song CX, Wu H, Dai Q, Irier H, Upadhyay AK, Gearing M, Levey AI, Vasanthakumar A, Godley LA, Chang Q, Cheng X, He C, Jin P: 5-hmC-mediated epigenetic dynamics during postnatal neurodevelopment and aging. Nat Neurosci. 2011, 14: 1607-1616. 10.1038/nn.2959.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Weber M, Davies JJ, Wittig D, Oakeley EJ, Haase M, Lam WL, Schübeler D: Chromosome-wide and promoter-specific analyses identify sites of differential DNA methylation in normal and transformed human cells. Nat Genet. 2005, 37: 853-862. 10.1038/ng1598.

    Article  CAS  PubMed  Google Scholar 

  38. Gilbert N, Thomson I, Boyle S, Allan J, Ramsahoye B, Bickmore WA: DNA methylation affects nuclear organization, histone modifications, and linker histone binding but not chromatin compaction. J Cell Biol. 2007, 177: 401-411. 10.1083/jcb.200607133.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  39. Branco MR, Ficz G, Reik W: Uncovering the role of 5-hydroxymethylcytosine in the epigenome. Nat Rev Genet. 2012, 13: 7-13.

    CAS  Google Scholar 

  40. Stroud H, Feng S, Morey Kinney S, Pradhan S, Jacobsen SE: 5-Hydroxymethylcytosine is associated with enhancers and gene bodies in human embryonic stem cells. Genome Biol. 2011, 12: R54-10.1186/gb-2011-12-6-r54.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  41. Ficz G, Branco MR, Seisenberger S, Santos F, Krueger F, Hore TA, Marques CJ, Andrews S, Reik W: Dynamic regulation of 5-hydroxymethylcytosine in mouse ES cells and during differentiation. Nature. 2011, 473: 398-402. 10.1038/nature10008.

    Article  CAS  PubMed  Google Scholar 

  42. Pastor WA, Pape UJ, Huang Y, Henderson HR, Lister R, Ko M, McLoughlin EM, Brudno Y, Mahapatra S, Kapranov P, Tahiliani M, Daley GQ, Liu XS, Ecker JR, Milos PM, Agarwal S, Rao A: Genome-wide mapping of 5-hydroxymethylcytosine in embryonic stem cells. Nature. 2011, 473: 394-397. 10.1038/nature10102.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  43. Nestor CE, Ottaviano R, Reddington J, Sproul D, Reinhardt D, Dunican D, Katz E, Dixon JM, Harrison DJ, Meehan RR: Tissue type is a major modifier of the 5-hydroxymethylcytosine content of human genes. Genome Res. 2012, 22: 467-477. 10.1101/gr.126417.111.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  44. Kinney SM, Chin HG, Vaisvila R, Bitinaite J, Zheng Y, Esteve PO, Feng S, Stroud H, Jacobsen SE, Pradhan S: Tissue-specific distribution and dynamic changes of 5-hydroxymethylcytosine in mammalian genomes. J Biol Chem. 2011, 286: 24685-24693. 10.1074/jbc.M110.217083.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  45. Kudo Y, Tateishi K, Yamamoto K, Yamamoto S, Asaoka Y, Ijichi H, Nagae G, Yoshida H, Aburatani H, Koike K: Loss of 5-hydroxymethylcytosine is accompanied with malignant cellular transformation. Cancer Sci. 2012, 103: 670-676. 10.1111/j.1349-7006.2012.02213.x.

    Article  CAS  PubMed  Google Scholar 

  46. Jin SG, Jiang Y, Qiu R, Rauch TA, Wang Y, Schackert G, Krex D, Lu Q, Pfeifer GP: 5-Hydroxymethylcytosine is strongly depleted in human cancers but its levels do not correlate with IDH1 mutations. Cancer Res. 2011, 71: 7360-7365. 10.1158/0008-5472.CAN-11-2023.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  47. Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, Franz M, Grouios C, Kazi F, Lopes CT, Maitland A, Mostafavi S, Montojo J, Shao Q, Wright G, Bader GD, Morris Q: The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010, 38: W214-W220. 10.1093/nar/gkq537.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  48. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.

    Article  CAS  PubMed  Google Scholar 

  49. Davidsson J, Andersson A, Paulsson K, Heidenblad M, Isaksson M, Borg Ã…, Heldrup J, Behrendtz M, Panagopoulos I, Fioretos T, Johansson B: Tiling resolution array comparative genomic hybridization, expression and methylation analyses of dup(1q) in Burkitt lymphomas and pediatric high hyperdiploid acute lymphoblastic leukemias reveal clustered near-centromeric breakpoints and overexpression of genes in 1q22-32.3. Hum Mol Genet. 2007, 16: 2215-2225. 10.1093/hmg/ddm173.

    Article  CAS  PubMed  Google Scholar 

  50. Scacheri PC, Crawford GE, Davis S: Statistics for ChIP-chip and DNase hypersensitivity experiments on NimbleGen arrays. Methods Enzymol. 2006, 411: 270-282.

    Article  CAS  PubMed  Google Scholar 

  51. Lingjaerde OC, Baumbusch LO, Liestøl K, Glad IK, Børresen-Dale A-L: CGH-Explorer: a program for analysis of array-CGH data. Bioinformatics. 2005, 21: 821-822. 10.1093/bioinformatics/bti113.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgments

This study was funded by grants from the Erik-Philip Sörensen Foundation, the Magnus Bergvalls Foundation, the Swedish Cancer Society, the Swedish Childhood Cancer Foundation, and the Swedish Research Council.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Josef Davidsson.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

JD and BJ participated in the concept and writing of the manuscript. JD executed specific experiments. JD, SV and BJ were involved in the analysis of the data. All authors have read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Table S1: Case reports on hematological disorders/malignancies in patients with CT8M. (DOC 42 KB)

13072_2013_289_MOESM2_ESM.doc

Additional file 2: Figure S1: Unsupervised HCA of gene and miRNA expression reveals accurate clustering of the trisomy 8 subgroup. (A) HCA of global gene expression in the disomy 8, trisomy 8, and reference cultures clustered trisomy 8 in the same branch. (B) A similar pattern was observed when co-analyzing gene and miRNA expression. (C) On the other hand, HCA of only global miRNA expression did not cluster the different subgroups accurately. (D) After long-term culturing of one trisomy 8 culture (6lt) and one disomy 8 culture (1lt), HCA of global gene and miRNA expression accurately grouped these together with their short-term cultured counterparts. (DOC 550 KB)

13072_2013_289_MOESM3_ESM.doc

Additional file 3: Figure S2: The majority of genes on chromosome 8 in the trisomy 8 cultures are upregulated compared with the disomy and reference cultures. Median-centered expression of genes on chromosome 8 demonstrates that the majority of these genes in the trisomy 8 cultures are upregulated compared with (A) the reference cultures as well as compared with both the disomy 8 and (B) reference cultures. (DOC 1 MB)

13072_2013_289_MOESM4_ESM.docx

Additional file 4: Figure S3: Validation of expression array results by real-time quantitative PCR (qPCR). Seven genes were selected from Table 1 and analyzed with qPCR using commercial Taq-Man probes. The expression levels are presented as a fold change for each culture using TBP as endogenous control. The qPCR analyses confirmed the over- or underexpression of the seven genes in the trisomy 8-positive cultures compared with the disomy 8 and reference cultures. (DOCX 1 MB)

13072_2013_289_MOESM5_ESM.doc

Additional file 5: Table S2: Ingenuity pathway analysis of 502 differentially expressed genes in the trisomy 8 and disomy 8 cultures. (DOC 37 KB)

13072_2013_289_MOESM6_ESM.doc

Additional file 6: Table S3: Hypermethylated and hyperhydroxymethylated promoters/CpG islands in relation to gene expression. (DOC 30 KB)

13072_2013_289_MOESM7_ESM.doc

Additional file 7: Figure S4: Promoter-specific methylation levels of autosomes do not differ between the trisomy 8 and disomy 8/refererence cultures. The levels of average promoter-specific methylation on chromosomes (A) 2, (B) 6, (C) 7, and (D) 8 in the trisomy 8 cultures and in the disomy 8 and reference cultures combined were log2-converted, median-centered, and plotted against genomic positions. No significant differences between the culture groups were observed. (DOC 1 MB)

13072_2013_289_MOESM8_ESM.doc

Additional file 8: Figure S5: Promoter-specific hydroxymethylation levels of chromosome 8, but not of the other autosomes, is significantly lower in the trisomy 8 compared with the disomy 8/reference cultures. The levels of average promoter-specific hydroxymethylation on chromosomes (A) 2, (B) 6, (C) 7, and (D) 8 in the trisomy 8 cultures and in the disomy 8 and references cultures combined cultures were log2-converted, median-centered, and plotted against genomic positions. Significantly lower levels (P < 0.05; t- test) were observed only for chromosome 8 in the trisomy 8 cultures. (DOC 3 MB)

13072_2013_289_MOESM9_ESM.doc

Additional file 9: Figure S6: The X chromosome is significantly less hydroxymethylated compared with all autosomes. Measuring and comparing the average intensity ratios of all autosomes combined with the X chromosome, the promoters/CpG islands on the latter chromosome were generally (irrespective of gender and culture group) significantly (P <0.05; t- test) less hydroxymethylated. (DOC 158 KB)

Additional file 10: Table S4: Candidate genes for clinical manifestations. (DOC 44 KB)

13072_2013_289_MOESM11_ESM.doc

Additional file 11: Figure S7: Functional association analyses of cardiovascular dysfunction, cancer, and neurological development and function reveal widespread interactions between the candidate genes. Analyses of association data including protein and genetic interactions, pathways, co-expression, co-localization, and protein domain similarity of candidate genes derived from Table 1 reveal widespread interactions between gene sets associated with (A) cardiovascular dysfunction, (B) cancer, and (C) neurological development and function. Circles in grey represent candidate genes (Table 1), whereas white circles represent associated genes identified by the GeneMANIA analyses. Since the analysis is three dimensional, the sizes of the circles vary in this two dimensional representation. The colors of the lines connecting the various genes denote the different types of functional associations, as indicated by the color boxes. (DOC 594 KB)

13072_2013_289_MOESM12_ESM.doc

Additional file 12: Figure S8: Schematic overview of single-cell cloning of cells disomic and trisomic for chromosome 8, respectively. From an original cell line from a patient with CT8M, with the karyotype 47,XY,+8[5]/46,XY[20], a total of three cultures with disomy 8 and three cultures with trisomy 8 were generated. DNA and RNA from these cultures were used for subsequent epigenetic analyses. In addition, two commercially available control cell lines, with normal male and female karyotypes, respectively, were included as references in all analyses. (DOC 78 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Davidsson, J., Veerla, S. & Johansson, B. Constitutional trisomy 8 mosaicism as a model for epigenetic studies of aneuploidy. Epigenetics & Chromatin 6, 18 (2013). https://doi.org/10.1186/1756-8935-6-18

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1756-8935-6-18

Keywords