DNA copy number analysis of fresh and formalin-fixed specimens by shallow whole-genome sequencing with identification and exclusion of problematic regions in the genome assembly.

Detection of DNA copy number aberrations by shallow whole-genome sequencing (WGS) faces many challenges, including lack of completion and errors in the human reference genome, repetitive sequences, polymorphisms, variable sample quality, and biases in the sequencing procedures. Formalin-fixed paraffin-embedded (FFPE) archival material, the analysis of which is important for studies of cancer, presents particular analytical difficulties due to degradation of the DNA and frequent lack of matched reference samples. We present a robust, cost-effective WGS method for DNA copy number analysis that addresses these challenges more successfully than currently available procedures. In practice, very useful profiles can be obtained with ∼0.1× genome coverage. We improve on previous methods by first implementing a combined correction for sequence mappability and GC content, and second, by applying this procedure to sequence data from the 1000 Genomes Project in order to develop a blacklist of problematic genome regions. A small subset of these blacklisted regions was previously identified by ENCODE, but the vast majority are novel unappreciated problematic regions. Our procedures are implemented in a pipeline called QDNAseq. We have analyzed over 1000 samples, most of which were obtained from the fixed tissue archives of more than 25 institutions. We demonstrate that for most samples our sequencing and analysis procedures yield genome profiles with noise levels near the statistical limit imposed by read counting. The described procedures also provide better correction of artifacts introduced by low DNA quality than prior approaches and better copy number data than high-resolution microarrays at a substantially lower cost.

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.

Regional copy number-independent deregulation of transcription in cancer.

Genetic and epigenetic alterations have been identified that lead to transcriptional deregulation in cancers. Genetic mechanisms may affect single genes or regions containing several neighboring genes, as has been shown for DNA copy number changes. It was recently reported that epigenetic suppression of gene expression can also extend to a whole region; this is known as long-range epigenetic silencing. Various techniques are available for identifying regional genetic alterations, but no large-scale analysis has yet been carried out to obtain an overview of regional epigenetic alterations. We carried out an exhaustive search for regions susceptible to such mechanisms using a combination of transcriptome correlation map analysis and array CGH data for a series of bladder carcinomas. We validated one candidate region experimentally, demonstrating histone methylation leading to the loss of expression of neighboring genes without DNA methylation.

A collection of breast cancer cell lines for the study of functionally distinct cancer subtypes.

Recent studies suggest that thousands of genes may contribute to breast cancer pathophysiologies when deregulated by genomic or epigenomic events. Here, we describe a model "system" to appraise the functional contributions of these genes to breast cancer subsets. In general, the recurrent genomic and transcriptional characteristics of 51 breast cancer cell lines mirror those of 145 primary breast tumors, although some significant differences are documented. The cell lines that comprise the system also exhibit the substantial genomic, transcriptional, and biological heterogeneity found in primary tumors. We show, using Trastuzumab (Herceptin) monotherapy as an example, that the system can be used to identify molecular features that predict or indicate response to targeted therapies or other physiological perturbations.

Genomic and transcriptional aberrations linked to breast cancer pathophysiologies.

This study explores the roles of genome copy number abnormalities (CNAs) in breast cancer pathophysiology by identifying associations between recurrent CNAs, gene expression, and clinical outcome in a set of aggressively treated early-stage breast tumors. It shows that the recurrent CNAs differ between tumor subtypes defined by expression pattern and that stratification of patients according to outcome can be improved by measuring both expression and copy number, especially high-level amplification. Sixty-six genes deregulated by the high-level amplifications are potential therapeutic targets. Nine of these (FGFR1, IKBKB, ERBB2, PROCC, ADAM9, FNTA, ACACA, PNMT, and NR1D1) are considered druggable. Low-level CNAs appear to contribute to cancer progression by altering RNA and cellular metabolism.

Array comparative genomic hybridization and its applications in cancer.

