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

Ultradense array CGH and discovery of micro-copy number alterations and gene fusions in the cancer genome

The characterization of molecular alterations specific to cancer facilitates the discovery of predictive and prognostic biomarkers important to targeted therapeutics. Alterations critical to cancer therapeutics include copy number alterations (CNAs) such as gene amplifications and deletions as well as genomic rearrangements resulting in gene fusions. There are two genome-wide technologies used to detect CNAs: next generation sequencing (NGS) and dense microarray based comparative genomic hybridization, termed array CGH (aCGH). aCGH is a mature robust technology of lower cost and more accessible than NGS. This chapter describes the protocol steps and analysis required to obtain reliable aCGH results from clinical samples. Technical options and various necessary compromises related to the nature of clinical material are considered and the consequences of these choices for data analysis and interpretation are discussed. The chapter includes brief description of the data analysis, even though analysis is often performed by bioinformaticians. Today's cancer research requires collaboration of clinicians, molecular biologists, and mathematicians. Acquaintance with the basic principles related to the extraction of the data from arrays, its normalization and the algorithms available for analysis provides a baseline for mutual understanding and communication.

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.

DNA amplification techniques in pharmacogenomics

The variable predisposition of patients, both to disease susceptibility and drug response, is well established. It is largely attributed to genetic, as well as epigenetic variations between individuals, which may be inherited or acquired. The most common variation in the human genome is the SNP, which occurs throughout the genome, both within coding and noncoding regions. Characterization of SNPs in the context of both inherited and acquired conditions, such as cancer, are a main focus of many genotyping procedures. The demand for identifying (diagnosing) targeted SNPs or other variations, as well as the application of genome-wide screens, is continuously directing the development of new technologies. In general, most methods require a DNA amplification step to provide the amounts of DNA needed for the SNP detection step. In addition, DNA amplification is an important step when investigating other types of genomic information, for instance when addressing repeat, deletion, copy number variation or epigenetic regulation by DNA methylation. Besides the widely used PCR technique, there are several alternative approaches for genomic DNA amplification suitable for supporting the detection of genomic variation. In this article, we describe and evaluate a number of techniques, and discuss possible future prospects of DNA amplification in the fields of pharmacogenetics and pharmacogenomics.

Future medical applications of single-cell sequencing in cancer

Advances in whole genome amplification and next-generation sequencing methods have enabled genomic analyses of single cells, and these techniques are now beginning to be used to detect genomic lesions in individual cancer cells. Previous approaches have been unable to resolve genomic differences in complex mixtures of cells, such as heterogeneous tumors, despite the importance of characterizing such tumors for cancer treatment. Sequencing of single cells is likely to improve several aspects of medicine, including the early detection of rare tumor cells, monitoring of circulating tumor cells (CTCs), measuring intratumor heterogeneity, and guiding chemotherapy. In this review we discuss the challenges and technical aspects of single-cell sequencing, with a strong focus on genomic copy number, and discuss how this information can be used to diagnose and treat cancer patients.

Laser capture microdissection (LCM) and whole genome amplification (WGA) of DNA from normal breast tissue --- optimization for genome wide array analyses

Background

Laser capture microdissection (LCM) can be applied to tissues where cells of interest are distinguishable from surrounding cell populations. Here, we have optimized LCM for fresh frozen normal breast tissue where large amounts of fat can cause problems during microdissection. Since the amount of DNA needed for genome wide analyses, such as single nucleotide polymorphism (SNP) arrays, is often greater than what can be obtained from the dissected tissue, we have compared three different whole genome amplification (WGA) kits for amplification of DNA from LCM material. In addition, the genome wide profiling methods commonly used today require extremely high DNA quality compared to PCR based techniques and DNA quality is thus critical for successful downstream analyses.

Findings

