A new look towards BAC-based array CGH through a comprehensive comparison with oligo-based array CGH

Genomic amplification of chromosome 20q13.33 is the early biomarker for the development of sporadic colorectal carcinoma

BACKGROUND

Colorectal carcinoma (CRC) is the third most common cancer in the world and also the third leading cause of cancer-related mortality in Taiwan. CRC tumorigenesis is a multistep process, starting from mutations causing loss of function of tumor suppressor genes, canonically demonstrated in adenomatous polyposis coli pathogenesis. Although many genes or chromosomal alterations have been shown to be involved in this process, there are still unrecognized molecular events within CRC tumorigenesis. Elucidating these mechanisms may help improve the management and treatment.

METHODS

In this study, we aimed to identify copy number alteration of the smallest chromosomal regions that is significantly associated with sporadic CRC tumorigenesis using high-resolution array-based Comparative Genomic Hybridization (aCGH) and quantitative Polymerase chain reaction (qPCR). In addition, microsatellite instability assay and sequencing-based mutation assay were performed to illustrate the initiation event of CRC tumorigenesis.

RESULTS

A total of 571 CRC patients were recruited and 377 paired CRC tissues from sporadic CRC cases were used to define the smallest regions with chromosome copy number changes. In addition, 198 colorectal polyps from 160 patients were also used to study the role of 20q13.33 gain in CRC tumorigenesis. We found that gain in 20q13.33 is the main chromosomal abnormalities in this patient population and counts 50.9 and 62.8% in CRC and colon polyps, respectively. Furthermore, APC and KRAS gene mutations were profiled simultaneously and co-analyzed with microsatellite instability and 20q13.33 gain in CRC patients. Our study showed that the frequency of 20q13.33 copy number gain was highest among all reported CRC mutations.

CONCLUSION

As APC or KRAS mutations are currently identified as the most important targets for CRC therapy, this study proposes that 20q13.33 copy number gain and the associated chromosomal genes function as promising biomarkers for both early stage detection and targeted therapy of sporadic CRCs in the future.

Challenges in identifying large germline structural variants for clinical use by long read sequencing

Genomic structural variations, previously considered rare events, are widely recognized as a major source of inter-individual variability and hence, a major hurdle in optimum patient stratification and disease management. Herein, we focus on large complex germline structural variations and present challenges towards target treatment via the synergy of state-of-the-art approaches and information technology tools. A complex structural variation detection remains challenging, as there is no gold standard for identifying such genomic variations with long reads, especially when the chromosomal rearrangement in question is a few Mb in length. A clinical case with a large complex chromosomal rearrangement serves as a paradigm. We feel that functional validation and data interpretation are of outmost importance for information growth to be translated into knowledge growth and hence, new working practices are highlighted.

Diagnostics based on nucleic acid sequence variant profiling: PCR, hybridization, and NGS approaches

Nucleic acid sequence variations have been implicated in many diseases, and reliable detection and quantitation of DNA/RNA biomarkers can inform effective therapeutic action, enabling precision medicine. Nucleic acid analysis technologies being translated into the clinic can broadly be classified into hybridization, PCR, and sequencing, as well as their combinations. Here we review the molecular mechanisms of popular commercial assays, and their progress in translation into in vitro diagnostics.

Comparison of the Glasgow Prognostic Score (GPS) and the modified Glasgow Prognostic Score (mGPS) in evaluating the prognosis of patients with operable and inoperable non-small cell lung cancer

PURPOSE

The Glasgow Prognostic Score (GPS) and modified Glasgow Prognostic Score (mGPS) are shown to be reliable prognostic indexes in patients with operable and inoperable non-small cell lung cancer (NSCLC). Considering the difference between the two indexes lies in whether hypoalbuminemia without an elevated C-reactive protein (CRP) is associated with worse survival, this study aims to evaluate the prognostic performance of hypoalbuminemia in patients without an elevated CRP and to compare the prognostic value of GPS and mGPS in patients with operable and inoperable NSCLC.

METHODS

The data of 2988 patients were retrospectively collected from the Shanghai Health Information Network. Univariate and multivariate Cox regression was performed to investigate the prognostic effect of albumin, CRP, GPS and mGPS. Restricted cubic spline was also performed to evaluate the relationship between albumin and hazard ratio. Kaplan-Meier survival curves were estimated and compared using the log-rank test. Additional discriminative ability of GPS and of mGPS was evaluated using the area under the curve and Harrell's concordance index.