Alteration in DNA copy number is one of the many ways in which gene expression and function may be modified. Some variations are found among normal individuals, others occur in the course of normal processes in some species and still others participate in causing various disease states. For example, many defects in human development are due to gains and losses of chromosomes and chromosomal segments that occur before or shortly after fertilization, and DNA dosage-alteration changes occurring in somatic cells are frequent contributors to cancer. Detecting these aberrations and interpreting them in the context of broader knowledge facilitates the identification of crucial genes and pathways involved in biological processes and disease. Over the past several years, array comparative genomic hybridization has proven its value for analyzing DNA copy-number variations. Here, we discuss the state of the art of array comparative genomic hybridization and its applications in cancer, emphasizing general concepts rather than specific results.

Epigenome analyses using BAC microarrays identify evolutionary conservation of tissue-specific methylation of SHANK3.

CpG islands are present in one-half of all human and mouse genes and typically overlap with promoters or exons. We developed a method for high-resolution analysis of the methylation status of CpG islands genome-wide, using arrays of BAC clones and the methylation-sensitive restriction enzyme NotI. Here we demonstrate the accuracy and specificity of the method. By computationally mapping all NotI sites, methylation events can be defined with single-nucleotide precision throughout the genome. We also demonstrate the unique expandability of the array method using a different methylation-sensitive restriction enzyme, BssHII. We identified and validated new CpG island loci that are methylated in a tissue-specific manner in normal human tissues. The methylation status of the CpG islands is associated with gene expression for several genes, including SHANK3, which encodes a structural protein in neuronal postsynaptic densities. Defects in SHANK3 seem to underlie human 22q13 deletion syndrome. Furthermore, these patterns for SHANK3 are conserved in mice and rats.

Fine-scale structural variation of the human genome.

Inversions, deletions and insertions are important mediators of disease and disease susceptibility. We systematically compared the human genome reference sequence with a second genome (represented by fosmid paired-end sequences) to detect intermediate-sized structural variants >8 kb in length. We identified 297 sites of structural variation: 139 insertions, 102 deletions and 56 inversion breakpoints. Using combined literature, sequence and experimental analyses, we validated 112 of the structural variants, including several that are of biomedical relevance. These data provide a fine-scale structural variation map of the human genome and the requisite sequence precision for subsequent genetic studies of human disease.

A tiling resolution DNA microarray with complete coverage of the human genome.

We constructed a tiling resolution array consisting of 32,433 overlapping BAC clones covering the entire human genome. This increases our ability to identify genetic alterations and their boundaries throughout the genome in a single comparative genomic hybridization (CGH) experiment. At this tiling resolution, we identified minute DNA alterations not previously reported. These alterations include microamplifications and deletions containing oncogenes, tumor-suppressor genes and new genes that may be associated with multiple tumor types. Our findings show the need to move beyond conventional marker-based genome comparison approaches, that rely on inference of continuity between interval markers. Our submegabase resolution tiling set for array CGH (SMRT array) allows comprehensive assessment of genomic integrity and thereby the identification of new genes associated with disease.

Integrated genomic and epigenomic analyses pinpoint biallelic gene inactivation in tumors.

Aberrant methylation of CpG islands and genomic deletion are two predominant mechanisms of gene inactivation in tumorigenesis, but the extent to which they interact is largely unknown. The lack of an integrated approach to study these mechanisms has limited the understanding of tumor genomes and cancer genes. Restriction landmark genomic scanning (RLGS; ref. 1) is useful for global analysis of aberrant methylation of CpG islands, but has not been amenable to alignment with deletion maps because the identity of most RLGS fragments is unknown. Here, we determined the nucleotide sequence and exact chromosomal position of RLGS fragments throughout the genome using the whole chromosome of origin of the fragments and in silico restriction digestion of the human genome sequence. To study the interaction of these gene-inactivation mechanisms in primary brain tumors, we integrated RLGS-based methylation analysis with high-resolution deletion maps from microarray-based comparative genomic hybridization (array CGH; ref. 3). Certain subsets of gene-associated CpG islands were preferentially affected by convergent methylation and deletion, including genes that exhibit tumor-suppressor activity, such as CISH1 (encoding SOCS1; ref. 4), as well as genes such as COE3 that have been missed by traditional non-integrated approaches. Our results show that most aberrant methylation events are focal and independent of deletions, and the rare convergence of these mechanisms can pinpoint biallelic gene inactivation without the use of positional cloning.

Deletion of chromosome 11q predicts response to anthracycline-based chemotherapy in early breast cancer.