We found that by using FrameSlides without glass backing for LCM and treating the slides with acetone after staining, the problems caused by excessive fat could be significantly decreased. The amount of DNA obtained after extraction from LCM tissue was not sufficient for direct SNP array analysis in our material. However, the two WGA kits based on Phi29 polymerase technology (Repli-g^® (Qiagen) and GenomiPhi (GE Healthcare)) gave relatively long amplification products, and amplified DNA from Repli-g^® gave call rates in the subsequent SNP analysis close to those from non-amplified DNA. Furthermore, the quality of the input DNA for WGA was found to be essential for successful SNP array results and initial DNA fragmentation problems could be reduced by switching from a regular halogen lamp to a VIS-LED lamp during LCM.

Conclusions

LCM must be optimized to work satisfactorily in difficult tissues. We describe a work flow for fresh frozen normal breast tissue where fat is inclined to cause problems if sample treatment is not adapted to this tissue. We also show that the Phi29-based Repli-g^® WGA kit (Qiagen) is a feasible approach to amplify DNA of high quality prior to genome wide analyses such as SNP profiling.

Tumour evolution inferred by single-cell sequencing

Genomic analysis provides insights into the role of copy number variation in disease, but most methods are not designed to resolve mixed populations of cells. In tumours, where genetic heterogeneity is common, very important information may be lost that would be useful for reconstructing evolutionary history. Here we show that with flow-sorted nuclei, whole genome amplification and next generation sequencing we can accurately quantify genomic copy number within an individual nucleus. We apply single-nucleus sequencing to investigate tumour population structure and evolution in two human breast cancer cases. Analysis of 100 single cells from a polygenomic tumour revealed three distinct clonal subpopulations that probably represent sequential clonal expansions. Additional analysis of 100 single cells from a monogenomic primary tumour and its liver metastasis indicated that a single clonal expansion formed the primary tumour and seeded the metastasis. In both primary tumours, we also identified an unexpectedly abundant subpopulation of genetically diverse 'pseudodiploid' cells that do not travel to the metastatic site. In contrast to gradual models of tumour progression, our data indicate that tumours grow by punctuated clonal expansions with few persistent intermediates.

The use of ultra-dense array CGH analysis for the discovery of micro-copy number alterations and gene fusions in the cancer genome

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.

Minimum altered regions in early prostate cancer progression identified by high resolution whole genome tiling path BAC array comparative hybridization

BACKGROUND

Carcinoma of the prostate (CaP) is a serious health problem. The altered molecular mechanisms that lead to this disease are poorly understood.

METHODS

Specimens from radical prostatectomies and blood were collected from 18 CaP surgery patients. For CGH studies, 20 CaP-related samples (16 Gleason grade 3, 3 higher grades, 1 BPH sample) and 18 samples of patient-matched normal epithelial cells were obtained by laser-assisted microdissection from frozen sections of the 18 prostatectomy specimens. High resolution SMRT aCGH was used to compare genomic profiles of prostatic samples to patient-matched blood and pooled female DNA. TMPRSS2-ERG fusion transcript analysis was performed by RT-PCR in relation to alterations detected at the TMPRSS2 locus.

RESULTS

Our comprehensive aCGH approach allowed us to define 35 regions of recurrent alterations while excluding germline copy number polymorphisms. Novel regions identified include 2q14.2, containing INHBB, and 17q21.31. The TMPRSS2 locus at 21q22.3 may be a hotspot for rearrangements with 75% of the alterations resulting in the expression of a TMPRSS2-ERG fusion transcript. Differences in fusion expression in different areas in an individual tumor focus and expression in adjacent normal epithelium supported intrafocal heterogeneity and field cancerization, respectively. Both features challenge our efforts to develop more objective markers for diagnosis and prediction of the severity of CaP.

CONCLUSION

The high-density array enabled precise mapping of genomic alterations and consequently definition of minimum altered regions smaller than previously reported thus facilitating identification of those genes that contribute to the cancer transformation process.

Genomic DNA amplification by the multiple displacement amplification (MDA) method

Large amounts of DNA are frequently required for use in detection assays and genomic analysis. The limited availability of DNA can be a critical obstacle to meeting research and clinical needs. DNA amplification methods are often required to generate sufficient material from small specimens or environmental samples with low DNA content. The MDA (multiple displacement amplification) reaction is increasingly the method of choice for many applications because of its extensive coverage of the genome, the generation of extremely long DNA products compared with older whole genome amplification methods and the high DNA yields, even from exceedingly low amounts of starting material. Remarkably, MDA enables genomic sequencing even from single microbial cells. Some of the uses of MDA and its strengths and limitations will be discussed.

