From microscopes to microarrays: dissecting recurrent chromosomal rearrangements

Submicroscopic chromosomal rearrangements that lead to copy-number changes have been shown to underlie distinctive and recognizable clinical phenotypes. The sensitivity to detect copy-number variation has escalated with the advent of array comparative genomic hybridization (CGH), including BAC and oligonucleotide-based platforms. Coupled with improved assemblies and annotation of genome sequence data, these technologies are facilitating the identification of new syndromes that are associated with submicroscopic genomic changes. Their characterization reveals the role of genome architecture in the aetiology of many clinical disorders. We review a group of genomic disorders that are mediated by segmental duplications, emphasizing the impact that high-throughput detection methods and the availability of the human genome sequence have had on their dissection and diagnosis.     

Segmental duplications and copy-number variation in the human genome

The human genome contains numerous blocks of highly homologous duplicated sequence. This higher-order architecture provides a substrate for recombination and recurrent chromosomal rearrangement associated with genomic disease. However, an assessment of the role of segmental duplications in normal variation has not yet been made. On the basis of the duplication architecture of the human genome, we defined a set of 130 potential rearrangement hotspots and constructed a targeted bacterial artificial chromosome (BAC) microarray (with 2,194 BACs) to assess copy-number variation in these regions by array comparative genomic hybridization. Using our segmental duplication BAC microarray, we screened a panel of 47 normal individuals, who represented populations from four continents, and we identified 119 regions of copy-number polymorphism (CNP), 73 of which were previously unreported. We observed an equal frequency of duplications and deletions, as well as a 4-fold enrichment of CNPs within hotspot regions, compared with control BACs (P < .000001), which suggests that segmental duplications are a major catalyst of large-scale variation in the human genome. Importantly, segmental duplications themselves were also significantly enriched >4-fold within regions of CNP. Almost without exception, CNPs were not confined to a single population, suggesting that these either are recurrent events, having occurred independently in multiple founders, or were present in early human populations. Our study demonstrates that segmental duplications define hotspots of chromosomal rearrangement, likely acting as mediators of normal variation as well as genomic disease, and it suggests that the consideration of genomic architecture can significantly improve the ascertainment of large-scale rearrangements. Our specialized segmental duplication BAC microarray and associated database of structural polymorphisms will provide an important resource for the future characterization of human genomic disorders.     

Development of a DNA-labeling system for array-based comparative genomic hybridization.

Chromosomal amplifications and deletions are critical components of tumorigenesis and DNA copy-number variations also correlate with changes in mRNA expression levels. Genome-wide microarray comparative genomic hybridization (CGH) has become an important method for detecting and mapping chromosomal changes in tumors. Thus, the ability to detect twofold differences in fluorescent intensity between samples on microarrays depends on the generation of high-quality labeled probes. To enhance array-based CGH analysis, a random prime genomic DNA labeling method optimized for improved sensitivity, signal-to-noise ratios, and reproducibility has been developed. The labeling system comprises formulated random primers, nucleotide mixtures, and notably a high concentration of the double mutant exo-large fragment of DNA polymerase I (exo-Klenow). Microarray analyses indicate that the genomic DNA-labeled templates yield hybridization signals with higher fluorescent intensities and greater signal-to-noise ratios and detect more positive features than the standard random prime and conventional nick translation methods. Also, templates generated by this system have detected twofold differences in gene copy number between male and female genomic DNA and identified amplification and deletions from the BT474 breast cancer cell line in microarray hybridizations. Moreover, alterations in gene copy number were routinely detected with 0.5 microg of genomic DNA starting sample. The method is flexible and performs efficiently with different fluorescently labeled nucleotides. Application of the optimized CGH labeling system may enhance the resolution and sensitivity of array-based CGH analysis in cancer and medical genetic studies.     

Array rank order regression analysis for the detection of gene copy-number changes in human cancer

cDNA microarray technology has been applied to the detection of DNA copy-number changes in malignant tumors. Test and control genomic DNA samples are differentially labeled and cohybridized to a spotted cDNA microarray. The ratio of test to control fluorescence intensities for each spot reflects relative gene copy number. The low signal-to-noise ratios of this assay and the variable levels of gene amplification and deletion among tumors hamper the detection of deviations from the diploid complement. We describe a regression-based statistical method to test for altered copy number on each gene and apply the technique to copy-number profiles in 10 thyroid tumors. We show that a novel transformation of fluorescence ratios into array rank order efficiently normalizes the heterogeneity among copy-number profiles and improves the reproducibility of the results. Array rank order regression analysis enhances the detection of consistent changes in gene copy number in solid tumors by cDNA microarray-based comparative genome hybridization.     