RESULTS

Hypoalbuminemia was associated with worse survival in both operable and inoperable patients without an elevated CRP. The Kaplan-Meier survival curve of hypoalbuminemic patients without an elevated CRP was more close to the curve of patients with an elevated CRP and a normal albumin than to the curve of patients with neither of these abnormalities. Multivariate analysis, AUC and C-index all indicated that GPS had a higher prognostic value than mGPS.

CONCLUSIONS

Hypoalbuminemia was associated with worse survival in patients with or without an elevated CRP. GPS was superior to mGPS in predicting survival in operable and inoperable NSCLC patients.

Robust BRCA1-like classification of copy number profiles of samples repeated across different datasets and platforms

Breast cancers with BRCA1 germline mutation have a characteristic DNA copy number (CN) pattern. We developed a test that assigns CN profiles to be 'BRCA1-like' or 'non-BRCA1-like', which refers to resembling a BRCA1-mutated tumor or resembling a tumor without a BRCA1 mutation, respectively. Approximately one third of the BRCA1-like breast cancers have a BRCA1 mutation, one third has hypermethylation of the BRCA1 promoter and one third has an unknown reason for being BRCA1-like. This classification is indicative of patients' response to high dose alkylating and platinum containing chemotherapy regimens, which targets the inability of BRCA1 deficient cells to repair DNA double strand breaks. We investigated whether this classification can be reliably obtained with next generation sequencing and copy number platforms other than the bacterial artificial chromosome (BAC) array Comparative Genomic Hybridization (aCGH) on which it was originally developed. We investigated samples from 230 breast cancer patients for which a CN profile had been generated on two to five platforms, comprising low coverage CN sequencing, CN extraction from targeted sequencing panels (CopywriteR), Affymetrix SNP6.0, 135K/720K oligonucleotide aCGH, Affymetrix Oncoscan FFPE (MIP) technology, 3K BAC and 32K BAC aCGH. Pairwise comparison of genomic position-mapped profiles from the original aCGH platform and other platforms revealed concordance. For most cases, biological differences between samples exceeded the differences between platforms within one sample. We observed the same classification across different platforms in over 80% of the patients and kappa values of at least 0.36. Differential classification could be attributed to CN profiles that were not strongly associated to one class. In conclusion, we have shown that the genomic regions that define our BRCA1-like classifier are robustly measured by different CN profiling technologies, providing the possibility to retro- and prospectively investigate BRCA1-like classification across a wide range of CN platforms.

Decoding of the spike timing of primary afferents during voluntary arm movements in monkeys

Understanding the mechanisms of encoding forelimb kinematics in the activity of peripheral afferents is essential for developing a somatosensory neuroprosthesis. To investigate whether the spike timing of dorsal root ganglion (DRG) neurons could be estimated from the forelimb kinematics of behaving monkeys, we implanted two multi-electrode arrays chronically in the DRGs at the level of the cervical segments in two monkeys. Neuronal activity during voluntary reach-to-grasp movements were recorded simultaneously with the trajectories of hand/arm movements, which were tracked in three-dimensional space using a motion capture system. Sixteen and 13 neurons, including muscle spindles, skin receptors, and tendon organ afferents, were recorded in the two monkeys, respectively. We were able to reconstruct forelimb joint kinematics from the temporal firing pattern of a subset of DRG neurons using sparse linear regression (SLiR) analysis, suggesting that DRG neuronal ensembles encoded information about joint kinematics. Furthermore, we estimated the spike timing of the DRG neuronal ensembles from joint kinematics using an integrate-and-fire model (IF) incorporating the SLiR algorithm. The temporal change of firing frequency of a subpopulation of neurons was reconstructed precisely from forelimb kinematics using the SLiR. The estimated firing pattern of the DRG neuronal ensembles encoded forelimb joint angles and velocities as precisely as the originally recorded neuronal activity. These results suggest that a simple model can be used to generate an accurate estimate of the spike timing of DRG neuronal ensembles from forelimb joint kinematics, and is useful for designing a proprioceptive decoder in a brain machine interface.

Evaluating the repair of DNA derived from formalin-fixed paraffin-embedded tissues prior to genomic profiling by SNP-CGH analysis

