Single-nucleotide polymorphism array coupled with multiple displacement amplification: accuracy and spatial resolution for analysis of chromosome copy numbers in few cells

When coupled with multiple displacement amplification (MDA), microarray-based comparative genomic intensity allows detection of chromosome copy number aberrations even in single or few cells, but the actual performance of the system and their influencing factors have not been well defined. Here, using single-nucleotide polymorphism (SNP) array, we analyzed copy number profiles from DNA amplified by MDA in 1-10 cells and estimated the accuracy and spatial resolution of the analysis. Based on the concordance of SNP copy numbers for DNA with and without MDA, the accuracy of the system can be significantly enhanced by using MDA-amplified DNA as reference and also by increasing the cell numbers. Analyses under different smoothing treatments revealed a practical resolution of 2?Mb for 10 cells and 10?Mb for a single cell. When both cells with known chromosomal duplication and deletion were analyzed, this platform detected a copy number "loss" more accurately than a "gain" (P < 0.01), particularly in single-cell MDA products. Together, we demonstrated that SNP array coupled with MDA is reliable and efficient for detection of copy number aberrations in a small number of cells, and its accuracy and resolution can both be significantly enhanced with increasing the number of cells as MDA template.     

Single-cell copy number variation detection

Detection of chromosomal aberrations from a single cell by array comparative genomic hybridization (single-cell array CGH), instead of from a population of cells, is an emerging technique. However, such detection is challenging because of the genome artifacts and the DNA amplification process inherent to the single cell approach. Current normalization algorithms result in inaccurate aberration detection for single-cell data. We propose a normalization method based on channel, genome composition and recurrent genome artifact corrections. We demonstrate that the proposed channel clone normalization significantly improves the copy number variation detection in both simulated and real single-cell array CGH data.     

Single-cell paired-end genome sequencing reveals structural variation per cell cycle

The nature and pace of genome mutation is largely unknown. Because standard methods sequence DNA from populations of cells, the genetic composition of individual cells is lost, de novo mutations in cells are concealed within the bulk signal and per cell cycle mutation rates and mechanisms remain elusive. Although single-cell genome analyses could resolve these problems, such analyses are error-prone because of whole-genome amplification (WGA) artefacts and are limited in the types of DNA mutation that can be discerned. We developed methods for paired-end sequence analysis of single-cell WGA products that enable (i) detecting multiple classes of DNA mutation, (ii) distinguishing DNA copy number changes from allelic WGA-amplification artefacts by the discovery of matching aberrantly mapping read pairs among the surfeit of paired-end WGA and mapping artefacts and (iii) delineating the break points and architecture of structural variants. By applying the methods, we capture DNA copy number changes acquired over one cell cycle in breast cancer cells and in blastomeres derived from a human zygote after in vitro fertilization. Furthermore, we were able to discover and fine-map a heritable inter-chromosomal rearrangement t(1;16)(p36;p12) by sequencing a single blastomere. The methods will expedite applications in basic genome research and provide a stepping stone to novel approaches for clinical genetic diagnosis.     

High frequency of common DNA copy number abnormalities detected by bacterial artificial chromosome array comparative genomic hybridization in 24 breast cancer cell lines

Breast cancer is a widespread disease in Japan and across the world. Breast cancer cells, as well as most other types of cancer cells, have diverse chromosomal aberrations. Clarifying the character of these chromosomal aberrations should contribute to the development of more suitable therapies, along with the predictions of metastasis and prognosis. Twenty-four breast cancer cell lines were analyzed by bacterial artificial chromosome (BAC) array comparative genomic hybridization (CGH). The array slide contained duplicate spots of 4030 BAC clone DNAs covering the entire human genome with 1 Mbp resolution. In all 24 breast cancer cell lines, frequent and significant amplifications as well as deletions were detected by BAC array CGH. Common DNA copy number gains, detected in 60% (above 15 cell lines) of the 24 breast cancer cell lines were found in 76 BAC clones, located at 1q, 5p, 8q, 9p, 16p, 17q, and 20q. Moreover, common DNA copy number loss was detected in 136 BAC clones, located at 1q, 2q, 3p, 4p, 6q, 8p, 9p, 11p, 13q, 17p, 18q, 19p, Xp, and Xq. The DNA copy number abnormalities found included abnormality of the well-known oncogene cMYC (8q24.21); however, most of them were not reported to relate to breast cancer. BAC array CGH has great potential to detect DNA copy number abnormalities, and has revealed that breast cancer cell lines have substantial heterogeneity.     