Despite the recent consensus on the eligibility of adjuvant systemic therapy in patients with lymph node-negative breast cancer (NNBC) based on clinicopathologic criteria, specific biological markers are needed to predict sensitivity to the different available therapeutic options. We examined the feasibility of developing a genomic predictor of chemotherapy response and recurrence risk in 185 patients with NNBC using assembled arrays containing 2,460 bacterial artificial chromosome clones for scanning the genome for DNA copy number changes. After surgery, 90 patients received anthracycline-based chemotherapy, whereas 95 did not. Tamoxifen was administered to patients with hormone receptor-positive tumors. The association of genomic and clinicopathologic data and outcome was computed using Cox proportional hazard models and multiple testing adjustment procedures. Analysis of NNBC genomes revealed a common genomic signature. Specific DNA copy number aberrations were associated with hormonal receptor status, but not with other clinicopathologic variables. In patients treated with chemotherapy, none of the genomic changes were significantly correlated with recurrence. In patients not receiving chemotherapy, deletion of eight bacterial artificial chromosome clones clustered to chromosome 11q was independently associated with relapse (disease-free survival at 10 years+/-SE, 40%+/-14% versus 86%+/-6%; P<0.0001). The 54 patients with deletion of 11q (29%) did not present more aggressive clinicopathologic features than those without 11q loss. The adverse influence of 11q deletion on clinical outcome was confirmed in an independent validation series of 88 patients with NNBC. Our data suggests that patients with NNBC with the 11q deletion might benefit from anthracycline-based chemotherapy despite other clinical, pathologic, or genetic features. However, these initial findings should be evaluated in randomized clinical trials.

Improving melanoma classification by integrating genetic and morphologic features.

BACKGROUND: In melanoma, morphology-based classification systems have not been able to provide relevant information for selecting treatments for patients whose tumors have metastasized. The recent identification of causative genetic alterations has revealed mutations in signaling pathways that offer targets for therapy. Identifying morphologic surrogates that can identify patients whose tumors express such alterations (or functionally equivalent alterations) would be clinically useful for therapy stratification and for retrospective analysis of clinical trial data. METHODOLOGY/PRINCIPAL FINDINGS: We defined and assessed a panel of histomorphologic measures and correlated them with the mutation status of the oncogenes BRAF and NRAS in a cohort of 302 archival tissues of primary cutaneous melanomas from an academic comprehensive cancer center. Melanomas with BRAF mutations showed distinct morphological features such as increased upward migration and nest formation of intraepidermal melanocytes, thickening of the involved epidermis, and sharper demarcation to the surrounding skin; and they had larger, rounder, and more pigmented tumor cells (all p-values below 0.0001). By contrast, melanomas with NRAS mutations could not be distinguished based on these morphological features. Using simple combinations of features, BRAF mutation status could be predicted with up to 90.8% accuracy in the entire cohort as well as within the categories of the current World Health Organization (WHO) classification. Among the variables routinely recorded in cancer registries, we identified age < 55 y as the single most predictive factor of BRAF mutation in our cohort. Using age < 55 y as a surrogate for BRAF mutation in an independent cohort of 4,785 patients of the Southern German Tumor Registry, we found a significant survival benefit (p < 0.0001) for patients who, based on their age, were predicted to have BRAF mutant melanomas in 69% of the cases. This group also showed a different pattern of metastasis, more frequently involving regional lymph nodes, compared to the patients predicted to have no BRAF mutation and who more frequently displayed satellite, in-transit metastasis, and visceral metastasis (p < 0.0001). CONCLUSIONS: Refined morphological classification of primary melanomas can be used to improve existing melanoma classifications by forming subgroups that are genetically more homogeneous and likely to differ in important clinical variables such as outcome and pattern of metastasis. We expect this information to improve classification and facilitate stratification for therapy as well as retrospective analysis of existing trial data.

High-efficiency microarray printer using fused-silica capillary tube printing pins.