Pathology archives contain vast resources of clinical material in the form of formalin-fixed paraffin-embedded (FFPE) tissue samples. Owing to the methods of tissue fixation and storage, the integrity of DNA and RNA available from FFPE tissue is compromized, which means obtaining informative data regarding epigenetic, genomic, and expression alterations can be challenging. Here, we have investigated the utility of repairing damaged DNA derived from FFPE tumors prior to single-nucleotide polymorphism (SNP) arrays for whole-genome DNA copy number analysis. DNA was extracted from FFPE samples spanning five decades, involving tumor material obtained from surgical specimens and postmortems. Various aspects of the protocol were assessed, including the method of DNA extraction, the role of Quality Control quantitative PCR (qPCR) in predicting sample success, and the effect of DNA restoration on assay performance, data quality, and the prediction of copy number aberrations (CNAs). DNA that had undergone the repair process yielded higher SNP call rates, reduced log R ratio variance, and improved calling of CNAs compared with matched FFPE DNA not subjected to repair. Reproducible mapping of genomic break points and detection of focal CNAs representing high-level gains and homozygous deletions (HD) were possible, even on autopsy material obtained in 1974. For example, DNA amplifications at the ERBB2 and EGFR gene loci and a HD mapping to 13q14.2 were validated using immunohistochemistry, in situ hybridization, and qPCR. The power of SNP arrays lies in the detection of allele-specific aberrations; however, this aspect of the analysis remains challenging, particularly in the distinction between loss of heterozygosity (LOH) and copy neutral LOH. In summary, attempting to repair DNA that is damaged during fixation and storage may be a useful pretreatment step for genomic studies of large archival FFPE cohorts with long-term follow-up or for understanding rare cancer types, where fresh frozen material is scarce.

Platform comparisons for identification of breast cancers with a BRCA-like copy number profile

Previously, we employed bacterial artificial chromosome (BAC) array comparative genomic hybridization (aCGH) profiles from BRCA1 and -2 mutation carriers and sporadic tumours to construct classifiers that identify tumour samples most likely to harbour BRCA1 and -2 mutations, designated 'BRCA1 and -2-like' tumours, respectively. The classifiers are used in clinical genetics to evaluate unclassified variants, and patients for which no good quality germline DNA is available. Furthermore, we have shown that breast cancer patients with BRCA-like tumour aCGH profiles benefit substantially from platinum-based chemotherapy, potentially due to their inability to repair DNA double strand breaks (DSB), providing a further important clinical application for the classifiers. The BAC array technology has been replaced with oligonucleotide arrays. To continue clinical use of existing classifiers, we mapped oligonucleotide aCGH data to the BAC domain, such that the oligonucleotide profiles can be employed as in the BAC classifier. We demonstrate that segmented profiles derived from oligonucleotide aCGH show high correlation with BAC aCGH profiles. Furthermore, we trained a support vector machine score to objectify aCGH profile quality. Using the mapped oligonucleotide aCGH data, we show equivalence in classification of biologically relevant cases between BAC and oligonucleotide data. Furthermore, the predicted benefit of DSB inducing chemotherapy due to a homologous recombination defect is retained. We conclude that oligonucleotide aCGH data can be mapped to and used in the previously developed and validated BAC aCGH classifiers. Our findings suggest that it is possible to map copy number data from any other technology in a similar way.

Exploring the genetic landscape in chronic lymphocytic leukemia using high-resolution technologies

Abstract During recent years, microarray-based technologies and next-generation sequencing (NGS) have been applied in chronic lymphocytic leukemia (CLL) in order to identify novel genomic aberrations that may contribute to the pathogenesis of the disease. Even though high-resolution microarray studies have confirmed the importance of the known recurrent aberrations, i.e. del(11q), trisomy 12, del(13q) and del(17p), and have more precisely delineated the genomic borders of these aberrations, only a few novel aberrations, found at a low frequency, have been detected with these techniques. In contrast to this, the application of NGS technology of the coding genome (exome sequencing) or the entire genome (whole-genome sequencing) has unveiled a number of novel recurrent mutations in e.g. the NOTCH1, SF3B1 and BIRC3 genes. Importantly, mutations in these latter genes were reported to be associated with a particularly poor outcome, similar to TP53 aberrations, and may play key roles in tumor development, treatment resistance and prognosis. In this review, we not only summarize the latest achievements using array-based or NGS technologies, but also point to new directions for research aiming to unravel the complex genetic "map" in CLL and its prognostic subsets.

