Acquired copy number alterations in adult acute myeloid leukemia genomes.

Cytogenetic analysis of acute myeloid leukemia (AML) cells has accelerated the identification of genes important for AML pathogenesis. To complement cytogenetic studies and to identify genes altered in AML genomes, we performed genome-wide copy number analysis with paired normal and tumor DNA obtained from 86 adult patients with de novo AML using 1.85 million feature SNP arrays. Acquired copy number alterations (CNAs) were confirmed using an ultra-dense array comparative genomic hybridization platform. A total of 201 somatic CNAs were found in the 86 AML genomes (mean, 2.34 CNAs per genome), with French-American-British system M6 and M7 genomes containing the most changes (10-29 CNAs per genome). Twenty-four percent of AML patients with normal cytogenetics had CNA, whereas 40% of patients with an abnormal karyotype had additional CNA detected by SNP array, and several CNA regions were recurrent. The mRNA expression levels of 57 genes were significantly altered in 27 of 50 recurrent CNA regions <5 megabases in size. A total of 8 uniparental disomy (UPD) segments were identified in the 86 genomes; 6 of 8 UPD calls occurred in samples with a normal karyotype. Collectively, 34 of 86 AML genomes (40%) contained alterations not found with cytogenetics, and 98% of these regions contained genes. Of 86 genomes, 43 (50%) had no CNA or UPD at this level of resolution. In this study of 86 adult AML genomes, the use of an unbiased high-resolution genomic screen identified many genes not previously implicated in AML that may be relevant for pathogenesis, along with many known oncogenes and tumor suppressor genes.

Dynamic histone H3 epigenome marking during the intraerythrocytic cycle of Plasmodium falciparum.

Epigenome profiling has led to the paradigm that promoters of active genes are decorated with H3K4me3 and H3K9ac marks. To explore the epigenome of Plasmodium falciparum asexual stages, we performed MS analysis of histone modifications and found a general preponderance of H3/H4 acetylation and H3K4me3. ChIP-on-chip profiling of H3, H3K4me3, H3K9me3, and H3K9ac from asynchronous parasites revealed an extensively euchromatic epigenome with heterochromatin restricted to variant surface antigen gene families (VSA) and a number of genes hitherto unlinked to VSA. Remarkably, the vast majority of the genome shows an unexpected pattern of enrichment of H3K4me3 and H3K9ac. Analysis of synchronized parasites revealed significant developmental stage specificity of the epigenome. In rings, H3K4me3 and H3K9ac are homogenous across the genes marking active and inactive genes equally, whereas in schizonts, they are enriched at the 5' end of active genes. This study reveals an unforeseen and unique plasticity in the use of the epigenetic marks and implies the presence of distinct epigenetic pathways in gene silencing/activation throughout the erythrocytic cycle.

High-resolution genome-wide cytosine methylation profiling with simultaneous copy number analysis and optimization for limited cell numbers.

Many genome-wide assays involve the generation of a subset (or representation) of the genome following restriction enzyme digestion. The use of enzymes sensitive to cytosine methylation allows high-throughput analysis of this epigenetic regulatory process. We show that the use of a dual-adapter approach allows us to generate genomic representations that includes fragments of <200 bp in size, previously not possible when using the standard approach of using a single adapter. By expanding the representation to smaller fragments using HpaII or MspI, we increase the representation by these isoschizomers to more than 1.32 million loci in the human genome, representing 98.5% of CpG islands and 91.1% of refSeq promoters. This advance allows the development of a new, high-resolution version of our HpaII-tiny fragment Enrichment by Ligation-mediated PCR (HELP) assay to study cytosine methylation. We also show that the MspI representation generates information about copy-number variation, that the assay can be used on as little as 10 ng of DNA and that massively parallel sequencing can be used as an alternative to microarrays to read the output of the assay, making this a powerful discovery platform for studies of genomic and epigenomic abnormalities.

An integrative genomic and epigenomic approach for the study of transcriptional regulation.

