Single-cell sequencing reveals karyotype heterogeneity in murine and human malignancies
- Bjorn Bakker†1,
- Aaron Taudt†1,
- Mirjam E. Belderbos1, 3,
- David Porubsky1,
- Diana C. J. Spierings1,
- Tristan V. de Jong1,
- Nancy Halsema1,
- Hinke G. Kazemier1,
- Karina Hoekstra-Wakker1,
- Allan Bradley2,
- Eveline S. J. M. de Bont3,
- Anke van den Berg4,
- Victor Guryev1,
- Peter M. Lansdorp1, 5, 6,
- Maria Colomé-Tatché1, 7Email author and
- Floris Foijer1Email author
© Bakker et al. 2016
Received: 29 December 2015
Accepted: 29 April 2016
Published: 31 May 2016
Chromosome instability leads to aneuploidy, a state in which cells have abnormal numbers of chromosomes, and is found in two out of three cancers. In a chromosomal instable p53 deficient mouse model with accelerated lymphomagenesis, we previously observed whole chromosome copy number changes affecting all lymphoma cells. This suggests that chromosome instability is somehow suppressed in the aneuploid lymphomas or that selection for frequently lost/gained chromosomes out-competes the CIN-imposed mis-segregation.
To distinguish between these explanations and to examine karyotype dynamics in chromosome instable lymphoma, we use a newly developed single-cell whole genome sequencing (scWGS) platform that provides a complete and unbiased overview of copy number variations (CNV) in individual cells. To analyse these scWGS data, we develop AneuFinder, which allows annotation of copy number changes in a fully automated fashion and quantification of CNV heterogeneity between cells. Single-cell sequencing and AneuFinder analysis reveals high levels of copy number heterogeneity in chromosome instability-driven murine T-cell lymphoma samples, indicating ongoing chromosome instability. Application of this technology to human B cell leukaemias reveals different levels of karyotype heterogeneity in these cancers.
Our data show that even though aneuploid tumours select for particular and recurring chromosome combinations, single-cell analysis using AneuFinder reveals copy number heterogeneity. This suggests ongoing chromosome instability that other platforms fail to detect. As chromosome instability might drive tumour evolution, karyotype analysis using single-cell sequencing technology could become an essential tool for cancer treatment stratification.
KeywordsAneuploidy Karyotype heterogeneity Single-cell sequencing Copy number detection Lymphoma Leukaemia
Chromosomal instability (CIN) is a process leading to structural and whole chromosome abnormalities and results in cells with an abnormal DNA content, a state defined as aneuploid. CIN and the resulting aneuploidy cause physiological stress and growth defects in yeast and primary mouse embryonic fibroblasts [1, 2]. Furthermore, some of the mouse models that were engineered to model CIN are characterised with a reduced lifespan, which can be rescued by reducing the levels of CIN [3–5]. Although aneuploidy has detrimental consequences for untransformed cells, more than two out of three cancers are aneuploid, suggesting a fundamental relationship between aneuploidy and tumourigenesis that so far remains poorly understood [6–8].
Various mouse models have been engineered to investigate the relationship between aneuploidy and tumourigenesis, typically by abrogating components of the spindle assembly checkpoint (SAC) [9–11]. The SAC monitors chromosome segregation by blocking mitotic cells in metaphase until proper kinetochore-microtubule attachment and tension have been established [12–14], and its (partial) abrogation will therefore lead to CIN. We have previously shown that CIN as provoked by truncating the SAC kinase Mps1 selectively in murine T-cells is not sufficient to trigger malignant proliferation. However, when combined with loss of p53, CIN significantly accelerated lymphomagenesis. When we analysed the average chromosome content in the emerging T cell acute lymphoblastic lymphomas (T-ALLs) by array comparative genomic hybridisation (aCGH), we found that all lymphomas displayed highly similar karyotypes, suggesting clonal selection for the recurring chromosomal abnormalities . This was surprising because we also showed that Mps1 truncation results in mitotic abnormalities in virtually each cell division in tissue culture cells, which would counteract clonal selection of any chromosomal abnormality in a tumour. One possible explanation is that aneuploid T-ALL cells can overcome the CIN phenotype and thus maintain a stable karyotype. Alternatively, selection forces for specific chromosome alterations outcompete the mis-segregation events. The resulting chromosomal imbalances are expected to have severe effects on gene expression and thus cellular fitness. Although most forms of aneuploidy are expected to decrease fitness, karyotype heterogeneity could result in selection of cells that have accumulated favourable copy numbers of chromosomes expressing genes important for tumour evolution and overall cellular fitness.
Indeed, karyotype heterogeneity occurs in human cancers [16, 17], has been associated with human tumour evolution, and might therefore impact therapeutic response to cancer therapy . However, most current cytogenetic and molecular techniques are limited in the number of karyotype alterations per cell they can detect, are biased towards dividing subpopulations, or can only measure the population-average chromosome copy number alterations [19, 20]. These shortcomings have precluded thorough analysis of intratumour chromosome copy number variations. Recent advances in single-cell genomics allow researchers to dissect the heterogeneity of the cancer genome with greater resolution than ever before [21, 22].
In this study we describe the application of single-cell whole genome sequencing (scWGS) to map and quantify karyotype heterogeneity in primary mouse lymphoma and human leukaemia samples. For this purpose, we have developed a new bioinformatics tool, AneuFinder, to identify chromosome copy numbers in scWGS data. Using these tools, we now show that besides recurrent chromosome copy number alterations, the aneuploid T-ALLs that arise in our mouse model display high-grade ongoing CIN as evidenced by severe intratumour karyotype heterogeneity. Importantly, analyses of a number of human paediatric B cell acute lymphoblastic leukaemia (B-ALL) samples using this platform revealed different grades of karyotype heterogeneity, demonstrating that CIN rates differ between these malignancies. Therefore, as ongoing CIN is an important hypothesised driver of cancer evolution , our platform might become an important tool to predict treatment outcome or even treatment stratification.
Chromosomal instable T-ALLs show recurrent chromosome copy number changes
Traditional methods to examine karyotypes depend on dividing cells (in case of normal and spectral karyotyping [SKY]), or are limited in the number of chromosomes that can be quantified per cell (in case of interphase FISH). An alternative to measure copy number alterations in a tumour is to measure the average DNA content (e.g. by aCGH) [19, 23], but this obscures intratumour heterogeneity. We therefore moved to single-cell sequencing as a method for karyotyping, making use of a modified scWGS protocol, described in more detail in van den Bos et al. . Briefly, this scWGS platform involves single-cell sorting of primary tumour cells as nuclei by flow cytometry, followed by automated DNA fragmentation, barcoded next generation sequencing library preparation and shallow multiplexed sequencing . To validate our platform, we first sequenced the genomes of 25 primary T-ALL cells isolated from an Mps1 f/f ; p53 f/f ; Lck-Cre + lymphoma that we had previously assessed  by aCGH-analysis (T-ALL 1, aCGH data in Fig. 1a). We first compared the single-cell sequencing data to the existing aCGH data by creating an artificial ‘bulk sequencing file’ that has the cumulative data of all individual single-cell sequencing libraries (Fig. 1b) to determine how representative the sampled cells are for the bulk tumour. Indeed, we found that the copy number changes in the bulk sequencing analysis were identical to those observed in the aCGH data (compare Fig. 1a and b). However, when we plotted the individual single-cell libraries, we detected many additional copy number changes: while some of the single cells showed the exact karyotype as found by aCGH analysis (Ts2, Ts4, Del7, Ts9, Ts14, Ts15; Fig. 1c; cell 1, more examples in Additional file 2: Figure S2), most cells displayed additional chromosomal aberrations (Fig. 1c; compare cells 2, 3 and 4 to cell 1, more single-cell libraries in Additional file 2: Figure S2), reflecting karyotype heterogeneity and suggestive of ongoing CIN. Indeed, when we manually annotated the individual karyotypes of all 25 cells, we found that 56 % of the cells had a unique karyotype (Additional file 3: Figure S3), further emphasising the heterogeneity that our earlier aCGH analysis had failed to detect.
AneuFinder: a tool to analyse high throughput single-cell sequencing data
scWGS and AneuFinder reveal karyotype heterogeneity in murine CIN T-ALLs
To determine aneuploidy at the single-cell level, we selected four Mps1 f/f ; p53 f/f ; Lck-Cre + T-ALLs, one Mps1 f/f ; p53 f/+ ; Lck-Cre + T-ALL (Fig. 3a, arrows) and one control thymus (Mps1 f/f ; p53 f/f ; Lck-Cre − ) for single-cell sequencing. For each sample, we sequenced 48 single-cell libraries. While AneuFinder did not detect any chromosome copy number alterations in the male Lck-Cre− control thymus (Fig. 3c), all five T-ALL samples revealed clonal copy number gains of chromosomes 4, 9, 14 and 15 (Fig. 3d–h), highly similar as identified in Mps1; p53; Lck-Cre T-ALLs assessed by aCGH analysis (Fig. 1a and ). However, more importantly, AneuFinder also detected many other chromosome copy number variations that were present in a minority of the lymphoma cells that we predict aCGH analysis would have failed to detect. Indeed, these additional chromosome copy number alterations disappeared when we analysed our single - cell sequencing in bulk using the ‘bulk analysis sequencing files’ similar to Fig. 1b (Additional file 8: Figure S5). Furthermore, AneuFinder also revealed heterogeneity in the copy number states for the frequently gained chromosomes that we failed to detect previously by aCGH . For instance, in tumours T158, T170 and T257 half of the cells had three copies of chromosome 14, while the other half had four or even more copies (compare Fig. 3e–g). All together, these data clearly show that Mps1 f/f ; p53 f/f ; Lck-Cre + T-ALLs are not only highly aneuploid, but also display high-grade karyotype heterogeneity, suggestive of ongoing CIN.
Quantifying karyotype heterogeneity in CIN T-ALLs
Since different chromosomes showed varying degrees of gain or loss in the T-ALLs (Fig. 3d–h and Additional file 9: Figure S6b, c), we wondered whether calculating the aneuploidy and heterogeneity scores for individual chromosomes would reveal whether specific chromosomes more often showed changes in copy number than others. To this end, we plotted both scores per chromosome for all samples that were analysed by single-cell sequencing. For the control thymus, all chromosomes clustered together in the bottom left, indicating that none of the cells displayed chromosome copy number alterations (Fig. 4b, control thymus). In contrast, in the tumours we identified three ‘types’ of chromosomes: (1) chromosomes that were (virtually) never lost or gained (Fig. 4b, green chromosomes in T260 and T158), presumably due to lethality associated with such gain/loss events; (2) chromosomes that show a high heterogeneity rate, but low aneuploidy rate, for which copy number changes are presumably not selected for but that are tolerated relatively well (Fig 4b, blue chromosomes); and (3) chromosomes with high aneuploidy, but low(er) heterogeneity scores, for which the observed copy number changes are selected for in the tumours (Fig 4b, red chromosomes). These results (additional plots in Additional file 10: Figure S7) indicate that while some chromosomes copy number changes are favoured and others are not tolerated in these tumours, there is a third group of chromosomes for which copy number changes contribute to heterogeneity, but that are not (yet) selected for.
Even though we observed high levels of karyotype heterogeneity, the single-cell sequencing experiments are endpoint measurements and we therefore formally could not rule out that the T-ALLs consisted of a large number of different aneuploid but chromosomal stable clones. To determine whether chromosome mis-segregation events still occurred in the primary malignancies, we derived primary lymphoma cell lines from three endpoint T-ALLs. This yielded three viable primary murine T-ALL lines (T302, T392 and T397) that displayed chromosome numbers similar as observed in primary lymphomas as assessed by metaphase spread-based chromosome counting (median 47.5–51 for the T-ALL lines [passage 8], compared to 46–50 for the primary T-ALLs [compare Fig. 4c to Fig. 3d–h]). Furthermore, the lines displayed variability in chromosome numbers, indicating karyotype heterogeneity (Fig. 4c). We then labelled primary T-ALL cell line 1 (T302) with an H2B-GFP fusion protein to monitor the actual rates of chromosome mis-segregation by time-lapse microscopy. Indeed, we found that half of the mitoses (n = 32) showed clear mitotic abnormalities, mostly lagging chromosomes and failed alignments, similar to those observed previously in Mps1 f/f ; p53 f/f mouse embryonic fibroblasts (Fig. 4d, e, Additional files 11, 12 and 13: Movies 1–3 ). We therefore conclude that the karyotype heterogeneity that we observed in developing and endpoint Mps1; p53; Lck-Cre lymphomas is the result of ongoing CIN.
However, as we observed roughly 10 % of the mitoses in our primary T-ALL lines to result in tetraploid cells (Fig. 4d, e) and as for single-cell sequencing purposes we typically sort near diploid (2n) cells, we next investigated whether this resulted in an underestimation of the heterogeneity in endpoint tumours. To test this, we first examined the fate of T-ALL cells undergoing binucleation or polyploidisation events by time-lapse imaging in primary murine T-ALL cultures. While all cells displaying polyploidisation or binucleation events died between 4 and 11 h after mitotic exit (n = 7, e.g. Additional file 14: Movie 4), only 20 % of the cells that showed no clear mitotic abnormalities died in this experiment (n = 10). Even though this indicates that tetraploidisation events are selected against in a tissue culture setting, this still could be different in vivo, given that whole genome duplication is known to occur in human cancers. However, as with our existing strains we could not distinguish G2 cells from G1 tetraploid cells by FACS, we therefore analysed aneuploid murine T-ALLs for cells with a DNA content larger than 4n as an indication for proliferating tetraploid cells. Indeed, two out of the four analysed tumours (T158 and T257) had cells with a DNA content larger than 4n (see arrows Additional file 15: Figure S8a) indicating that some of the cells with a 4n DNA content had to be tetraploid. We then sorted individual G1 and G2-diploid/G1-(near-)tetraploid cells followed by single-cell sequencing from T158 (Fig. 4f, g) and T257 (Fig. 4h, i), respectively. Even though in this experiment we could not distinguish G2 cells from (near) tetraploid cells, the data revealed that near tetraploid cells were present as some of the cells had odd copy numbers of chromosomes and G2 cells by definition must have an even number of chromosome copies. Importantly, this allowed us to confirm that AneuFinder can detect such events: when we forced AneuFinder to fit the 4n population as diploid, the resulting fit was poor for cells with odd numbers (and thus true near-tetraploid cells; compare left and right hand panels in Additional file 15: Figure S8b, c) indicating that AneuFinder can detect tetraploid events even without a priori knowledge about the most prevalent state (diploid/tetraploid), similar to Ginkgo.
Further analysis of the 4n population of T158 revealed that 12 out of 37 analysed cells had an odd copy number for at least one of the chromosomes, indicating that at least 32 % of the 4n cells were true near-tetraploid and not G2 cells. The heterogeneity score for these near-tetraploid cells (0.6503) was much larger than for the 2n population (0.3099), indicating that the near-tetraploid cells can generate greater karyotype diversity than near-diploid cells, presumably because of the larger number of copy number states available. For the other 4n cells (68 %) we could not discriminate whether these were G2 cells or G1 near-tetraploid cells in this experimental setup. However, the 4n population only represented 4.1 % of all the live tumour cells, and therefore only between 1.3 % (32 % of 4.1 % near-tetraploid cells with odd numbers of chromosomes) and 4.1 % (all cells with near 4n DNA content) in tumour T158 were near-tetraploid cells. The contribution of these cells to intratumour heterogeneity is therefore limited and in agreement with the observed cell death events in the time-lapse experiments in primary aneuploid T-ALL cultures (Additional file 14: Movie 4). In the T257 4n population (5.3 % of all tumour cells) we found two out of 34 analysed cells to have odd chromosome numbers, representing ~0.3 % of the total tumour. Therefore, the contribution of genuine near-tetraploid cells to tumour T257 lies between 0.3 % and 5.3 %. We conclude that while the near-tetraploid cells contribute to intratumour heterogeneity, this contribution is limited, presumably due to cell death events after polyploidisation leading to low fractions of near-tetraploid cells in the primary aneuploid T-ALLs.
Karyotype heterogeneity in human B cell acute lymphoblastic leukaemia
To assess whether the observed karyotype heterogeneity in the primary human B-ALL samples could be indicative for ongoing chromosomal instability, we next engrafted B-ALL B into sub-lethally irradiated immune-deficient mice. Mice were sampled every four weeks to assess the levels of chimerism (the fraction of human-specific CD45+ cells in the peripheral blood; Additional file 17: Figure S10), and were sacrificed when these levels became larger than 90 %. Interestingly, when we analysed the bone marrow of these mice (B-ALL B-1 and B-2) using scWGS and AneuFinder, we found that leukaemic cell karyotypes had changed after transplantation. In both transplanted samples, a fraction of cells gained chromosome 2 while none of the primary B-ALL B cells displayed chromosome 2 trisomy, and one cell in transplant B2 showed unique gain of chromosome 10 (Fig. 5c). In addition, the cells that had gained chromosome 2 had lost the extra copy of chromosome 9. This could be the result of a copy neutral loss of heterozygosity (cnLOH) event. To address this, we extracted all single nucleotide polymorphisms (SNPs) on chromosome 9 that were identified by our single-cell sequencing analysis in the cells harbouring two copies of chromosome 9. Unfortunately, we found that the coverage of the single-cell libraries was too low to sufficiently sample alternative alleles at defined SNP positions. To detect such events deeper sequencing is required (which increases cost price).
Heterogeneity scores for both transplants revealed an increased heterogeneity score for both transplanted leukaemias B1 and B2 (Fig. 5d, an increase from 0.0176 [black diamond] to 0.0565 and 0.0398 [blue triangles], respectively). Even though we formally cannot rule out that the cells with ‘new’ karyotype patterns were pre-existing in the primary leukaemia as low-abundant clones, these data strongly suggest that human tumours also exhibit ongoing CIN that can facilitate tumour evolution when tumour cells are exposed to stress, in this case survival and propagation in another niche (i.e. the mouse hematopoietic system).
Discussion and conclusions
Aneuploidy is a hallmark of cancer cells and results from chromosome missegregation events during mitosis, a process called chromosomal instability (CIN). Indeed, most cancer cell lines exhibit CIN with ~10 % to 60 % of the cells displaying lagging chromosomes in mitosis, resulting in a distribution of cells with various aneuploidies [28–30]. However, while measuring CIN rates and the resulting karyotype variation in cancer cell lines is straight-forward, assessing karyotype heterogeneity in primary tumours is technically challenging . Only recently, the first platforms to measure full karyotypes of single non-dividing primary cells have been reported [22, 23, 31].
In this study we describe a whole genome single-cell sequencing pipeline to quantify aneuploidy in primary (tumour) cells in an unbiased and high throughput fashion. We use this pipeline to measure karyotype heterogeneity in aneuploid murine T-ALLs that arise in our earlier published Mps1; p53; Lck-Cre mouse model as a measure for ongoing CIN . We showed that our pipeline confirms our earlier observations that Mps1; p53; Lck-Cre lymphomas display recurrent chromosome copy number alterations for some chromosomes. In addition, we observed groups of chromosomes to show distinct degrees of aneuploidy and heterogeneity, with specific chromosomes favouring either aneuploidy or euploidy, and a third group of chromosomes showing apparent random copy number changes. This latter group of chromosomes could be important for tumour evolution as CIN events for these chromosomes are not selected against in the primary tumours. However, these copy number changes might become beneficial when the tumour environment changes, for instance during therapy and metastasis and thus lead to tumour recurrence or disease progression. It will therefore be extremely interesting to measure chromosome dynamics in primary tumours, metastases and recurring tumours in future experiments.
Using AneuFinder, we were able to detect high-grade karyotype heterogeneity that would likely have been obscured using other platforms such as array CGH. This is important for two reasons: (1) it confirms that chromosome mis-segregation is ongoing in the T-ALLs that arise in our Mps1; p53; Lck-Cre mice despite clonal selection for some chromosome number alterations and therefore that the selection forces outcompete mis-segregation rates. Indeed, we confirm ongoing CIN by showing that 50 % of the primary T-ALL cell divisions display mitotic abnormalities; and (2) it shows that karyotype heterogeneity as a result from high grade CIN will go unnoticed when using one of the most commonly employed platforms to assess (aCGH). Therefore, karyotype heterogeneity might be a gravely underestimated feature of cancer. Indeed, when we assessed karyotype heterogeneity in three paediatric B-ALL samples we found that heterogeneity greatly varied between the three assessed samples. The heterogeneity was not directly predictive of tumour outcome, as the two tumours with the lowest and highest heterogeneity rate had a favourable outcome while the tumour with medium CIN rates had a poor outcome. However, even though our sample size is too small to draw any conclusions, it is tempting to hypothesise that tumours benefit most from medium CIN rates as low CIN rates do not allow for any karyotype evolution and too high CIN rates prevent clonal selection of chromosomes when the environment has changed. This hypothesis appears to hold true in the mouse , but further studies comparing human tumours with known outcome are required to test this. Importantly, when we xenografted the human aneuploid B-ALL with a medium CIN rate into mice, we found that its propagation led to additional karyotypic changes, presumably the result of propagation in another niche in line with our hypothesis. Therefore, we argue that ongoing CIN and the resulting karyotypic heterogeneity are important factors to consider when predicting disease progression/outcome and potentially treatment response.
Genomic heterogeneity is a common feature of many cancers, but because of technical limitations, it is rarely properly quantified. Instead, primary tumour cell population averages are used for research and diagnostics purposes that may not entirely represent the diversity in the primary tumour . The use of single-cell sequencing technology allows us to better estimate this level of tumour heterogeneity [22, 33, 34]. High resolution quantification of this heterogeneity is of the utmost importance as karyotype heterogeneity will drive tumour evolution and ultimately affect response to treatment [34–37]. While the mutational landscape of many tumour types has been characterised in detail, we have only begun to understand the role of karyotype heterogeneity in tumour evolution [18, 33, 38], a role that can be further examined using the scWGS-AneuFinder pipeline.
Animal housing and experiments
Animals used in this study had mixed C57BL/6 genetic backgrounds. All mice were bred in the Central Animal Facility (University Medical Centre Groningen [UMCG], Groningen, The Netherlands). The conditional deletion of Mps1 was described before . For survival studies, mice were monitored for tumour development starting at the age of 2.5 months by looking for signs of dyspnoea (a consequence of thymic hypertrophy) and general animal wellbeing.
Animal protocols were approved by the UMCG Committee on Animal Care (DEC) (RUG-DEC-6369). Human samples (bone marrow isolated from patients with acute lymphoblastic leukaemia) were collected as part of routine diagnostic procedures and leftover cells were used in this study (METc 2013.281). The use of these leftover samples was approved by the Medical Ethical Committee of the University Medical Center Groningen when patients (or guardians) had signed an informed consent letter. All patients and their guardians provided this written informed consent. All experimental methods comply with the Helsinki declaration for medical research involving human subjects.
Array CGH data analysis
For the aCGH data analysis, probe signals were divided by corresponding signals of a euploid control reference and log2-transformed. The data were binned into 1 Mb segments in steps of 500 Kb and plotted. Array CGH data have previously been deposited in the National Center for Biotechnology Information Gene Expression Omnibus (NCBI GEO) database under accession no. GSE57334 .
Single-cell whole genome sequencing and data processing
Single cells for the purpose of single-cell sequencing were isolated from thymuses by dissecting the organ and mincing the tissue through a 100 μm cell strainer. Cells were washed in PBS twice before stored in FBS with 10 % DMSO at −80 °C until sorted. Prior to sorting, cells were thawed, and incubated for 10 min on ice in the dark in a DNA-staining buffer containing 10 μg/mL propidium iodide and 10 μg/mL Hoechst, as well as 0.1 % Nonidet P-40 to dissociate the cytoplasm. The G1 population was sorted as individual nuclei into a 96-well plate format. Per sample, 10 nuclei were sorted into a single well to serve as a positive control, as well as an empty well acting as a negative control. Library preparation was performed using a modified scWGS protocol [24, 31] using a Bravo Automated Liquid Handling Platform (Agilent Technologies, Santa Clara, CA, USA). Clusters for sequencing were generated on the cBot (Illumina). Single-end 50 bp sequencing was performed on an Illumina HiSeq 2500 at ERIBA (Illumina, San Diego, CA, USA). Raw sequencing data were demultiplexed based on library-specific barcodes and converted to fastq format using standard Illumina software (bcl2fastq version 1.8.4). The resulting reads were mapped to mouse (mm10) or human (hg19) reference genome using Bowtie2 (version 2.2.4). Duplicate reads were marked using BamUtil (version 1.0.3). All single-cell sequencing data have been deposited at ArrayExpress under accession no. E-MTAB-4183.
Copy number analysis using AneuFinder
Analysis of single-cell sequencing data was performed with the AneuFinder pipeline, using a blacklist strategy to exclude reads from artefact-prone regions. For analyses, we divided genomes into non-overlapping bins of variable size based on mappability with a mean of 1 Mb (for more details, see Additional file 5: Methods). Quality control was performed using a clustering approach (minimum of 1 and maximum of 3 clusters) based on the quality parameters described in the main text.
T-ALL cell culture and live-cell time-lapse imaging
Mouse T-ALL cell lines were derived from primary mouse lymphomas and cultured as described by Jinadasa et al. . Cells were incubated at 37 °C, 5 % CO2 and low oxygen (3 % O2). Cells were retrovirally transduced with a constitutive H2B-GFP construct (pSF91.2 H2B-GFP) by spinfection for 90 min at 1000 × g and 32 °C in retroviral supernatant containing polybrene 4 μg/mL (Sigma). Transduced cells were harvested, washed and seeded onto Glass Bottom Dishes (Greiner) 1 h prior to imaging. Mitoses were captured at 2-min intervals over the course of 16 h on a DeltaVision microscope (GE healthcare) at 100× magnification using SoftWorx software (GE Healthcare).
In vivo leukaemia transplantation
Diagnostic bone marrow cells of patients with acute lymphoblastic leukaemia were thawed and transplanted by tail vein injection into sub-lethally irradiated immune deficient NOD/SCID/IL2Ry−/− (NSG) mice, at a dose of 2.5 × 106 cells/mouse. At 4-weekly intervals, blood was drawn from the retro-orbital plexus and leukaemia development was followed by flow cytometric analysis using the following markers against human antigens: CD19-PE, CD45-PECy7, CD34-APC, CD10-BV410 and CD20-BV605. Propidium Iodide was used to exclude dead cells. Human chimerism was defined as >1 % hCD45+ cells among total peripheral blood mononuclear cells (PBMC).
We are grateful to the other members of the Foijer and Lansdorp groups for their technical support and fruitful discussions. We thank Maria Sarkis Azkanaz for her assistance with analysing time-lapse microscopy data. Furthermore, we thank the operators at the Central Flow Cytometry Unit (UMCG), Geert Mesander, Henk Moes and Roelof Jan van der Lei, for their assistance with sorting, and the animal caretakers at the Central Animal Facility (UMCG). The authors are indebted to all patients who provided consent for the use of their material in this study.
This work was supported by a Dutch Cancer Society grant (KWF grant RUG-2012-5549) and a Groningen Foundation for Paediatric Oncology (SKOG) grant to FF and PML, and by an NWO (The Netherlands Organisation for Scientific Research) MEERVOUD grant and a University of Groningen Rosalind Franklin Fellowship grant to MCT.
Availability of data and materials
All array CGH and sequencing data have been deposited at public databases. aCGH data can be accessed from the Gene Expression Omnibus (NCBI GEO) database (http://www.ncbi.nlm.nih.gov/geo/) under accession no. GSE57334 . Sequencing data can be downloaded from ArrayExpress (https://www.ebi.ac.uk/arrayexpress/) under accession no. E-MTAB-4183. AneuFinder is available through BioConductor at: http://bioconductor.org/packages/AneuFinder/ and is licensed under Artistic-2.0
BB, AT, MCT and FF conceived the study. BB and MEB performed experiments and ESJMB provided human B-ALL samples. AT, BB, MCT and DP designed and trained the AneuFinder pipeline. BB, AT, MCT and FF analysed data. DCJS, NH, HGK and KHW performed single-cell library preparation and provided sequencing support under the supervision of DCJS and PML. AT, DP, TVJ, VG and MCT provided bioinformatics support. PML, MCT and FF provided financial support. BB, AT, MCT and FF prepared the manuscript with input from MEB, DP, DCJS, PML and TVJ. All authors read and approved the manuscript.
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Williams BR, Prabhu VR, Hunter KE, Glazier CM, Whittaker CA, Housman DE, et al. Aneuploidy affects proliferation and spontaneous immortalization in mammalian cells. Science. 2008;322:703–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Torres EM, Sokolsky T, Tucker CM, Chan LY, Boselli M, Dunham MJ, et al. Effects of aneuploidy on cellular physiology and cell division in haploid yeast. Science. 2007;317:916–24.View ArticlePubMedGoogle Scholar
- Baker DJ, Jeganathan KB, Cameron JD, Thompson M, Juneja S, Kopecka A, et al. BubR1 insufficiency causes early onset of aging-associated phenotypes and infertility in mice. Nat Genet. 2004;36:744–9.View ArticlePubMedGoogle Scholar
- Baker DJ, Jeganathan KB, Malureanu L, Perez-Terzic C, Terzic A, van Deursen JM. Early aging-associated phenotypes in Bub3/Rae1 haploinsufficient mice. J Cell Biol. 2006;172:529–40.View ArticlePubMedPubMed CentralGoogle Scholar
- Baker DJ, Dawlaty MM, Wijshake T, Jeganathan KB, Malureanu L, van Ree JH, et al. Increased expression of BubR1 protects against aneuploidy and cancer and extends healthy lifespan. Nat Cell Biol. 2013;15:96–102.View ArticlePubMedPubMed CentralGoogle Scholar
- Sheltzer JM, Amon A. The aneuploidy paradox: costs and benefits of an incorrect karyotype. Trends Genet. 2011;27:446–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Giam M, Rancati G. Aneuploidy and chromosomal instability in cancer: a jackpot to chaos. Cell Div. 2015;10:3.View ArticlePubMedPubMed CentralGoogle Scholar
- Siegel JJ, Amon A. New insights into the troubles of aneuploidy. Annu Rev Cell Dev Biol. 2012;28:189–214.View ArticlePubMedPubMed CentralGoogle Scholar
- Foijer F. Mouse models for chromosomal instability. In: Storchova Z, editor. Aneuploidy in health and disease. Rijeka: Intech; 2012. p. 59–78.Google Scholar
- Foijer F. CINister thoughts. Biochim Biophys Acta. 2010;38:1–7.Google Scholar
- Simon J, Bakker B, Foijer F. CINcere modelling: What have mouse models for chromosome instability taught us? Results Cancer Res. 2015;200:39–60.View ArticleGoogle Scholar
- Musacchio A, Salmon ED. The spindle-assembly checkpoint in space and time. Mol Cell. 2007;8:379–93.Google Scholar
- Kops GJPL, Shah JV. Connecting up and clearing out: how kinetochore attachment silences the spindle assembly checkpoint. Chromosoma. 2012;121:509–25.View ArticlePubMedGoogle Scholar
- Vleugel M, Hoogendoorn E, Snel B, Kops GJPL. Evolution and function of the mitotic checkpoint. Dev Cell. 2012;23:239–50.View ArticlePubMedGoogle Scholar
- Foijer F, Xie SZ, Simon JE, Bakker PL, Conte N, Davis SH, et al. Chromosome instability induced by Mps1 and p53 mutation generates aggressive lymphomas exhibiting aneuploidy-induced stress. Proc Natl Acad Sci. 2014;111:13427–32.View ArticlePubMedPubMed CentralGoogle Scholar
- Nicholson JM, Cimini D. Cancer karyotypes: survival of the fittest. Front Oncol. 2013;3:148.View ArticlePubMedPubMed CentralGoogle Scholar
- Gerlinger M, Rowan AJ, Horswell S, Larkin J, Endesfelder D, Gronroos E, et al. Intratumor heterogeneity and branched evolution revealed by multiregion sequencing. N Engl J Med. 2012;366:883–92.View ArticlePubMedPubMed CentralGoogle Scholar
- Heng HHQ, Liu G, Stevens JB, Abdallah BY, Horne SD, Ye KJ, et al. Karyotype heterogeneity and unclassified chromosomal abnormalities. Cytogenet Genome Res. 2013;139:144–57.View ArticlePubMedGoogle Scholar
- Bakker B, van den Bos H, Lansdorp PM, Foijer F. How to count chromosomes in a cell: An overview of current and novel technologies. BioEssays. 2015;37:570–7.View ArticlePubMedGoogle Scholar
- McGranahan N, Burrell RA, Endesfelder D, Novelli MR, Swanton C. Cancer chromosomal instability: therapeutic and diagnostic challenges. EMBO Rep. 2012;13:528–38.View ArticlePubMedPubMed CentralGoogle Scholar
- Macaulay IC, Voet T. Single cell genomics: advances and future perspectives. PLoS Genetics. 2014;10:e1004126.View ArticlePubMedPubMed CentralGoogle Scholar
- Navin N, Kendall J, Troge J, Andrews P, Rodgers L, McIndoo J, et al. Tumour evolution inferred by single-cell sequencing. Nature. 2011;472:90–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Knouse K, Wu J, Whittaker C, Amon A. Single cell sequencing reveals low levels of aneuploidy across mammalian tissues. Proc Natl Acad Sci. 2014;111:13409–14.View ArticlePubMedPubMed CentralGoogle Scholar
- van den Bos H, Spierings DCJ, Taudt AS, Bakker B, Porubský D, Falconer E, et al. Single-cell whole genome sequencing reveals no evidence for common aneuploidy in normal and Alzheimer's disease neurons. Genome Biol. 2016. doi:10.1186/s13059-016-0976-2.
- Baum L, Petrie T, Soules G, Weiss N. A Maximization technique occurring in the statistical analysis of probabilistic functions of Markov chains. Ann Math Stat. 1970;41:164–71.View ArticleGoogle Scholar
- Bhattacharyya A. On a measure of divergence between two statistical populations defined by their probability distributions. Bull Calcutta Math Soc. 1943;35:239–58.Google Scholar
- Garvin T, Aboukhalil R, Kendall J, Baslan T, Atwal GS, Hicks J, et al. Interactive analysis and assessment of single-cell copy-number variations. Nat Methods. 2015;12:1058–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Thompson SL, Compton DA. Examining the link between chromosomal instability and aneuploidy in human cells. J Cell Biol. 2008;180:665–72.View ArticlePubMedPubMed CentralGoogle Scholar
- Ganem NJ, Godinho SA, Pellman D. A mechanism linking extra centrosomes to chromosomal instability. Nature. 2009;460:278–82.View ArticlePubMedPubMed CentralGoogle Scholar
- Silkworth WT, Nardi IK, Scholl LM, Cimini D. Multipolar spindle pole coalescence is a major source of kinetochore mis-attachment and chromosome mis-segregation in cancer cells. PLoS One. 2009;4:e6564.View ArticlePubMedPubMed CentralGoogle Scholar
- Falconer E, Hills M, Naumann U, Poon SSS, Chavez EA, Sanders AD, et al. DNA template strand sequencing of single-cells maps genomic rearrangements at high resolution. Nat Methods. 2012;9:1107–12.View ArticlePubMedPubMed CentralGoogle Scholar
- Silk AD, Zasadil LM, Holland AJ, Vitre B, Cleveland DW, Weaver BA. Chromosome missegregation rate predicts whether aneuploidy will promote or suppress tumors. Proc Natl Acad Sci U S A. 2013;110:E4134–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Abdallah BY, Horne SD, Stevens JB, Liu G, Ying AY, Vanderhyden B, et al. Single cell heterogeneity: Why unstable genomes are incompatible with average profiles. Cell Cycle. 2013;12:3640–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Navin N, Hicks J. Future medical applications of single-cell sequencing in cancer. Genome Med. 2011;3:31.View ArticlePubMedPubMed CentralGoogle Scholar
- Alizadeh AA, Aranda V, Bardelli A, Blanpain C, Bock C, Borowski C, et al. Toward understanding and exploiting tumor heterogeneity. Nat Med. 2015;21:1–8.View ArticleGoogle Scholar
- Sotillo R, Hernando E, Díaz-Rodríguez E, Teruya-Feldstein J, Cordón-Cardo C, Lowe SW, et al. Mad2 overexpression promotes aneuploidy and tumorigenesis in mice. Cancer Cell. 2007;11:9–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Laughney AM, Elizalde S, Genovese G, Bakhoum SF. Dynamics of tumor heterogeneity derived from clonal karyotypic evolution. Cell Rep. 2015;12:1–12.View ArticleGoogle Scholar
- Saadatpour A, Lai S, Guo G, Yuan G-C. Single-cell analysis in cancer genomics. Trends Genet. 2015;31:576–86.View ArticlePubMedGoogle Scholar
- Jinadasa R, Balmus G, Gerwitz L, Roden J, Weiss R, Duhamel G. Derivation of thymic lymphoma T-cell lines from Atm(−/−) and p53(−/−) mice. J Vis Exp. 2011;50:6–8.Google Scholar