Use of BAC array CGH for evaluation of chromosomal stability of clinically used human mesenchymal stem cells and of cancer cell lines

Array-based comparative genomic hybridization (aCGH) using bacterial artificial chromosomes (BAC) is a powerful method to analyze DNA copy number aberrations of the entire human genome. In fact, CGH and aCGH have revealed various DNA copy number aberrations in numerous cancer cells and cancer cell lines examined so far. In this report, BAC aCGH was applied to evaluate the stability or instability of cell lines. Established cell lines have greatly contributed to advancements in not only biology but also medical science. However, cell lines have serious problems, such as alteration of biological properties during long-term cultivation. Firstly, we investigated two cancer cell lines, HeLa and Caco-2. HeLa cells, established from a cervical cancer, showed significantly increased DNA copy number alterations with passage time. Caco-2 cells, established from a colon cancer, showed no remarkable differences under various culture conditions. These results indicate that BAC aCGH can be used for the evaluation and validation of genomic stability of cultured cells. Secondly, BAC aCGH was applied to evaluate and validate the genomic stabilities of three patient's mesenchymal stem cells (MSCs), which were already used for their treatments. These three MSCs showed no significant differences in DNA copy number aberrations over their entire chromosomal regions. Therefore, BAC aCGH is highly recommended for use for a quality check of various cells before using them for any kind of biological investigation or clinical application.     

Whole-genome multiple displacement amplification from single cells

Multiple displacement amplification (MDA) is a recently described method of whole-genome amplification (WGA) that has proven efficient in the amplification of small amounts of DNA, including DNA from single cells. Compared with PCR-based WGA methods, MDA generates DNA with a higher molecular weight and shows better genome coverage. This protocol was developed for preimplantation genetic diagnosis, and details a method for performing single-cell MDA using the phi29 DNA polymerase. It can also be useful for the amplification of other minute quantities of DNA, such as from forensic material or microdissected tissue. The protocol includes the collection and lysis of single cells, and all materials and steps involved in the MDA reaction. The whole procedure takes 3 h and generates 1-2 microg of DNA from a single cell, which is suitable for multiple downstream applications, such as sequencing, short tandem repeat analysis or array comparative genomic hybridization.     

Genome-wide copy number analysis of single cells

Copy number variation (CNV) is increasingly recognized as an important contributor to phenotypic variation in health and disease. Most methods for determining CNV rely on admixtures of cells in which information regarding genetic heterogeneity is lost. Here we present a protocol that allows for the genome-wide copy number analysis of single nuclei isolated from mixed populations of cells. Single-nucleus sequencing (SNS), combines flow sorting of single nuclei on the basis of DNA content and whole-genome amplification (WGA); this is followed by next-generation sequencing to quantize genomic intervals in a genome-wide manner. Multiplexing of single cells is discussed. In addition, we outline informatic approaches that correct for biases inherent in the WGA procedure and allow for accurate determination of copy number profiles. All together, the protocol takes ～3 d from flow cytometry to sequence-ready DNA libraries.     

Genome-wide copy number profiling of single cells in S-phase reveals DNA-replication domains