De Novo Identification of Single Nucleotide Mutations in Caenorhabditis elegans Using Array Comparative Genomic Hybridization

Array comparative genomic hybridization (aCGH) has been used primarily to detect copy-number variants between two genomes. Here we report using aCGH to detect single nucleotide mutations on oligonucleotide microarrays with overlapping 50-mer probes. This technique represents a powerful method for rapidly detecting novel homozygous single nucleotide mutations in any organism with a sequenced reference genome.     

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, 29 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.     

Single-cell chromosomal imbalances detection by array CGH

Genomic imbalances are a major cause of constitutional and acquired disorders. Therefore, aneuploidy screening has become the cornerstone of preimplantation, prenatal and postnatal genetic diagnosis, as well as a routine aspect of the diagnostic workup of many acquired disorders. Recently, array comparative genomic hybridization (array CGH) has been introduced as a rapid and high-resolution method for the detection of both benign and disease-causing genomic copy-number variations. Until now, array CGH has been performed using a significant quantity of DNA derived from a pool of cells. Here, we present an array CGH method that accurately detects chromosomal imbalances from a single lymphoblast, fibroblast and blastomere within a single day. Trisomy 13, 18, 21 and monosomy X, as well as normal ploidy levels of all other chromosomes, were accurately determined from single fibroblasts. Moreover, we showed that a segmental deletion as small as 34 Mb could be detected. Finally, we demonstrated the possibility to detect aneuploidies in single blastomeres derived from preimplantation embryos. This technique offers new possibilities for genetic analysis of single cells in general and opens the route towards aneuploidy screening and detection of unbalanced translocations in preimplantation embryos in particular.     

Array-MAPH: a methodology for the detection of locus copy-number changes in complex genomes

High-throughput genome-wide screening methods to detect subtle genomic imbalances are extremely important for diagnostic genetics and genomics. Here, we provide a detailed protocol for a microarray-based technique, applying the principle of multiplex amplifiable probe hybridization (MAPH). Methodology and software have been developed for designing unique PCR-amplifiable sequences (400-600 bp) covering any genomic region of interest. These sequences are amplified, cloned and spotted onto arrays (targets). A mixture of the same sequences (probes) is hybridized to genomic DNA immobilized on a membrane. Bound probes are recovered and quantitatively amplified by PCR, labeled and hybridized to the array. The procedure can be completed in 4-5 working days, excluding microarray preparation. Unlike array-comparative genomic hybridization (array-CGH), test DNA of specifically reduced complexity is hybridized to an array of identical small amplifiable target sequences, resulting in increased hybridization specificity and higher potential for increasing resolution. Array-MAPH can be used for detection of small-scale copy-number changes in complex genomes, leading to genotype-phenotype correlations and the discovery of new genes.     

Use of array CGH in the evaluation of dysmorphology, malformations, developmental delay, and idiopathic mental retardation

The clinical implementation of array comparative genomic hybridization has revolutionized the diagnosis of patients with syndromic or nonsyndromic mental retardation. Multiple studies of hundreds of patients with idiopathic mental retardation, and normal karyotype and/or subtelomeric testing using genome-wide microarray platforms with approximately 2000 to >30,000 (tiling-path) interrogating BAC/PAC probes have detected chromosome abnormalities in up to 17% of cases. Surprisingly, some of the pathogenic changes are mosaic and not detectable in conventional karyotyping. Commercially available genome-wide microarrays with >300,000 synthesized oligonucleotide probes enable higher resolution and sensitivity and will probably replace the BAC/PAC arrays in clinical laboratories.     

A stepwise framework for the normalization of array CGH data

Background  

In two-channel competitive genomic hybridization microarray experiments, the ratio of the two fluorescent signal intensities at each spot on the microarray is commonly used to infer the relative amounts of the test and reference sample DNA levels. This ratio may be influenced by systematic measurement effects from non-biological sources that can introduce biases in the estimated ratios. These biases should be removed before drawing conclusions about the relative levels of DNA. The performance of existing gene expression microarray normalization strategies has not been evaluated for removing systematic biases encountered in array-based comparative genomic hybridization (CGH), which aims to detect single copy gains and losses typically in samples with heterogeneous cell populations resulting in only slight shifts in signal ratios. The purpose of this work is to establish a framework for correcting the systematic sources of variation in high density CGH array images, while maintaining the true biological variations.     

Progress in the detection of human genome structural variations