Comprehensive DNA copy number profile and BAC library construction of an Indian individual

Bacterial artificial chromosomes (BACs) are used in genomic variation studies due to their capacity to carry a large insert, their high clonal stability, low rate of chimerism and ease of manipulation. In the present study, an attempt was made to create the first genomic BAC library of an anonymous Indian male (IMBL4) consisting of 100,224 clones covering the human genome more than three times. Restriction mapping of 255 BAC clones by pulse field gel electrophoresis confirmed an average insert size of 120 kb. The library was screened by PCR using SHANK3 (SH3 and multiple ankyrin repeat domains 3) and OLFM3 (olfactomedin 3) specific primers. A selection of clones was analyzed by fluorescent in situ hybridization (FISH) and sequencing. Fine mapping of copy number variable regions by array based comparative genomic hybridization identified 467 CNVRs in the IMBL4 genome. The IMBL4 BAC library represents the first cataloged Indian genome resource for applications in basic and clinical research.

Comprehensive assessment of array-based platforms and calling algorithms for detection of copy number variants

We have systematically compared copy number variant (CNV) detection on eleven microarrays to evaluate data quality and CNV calling, reproducibility, concordance across array platforms and laboratory sites, breakpoint accuracy and analysis tool variability. Different analytic tools applied to the same raw data typically yield CNV calls with <50% concordance. Moreover, reproducibility in replicate experiments is <70% for most platforms. Nevertheless, these findings should not preclude detection of large CNVs for clinical diagnostic purposes because large CNVs with poor reproducibility are found primarily in complex genomic regions and would typically be removed by standard clinical data curation. The striking differences between CNV calls from different platforms and analytic tools highlight the importance of careful assessment of experimental design in discovery and association studies and of strict data curation and filtering in diagnostics. The CNV resource presented here allows independent data evaluation and provides a means to benchmark new algorithms.

Array comparative genomic hybridization analysis identifies recurrent gain of chromosome 2p25.3 involving the ACP1 and MYCN genes in chronic lymphocytic leukemia

Chromosomal aberrations are independent prognostic markers in chronic lymphocytic leukemia (CLL). Recent studies using genomic arrays have shown recurrent gains of the short arm of chromosome 2 (2p) in a subset of CLL. We evaluated 178 CLL cases for 2p gains using custom-designed oligonucleotide array-based comparative genomic hybridization (aCGH). A high frequency of 2p gains was observed in 53 of 178 (30%) cases, which ranged from a small 29-kb region to large segments involving the entire short arm. Besides several common chromosomal aberrations associated with 2p gain, we demonstrated a novel observation that gain of the telomeric region 2p25.3 harboring the ACP1 gene is common in CLL (25%, 44 of 178 cases). The ACP1 gene has been previously shown to regulate T-cell receptor signaling through ZAP-70, and both genes are unfavorable clinical markers for CLL. Quantitative polymerase chain reaction (qPCR) confirmed the presence of 3-6 copies of ACP1 in 35 of 40 (88%) of these cases. Interestingly, none of the aCGH diploid CLL cases showed gain of ACP1. Assessment of 73 healthy individuals by qPCR revealed ACP1 copy number gain in only two cases (2.7%). Gain of 2p25.3 was associated with ZAP-70 expression (P < .002) and unmutated immunoglobulin heavy chain variable (IGHV) gene mutation (P < .0001). A high frequency of MYCN co-amplication with ACP1 was observed (14 of 40 cases, 35%). The frequent 2p25.3 gain involving the ACP1 and MYCN genes may help define the critical region of 2p that contributes to pathogenesis of CLL together with other chromosomal abnormalities.

Comparison of genome-wide array genomic hybridization platforms for the detection of copy number variants in idiopathic mental retardation

BACKGROUND

Clinical laboratories are adopting array genomic hybridization as a standard clinical test. A number of whole genome array genomic hybridization platforms are available, but little is known about their comparative performance in a clinical context.

METHODS