Whole genome amplification and its impact on CGH array profiles

BACKGROUND

Some array comparative genomic hybridisation (array CGH) platforms require a minimum of micrograms of DNA for the generation of reliable and reproducible data. For studies where there are limited amounts of genetic material, whole genome amplification (WGA) is an attractive method for generating sufficient quantities of genomic material from miniscule amounts of starting material. A range of WGA methods are available and the multiple displacement amplification (MDA) approach has been shown to be highly accurate, although amplification bias has been reported. In the current study, WGA was used to amplify DNA extracted from whole blood. In total, six array CGH experiments were performed to investigate whether the use of whole genome amplified DNA (wgaDNA) produces reliable and reproducible results. Four experiments were conducted on amplified DNA compared to unamplified DNA and two experiments on unamplified DNA compared to unamplified DNA.

FINDINGS

All the experiments involving wgaDNA resulted in a high proportion of losses and gains of genomic material. Previously, amplification bias has been overcome by using amplified DNA in both the test and reference DNA. Our data suggests that this approach may not be effective, as the gains and losses introduced by WGA appears to be random and are not reproducible between different experiments using the same DNA.

CONCLUSION

In light of these findings, the use of both amplified test and reference DNA on CGH arrays may not provide an accurate representation of copy number variation in the DNA.

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.

Large fragment Bst DNA polymerase for whole genome amplification of DNA from formalin-fixed paraffin-embedded tissues

BACKGROUND

Formalin-fixed paraffin-embedded (FFPE) tissues represent the largest source of archival biological material available for genomic studies of human cancer. Therefore, it is desirable to develop methods that enable whole genome amplification (WGA) using DNA extracted from FFPE tissues. Multiple-strand Displacement Amplification (MDA) is an isothermal method for WGA that uses the large fragment of Bst DNA polymerase. To date, MDA has been feasible only for genomic DNA isolated from fresh or snap-frozen tissue, and yields a representational distortion of less than threefold.

RESULTS

We amplified genomic DNA of five FFPE samples of normal human lung tissue with the large fragment of Bst DNA polymerase. Using quantitative PCR, the copy number of 7 genes was evaluated in both amplified and original DNA samples. Four neuroblastoma xenograft samples derived from cell lines with known N-myc gene copy number were also evaluated, as were 7 samples of non-small cell lung cancer (NSCLC) tumors with known Skp2 gene amplification. In addition, we compared the array comparative genomic hybridization (CGH)-based genome profiles of two NSCLC samples before and after Bst MDA. A median 990-fold amplification of DNA was achieved. The DNA amplification products had a very high molecular weight (> 23 Kb). When the gene content of the amplified samples was compared to that of the original samples, the representational distortion was limited to threefold. Array CGH genome profiles of amplified and non-amplified FFPE DNA were similar.

CONCLUSION

Large fragment Bst DNA polymerase is suitable for WGA of DNA extracted from FFPE tissues, with an expected maximal representational distortion of threefold. Amplified DNA may be used for the detection of gene copy number changes by quantitative realtime PCR and genome profiling by array CGH.

Cytogenetically balanced translocations are associated with focal copy number alterations

