Detailed chromosomal and molecular genetic analysis of single cells by whole genome amplification and comparative genomic hybridisation.

Molecular genetic analysis of isolated single cells and other minute DNA samples is limited because there is insufficient DNA to perform more than one independent PCR amplification. One solution to this problem is to first amplify the entire genome, thus providing enough DNA for numerous subsequent PCRs. In this study we have investigated four different methods of whole genome amplification performed on single cells, and have identified a protocol that generates sufficient quantities of DNA for comparative genomic hybridisation (CGH) as well as more than 90 independent amplification reactions. Thus, numerous specific loci and the copy number of every chromosome can be assessed in a single cell. We report here the first reliable application of CGH to single cells from human preimplantation embryos (blastomeres) and to single fibroblasts, buccal cells and amniocytes.     

Multiple displacement amplification for malaria parasite DNA

Multiple displacement amplification (MDA) using Phi29 has proved to be an efficient, high-fidelity method for whole genome amplification in many organisms. This project was designed to evaluate this approach for use with the malaria parasite Plasmodium falciparum. In particular, we were concerned that the AT richness and presence of contaminating human DNA could limit efficiency of MDA in this system. We amplified 60 DNA samples using phi29 and scored 14 microsatellites, 9 single-nucleotide polymorphisms (SNPs), and gene copy number at GTP-cyclohydrolase I both before and after MDA. We observed 100% concordance in 829 microsatellite genotypes and in 499 SNP genotypes. Furthermore, copy number estimates for the GTP-cyclohydrolase I gene were correlated (r(2) = 0.67) in pre- and postamplification samples. These data confirm that MDA permits scoring of a range of different types of polymorphisms in P. falciparum malaria and can be used to extend the life of valuable DNA stocks.     

Allele-specific disparity in breast cancer

Background

Molecular alterations critical to development of cancer include mutations, copy number alterations (amplifications and deletions) as well as genomic rearrangements resulting in gene fusions. Massively parallel next generation sequencing, which enables the discovery of such changes, uses considerable quantities of genomic DNA (> 5 ug), a serious limitation in ever smaller clinical samples. However, a commonly available microarray platforms such as array comparative genomic hybridization (array CGH) allows the characterization of gene copy number at a single gene resolution using much smaller amounts of genomic DNA. In this study we evaluate the sensitivity of ultra-dense array CGH platforms developed by Agilent, especially that of the 1 million probe array (1 M array), and their application when whole genome amplification is required because of limited sample quantities.

Methods

We performed array CGH on whole genome amplified and not amplified genomic DNA from MCF-7 breast cancer cells, using 244 K and 1 M Agilent arrays. The ADM-2 algorithm was used to identify micro-copy number alterations that measured less than 1 Mb in genomic length.

Results

DNA from MCF-7 breast cancer cells was analyzed for micro-copy number alterations, defined as measuring less than 1 Mb in genomic length. The 4-fold extra resolution of the 1 M array platform relative to the less dense 244 K array platform, led to the improved detection of copy number variations (CNVs) and micro-CNAs. The identification of intra-genic breakpoints in areas of DNA copy number gain signaled the possible presence of gene fusion events. However, the ultra-dense platforms, especially the densest 1 M array, detect artifacts inherent to whole genome amplification and should be used only with non-amplified DNA samples.

Conclusions

This is a first report using 1 M array CGH for the discovery of cancer genes and biomarkers. We show the remarkable capacity of this technology to discover CNVs, micro-copy number alterations and even gene fusions. However, these platforms require excellent genomic DNA quality and do not tolerate relatively small imperfections related to the whole genome amplification.     

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.     

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.     

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.     

Identification of Nine Genomic Regions of Amplification in Urothelial Carcinoma,Correlation with Stage,and Potential Prognostic and Therapeutic Value