The emerging of high-throughput and high-resolution genomic technologies led to the detection of submicroscopic variants ranging from 1 kb to 3 Mb in the human genome. These variants include copy number variations (CNVs), inversions, insertions, deletions and other complex rearrangements of DNA sequences. This paper briefly reviews the commonly used technologies to discover both genomic structural variants and their potential influences. Particularly, we highlight the array-based, PCR-based and sequencing-based assays, including array-based comparative genomic hybridization (aCGH), representational oligonucleotide microarray analysis (ROMA), multiplex amplifiable probe hybridization (MAPH), multiplex ligation-dependent probe amplification (MLPA), paired-end mapping (PEM), and next-generation DNA sequencing technologies. Furthermore, we discuss the limitations and challenges of current assays and give advices on how to make the database of genomic variations more reliable. Supported by the National High Technology Research and Development Program of China (Grant No. 2006AA020704).     

Human BAC library: construction and rapid screening

We have constructed a human genomic bacterial artificial chromosome (BAC) library using high molecular weight DNA from a pre-pro-B cell line, FLEB14-14, with a normal male diploid karyotype. This BAC library consists of 96 000 clones with an average DNA insert size of 110 kb, covering the human genome approximately 3 times. The library can be screened by three different methods. (1) Probe hybridization to 31 high-density replica (HDR) filters: each filter contains 3072 BAC clones which were gridded in a 6×6 pattern. (2) Probe hybridization to two Southern blot filters to which 31 HindIII digests of the pooled 3072 BAC clones were loaded. This identifies a particular HDR filter for which further probe hybridization is performed to identify a particular clone(s). (3) Two-step polymerase chain reaction (PCR). First, PCR is applied to DNA samples prepared from ten superpools of 9600 BAC clones each to identify a particular superpool and the second PCR is applied to 40 unique DNA samples prepared from the four-dimensionally assigned BAC clones of the particular superpool. We present typical examples of the library screening using these three methods. The two-step PCR screening is particularly powerful since it allows us to isolate a desired BAC clone(s) within a day or so. The theoretical consideration of the advantage of this method is presented. Furthermore, we have adapted Vectorette method to our BAC library for the isolation of terminal sequences of the BAC DNA insert to facilitate contig formation by BAC walking.     

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.     

Microarray detection of food-borne pathogens using specific probes prepared by comparative genomics

In order to design a method for the accurate detection and identification of food-borne pathogens, we used comparative genomics to select 70-mer oligonucleotide probes specific for 11 major food-borne pathogens (10 overlapping probes per pathogen) for use in microarray analysis. We analyzed the hybridization pattern of this constructed microarray with the Cy3-labeled genomic DNA of various food-borne pathogens and other bacteria. Our microarray showed a highly specific hybridization pattern with the genomic DNA of each food-borne pathogen; little unexpected cross-hybridization was observed. Microarray data were analyzed and clustered using the GenePix Pro 6.0 and GeneSpring GX 7.3.1 programs. The analyzed dendrogram revealed the discriminating power of constructed microarray. Each food-borne pathogen clustered according to its hybridization specificity and non-pathogenic species were discriminated from pathogenic species. Our method can be applied to the rapid and accurate detection and identification of food-borne pathogens in the food industry. In addition, this study demonstrates that genome sequence comparison and DNA microarray analysis have a powerful application in epidemiologic and taxonomic studies, as well as in the food safety and biodefense fields.     

Diagnostic genome profiling in mental retardation

Mental retardation (MR) occurs in 2%-3% of the general population. Conventional karyotyping has a resolution of 5-10 million bases and detects chromosomal alterations in approximately 5% of individuals with unexplained MR. The frequency of smaller submicroscopic chromosomal alterations in these patients is unknown. Novel molecular karyotyping methods, such as array-based comparative genomic hybridization (array CGH), can detect submicroscopic chromosome alterations at a resolution of 100 kb. In this study, 100 patients with unexplained MR were analyzed using array CGH for DNA copy-number changes by use of a novel tiling-resolution genomewide microarray containing 32,447 bacterial artificial clones. Alterations were validated by fluorescence in situ hybridization and/or multiplex ligation-dependent probe amplification, and parents were tested to determine de novo occurrence. Reproducible DNA copy-number changes were present in 97% of patients. The majority of these alterations were inherited from phenotypically normal parents, which reflects normal large-scale copy-number variation. In 10% of the patients, de novo alterations considered to be clinically relevant were found: seven deletions and three duplications. These alterations varied in size from 540 kb to 12 Mb and were scattered throughout the genome. Our results indicate that the diagnostic yield of this approach in the general population of patients with MR is at least twice as high as that of standard GTG-banded karyotyping.     

Genome-wide examination of chromosomal aberrations in neuroblastoma SH-SY5Y cells by array-based comparative genomic hybridization