Current cytogenetic methods (e.g., G-banding and multicolor chromosomal painting) allow detection of translocation events but lack the resolution to (a) locate the breakpoints precisely at the chromosome band level or (b) discriminate balanced translocations from translocations with copy number alterations not previously reported, or imperfectly balanced translocations. In this study, we demonstrate that cytogenetically balanced translocations are in fact frequently associated with segmental gain or loss of DNA. The recent development of a whole genome tiling path BAC array has enabled tiling resolution analysis of genomic segmental copy number status. Combining tiling resolution BAC array comparative genomic hybridization (array CGH) with G-Banding analysis and multicolor chromosomal painting approaches such as spectral karyotyping (SKY) facilitates high-resolution mapping of genomic alterations associated with imperfectly balanced translocations. Using a refined version of our CGH array we have deduced the copy number status throughout the genomes of three cytogenetically well-characterized prostate cancer cell lines (PC3, DU145, LNCaP) to determine whether translocations are associated with focal gains and losses of DNA. At 78 kb tiling resolution we identified the boundaries of 170, 80, and 34 known and novel copy number alterations (CNA) in these cell line genomes, respectively. Thirty-three of the 36 known translocations (92%, P < 0.001) in DU145 were associated with segmental CNA. Likewise, 80% (P < 0.001) of the known translocations showed association in LNCaP. Although many translocation breakpoints exhibit segmental alteration in PC3, the pattern of chromosomal rearrangements is too complex for use in comprehensive association with CNA boundaries. Our results reveal that imperfectly balanced translocations in tumor genomes are a phenomenon that occurs at frequencies much higher than previously demonstrated.

Metaphase and array comparative genomic hybridization: unique copy number changes and gene amplification of medulloblastomas in South America

Tumors of the central nervous system are the second most frequent malignancy of childhood, accounting for the majority of cancer-related deaths in this age group. Among these tumors, medulloblastomas (MB) remain in need of further genomic characterization toward understanding of pathogenesis and outcome predictors. Eight pediatric embryonal brain tumors were analyzed: five MB (one being desmoplastic), one PNET, one medulloepithelioma, and one ependymoblastoma. Analyses identified genomic imbalances, including the gain of 16p and the nonsyntenic coamplification of MYCN and TERT loci. More detailed FISH analysis showed that coamplification of MYCN and TERT in one of the MBs manifested as dispersed nuclear speckling, consistent with the presence of double minute chromosomes. There was considerable cell-to-cell copy number heterogeneity present, but it was clear that both genes were amplified concordantly. The amplification of oncogenes seems to play an important role in the pathogenesis of MB, and the association between MYCN and TERT amplifications and poor prognosis has not been well recognized. The uncharacteristic pattern of genomic imbalances detected in MB tumors may be a reflection of the characteristics of these tumors occurring in South America.

The use of whole genome amplification to study chromosomal changes in prostate cancer: insights into genome-wide signature of preneoplasia associated with cancer progression

Background

Prostate cancer (CaP) is a disease with multifactorial etiology that includes both genetic and environmental components. The knowledge of the genetic basis of CaP has increased over the past years, mainly in the pathways that underlie tumourigenesis, progression and drug resistance. The vast majority of cases of CaP are adenocarcinomas that likely develop through a pre-malignant lesion and high-grade prostatic intraepithelial neoplasia (HPIN). Histologically, CaP is a heterogeneous disease consisting of multiple, discrete foci of invasive carcinoma and HPIN that are commonly interspersed with benign glands and stroma. This admixture with benign tissue can complicate genomic analyses in CaP. Specifically, when DNA is bulk-extracted the genetic information obtained represents an average for all of the cells within the sample.

Results

To minimize this problem, we obtained DNA from individual foci of HPIN and CaP by laser capture microdissection (LCM). The small quantities of DNA thus obtained were then amplified by means of multiple-displacement amplification (MDA), for use in genomic DNA array comparative genomic hybridisation (gaCGH). Recurrent chromosome copy number abnormalities (CNAs) were observed in both HPIN and CaP. In HPIN, chromosomal imbalances involving chromosome 8 where common, whilst in CaP additional chromosomal changes involving chromosomes 6, 10, 13 and 16 where also frequently observed.

Conclusion

An overall increase in chromosomal changes was seen in CaP compared to HPIN, suggesting a universal breakdown in chromosomal stability. The accumulation of CNAs, which occurs during this process is non-random and may indicate chromosomal regions important in tumourigenesis. It is therefore likely that the alterations in copy number are part of a programmed cycle of events that promote tumour development, progression and survival. The combination of LCM, MDA and gaCGH is ideally suited for the identification of CNAs from small cell clusters and may assist in the discovery of potential genomic markers for early diagnosis, or identify the location of tumour suppressor genes (TSG) or oncogenes previously unreported in HPIN and CaP.