We describe a contact printing approach for microarrays that uses fused-silica capillary tubes with tapered tips for printing pins and a pressure/vacuum system to control pin loading, printing, and cleaning. The printing process is insensitive to variable environmental factors such as humidity, and the small diameter of the pins allows routine printing from 1536 well source plates. Pin load capacity, 0.2 microL in the current system, is adjustable by controlling pin length. More than 2000 spots can be printed per 0.2-microL pin load (<100 pl/spot), and densities of >12,000 spots/cm(2) are readily achievable. Solutions with a wide range of viscosities and chemical properties can be printed. The system can print tens of thousands of different solutions at high speed, due to the ability to use large numbers of pins simultaneously, and can produce a large number of replicate arrays since all of the solution picked up by the pins is available for deposition.

Analytical description of mutational effects in competing asexual populations.

The adaptation of a population to a new environment is a result of selection operating on a suite of stochastically occurring mutations. This article presents an analytical approach to understanding the population dynamics during adaptation, specifically addressing a system in which periods of growth are separated by selection in bottlenecks. The analysis derives simple expressions for the average properties of the evolving population, including a quantitative description of progressive narrowing of the range of selection coefficients of the predominant mutant cells and of the proportion of mutant cells as a function of time. A complete statistical description of the bottlenecks is also presented, leading to a description of the stochastic behavior of the population in terms of effective mutation times. The effective mutation times are related to the actual mutation times by calculable probability distributions, similar to the selection coefficients being highly restricted in their probable values. This analytical approach is used to model recently published experimental data from a bacterial coculture experiment, and the results are compared to those of a numerical model published in conjunction with the data. Finally, experimental designs that may improve measurements of fitness distributions are suggested.

Roadmap for new opportunities in melanoma research.

Investigators representing all major melanoma research areas present an overview of the most important challenges for the field. Four major research areas are covered plus the training of new investigators. For each area we first describe the present status, its strengths and weaknesses, and then outline specific recommendations. In basic research of melanoma, we outline the pertinent issues for melanoma classification, understanding melanocyte development and transformation, melanoma resistance, tumor microenvironment, metastasis, animal models, immune response, and blood and tissue diagnostics. In clinical research we provide an overview of the current challenges and the strategies for characterization, monitoring, and therapy. It will be important to develop strong research and clinical infrastructures by establishing tumor banks, identifying and validating biomarkers, developing new imaging techniques, and increasing multidisciplinary collaboration and communication. To strengthen the field we need to recruit both young and established investigators and foster career development plans that cover all disciplines. Recent research advances provide significant opportunities to have a major impact on this devastating disease. This group provides recommendations for both short- and long-term strategies that build on research strengths and opportunities established by the many members of the research community.

Centromeric chromosomal translocations show tissue-specific differences between squamous cell carcinomas and adenocarcinomas.

Structural chromosomal aberrations are common in epithelial tumors. Here, we compared the location of centromeric breaks associated with whole arm translocations in seven adenocarcinoma cell lines and nine squamous cell carcinoma cell lines using SKY, microarray-based comparative genomic hybridization (array CGH) and fluorescence in situ hybridization (FISH). Whole arm translocations were more frequent in squamous cell carcinomas (112 in nine cell lines and nine in one short-term culture) than in adenocarcinomas (13 in seven cases) and most often resulted in copy number alterations. Array CGH analysis demonstrated that in all squamous cell carcinomas and in most adenocarcinomas, the breakpoints of unbalanced whole arm translocations occurred between the two clones on the array flanking the centromeres. However, FISH with centromeric probes revealed that in squamous cell carcinomas, the marker chromosomes with whole arm translocations contained centromeres comprised of material from both participating chromosomes, while in adenocarcinomas centromeric material from only one of the chromosomes was present. These observations suggest that different mechanisms of centromeric instability underlie the formation of chromosomal aberrations in adenocarcinomas and squamous cell carcinomas.

Rare amplicons implicate frequent deregulation of cell fate specification pathways in oral squamous cell carcinoma.

Genomes of solid tumors are characterized by gains and losses of regions, which may contribute to tumorigenesis by altering gene expression. Often the aberrations are extensive, encompassing whole chromosome arms, which makes identification of candidate genes in these regions difficult. Here, we focused on narrow regions of gene amplification to facilitate identification of genetic pathways important in oral squamous cell carcinoma (SCC) development. We used array comparative genomic hybridization (array CGH) to define minimum common amplified regions and then used expression analysis to identify candidate driver genes in amplicons that spanned <3 Mb. We found genes involved in integrin signaling (TLN1), survival (YAP1, BIRC2), and adhesion and migration (TLN1, LAMA3, MMP7), as well as members of the hedgehog (GLI2) and notch (JAG1, RBPSUH, FJX1) pathways to be amplified and overexpressed. Deregulation of these and other members of the hedgehog and notch pathways (HHIP, SMO, DLL1, NOTCH4) implicates deregulation of developmental and differentiation pathways, cell fate misspecification, in oral SCC development.