Single-cell genomics is revolutionizing basic genome research and clinical genetic diagnosis. However, none of the current research or clinical methods for single-cell analysis distinguishes between the analysis of a cell in G1-, S- or G2/M-phase of the cell cycle. Here, we demonstrate by means of array comparative genomic hybridization that charting the DNA copy number landscape of a cell in S-phase requires conceptually different approaches to that of a cell in G1- or G2/M-phase. Remarkably, despite single-cell whole-genome amplification artifacts, the log2 intensity ratios of single S-phase cells oscillate according to early and late replication domains, which in turn leads to the detection of significantly more DNA imbalances when compared with a cell in G1- or G2/M-phase. Although these DNA imbalances may, on the one hand, be falsely interpreted as genuine structural aberrations in the S-phase cell’s copy number profile and hence lead to misdiagnosis, on the other hand, the ability to detect replication domains genome wide in one cell has important applications in DNA-replication research. Genome-wide cell-type-specific early and late replicating domains have been identified by analyses of DNA from populations of cells, but cell-to-cell differences in DNA replication may be important in genome stability, disease aetiology and various other cellular processes.     

Concurrent Whole-Genome Haplotyping and Copy-Number Profiling of Single Cells

Methods for haplotyping and DNA copy-number typing of single cells are paramount for studying genomic heterogeneity and enabling genetic diagnosis. Before analyzing the DNA of a single cell by microarray or next-generation sequencing, a whole-genome amplification (WGA) process is required, but it substantially distorts the frequency and composition of the cell’s alleles. As a consequence, haplotyping methods suffer from error-prone discrete SNP genotypes (AA, AB, BB) and DNA copy-number profiling remains difficult because true DNA copy-number aberrations have to be discriminated from WGA artifacts. Here, we developed a single-cell genome analysis method that reconstructs genome-wide haplotype architectures as well as the copy-number and segregational origin of those haplotypes by employing phased parental genotypes and deciphering WGA-distorted SNP B-allele fractions via a process we coin haplarithmisis. We demonstrate that the method can be applied as a generic method for preimplantation genetic diagnosis on single cells biopsied from human embryos, enabling diagnosis of disease alleles genome wide as well as numerical and structural chromosomal anomalies. Moreover, meiotic segregation errors can be distinguished from mitotic ones.     

Single-cell isolation from cell suspensions and whole genome amplification from single cells to provide templates for CGH analysis

A comprehensive genomic analysis of single cells is instrumental for numerous applications in tumor genetics, clinical diagnostics and forensic analyses. Here, we provide a protocol for single-cell isolation and whole genome amplification, which includes the following stages: preparation of single-cell suspensions from blood or bone marrow samples and cancer cell lines; their characterization on the basis of morphology, interphase fluorescent in situ hybridization pattern and antibody staining; isolation of single cells by either laser microdissection or micromanipulation; and unbiased amplification of single-cell genomes by either linker-adaptor PCR or GenomePlex library technology. This protocol provides a suitable template to screen for chromosomal copy number changes by conventional comparative genomic hybridization (CGH) or array CGH. Expected results include the generation of several micrograms of DNA from single cells, which can be used for CGH or other analyses, such as sequencing. Using linker-adaptor PCR or GenomePlex library technology, the protocol takes 72 or 30 h, respectively.     

Use of whole genome amplification and comparative genomic hybridisation to detect chromosomal copy number alterations in cell line material and tumour tissue

We have established that whole genome amplification (WGA), in conjunction with genomic DNA array comparative genomic hybridisation (gaCGH) allows for the identification of genome-wide copy number abnormalities (CNAs) in DNA extracted from both cell line and patient material. To determine the fidelity and reproducibility of WGA to detect copy number imbalances using gaCGH, well characterized cell line genomic DNA was analysed. The gaCGH data obtained from non-amplified DNA and amplified DNA for the neuroblastoma cell line NUB7 and a paediatric medulloblastoma patient was almost identical. In addition, laser capture microdissection (LCM) of prostate tumour cells and subsequent WGA allowed for the detection of a number of CNAs that may not have been identified if DNA had been extracted in bulk from heterogeneous tissue. The results presented here demonstrate the use of WGA for generating sufficient DNA for gaCGH analysis without the introduction of significant sequence representation bias. The combination of amplification and gaCGH using DNA extracted from archival patient material has the potential for permitting the studying of DNA from small cancerous or pre-cancerous foci, which may help to identify potential genomic markers for early diagnosis.     