We studied 30 children with idiopathic MR and both unaffected parents of each child using Affymetrix 500 K GeneChip SNP arrays, Agilent Human Genome 244 K oligonucleotide arrays and NimbleGen 385 K Whole-Genome oligonucleotide arrays. We also determined whether CNVs called on these platforms were detected by Illumina Hap550 beadchips or SMRT 32 K BAC whole genome tiling arrays and tested 15 of the 30 trios on Affymetrix 6.0 SNP arrays.

RESULTS

The Affymetrix 500 K, Agilent and NimbleGen platforms identified 3061 autosomal and 117 X chromosomal CNVs in the 30 trios. 147 of these CNVs appeared to be de novo, but only 34 (22%) were found on more than one platform. Performing genotype-phenotype correlations, we identified 7 most likely pathogenic and 2 possibly pathogenic CNVs for MR. All 9 of these putatively pathogenic CNVs were detected by the Affymetrix 500 K, Agilent, NimbleGen and the Illumina arrays, and 5 were found by the SMRT BAC array. Both putatively pathogenic CNVs identified in the 15 trios tested with the Affymetrix 6.0 were identified by this platform.

CONCLUSIONS

Our findings demonstrate that different results are obtained with different platforms and illustrate the trade-off that exists between sensitivity and specificity. The large number of apparently false positive CNV calls on each of the platforms supports the need for validating clinically important findings with a different technology.

Genome structural variation discovery and genotyping

Comparisons of human genomes show that more base pairs are altered as a result of structural variation - including copy number variation - than as a result of point mutations. Here we review advances and challenges in the discovery and genotyping of structural variation. The recent application of massively parallel sequencing methods has complemented microarray-based methods and has led to an exponential increase in the discovery of smaller structural-variation events. Some global discovery biases remain, but the integration of experimental and computational approaches is proving fruitful for accurate characterization of the copy, content and structure of variable regions. We argue that the long-term goal should be routine, cost-effective and high quality de novo assembly of human genomes to comprehensively assess all classes of structural variation.

Evaluation of high-resolution microarray platforms for genomic profiling of bone tumours

Background

Several high-density oligonucleotide microarray platforms are available for genome-wide single nucleotide polymorphism (SNP) detection and microarray-based comparative genomic hybridisation (array CGH), which may be used to detect copy number aberrations in human tumours. As part of the EuroBoNeT network of excellence for research on bone tumours (eurobonet.eu), we have evaluated four different commercial high-resolution microarray platforms in order to identify the most appropriate technology for mapping DNA copy number aberrations in such tumours.

Findings

DNA from two different cytogenetically well-characterized bone sarcoma cell lines, representing a simple and a complex karyotype, respectively, was tested in duplicate on four high-resolution microarray platforms; Affymetrix Genome-Wide Human SNP Array 6.0, Agilent Human Genome CGH 244A, Illumina HumanExon510s-duo and Nimblegen HG18 CGH 385 k WG tiling v1.0. The data was analysed using the platform-specific analysis software, as well as a platform-independent analysis algorithm. DNA copy number was measured at six specific chromosomes or chromosomal regions, and compared with the expected ratio based on available cytogenetic information. All platforms performed well in terms of reproducibility and were able to delimit and score small amplifications and deletions at similar resolution, but Agilent microarrays showed better linearity and dynamic range. The platform-specific analysis software provided with each platform identified in general correct copy numbers, whereas using a platform-independent analysis algorithm, correct copy numbers were determined mainly for Agilent and Affymetrix microarrays.

Conclusions

All platforms performed reasonably well, but Agilent microarrays showed better dynamic range, and like Affymetrix microarrays performed well with the platform-independent analysis software, implying more robust data. Bone tumours like osteosarcomas are heterogeneous tumours with complex karyotypes that may be difficult to interpret, and it is of importance to be able to well separate the copy number levels and detect copy number changes in subpopulations. Taking all this into consideration, the Agilent and Affymetrix microarray platforms were found to be a better choice for mapping DNA copy numbers in bone tumours, the latter having the advantage of also providing heterozygosity information.

The pitfalls of platform comparison: DNA copy number array technologies assessed

BACKGROUND