The molecular heterogeneity of acute leukemias and other tumors constitutes a major obstacle towards understanding disease pathogenesis and developing new targeted-therapies. Aberrant gene regulation is a hallmark of cancer and plays a central role in determining tumor phenotype. We predicted that integration of different genome-wide epigenetic regulatory marks along with gene expression levels would provide greater power in capturing biological differences between leukemia subtypes. Gene expression, cytosine methylation and histone H3 lysine 9 (H3K9) acetylation were measured using high-density oligonucleotide microarrays in primary human acute myeloid leukemia (AML) and acute lymphocytic leukemia (ALL) specimens. We found that DNA methylation and H3K9 acetylation distinguished these leukemias of distinct cell lineage, as expected, but that an integrative analysis combining the information from each platform revealed hundreds of additional differentially expressed genes that were missed by gene expression arrays alone. This integrated analysis also enhanced the detection and statistical significance of biological pathways dysregulated in AML and ALL. Integrative epigenomic studies are thus feasible using clinical samples and provide superior detection of aberrant transcriptional programming than single-platform microarray studies.

wuHMM: a robust algorithm to detect DNA copy number variation using long oligonucleotide microarray data.

Copy number variants (CNVs) are currently defined as genomic sequences that are polymorphic in copy number and range in length from 1000 to several million base pairs. Among current array-based CNV detection platforms, long-oligonucleotide arrays promise the highest resolution. However, the performance of currently available analytical tools suffers when applied to these data because of the lower signal:noise ratio inherent in oligonucleotide-based hybridization assays. We have developed wuHMM, an algorithm for mapping CNVs from array comparative genomic hybridization (aCGH) platforms comprised of 385 000 to more than 3 million probes. wuHMM is unique in that it can utilize sequence divergence information to reduce the false positive rate (FPR). We apply wuHMM to 385K-aCGH, 2.1M-aCGH and 3.1M-aCGH experiments comparing the 129X1/SvJ and C57BL/6J inbred mouse genomes. We assess wuHMM's performance on the 385K platform by comparison to the higher resolution platforms and we independently validate 10 CNVs. The method requires no training data and is robust with respect to changes in algorithm parameters. At a FPR of <10%, the algorithm can detect CNVs with five probes on the 385K platform and three on the 2.1M and 3.1M platforms, resulting in effective resolutions of 24 kb, 2-5 kb and 1 kb, respectively.

Methods in high-resolution, array-based comparative genomic hybridization.

A method of high resolution, array-based comparative genomic hybridization is described for the mapping of copy-number changes associated with chromosomal amplifications, deletions, and translocations. The method involves the design of whole-genome or targeted, fine-tiling arrays for synthesis on a high-density digital microarray-synthesis platform. The arrays can span entire eukaryotic genomes or be targeted to specific chromosomal regions for high-resolution identification of copy-number changes and the corresponding breakpoint locations. The methods described include the bioinformatics required for array design, and the protocols for DNA fragmentation, dual-color labeling, microarray hybridization, and array scanning. The processes for data extraction, normalization, and segmentation analysis are also described.

Discovery of a previously unrecognized microdeletion syndrome of 16p11.2-p12.2.

We have identified a recurrent de novo pericentromeric deletion in 16p11.2-p12.2 in four individuals with developmental disabilities by microarray-based comparative genomic hybridization analysis. The identification of common clinical features in these four individuals along with the characterization of complex segmental duplications flanking the deletion regions suggests that nonallelic homologous recombination mediated these rearrangements and that deletions in 16p11.2-p12.2 constitute a previously undescribed syndrome.

Characterization of a recurrent 15q24 microdeletion syndrome.

We describe multiple individuals with mental retardation and overlapping de novo submicroscopic deletions of 15q24 (1.7-3.9 Mb in size). High-resolution analysis showed that in three patients both proximal and distal breakpoints co-localized to highly identical segmental duplications (>51 kb in length, > 94% identity), suggesting non-allelic homologous recombination as the likely mechanism of origin. Sequencing studies in a fourth individual provided base pair resolution and showed that both breakpoints in this case were located in unique sequence. Despite the differences in the size and location of the deletions, all four individuals share several major features (growth retardation, microcephaly, digital abnormalities, hypospadias and loose connective tissue) and resemble one another facially (high anterior hair line, broad medial eyebrows, hypertelorism, downslanted palpebral fissures, broad nasal base, long smooth philtrum and full lower lip), indicating that this represents a novel syndrome caused by haploinsufficiency of one or more dosage-sensitive genes in the minimal deletion region. Our results define microdeletion of 15q24 as a novel recurrent genomic disorder.

Efficient high-resolution deletion discovery in Caenorhabditis elegans by array comparative genomic hybridization.

