Genome-Wide Copy Number Analysis in Esophageal Adenocarcinoma Using High-Density Single-Nucleotide Polymorphism Arrays

We applied whole-genome single-nucleotide polymorphism arrays to define a comprehensive genetic profile of 23 esophageal adenocarcinoma (EAC) primary tumor biopsies based on loss of heterozygosity (LOH) and DNA copy number changes. Alterations were common, averaging 97 (range, 23-208) per tumor. LOH and gains averaged 33 (range, 3-83) and 31 (range, 11-73) per tumor, respectively. Copy neutral LOH events averaged 27 (range, 7-57) per EAC. We noted 126 homozygous deletions (HD) across the EAC panel (range, 0-11 in individual tumors). Frequent HDs within FHIT (17 of 23), WWOX (8 of 23), and DMD (6 of 23) suggest a role for common fragile sites or genomic instability in EAC etiology. HDs were also noted for known tumor suppressor genes (TSG), including CDKN2A, CDKN2B, SMAD4, and GALR1, and identified PDE4D and MGC48628 as potentially novel TSGs. All tumors showed LOH for most of chromosome 17p, suggesting that TSGs other than TP53 may be targeted. Frequent gains were noted around MYC (13 of 23), BCL9 (12 of 23), CTAGE1 (14 of 23), and ZNF217 (12 of 23). Thus, we have confirmed previous reports indicating frequent changes to FHIT, CDKN2A, TP53, and MYC in EAC and identified additional genes of interest. Meta-analysis of previous genome-wide EAC studies together with the data presented here highlighted consistent regions of gain on 8q, 18q, and 20q and multiple LOH regions on 4q, 5q, 17p, and 18q, suggesting that more than one gene may be targeted on each of these chromosome arms. The focal gains and deletions documented here are a step toward identifying the key genes involved in EAC development

SiDCoN: A Tool to Aid Scoring of DNA Copy Number Changes in SNP Chip Data

The recent application of genome-wide, single nucleotide polymorphism (SNP) microarrays to investigate DNA copy number aberrations in cancer has provided unparalleled sensitivity for identifying genomic changes. In some instances the complexity of these changes makes them difficult to interpret, particularly when tumour samples are contaminated with normal (stromal) tissue. Current automated scoring algorithms require considerable manual data checking and correction, especially when assessing uncultured tumour specimens. To address these limitations we have developed a visual tool to aid in the analysis of DNA copy number data. Simulated DNA Copy Number (SiDCoN) is a spreadsheet-based application designed to simulate the appearance of B-allele and logR plots for all known types of tumour DNA copy number changes, in the presence or absence of stromal contamination. The system allows the user to determine the level of stromal contamination, as well as specify up to 3 different DNA copy number aberrations for up to 5000 data points (representing individual SNPs). This allows users great flexibility to assess simple or complex DNA copy number combinations. We demonstrate how this utility can be used to estimate the level of stromal contamination within tumour samples and its application in deciphering the complex heterogeneous copy number changes we have observed in a series of tumours. We believe this tool will prove useful to others working in the area, both as a training tool, and to aid in the interpretation of complex copy number changes

Similarity of aberrant DNA methylation in Barrett's esophagus and esophageal adenocarcinoma

Background Barrett's esophagus (BE) is the metaplastic replacement of squamous with columnar epithelium in the esophagus, as a result of reflux. It is the major risk factor for the development of esophageal adenocarcinoma (EAC). Methylation of CpG dinucleotides of normally unmethylated genes is associated with silencing of their expression, and is common in EAC. This study was designed to determine at what stage, in the progression from BE to EAC, methylation of key genes occurs. Results We examined nine genes (APC, CDKN2A, ID4, MGMT, RBP1, RUNX3, SFRP1, TIMP3, and TMEFF2), frequently methylated in multiple cancer types, in a panel of squamous (19 biopsies from patients without BE or EAC, 16 from patients with BE, 21 from patients with EAC), BE (40 metaplastic, seven high grade dysplastic) and 37 EAC tissues. The methylation frequency, the percentage of samples that had any extent of methylation, for each of the nine genes in the EAC (95%, 59%, 76%, 57%, 70%, 73%, 95%, 74% and 83% respectively) was significantly higher than in any of the squamous groups. The methylation frequency for each of the nine genes in the metaplastic BE (95%, 28%, 78%, 48%, 58%, 48%, 93%, 88% and 75% respectively) was significantly higher than in the squamous samples except for CDKN2A and RBP1. The methylation frequency did not differ between BE and EAC samples, except for CDKN2A and RUNX3 which were significantly higher in EAC. The methylation extent was an estimate of both the number of methylated alleles and the density of methylation on these alleles. This was significantly greater in EAC than in metaplastic BE for all genes except APC, MGMT and TIMP3. There was no significant difference in methylation extent for any gene between high grade dysplastic BE and EAC. Conclusion We found significant methylation in metaplastic BE, which for seven of the nine genes studied did not differ in frequency from that found in EAC. This is also the first report of gene silencing by methylation of ID4 in BE or EAC. This study suggests that metaplastic BE is a highly abnormal tissue, more similar to cancer tissue than to normal epithelium.Eric Smith, Neville J De Young, Sandra J Pavey, Nicholas K Hayward, Derek J Nancarrow, David C Whiteman, B Mark Smithers, Andrew R Ruszkiewicz, Andrew D Clouston, David C Gotley, Peter G Devitt, Glyn G Jamieson and Paul A Dre