Most neuroblastoma cells have chromosomal aberrations such as gains, losses, amplifications and deletions of DNA. Conventional approaches like fluorescence in situ hybridization (FISH) or metaphase comparative genomic hybridization (CGH) can detect chromosomal aberrations, but their resolution is low. In this study we used array-based comparative genomic hybridization to identify the chromosomal aberrations in human neuroblastoma SH-SY5Y cells. The DNA microarray consisting of 4000 bacterial artificial chromosome (BAC) clones was able to detect chromosomal regions with aberrations. The SH-SY5Y cells showed chromosomal gains in 1q12 approximately q44 (Chr1:142188905-246084832), 7 (over the whole chromosome), 2p25.3 approximately p16.3 (Chr2:18179-47899074), and 17q 21.32 approximately q25.3 (Chr17:42153031-78607159), while chromosomal losses detected were the distal deletion of 1p36.33 (Chr1:552910-563807), 14q21.1 approximately q21.3 (Chr14:37666271- 47282550), and 22q13.1 approximately q13.2 (Chr22:36885764-4190 7123). Except for the gain in 17q21 and the loss in 1p36, the other regions of gain or loss in SH-SY5Y cells were newly identified.     

Mapping tumor-suppressor genes with multipoint statistics from copy-number-variation data

Array-based comparative genomic hybridization (arrayCGH) is a microarray-based comparative genomic hybridization technique that has been used to compare tumor genomes with normal genomes, thus providing rapid genomic assays of tumor genomes in terms of copy-number variations of those chromosomal segments that have been gained or lost. When properly interpreted, these assays are likely to shed important light on genes and mechanisms involved in the initiation and progression of cancer. Specifically, chromosomal segments, deleted in one or both copies of the diploid genomes of a group of patients with cancer, point to locations of tumor-suppressor genes (TSGs) implicated in the cancer. In this study, we focused on automatic methods for reliable detection of such genes and their locations, and we devised an efficient statistical algorithm to map TSGs, using a novel multipoint statistical score function. The proposed algorithm estimates the location of TSGs by analyzing segmental deletions (hemi- or homozygous) in the genomes of patients with cancer and the spatial relation of the deleted segments to any specific genomic interval. The algorithm assigns, to an interval of consecutive probes, a multipoint score that parsimoniously captures the underlying biology. It also computes a P value for every putative TSG by using concepts from the theory of scan statistics. Furthermore, it can identify smaller sets of predictive probes that can be used as biomarkers for diagnosis and therapeutics. We validated our method using different simulated artificial data sets and one real data set, and we report encouraging results. We discuss how, with suitable modifications to the underlying statistical model, this algorithm can be applied generally to a wider class of problems (e.g., detection of oncogenes).     

Non-Denaturing Fluorescence in Situ Hybridization to Find Replication Origins in a Specific Genome Region on the DNA Fiber

Fluorescence in situ hybridization (FISH) is a useful method of determining the replication timing of specific genomic loci in mammals and of delineating replicon structures on DNA fibers in combination with in vivo replication labeling. In the case of simultaneous detection of a FISH probe and replicated forks, however, the DNA fibers are damaged by the DNA denaturation step for FISH detection, and the resulting fragmented fluorescence signals prevent analysis at high resolution. Here we found that hybridization of the probe to the genomic DNA was possible even under non-denaturing condition, but only at the time its genomic region replicated. Using the method designated non-denaturing FISH, we determined the replication timing of a specific BAC clone and the standard clones, and found that at least one replication origin exists within the genomic region covered by its BAC clone as an example.     

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.     

DNA Microarray for Rapid Detection of Mitochondrial DNA Haplotypes of Chum Salmon

For use in genetic stock identification, we developed an oligonucleotide (DNA) microarray hybridization method for rapid and accurate detection of nucleotide sequence variations in 20 previously identified variable nucleotide sites in about 500 bp within the 5 half of the control region of mitochondrial DNA of chum salmon (Oncorhynchus keta). The method includes immobilization of synthesized oligonucleotides containing respective polymorphic sites on a glass slide precoated with polycarbodiimide resin, a 2-hour hybridization with DNA microarray of biotinylated polymerase chain reaction fragments spanning the 5 variable portion followed by short washing, and visualization of hybridization signals by conventional ABC method and scanner-assisted computation of signal intensity on a computer. The entire process of hybridization and detection was completed within 4 hours. The resulting DNA microarray could detect all of the single nucleotide mutations and therefore could be used to identity the sequence variations defining 30 mtDNA haplotypes of chum salmon as revealed previously by nucleotide sequence analysis.