We have developed array Comparative Genomic Hybridization for Caenorhabditis elegans as a means of screening for novel induced deletions in this organism. We designed three microarrays consisting of overlapping 50-mer probes to annotated exons and micro-RNAs, the first with probes to chromosomes X and II, the second with probes to chromosome II alone, and a third to the entire genome. These arrays were used to reliably detect both a large (50 kb) multigene deletion and a small (1 kb) single-gene deletion in homozygous and heterozygous samples. In one case, a deletion breakpoint was resolved to fewer than 50 bp. In an experiment designed to identify new mutations we used the X:II and II arrays to detect deletions associated with lethal mutants on chromosome II. One is an 8-kb deletion targeting the ast-1 gene on chromosome II and another is a 141-bp deletion in the gene C06A8.1. Others span large sections of the chromosome, up to >750 kb. As a further application of array Comparative Genomic Hybridization in C. elegans we used the whole-genome array to detect the extensive natural gene content variation (almost 2%) between the N2 Bristol strain and the strain CB4856, a strain isolated in Hawaii and JU258, a strain isolated in Madeira.

A high-resolution map of segmental DNA copy number variation in the mouse genome.

Submicroscopic (less than 2 Mb) segmental DNA copy number changes are a recently recognized source of genetic variability between individuals. The biological consequences of copy number variants (CNVs) are largely undefined. In some cases, CNVs that cause gene dosage effects have been implicated in phenotypic variation. CNVs have been detected in diverse species, including mice and humans. Published studies in mice have been limited by resolution and strain selection. We chose to study 21 well-characterized inbred mouse strains that are the focus of an international effort to measure, catalog, and disseminate phenotype data. We performed comparative genomic hybridization using long oligomer arrays to characterize CNVs in these strains. This technique increased the resolution of CNV detection by more than an order of magnitude over previous methodologies. The CNVs range in size from 21 to 2,002 kb. Clustering strains by CNV profile recapitulates aspects of the known ancestry of these strains. Most of the CNVs (77.5%) contain annotated genes, and many (47.5%) colocalize with previously mapped segmental duplications in the mouse genome. We demonstrate that this technique can identify copy number differences associated with known polymorphic traits. The phenotype of previously uncharacterized strains can be predicted based on their copy number at these loci. Annotation of CNVs in the mouse genome combined with sequence-based analysis provides an important resource that will help define the genetic basis of complex traits.

Chromosome Conformation Capture Carbon Copy (5C): a massively parallel solution for mapping interactions between genomic elements.

Physical interactions between genetic elements located throughout the genome play important roles in gene regulation and can be identified with the Chromosome Conformation Capture (3C) methodology. 3C converts physical chromatin interactions into specific ligation products, which are quantified individually by PCR. Here we present a high-throughput 3C approach, 3C-Carbon Copy (5C), that employs microarrays or quantitative DNA sequencing using 454-technology as detection methods. We applied 5C to analyze a 400-kb region containing the human beta-globin locus and a 100-kb conserved gene desert region. We validated 5C by detection of several previously identified looping interactions in the beta-globin locus. We also identified a new looping interaction in K562 cells between the beta-globin Locus Control Region and the gamma-beta-globin intergenic region. Interestingly, this region has been implicated in the control of developmental globin gene switching. 5C should be widely applicable for large-scale mapping of cis- and trans- interaction networks of genomic elements and for the study of higher-order chromosome structure.

Discovery of previously unidentified genomic disorders from the duplication architecture of the human genome.

Genomic disorders are characterized by the presence of flanking segmental duplications that predispose these regions to recurrent rearrangement. Based on the duplication architecture of the genome, we investigated 130 regions that we hypothesized as candidates for previously undescribed genomic disorders. We tested 290 individuals with mental retardation by BAC array comparative genomic hybridization and identified 16 pathogenic rearrangements, including de novo microdeletions of 17q21.31 found in four individuals. Using oligonucleotide arrays, we refined the breakpoints of this microdeletion, defining a 478-kb critical region containing six genes that were deleted in all four individuals. We mapped the breakpoints of this deletion and of four other pathogenic rearrangements in 1q21.1, 15q13, 15q24 and 17q12 to flanking segmental duplications, suggesting that these are also sites of recurrent rearrangement. In common with the 17q21.31 deletion, these breakpoint regions are sites of copy number polymorphism in controls, indicating that these may be inherently unstable genomic regions.

Comparative isoschizomer profiling of cytosine methylation: the HELP assay.