We performed a genome wide analysis of 164 urothelial carcinoma samples and 27 bladder cancer cell lines to identify copy number changes associated with disease characteristics, and examined the association of amplification events with stage and grade of disease. Multiplex inversion probe (MIP) analysis, a recently developed genomic technique, was used to study 80 urothelial carcinomas to identify mutations and copy number changes. Selected amplification events were then analyzed in a validation cohort of 84 bladder cancers by multiplex ligation-dependent probe assay (MLPA). In the MIP analysis, 44 regions of significant copy number change were identified using GISTIC. Nine gene-containing regions of amplification were selected for validation in the second cohort by MLPA. Amplification events at these 9 genomic regions were found to correlate strongly with stage, being seen in only 2 of 23 (9%) Ta grade 1 or 1–2 cancers, in contrast to 31 of 61 (51%) Ta grade 3 and T2 grade 2 cancers, p<0.001. These observations suggest that analysis of genomic amplification of these 9 regions might help distinguish non-invasive from invasive urothelial carcinoma, although further study is required. Both MIP and MLPA methods perform well on formalin-fixed paraffin-embedded DNA, enhancing their potential clinical use. Furthermore several of the amplified genes identified here (ERBB2, MDM2, CCND1) are potential therapeutic targets.     

Sequence quality is maintained after multiple displacement amplification of non-invasively obtained macaque semen DNA

Whole genome amplification protocols are revolutionizing the fields of molecular and conservation biology as they open the possibility of obtaining a large number of copies of a complete genome from minute amounts of sample. Multiple displacement amplification (MDA) is a whole genome amplification technique based on the properties of the phi29 DNA polymerase, which leads to a uniform representation of the genome with very low error rates. In this study we performed MDA on 28 macaque DNA samples extracted from blood or non-invasively collected semen from which we obtained mitochondrial control region sequences both before and after MDA. The length of the readable sequences was longer for the original samples than for the MDA products, but the number of unresolved positions was comparable both before and after MDA. We conclude that the MDA technique is useful for increasing the amount of DNA for sequencing mitochondrial regions in the case of non-invasively collected semen samples.     

Balanced-PCR amplification allows unbiased identification of genomic copy changes in minute cell and tissue samples

Analysis of genomic DNA derived from cells and fresh or fixed tissues often requires whole genome amplification prior to microarray screening. Technical hurdles to this process are the introduction of amplification bias and/or the inhibitory effects of formalin fixation on DNA amplification. Here we demonstrate a balanced-PCR procedure that allows unbiased amplification of genomic DNA from fresh or modestly degraded paraffin-embedded DNA samples. Following digestion and ligation of a target and a control genome with distinct linkers, the two are mixed and amplified in a single PCR, thereby avoiding biases associated with PCR saturation and impurities. We demonstrate genome-wide retention of allelic differences following balanced-PCR amplification of DNA from breast cancer and normal human cells and genomic profiling by array-CGH (cDNA arrays, 100 kb resolution) and by real-time PCR (single gene resolution). Comparison of balanced-PCR with multiple displacement amplification (MDA) demonstrates equivalent performance between the two when intact genomic DNA is used. When DNA from paraffin-embedded samples is used, balanced PCR overcomes problems associated with modest DNA degradation and produces unbiased amplification whereas MDA does not. Balanced-PCR allows amplification and recovery of modestly degraded genomic DNA for subsequent retrospective analysis of human tumors with known outcomes.     

Genome coverage and sequence fidelity of phi29 polymerase-based multiple strand displacement whole genome amplification

Major efforts are underway to systematically define the somatic and germline genetic variations causally associated with disease. Genome-wide genetic analysis of actual clinical samples is, however, limited by the paucity of genomic DNA available. Here we have tested the fidelity and genome representation of phi29 polymerase-based genome amplification (phi29MDA) using direct sequencing and high density oligonucleotide arrays probing >10,000 SNP alleles. Genome representation was comprehensive and estimated to be 99.82% complete, although six regions encompassing a maximum of 5.62 Mb failed to amplify. There was no degradation in the accuracy of SNP genotyping and, in direct sequencing experiments sampling 500,000 bp, the estimated error rate (9.5 x 10(-6)) was the same as in paired unamplified samples. The detection of cancer-associated loss of heterozygosity and copy number changes, including homozygous deletion and gene amplification, were similarly robust. These results suggest that phi29MDA yields high fidelity, near-complete genome representation suitable for high resolution genetic analysis.     