High-resolution array comparative genomic hybridization of single micrometastatic tumor cells

Only few selected cancer cells drive tumor progression and are responsible for therapy resistance. Their specific genomic characteristics, however, are largely unknown because high-resolution genome analysis is currently limited to DNA pooled from many cells. Here, we describe a protocol for array comparative genomic hybridization (array CGH), which enables the detection of DNA copy number changes in single cells. Combining a PCR-based whole genome amplification method with arrays of highly purified BAC clones we could accurately determine known chromosomal changes such as trisomy 21 in single leukocytes as well as complex genomic imbalances of single cell line cells. In single T47D cells aberrant regions as small as 1–2 Mb were identified in most cases when compared to non-amplified DNA from 10⁶ cells. Most importantly, in single micrometastatic cancer cells isolated from bone marrow of breast cancer patients, we retrieved and confirmed amplifications as small as 4.4 and 5 Mb. Thus, high-resolution genome analysis of single metastatic precursor cells is now possible and may be used for the identification of novel therapy target genes.     

Human and mouse oligonucleotide-based array CGH

Array-based comparative genomic hybridization is a high resolution method for measuring chromosomal copy number changes. Here we present a validated protocol using in-house spotted oligonucleotide libraries for array comparative genomic hybridization (CGH). This oligo array CGH platform yields reproducible results and is capable of detecting single copy gains, multi-copy amplifications as well as homozygous and heterozygous deletions as small as 100 kb with high resolution. A human oligonucleotide library was printed on amine binding slides. Arrays were hybridized using a hybstation and analysed using BlueFuse feature extraction software, with >95% of spots passing quality control. The protocol allows as little as 300 ng of input DNA and a 90% reduction of Cot-1 DNA without compromising quality. High quality results have also been obtained with DNA from archival tissue. Finally, in addition to human oligo arrays, we have applied the protocol successfully to mouse oligo arrays. We believe that this oligo-based platform using ‘off-the-shelf’ oligo libraries provides an easy accessible alternative to BAC arrays for CGH, which is cost-effective, available at high resolution and easily implemented for any sequenced organism without compromising the quality of the results.     

不同单细胞全基因组扩增方法的比较及MALBAC在辅助生殖中的应用

单细胞全基因组扩增(whole genome amplification, WGA)是指在单细胞水平对全基因组进行扩增的新技术,其原理是将分离的单个细胞的微量全基因组DNA进行扩增,获得高覆盖率的完整的基因组后进行高通量测序,用于揭示细胞异质性。目前,WGA方法主要包括引物延伸预扩增(primer extension preamplification PCR, PEP-PCR)、简并寡核苷酸引物PCR (degenerate oligonucleotide primed PCR, DOP-PCR)、多重置换扩增(multiple displacement amplification, MDA)、多次退火环状循环扩增(multiple annealing and looping-based amplification cycles, MALBAC)等。本文对不同的单细胞WGA方法的原理及应用情况分别进行了阐述,并对其扩增效率进行评价和比较,包括基因组覆盖度、均一性、重现性、SNV (single-nucleotide variants)和CNV (copy number variants)检测力等。综合对比不同单细胞WGA方法后发现,MALBAC的扩增均一性最高、等位基因脱扣率最低、重现性最好,且对于CNV和SNV的检测效果最好。本文还阐述了MALBAC技术在人类单精子减数重组、非整倍体分析以及人类卵细胞基因组研究中的应用。     

Improved multiple displacement amplification with phi29 DNA polymerase for genotyping of single human cells