Whole genome amplification of Chelex-extracted DNA from a single mite: a method for studying genetics of the predatory mite Phytoseiulus persimilis

We developed and optimized a method using Chelex DNA extraction followed by whole genome amplification (WGA) to overcome problems conducting molecular genetic studies due to the limited amount of DNA obtainable from individual small organisms such as predatory mites. The DNA from a single mite, Phytoseiulus persimilis Athias-Henrot (Acari: Phytoseiidae), isolated in Chelex suspension was subjected to WGA. More than 1000-fold amplification of the DNA was achieved using as little as 0.03 ng genomic DNA template. The DNA obtained by the WGA was used for polymerase chain reaction followed by direct sequencing. From WGA DNA, nuclear DNA intergenic spacers ITS1 and ITS2 and a mitochondrial DNA 12S marker were tested in three different geographical populations of the predatory mite: California, the Netherlands, and Sicily. We found a total of four different alleles of the 12S in the Sicilian population, but no polymorphism was identified in the ITS marker. The combination of Chelex DNA extraction and WGA is thus shown to be a simple and robust technique for examining molecular markers for multiple loci by using individual mites. We conclude that the methods, Chelex extraction of DNA followed by WGA, provide a large quantity of DNA template that can be used for multiple PCR reactions useful for genetic studies requiring the genotypes of individual mites.

Analysis by comparative genomic hybridization of gastric cancer with peritoneal dissemination and/or positive peritoneal cytology

Peritoneal metastasis is an important prognostic factor in cases of gastric cancer. Although studies on comparative genomic hybridization (CGH) in gastric cancer have been reported, there are few reports on the peritoneal metastasis (P) and peritoneal cytology (CY) factors in this cancer. In this study, we analyzed the chromosomal changes in the primary tumor with a combination of laser microdissection analysis and CGH in an attempt to detect the unknown abnormal chromosomal regions. We analyzed 34 primary tumors, including 13 primary tumors with peritoneal metastasis (P1) and/or positive peritoneal cytology (CY1) using a combination of laser microdissection and CGH. The minimal overlapping regions in gains were assigned to 5p14 (46.2%), 7q21.3 (61.5%), 7q31 (46.2%), 7q36 (46.2%), 8q23 (53.8%), 15q26 (46.2%), 20q12 (61.5%), 20q13.1 (53.8%), and 20q13.2 (53.8%) in primary tumors with P1 and/or CY1. The minimal regions of losses that occurred most frequently were 4q34-q35 (23.1%) and 22q11.2 (23.1%). There were significant differences in the minimal regions of 5p14 (P=0.033), 7q21.3 (P < 0.0001), 7q31 (P=0.013), 7q36 (P=0.033), and 22q11.2 (P=0.048) between primary tumors with and without P1 and/or CY1. In this study, gain/amplification of 5p14, 7q21.3, 7q31, and 7q36, and loss of 22q11.2 were significant in gastric cancer cases with peritoneal dissemination and/or positive peritoneal cytology.

Genomic profiling by DNA amplification of laser capture microdissected tissues and array CGH

Comparative genomic hybridization by means of BAC microarrays (array CGH) allows high-resolution profiling of copy-number aberrations in tumor DNA. However, specific genetic lesions associated with small but clinically relevant tumor areas may pass undetected due to intra-tumor heterogeneity and/or the presence of contaminating normal cells. Here, we show that the combination of laser capture microdissection, phi29 DNA polymerase-mediated isothermal genomic DNA amplification, and array CGH allows genomic profiling of very limited numbers of cells. Moreover, by means of simple statistical models, we were able to bypass the exclusion of amplification distortions and variability prone areas, and to detect tumor-specific chromosomal gains and losses. We applied this new combined experimental and analytical approach to the genomic profiling of colorectal adenomatous polyps and demonstrated our ability to accurately detect single copy gains and losses affecting either whole chromosomes or small genomic regions from as little as 2 ng of DNA or 1000 microdissected cells.