Breast tumor copy number aberration phenotypes and genomic instability.

BACKGROUND: Genomic DNA copy number aberrations are frequent in solid tumors, although the underlying causes of chromosomal instability in tumors remain obscure. Genes likely to have genomic instability phenotypes when mutated (e.g. those involved in mitosis, replication, repair, and telomeres) are rarely mutated in chromosomally unstable sporadic tumors, even though such mutations are associated with some heritable cancer prone syndromes. METHODS: We applied array comparative genomic hybridization (CGH) to the analysis of breast tumors. The variation in the levels of genomic instability amongst tumors prompted us to investigate whether alterations in processes/genes involved in maintenance and/or manipulation of the genome were associated with particular types of genomic instability. RESULTS: We discriminated three breast tumor subtypes based on genomic DNA copy number alterations. The subtypes varied with respect to level of genomic instability. We find that shorter telomeres and altered telomere related gene expression are associated with amplification, implicating telomere attrition as a promoter of this type of aberration in breast cancer. On the other hand, the numbers of chromosomal alterations, particularly low level changes, are associated with altered expression of genes in other functional classes (mitosis, cell cycle, DNA replication and repair). Further, although loss of function instability phenotypes have been demonstrated for many of the genes in model systems, we observed enhanced expression of most genes in tumors, indicating that over expression, rather than deficiency underlies instability. CONCLUSION: Many of the genes associated with higher frequency of copy number aberrations are direct targets of E2F, supporting the hypothesis that deregulation of the Rb pathway is a major contributor to chromosomal instability in breast tumors. These observations are consistent with failure to find mutations in sporadic tumors in genes that have roles in maintenance or manipulation of the genome.

Array comparative genomic hybridization identifies genetic subgroups in grade 4 human astrocytoma.

Alterations of DNA copy number are believed to be important indicators of tumor progression in human astrocytoma. We used an array of bacterial artificial chromosomes to map relative DNA copy number in 50 primary glioblastoma multiforme tumors at approximately 1.4-Mb resolution. We identified 33 candidate sites for amplification and homozygous deletion in these tumors. We identified three major genetic subgroups within these glioblastoma multiforme tumors: tumors with chromosome 7 gain and chromosome 10 loss, tumors with only chromosome 10 loss in the absence of chromosome 7 gain, and tumors without copy number change in chromosomes 7 or 10. The significance of these genetic groups to therapeutics needs further study.

Duplication of distal 20q: clinical, cytogenetic and array CGH. Characterization of a new case.

Trisomy of the long arm of chromosome 20 is rare. We describe an 18-month-old male who was born at 36 weeks via Caesarian section after an uneventful pregnancy. During the newborn period he was found to have a right-sided cleft lip and cleft palate, hypertelorism, strabismus and mildly over-folded ears with cupping. Cardiovascular examination was consistent with the diagnosis of severe aortic coarctation, which was confirmed by echocardiogram. Additionally, hypothyroidism was diagnosed. Neurological evaluation at 18 months revealed a hypotonic infant with delayed acquisition of motor milestones. Cytogenetic analysis showed additional material on the long arm of chromosome 20, confirmed by fluorescence in situ hybridization (FISH) analysis as being of chromosome 20 origin. Because of the indistinct GTG-banding pattern it was not possible to distinguish between a proximal [dup(20)(q11.2q13.1)] or distal duplication [dup(20)(q13.1q13.3)]. To further define the duplication we used array comparative genomic hybridization (CGH) which demonstrated a 7.8 Mb interstitial duplication in distal 20q. Thus, the proband's karyotype was interpreted as 46,XY,dup(20)(q13.2q13.2). The proband is the first reported case of a pure duplication of this region. This case further highlights the utility of array CGH in characterizing aneusomies and, in particular, for accurate breakpoint designation and quantitation of ambiguous rearrangements.