The ability to genotype multiple loci of single cells would be of significant benefit to investigations of cellular processes such as oncogenesis, meiosis, fertilization, and embryogenesis. We report a simple two-step, single-tube protocol for whole-genome amplification (WGA) from single human cells using components of the GenomiPhi V2 DNA Amplification kit. For the first time, we demonstrate reliable generation of 4-7 microg amplified DNA from a single human cell within 4 h with a minimum amount of artifactual DNA synthesis. DNA amplified from single cells was genotyped for 13 heterozygous short tandem repeats (STRs) and 7 heterozygous single nucleotide polymorphisms (SNPs), and the genotyping results were compared with purified genomic DNA. Accuracy of genotyping (percent of single-cell amplifications genotyped accurately for any particular STR or SNP) varied from 37% to 100% (with an average of 80%) for STRs and from 89% to 100% (averaging 94%) for SNPs. We suggest that the method described in this report is suitable for WGA from single cells, the product of which can be subsequently used for many applications, such as preimplantation genetic analysis (PGD).     

High fidelity of whole-genome amplified DNA on high-density single nucleotide polymorphism arrays

Current microarray technology allows researchers to genotype a large number of SNPs with relatively small amounts of DNA. Nevertheless, researchers and clinicians still frequently face the problem of acquiring enough high-quality DNA for analysis. Whole-genome amplification (WGA) methods offer a solution for this problem, and earlier studies have shown that WGA samples perform reasonably well in small-scale genetic analyses (e.g. Affymetrix 10K array). To determine the performance of WGA products on a large-scale genotyping array, we compared the Affymetrix 250K array genotyping results of genomic DNA and their WGA products from four individuals. Our results indicate that WGA product performs well on the 250K array compared to genomic DNA, especially when using the BRLMM calling algorithm. WGA samples have high call rates (97.5% on average, compared to 99.4% for genomic DNA) and excellent concordance rates with their corresponding genomic DNA samples (98.7% on average). In addition, no apparent systematic genomic amplification bias can be detected. This study demonstrates that, although there is a slight decrease in the total call rates, WGA methods provide a reliable approach for increasing the amount of DNA samples for use with a common SNP genotyping array.     

A Robust Method to Analyze Copy Number Alterations of Less than 100 kb in Single Cells Using Oligonucleotide Array CGH

Comprehensive genome wide analyses of single cells became increasingly important in cancer research, but remain to be a technically challenging task. Here, we provide a protocol for array comparative genomic hybridization (aCGH) of single cells. The protocol is based on an established adapter-linker PCR (WGAM) and allowed us to detect copy number alterations as small as 56 kb in single cells. In addition we report on factors influencing the success of single cell aCGH downstream of the amplification method, including the characteristics of the reference DNA, the labeling technique, the amount of input DNA, reamplification, the aCGH resolution, and data analysis. In comparison with two other commercially available non-linear single cell amplification methods, WGAM showed a very good performance in aCGH experiments. Finally, we demonstrate that cancer cells that were processed and identified by the CellSearch® System and that were subsequently isolated from the CellSearch® cartridge as single cells by fluorescence activated cell sorting (FACS) could be successfully analyzed using our WGAM-aCGH protocol. We believe that even in the era of next-generation sequencing, our single cell aCGH protocol will be a useful and (cost-) effective approach to study copy number alterations in single cells at resolution comparable to those reported currently for single cell digital karyotyping based on next generation sequencing data.     

Single-Cell Genetic Analysis Using Automated Microfluidics to Resolve Somatic Mosaicism