The accurate and high resolution mapping of DNA copy number aberrations has become an important tool by which to gain insight into the mechanisms of tumourigenesis. There are various commercially available platforms for such studies, but there remains no general consensus as to the optimal platform. There have been several previous platform comparison studies, but they have either described older technologies, used less-complex samples, or have not addressed the issue of the inherent biases in such comparisons. Here we describe a systematic comparison of data from four leading microarray technologies (the Affymetrix Genome-wide SNP 5.0 array, Agilent High-Density CGH Human 244A array, Illumina HumanCNV370-Duo DNA Analysis BeadChip, and the Nimblegen 385 K oligonucleotide array). We compare samples derived from primary breast tumours and their corresponding matched normals, well-established cancer cell lines, and HapMap individuals. By careful consideration and avoidance of potential sources of bias, we aim to provide a fair assessment of platform performance.

RESULTS

By performing a theoretical assessment of the reproducibility, noise, and sensitivity of each platform, notable differences were revealed. Nimblegen exhibited between-replicate array variances an order of magnitude greater than the other three platforms, with Agilent slightly outperforming the others, and a comparison of self-self hybridizations revealed similar patterns. An assessment of the single probe power revealed that Agilent exhibits the highest sensitivity. Additionally, we performed an in-depth visual assessment of the ability of each platform to detect aberrations of varying sizes. As expected, all platforms were able to identify large aberrations in a robust manner. However, some focal amplifications and deletions were only detected in a subset of the platforms.

CONCLUSION

Although there are substantial differences in the design, density, and number of replicate probes, the comparison indicates a generally high level of concordance between platforms, despite differences in the reproducibility, noise, and sensitivity. In general, Agilent tended to be the best aCGH platform and Affymetrix, the superior SNP-CGH platform, but for specific decisions the results described herein provide a guide for platform selection and study design, and the dataset a resource for more tailored comparisons.

Comparison of comparative genomic hybridization technologies across microarray platforms

In the 2007 Association of Biomolecular Resource Facilities Microarray Research Group project, we analyzed HL-60 DNA with five platforms: Agilent, Affymetrix 500K, Affymetrix U133 Plus 2.0, Illumina, and RPCI 19K BAC arrays. Copy number variation was analyzed using circular binary segmentation (CBS) analysis of log ratio scores from four independently assessed hybridizations of each platform. Data obtained from these platforms were assessed for reproducibility and the ability to detect formerly reported copy number variations in HL-60. In HL-60, all of the tested platforms detected genomic DNA amplification of the 8q24 locus, trisomy 18, and monosomy X; and deletions at loci 5q11.2~q31, 9p21.3~p22, 10p12~p15, 14q22~q31, and 17p12~p13.3. In the HL-60 genome, at least two of the five platforms detected five novel losses and five novel gains. This report provides guidance in the selection of platforms based on this wide-ranging evaluation of available CGH platforms.

Bacterial artificial chromosome array-based comparative genomic hybridization using paired formalin-fixed, paraffin-embedded and fresh frozen tissue specimens in multiple myeloma

BACKGROUND

Multiple myeloma (MM) is a neoplasm of malignant plasma cells that often harbors many chromosomal aberrations. Currently, fresh frozen tissues (FT) are considered the most reliable for molecular genetic analysis; however, formalin-fixed, paraffin-embedded (FFPE) tissues are easily retrievable. Compared with conventional cytogenetics, bacterial artificial chromosome (BAC) array-comparative genomic hybridization (CGH) allows more sensitive detection of chromosomal abnormalities.

METHODS

The authors analyzed 7 paired FT and FFPE samples of bone marrow aspirate materials obtained from patients with MM in parallel to determine the efficacy of BAC array-CGH using FFPE.

RESULTS

Thirty-four aberrations were identified, including 29 that were observed in both sample types, yielding 85% concordance. Nonrandom anomalies, including gains on 7q, 9q, 15q, and 19p and losses on 8p and 13q, were observed in paired samples from at least 2 patients. To verify these results, fluorescence in situ hybridization (FISH) was performed using probes specific for 7q and 15q, and gains were observed in the 4 samples that were examined. Furthermore, 1 of 3 samples from patients who had monoclonal gammopathy of undetermined significance that were tested also carried gain on 7q, suggesting that this aberration may be an early transforming event.

CONCLUSIONS

The current results indicated that BAC array-CGH can be effective using FFPE samples and is a sensitive method for the identification of nonrandom chromosomal aberrations in MM.