Improved prediction of radiation pneumonitis by combining biological and radiobiological parameters using a data-driven Bayesian network analysis

Grade 2 and higher radiation pneumonitis (RP2) is a potentially fatal toxicity that limits efficacy of radiation therapy (RT). We wished to identify a combined biomarker signature of circulating miRNAs and cytokines which, along with radiobiological and clinical parameters, may better predict a targetable RP2 pathway. In a prospective clinical trial of response-adapted RT for patients (n = 39) with locally advanced non-small cell lung cancer, we analyzed patients\u27 plasma, collected pre- and during RT, for microRNAs (miRNAs) and cytokines using array and multiplex enzyme linked immunosorbent assay (ELISA), respectively. Interactions between candidate biomarkers, radiobiological, and clinical parameters were analyzed using data-driven Bayesian network (DD-BN) analysis. We identified alterations in specific miRNAs (miR-532, -99b and -495, let-7c, -451 and -139-3p) correlating with lung toxicity. High levels of soluble tumor necrosis factor alpha receptor 1 (sTNFR1) were detected in a majority of lung cancer patients. However, among RP patients, within 2 weeks of RT initiation, we noted a trend of temporary decline in sTNFR1 (a physiological scavenger of TNFα) and ADAM17 (a shedding protease that cleaves both membrane-bound TNFα and TNFR1) levels. Cytokine signature identified activation of inflammatory pathway. Using DD-BN we combined miRNA and cytokine data along with generalized equivalent uniform dose (gEUD) to identify pathways with better accuracy of predicting RP2 as compared to either miRNA or cytokines alone. This signature suggests that activation of the TNFα-NFκB inflammatory pathway plays a key role in RP which could be specifically ameliorated by etanercept rather than current therapy of non-specific leukotoxic corticosteroids

Whole Genome Expression Array Profiling Highlights Differences in Mucosal Defense Genes in Barrett's Esophagus and Esophageal Adenocarcinoma

Esophageal adenocarcinoma (EAC) has become a major concern in Western countries due to rapid rises in incidence coupled with very poor survival rates. One of the key risk factors for the development of this cancer is the presence of Barrett's esophagus (BE), which is believed to form in response to repeated gastro-esophageal reflux. In this study we performed comparative, genome-wide expression profiling (using Illumina whole-genome Beadarrays) on total RNA extracted from esophageal biopsy tissues from individuals with EAC, BE (in the absence of EAC) and those with normal squamous epithelium. We combined these data with publically accessible raw data from three similar studies to investigate key gene and ontology differences between these three tissue states. The results support the deduction that BE is a tissue with enhanced glycoprotein synthesis machinery (DPP4, ATP2A3, AGR2) designed to provide strong mucosal defenses aimed at resisting gastro-esophageal reflux. EAC exhibits the enhanced extracellular matrix remodeling (collagens, IGFBP7, PLAU) effects expected in an aggressive form of cancer, as well as evidence of reduced expression of genes associated with mucosal (MUC6, CA2, TFF1) and xenobiotic (AKR1C2, AKR1B10) defenses. When our results are compared to previous whole-genome expression profiling studies keratin, mucin, annexin and trefoil factor gene groups are the most frequently represented differentially expressed gene families. Eleven genes identified here are also represented in at least 3 other profiling studies. We used these genes to discriminate between squamous epithelium, BE and EAC within the two largest cohorts using a support vector machine leave one out cross validation (LOOCV) analysis. While this method was satisfactory for discriminating squamous epithelium and BE, it demonstrates the need for more detailed investigations into profiling changes between BE and EAC