Somatic mosaicism occurs throughout normal development and contributes to numerous disease etiologies, including tumorigenesis and neurological disorders. Intratumor genetic heterogeneity is inherent to many cancers, creating challenges for effective treatments. Unfortunately, analysis of bulk DNA masks subclonal phylogenetic architectures created by the acquisition and distribution of somatic mutations amongst cells. As a result, single-cell genetic analysis is becoming recognized as vital for accurately characterizing cancers. Despite this, methods for single-cell genetics are lacking. Here we present an automated microfluidic workflow enabling efficient cell capture, lysis, and whole genome amplification (WGA). We find that ~90% of the genome is accessible in single cells with improved uniformity relative to current single-cell WGA methods. Allelic dropout (ADO) rates were limited to 13.75% and variant false discovery rates (SNV FDR) were 4.11x10^-6, on average. Application to ER-/PR-/HER2+ breast cancer cells and matched normal controls identified novel mutations that arose in a subpopulation of cells and effectively resolved the segregation of known cancer-related mutations with single-cell resolution. Finally, we demonstrate effective cell classification using mutation profiles with 10X average exome coverage depth per cell. Our data demonstrate an efficient automated microfluidic platform for single-cell WGA that enables the resolution of somatic mutation patterns in single cells.     

Detection of chromosomal structural alterations in single cells by SNP arrays: a systematic survey of amplification bias and optimized workflow

Background

In single-cell human genome analysis using whole-genome amplified product, a strong amplification bias involving allele dropout and preferential amplification hampers the quality of results. Using an oligonucleotide single nucleotide polymorphism (SNP) array, we systematically examined the nature of this amplification bias, including frequency, degree, and preference for genomic location, and we assessed the effects of this amplification bias on subsequent genotype and chromosomal copy number analyses.

Methodology/Principal Findings

We found a large variability in amplification bias among the amplified products obtained by multiple displacement amplification (MDA), and this bias had a severe effect on the genotype and chromosomal copy number analyses. We established optimal experimental conditions for pre-screening for high-quality amplified products, processing array data, and analyzing chromosomal structural alterations. Using this optimized protocol, we successfully detected previously unidentified chromosomal structural alterations in single cells from a lymphoblastoid cell line. These alterations were subsequently confirmed by karyotype analysis. In addition, we successfully obtained reproducible chromosomal copy number profiles of single cells from the cell line with a complex karyotype, indicating the applicability and potential of our optimized workflow.

Conclusions/Significance

Our results suggest that the quality of amplification products should be critically assessed before using them for genomic analyses. The method of MDA-based whole-genome amplification followed by SNP array analysis described here will be useful for exploring chromosomal alterations in single cells.     

Assessing the significance of conserved genomic aberrations using high resolution genomic microarrays

Genomic aberrations recurrent in a particular cancer type can be important prognostic markers for tumor progression. Typically in early tumorigenesis, cells incur a breakdown of the DNA replication machinery that results in an accumulation of genomic aberrations in the form of duplications, deletions, translocations, and other genomic alterations. Microarray methods allow for finer mapping of these aberrations than has previously been possible; however, data processing and analysis methods have not taken full advantage of this higher resolution. Attention has primarily been given to analysis on the single sample level, where multiple adjacent probes are necessarily used as replicates for the local region containing their target sequences. However, regions of concordant aberration can be short enough to be detected by only one, or very few, array elements. We describe a method called Multiple Sample Analysis for assessing the significance of concordant genomic aberrations across multiple experiments that does not require a-priori definition of aberration calls for each sample. If there are multiple samples, representing a class, then by exploiting the replication across samples our method can detect concordant aberrations at much higher resolution than can be derived from current single sample approaches. Additionally, this method provides a meaningful approach to addressing population-based questions such as determining important regions for a cancer subtype of interest or determining regions of copy number variation in a population. Multiple Sample Analysis also provides single sample aberration calls in the locations of significant concordance, producing high resolution calls per sample, in concordant regions. The approach is demonstrated on a dataset representing a challenging but important resource: breast tumors that have been formalin-fixed, paraffin-embedded, archived, and subsequently UV-laser capture microdissected and hybridized to two-channel BAC arrays using an amplification protocol. We demonstrate the accurate detection on simulated data, and on real datasets involving known regions of aberration within subtypes of breast cancer at a resolution consistent with that of the array. Similarly, we apply our method to previously published datasets, including a 250K SNP array, and verify known results as well as detect novel regions of concordant aberration. The algorithm has been fully implemented and tested and is freely available as a Java application at http://www.cbil.upenn.edu/MSA.