Screening for copy-number alterations and loss of heterozygosity in chronic lymphocytic leukemia--a comparative study of four differently designed, high resolution microarray platforms

Screening for gene copy-number alterations (CNAs) has improved by applying genome-wide microarrays, where SNP arrays also allow analysis of loss of heterozygozity (LOH). We here analyzed 10 chronic lymphocytic leukemia (CLL) samples using four different high-resolution platforms: BAC arrays (32K), oligonucleotide arrays (185K, Agilent), and two SNP arrays (250K, Affymetrix and 317K, Illumina). Cross-platform comparison revealed 29 concordantly detected CNAs, including known recurrent alterations, which confirmed that all platforms are powerful tools when screening for large aberrations. However, detection of 32 additional regions present in 2-3 platforms illustrated a discrepancy in detection of small CNAs, which often involved reported copy-number variations. LOH analysis using dChip revealed concordance of mainly large regions, but showed numerous, small nonoverlapping regions and LOH escaping detection. Evaluation of baseline variation and copy-number ratio response showed the best performance for the Agilent platform and confirmed the robustness of BAC arrays. Accordingly, these platforms demonstrated a higher degree of platform-specific CNAs. The SNP arrays displayed higher technical variation, although this was compensated by high density of elements. Affymetrix detected a higher degree of CNAs compared to Illumina, while the latter showed a lower noise level and higher detection rate in the LOH analysis. Large-scale studies of genomic aberrations are now feasible, but new tools for LOH analysis are requested.

Bacterial artificial chromosome-emulation oligonucleotide arrays for targeted clinical array-comparative genomic hybridization analyses

PURPOSE

The goal of this work was to test the ability of oligonucleotide-based arrays to reproduce the results of focused bacterial artificial chromosome (BAC)-based arrays used clinically in comparative genomic hybridization experiments to detect constitutional copy number changes in genomic DNA.

METHODS

Custom oligonucleotide (oligo) arrays were designed using the Agilent Technologies platform to give high-resolution coverage of regions within the genome sequence coordinates of BAC/P1 artificial chromosome (PAC) clones that had already been validated for use in previous versions of clone arrays used in clinical practice. Standard array-comparative genomic hybridization experiments, including a simultaneous blind analysis of a set of clinical samples, were conducted on both array platforms to identify copy number differences between patient samples and normal reference controls.

RESULTS

Initial experiments successfully demonstrated the capacity of oligo arrays to emulate BAC data without the need for dye-reversal comparisons. Empirical data and computational analyses of oligo response and distribution from a pilot array were used to design an optimized array of 44,000 oligos (44K). This custom 44K oligo array consists of probes localized to the genomic positions of >1400 fluorescence in situ hybridization-verified BAC/PAC clones covering more than 140 regions implicated in genetic diseases, as well as all clinically relevant subtelomeric and pericentromeric regions.

CONCLUSIONS

Our data demonstrate that oligo-based arrays offer a valid alternative for focused BAC arrays. Furthermore, they have significant advantages, including better design flexibility, avoidance of repetitive sequences, manufacturing processes amenable to good manufacturing practice standards in the future, increased robustness because of an enhanced dynamic range (signal to background), and increased resolution that allows for detection of smaller regions of change.

High-resolution aCGH and expression profiling identifies a novel genomic subtype of ER negative breast cancer

High resolution array-CGH and expression profiling identifies a novel genomic subtype of ER negative breast cancer, and provides a genome-wide list of common copy number alterations associated with aberrant expression and poor prognosis.

Background

The characterization of copy number alteration patterns in breast cancer requires high-resolution genome-wide profiling of a large panel of tumor specimens. To date, most genome-wide array comparative genomic hybridization studies have used tumor panels of relatively large tumor size and high Nottingham Prognostic Index (NPI) that are not as representative of breast cancer demographics.

Results