Cross-Platform Array Screening Identifies COL1A2, THBS1, TNFRSF10D and UCHL1 as Genes Frequently Silenced by Methylation in Melanoma

Epigenetic regulation of tumor suppressor genes (TSGs) has been shown to play a central role in melanomagenesis. By integrating gene expression and methylation array analysis we identified novel candidate genes frequently methylated in melanoma. We validated the methylation status of the most promising genes using highly sensitive Sequenom Epityper assays in a large panel of melanoma cell lines and resected melanomas, and compared the findings with those from cultured melanocytes. We found transcript levels of UCHL1, COL1A2, THBS1 and TNFRSF10D were inversely correlated with promoter methylation. For THBS1 and UCHL1 the effect of this methylation on expression was confirmed at the protein level. Identification of these candidate TSGs and future research designed to understand how their silencing is related to melanoma development will increase our understanding of the etiology of this cancer and may provide tools for its early diagnosis

Linkage to chromosome 2q32.2-q33.3 in familial serrated neoplasia (Jass syndrome)

Causative genetic variants have to date been identified for only a small proportion of familial colorectal cancer (CRC). While conditions such as Familial Adenomatous Polyposis and Lynch syndrome have well defined genetic causes, the search for variants underlying the remainder of familial CRC is plagued by genetic heterogeneity. The recent identification of families with a heritable predisposition to malignancies arising through the serrated pathway (familial serrated neoplasia or Jass syndrome) provides an opportunity to study a subset of familial CRC in which heterogeneity may be greatly reduced. A genome-wide linkage screen was performed on a large family displaying a dominantly-inherited predisposition to serrated neoplasia genotyped using the Affymetrix GeneChip Human Mapping 10 K SNP Array. Parametric and nonparametric analyses were performed and resulting regions of interest, as well as previously reported CRC susceptibility loci at 3q22, 7q31 and 9q22, were followed up by finemapping in 10 serrated neoplasia families. Genome-wide linkage analysis revealed regions of interest at 2p25.2-p25.1, 2q24.3-q37.1 and 8p21.2-q12.1. Finemapping linkage and haplotype analyses identified 2q32.2-q33.3 as the region most likely to harbour linkage, with heterogeneity logarithm of the odds (HLOD) 2.09 and nonparametric linkage (NPL) score 2.36 (P = 0.004). Five primary candidate genes (CFLAR, CASP10, CASP8, FZD7 and BMPR2) were sequenced and no segregating variants identified. There was no evidence of linkage to previously reported loci on chromosomes 3, 7 and 9

WebFOG: A web tool to map genomic features onto genes

A large number of new genomic features are being discovered using high throughput techniques. The next challenge is to automatically map them to the reference genome for further analysis and functional annotation. We have developed a tool that can be used to map important genomic features to the latest version of the human genome and also to annotate new features. These genomic features could be of many different source types, including miRNAs, microarray primers or probes, Chip-on-Chip data, CpG islands and SNPs to name a few. A standalone version and web interface for the tool can be accessed through: http://populationhealth.qimr.edu.au/cgi-bin/webFOG/index.cgi. The project details and source code is also available at http://www.bioinformatics.org/webfog. © 2010 Elsevier Inc

Microsatellite stable colorectal cancers stratified by the BRAF V600E mutation show distinct patterns of chromosomal instability

The BRAF (V600E) mutation in colorectal cancers that are microsatellite stable (MSS) confers a poor patient prognosis, whereas BRAF mutant microsatellite-unstable (MSI) colorectal cancers have an excellent prognosis. BRAF wild type cancers are typically MSS and display chromosomal instability (CIN). CIN has not been extensively studied on a genome-wide basis in relation to BRAF mutational status in colorectal cancer. BRAF mutant/MSS ( BRAFmut/MSS) cancers (n = 33) and BRAF mutant/MSI (BRAFmut/MSI) cancers (n = 30) were compared for presence of copy number aberrations (CNAs) indicative of CIN, with BRAF wild type/MSS (BRAFwt/MSS) cancers (n = 18) using Illumina CytoSNP-12 arrays. BRAFmut/MSS and BRAFwt/ MSS cancers showed comparable numbers of CNAs/cancer at 32.8 and 29.8 respectively. However, there were differences in patterns of CNA length between MSS cohorts, with BRAFmut/MSS cancers having significantly greater proportions of focal CNAs compared to BRAFwt/MSS cancers (