Multiple displacement amplification prior to single nucleotide polymorphism genotyping in epidemiologic studies

We assessed the whole genome amplification strategy, known as multiple displacement amplification (MDA), for use with the TaqMan genotyping platform for DNA samples derived from two case-control studies nested in the Nurses' Health Study and the Physicians' Health Study. Our objectives were to (1) quantify DNA yield from samples of varying starting concentrations and (2) assess whether MDA products give an accurate representation of the original genomic sequence. Multiple displacement amplification yielded a mean 23000-fold increase in DNA quantity and genotyping results demonstrate 99.95% accuracy across six SNPs from four genes for 352 samples included in this study. These results suggest that MDA will provide a sufficiently robust amplification of limiting samples of genomic DNA that can be used for SNP genotyping in large case-control studies of complex diseases.     

Whole genome amplification strategy for forensic genetic analysis using single or few cell equivalents of genomic DNA

Evidentiary items sometimes contain an insufficient quantity of DNA for routine forensic genetic analysis. These so-called low copy number DNA samples (< 100 pg of genomic DNA) often fall below the sensitivity limitations of routine DNA analysis methods. Theoretically, one way of making such intractable samples amenable to analysis would be to increase the number of starting genomes available for subsequent STR (short tandem repeat) analysis by a whole genome amplification strategy (WGA). Although numerous studies employing WGA have focused primarily on clinical applications, few in-depth studies have been conducted to evaluate the potential usefulness of these methods in forensic casework. After an initial evaluation of existing methods, a modified WGA strategy was developed that appears to have utility for low copy number forensic casework specimens. The method employs a slight, but important, modification of the "improved primer extension preamplification PCR" method (I-PEP-PCR), which we term mIPEP (modified-I-PEP-PCR). Complete autosomal STR and Y-STR (Y chromosome short tandem repeat) profiles were routinely obtained with 5 pg of template DNA, which is equivalent to 1-2 diploid cells. Remarkably, partial Y- and autosomal STR profiles were obtained from mIPEP-treated DNA recovered from bloodstains exposed to the outside environment for 1 year whereas non-mIPEP-treated samples did not produce profiles. STR profiles were obtained from contact DNA from single dermal ridge fingerprints when the DNA was subjected to prior mIPEP amplification but not when the mIPEP step was omitted.     

不同单细胞全基因组扩增方法的比较及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技术在人类单精子减数重组、非整倍体分析以及人类卵细胞基因组研究中的应用。     

Allele-specific amplification in cancer revealed by SNP array analysis

Amplification, deletion, and loss of heterozygosity of genomic DNA are hallmarks of cancer. In recent years a variety of studies have emerged measuring total chromosomal copy number at increasingly high resolution. Similarly, loss-of-heterozygosity events have been finely mapped using high-throughput genotyping technologies. We have developed a probe-level allele-specific quantitation procedure that extracts both copy number and allelotype information from single nucleotide polymorphism (SNP) array data to arrive at allele-specific copy number across the genome. Our approach applies an expectation-maximization algorithm to a model derived from a novel classification of SNP array probes. This method is the first to our knowledge that is able to (a) determine the generalized genotype of aberrant samples at each SNP site (e.g., CCCCT at an amplified site), and (b) infer the copy number of each parental chromosome across the genome. With this method, we are able to determine not just where amplifications and deletions occur, but also the haplotype of the region being amplified or deleted. The merit of our model and general approach is demonstrated by very precise genotyping of normal samples, and our allele-specific copy number inferences are validated using PCR experiments. Applying our method to a collection of lung cancer samples, we are able to conclude that amplification is essentially monoallelic, as would be expected under the mechanisms currently believed responsible for gene amplification. This suggests that a specific parental chromosome may be targeted for amplification, whether because of germ line or somatic variation. An R software package containing the methods described in this paper is freely available at http://genome.dfci.harvard.edu/~tlaframb/PLASQ.     

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.     