We performed an oligo-array-based high-resolution analysis of copy number alterations in 171 primary breast tumors of relatively small size and low NPI, which was therefore more representative of breast cancer demographics. Hierarchical clustering over the common regions of alteration identified a novel subtype of high-grade estrogen receptor (ER)-negative breast cancer, characterized by a low genomic instability index. We were able to validate the existence of this genomic subtype in one external breast cancer cohort. Using matched array expression data we also identified the genomic regions showing the strongest coordinate expression changes ('hotspots'). We show that several of these hotspots are located in the phosphatome, kinome and chromatinome, and harbor members of the 122-breast cancer CAN-list. Furthermore, we identify frequently amplified hotspots on 8q22.3 (EDD1, WDSOF1), 8q24.11-13 (THRAP6, DCC1, SQLE, SPG8) and 11q14.1 (NDUFC2, ALG8, USP35) associated with significantly worse prognosis. Amplification of any of these regions identified 37 samples with significantly worse overall survival (hazard ratio (HR) = 2.3 (1.3-1.4) p = 0.003) and time to distant metastasis (HR = 2.6 (1.4-5.1) p = 0.004) independently of NPI.

Conclusion

We present strong evidence for the existence of a novel subtype of high-grade ER-negative tumors that is characterized by a low genomic instability index. We also provide a genome-wide list of common copy number alteration regions in breast cancer that show strong coordinate aberrant expression, and further identify novel frequently amplified regions that correlate with poor prognosis. Many of the genes associated with these regions represent likely novel oncogenes or tumor suppressors.

Promoter methylation correlates with reduced Smad4 expression in advanced prostate cancer

BACKGROUND

Transforming growth factor-beta (TGF-beta) is a potent growth inhibitor in a wide range of cell types. A transducer of TGF-beta signaling known as Mothers against decapentaplegic homologue 4 (Smad4) is a known tumor suppressor found on chromosome 18q21.1 and is typically inactivated by deletion or mutation in pancreatic and colorectal cancers. The purpose of the article is to investigate Smad4 expression, gene copy number and methylation status in advanced cases of prostate cancer.

METHODS

We have employed Methylation Specific PCR (MSP) to identify methylation sites within the Smad4 promoter and combined this with quantitative real-time PCR to look for correlates between methylation status and Smad4 expression and to examine androgen receptor (AR) expression. Bacterial artificial chromosome-comparative genomic hybridization (BAC-CGH) has been used to look for genomic amplifications and deletions which may also contribute to expression changes.

RESULTS

We fail to find evidence of genomic deletions or amplifications affecting the Smad4 locus on chromosome 18 but show a correlation between promoter methylation and the loss of Smad4 expression in the same material. We confirm that the AR locus on the X chromosome is amplified in 30% of the advanced clinical samples and that this correlates with increased transcript levels as previously reported by other groups.

CONCLUSION

This indicates that epigenetic changes affect the expression of the Smad4 protein in prostate cancer and points to methylation of the promoter as a novel marker of and contributor to the disease warranting further study.

High-resolution array comparative genomic hybridization analysis of human bronchial and salivary adenoid cystic carcinoma

Adenoid cystic carcinoma (ACC) is a rare but distinctive tumor. Oligonucleotide array comparative genomic hybridization has been applied for cataloging genomic copy number alterations (CNAs) in 17 frozen salivary or bronchial tumors. Only four whole chromosome CNAs were found, and most cases had 2-4 segmental CNAs. No high level amplification was observed. There were recurrent gains at 7p15.2, 17q21-25, and 22q11-13, and recurrent losses at 1p35, 6q22-25, 8q12-13, 9p21, 12q12-13, and 17p11-13. The minimal region of gain at 7p15.2 contained the HOXA cluster. The minimal common regions of deletions contained the CDKN2A/CDKN2B, TP53, and LIMA1 tumor suppressor genes. The recurrent deletion at 8q12.3-13.1 contained no straightforward tumor suppressor gene, but the MIRN124A2 microRNA gene, whose product regulates MMP2 and CDK6. Among unique CNAs, gains harbored CCND1, KIT/PDGFRA/KDR, MDM2, and JAK2. The CNAs involving CCND1, MDM2, KIT, CDKN2A/2B, and TP53 were validated by FISH and/or multiplex ligation-dependent probe amplification. Although most tumors overexpressed cyclin D1 compared with surrounding glands, the only case to overexpress MDM2 had the corresponding CNA. In conclusion, our report suggests that ACC is characterized by a relatively low level of structural complexity. Array CGH and immunohistochemical data implicate MDM2 as the oncogene targeted at 12q15. The gain at 4q12 warrants further exploration as it contains a cluster of receptor kinase genes (KIT/PDGFRA/KDR), whose products can be responsive to specific therapies.

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.