The distribution of cytosine methylation in 6.2 Mb of the mouse genome was tested using cohybridization of genomic representations from a methylation-sensitive restriction enzyme and its methylation-insensitive isoschizomer. This assay, termed HELP (HpaII tiny fragment Enrichment by Ligation-mediated PCR), allows both intragenomic profiling and intergenomic comparisons of cytosine methylation. The intragenomic profile shows most of the genome to be contiguous methylated sequence with occasional clusters of hypomethylated loci, usually but not exclusively at promoters and CpG islands. Intergenomic comparison found marked differences in cytosine methylation between spermatogenic and brain cells, identifying 223 new candidate tissue-specific differentially methylated regions (T-DMRs). Bisulfite pyrosequencing confirmed the four candidates tested to be T-DMRs, while quantitative RT-PCR for two genes with T-DMRs located at their promoters showed the HELP data to be correlated with gene activity at these loci. The HELP assay is robust, quantitative, and accurate and is providing new insights into the distribution and dynamic nature of cytosine methylation in the genome.

Analysis of chromosome breakpoints in neuroblastoma at sub-kilobase resolution using fine-tiling oligonucleotide array CGH.

Understanding the genes and genetic pathways targeted by recurrent chromosomal imbalances in malignancy, along with the molecular mechanisms that generate the imbalances, are important problems in cancer biology. In this report, we demonstrate that oligonucleotide array CGH (oaCGH) analysis can routinely map chromosomal imbalance breakpoints at exon-level resolution, including imbalances that are single copy number genomic alterations. Different tiling-path array designs were used in this study: a whole-genome array with a 6-kb median probe spacing and fine-tiling arrays for selected genomic regions with either 50- or 140-bp median probe spacing. In both array formats, oligonucleotide probes were of isothermal design and were tiled through genic and inter-genic regions. Whole-genome oaCGH analysis of two neuroblastoma cell lines and three primary tumors led to the identification of 58 chromosomal breakpoints that generated 45 large-scale partial chromosomal imbalances (> 2 Mb). An unexpectedly high proportion (34%) of these breakpoint intervals mapped to regions containing segmental duplications. In addition, 88 smaller-sized regions (< 2 Mb) of imbalance were detected, the majority of which mapped to segmentally duplicated regions and may reflect constitutional copy number polymorphisms. The chromosomal breakpoints for 12 recurrent abnormalities exhibited in neuroblastoma tumors and cell lines, including MYCN amplicon boundaries, loss of 3p, loss of 11q, and gain of 17q, could be mapped to intervals ranging from 50 bp to 10 kb in size using high-density fine-tiling oligonucleotide microarrays. Fine-tiling oaCGH analysis provides an unprecedented level of resolution, allowing detailed mapping of recurrent unbalanced chromosomal abnormalities. Supplementary material for this article can be found on the Genes, Chromosomes, and Cancer website at http://www.interscience.wiley.com/jpages/1045-2257/suppmat/index.html.

Differential requirement for the ATPase domain of the Cockayne syndrome group B gene in the processing of UV-induced DNA damage and 8-oxoguanine lesions in human cells.

Cockayne syndrome (CS) is a rare inherited human genetic disorder characterized by UV sensitivity, developmental abnormalities and premature aging. The cellular and molecular phenotypes of CS include increased sensitivity to oxidative and UV-induced DNA lesions. The CSB protein is thought to play a pivotal role in transcription-coupled repair and CS-B cells are defective in the repair of the transcribed strand of active genes, both after exposure to UV and in the presence of oxidative DNA lesions. A previous study has indicated that a conserved helicase ATPase motif II residue is essential for the function of the CSB protein in responding to UV-induced DNA damage in a hamster cell line. Due to the limitations in studying a complex human disorder in another species, this study introduced the site-directed mutation of the ATPase motif II in the human CSB gene in an isogenic human cell line. The CSB mutant allele was tested for genetic complementation of UV-sensitive phenotypes in the human CS-B cell line CS1AN.S3.G2. In addition, the incision of an 8-oxoguanine lesion by extracts of the CS-B cell lines stably transfected with the wild-type or ATPase mutant CSB gene has been investigated. The ATPase motif II point mutation (E646Q) abolished the function of the CSB protein to complement the UV-sensitive phenotypes of survival, RNA synthesis recovery and apoptosis. Interestingly, whole-cell extract prepared from these mutant cells retained wild-type incision activity on an oligonucleotide containing a single 8-oxoguanine lesion, whereas the absence of the CSB gene altogether resulted in reduced incision activity relative to wild-type. These results suggest damage-specific functional requirements for CSB in the repair of UV-induced and oxidative lesions in human cells. The transfection of the mutant or wild-type CSB gene into the CS1AN.S3.G2 cells did not alter the expression of the subset of genes examined by cDNA array analysis.