Whole genome sequencing of environmental Vibrio cholerae O1 from 10 nanograms of DNA using short reads

Multiple Displacement Amplification (MDA) of DNA using φ29 (phi29) DNA polymerase amplifies DNA several billion-fold, which has proved to be potentially very useful for evaluating genome information in a culture-independent manner. Whole genome sequencing using DNA from a single prokaryotic genome copy amplified by MDA has not yet been achieved due to the formation of chimeras and skewed amplification of genomic regions during the MDA step, which then precludes genome assembly. We have hereby addressed the issue by using 10 ng of genomic Vibrio cholerae DNA extracted within an agarose plug to ensure circularity as a starting point for MDA and then sequencing the amplified yield using the SOLiD platform. We successfully managed to assemble the entire genome of V. cholerae strain LMA3984-4 (environmental O1 strain isolated in urban Amazonia) using a hybrid de novo assembly strategy. Using our method, only 178 out of 16,713 (1%) of contigs were not able to be inserted into either chromosome scaffold, and out of these 178, only 3 appeared to be chimeras. The other contigs seem to be the result of template-independent non-specific amplification during MDA, yielding spurious reads. Extraction of genomic DNA within an agarose plug in order to ensure circularity of the extracted genome might be key to minimizing amplification bias by MDA for WGS.     

The Performance of Whole Genome Amplification Methods and Next-Generation Sequencing for Pre-Implantation Genetic Diagnosis of Chromosomal Abnormalities

Reliable and accurate pre-implantation genetic diagnosis(PGD) of patient’s embryos by next-generation sequencing(NGS) is dependent on efficient whole genome amplification(WGA) of a representative biopsy sample. However, the performance of the current state of the art WGA methods has not been evaluated for sequencing. Using low template DNA(15 pg) and single cells, we showed that the two PCR-based WGA systems Sure Plex and MALBAC are superior to the REPLI-g WGA multiple displacement amplification(MDA) system in terms of consistent and reproducible genome coverage and sequence bias across the 24 chromosomes, allowing better normalization of test to reference sequencing data. When copy number variation sequencing(CNV-Seq) was applied to single cell WGA products derived by either Sure Plex or MALBAC amplification, we showed that known disease CNVs in the range of 3e15 Mb could be reliably and accurately detected at the correct genomic positions. These findings indicate that our CNV-Seq pipeline incorporating either Sure Plex or MALBAC as the key initial WGA step is a powerful methodology for clinical PGD to identify euploid embryos in a patient’s cohort for uterine transplantation.     

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.     

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.     

Assessment of multiple displacement amplification in molecular epidemiology

Well-characterized epidemiological resources are generated with great effort, yet associated patient DNA samples can be limiting. The efficacy of the whole genome amplification (WGA) method, termed multiple displacement amplification (MDA), was assessed for detecting heterozygous sequence variants, mutation scanning, and PCR for challenging segments. Fifteen common polymorphisms from 10 genes located on 8 chromosomes were genotyped by direct sequencing of 300 PCR products from 115 high-quality MDA-amplified DNA samples extracted by different methods. The GC content of these analyzed segments ranges from 30% to 69%. Genotyping results demonstrate 100% accuracy. For heterozygotes, the relative intensity of peaks generated by the two alleles is highly similar for genomic and MDA-amplified genomic DNA, independent of GC content. In contrast, one of four heterozygous loci was mistyped when lower quality MDA-amplified DNA samples were used. The results of single-stranded conformation polymorphism (SSCP)-type of mutation scanningfor seven MDA-amplified DNA samples in four genes were concordant with the genomic DNA samples. PCR on MDA-amplified DNA was routinely successful for challenging 10- and 12-kb segments with GC content ranging from 30% to 80%, demonstrating that rather long segments, which are difficult to amplify with PCR, are amplified well with MDA. These results suggest that MDA is an effective method of WGA with utility in molecular epidemiology. Quality control of the MDA-amplified DNA is critical for high performance.