An Integrative Analysis to Identify Driver Genes in Esophageal Squamous Cell Carcinoma

BackgroundFew driver genes have been well established in esophageal squamous cell carcinoma (ESCC). Identification of the genomic aberrations that contribute to changes in gene expression profiles can be used to predict driver genes.MethodsWe searched for driver genes in ESCC by integrative analysis of gene expression microarray profiles and copy number data. To narrow down candidate genes, we performed survival analysis on expression data and tested the genetic vulnerability of each genes using public RNAi screening data. We confirmed the results by performing RNAi experiments and evaluating the clinical relevance of candidate genes in an independent ESCC cohort.ResultsWe found 10 significantly recurrent copy number alterations accompanying gene expression changes, including loci 11q13.2, 7p11.2, 3q26.33, and 17q12, which harbored CCND1, EGFR, SOX2, and ERBB2, respectively. Analysis of survival data and RNAi screening data suggested that GRB7, located on 17q12, was a driver gene in ESCC. In ESCC cell lines harboring 17q12 amplification, knockdown of GRB7 reduced the proliferation, migration, and invasion capacities of cells. Moreover, siRNA targeting GRB7 had a synergistic inhibitory effect when combined with trastuzumab, an anti-ERBB2 antibody. Survival analysis of the independent cohort also showed that high GRB7 expression was associated with poor prognosis in ESCC.ConclusionOur integrative analysis provided important insights into ESCC pathogenesis. We identified GRB7 as a novel ESCC driver gene and potential new therapeutic target.     

Identification of Candidate Driver Genes in Common Focal Chromosomal Aberrations of Microsatellite Stable Colorectal Cancer

Colorectal cancer (CRC) is a leading cause of cancer deaths worldwide. Chromosomal instability (CIN) is a major driving force of microsatellite stable (MSS) sporadic CRC. CIN tumours are characterised by a large number of somatic chromosomal copy number aberrations (SCNA) that frequently affect oncogenes and tumour suppressor genes. The main aim of this work was to identify novel candidate CRC driver genes affected by recurrent and focal SCNA. High resolution genome-wide comparative genome hybridisation (CGH) arrays were used to compare tumour and normal DNA for 53 sporadic CRC cases. Context corrected common aberration (COCA) analysis and custom algorithms identified 64 deletions and 32 gains of focal minimal common regions (FMCR) at high frequency (>10%). Comparison of these FMCR with published genomic profiles from CRC revealed common overlap (42.2% of deletions and 34.4% of copy gains). Pathway analysis showed that apoptosis and p53 signalling pathways were commonly affected by deleted FMCR, and MAPK and potassium channel pathways by gains of FMCR. Candidate tumour suppressor genes in deleted FMCR included RASSF3, IFNAR1, IFNAR2 and NFKBIA and candidate oncogenes in gained FMCR included PRDM16, TNS1, RPA3 and KCNMA1. In conclusion, this study confirms some previously identified aberrations in MSS CRC and provides in silico evidence for some novel candidate driver genes.     

Integration of DNA copy number alterations and transcriptional expression analysis in human gastric cancer

Background

Genomic instability with frequent DNA copy number alterations is one of the key hallmarks of carcinogenesis. The chromosomal regions with frequent DNA copy number gain and loss in human gastric cancer are still poorly defined. It remains unknown how the DNA copy number variations contributes to the changes of gene expression profiles, especially on the global level.

Principal Findings

We analyzed DNA copy number alterations in 64 human gastric cancer samples and 8 gastric cancer cell lines using bacterial artificial chromosome (BAC) arrays based comparative genomic hybridization (aCGH). Statistical analysis was applied to correlate previously published gene expression data obtained from cDNA microarrays with corresponding DNA copy number variation data to identify candidate oncogenes and tumor suppressor genes. We found that gastric cancer samples showed recurrent DNA copy number variations, including gains at 5p, 8q, 20p, 20q, and losses at 4q, 9p, 18q, 21q. The most frequent regions of amplification were 20q12 (7/72), 20q12–20q13.1 (12/72), 20q13.1–20q13.2 (11/72) and 20q13.2–20q13.3 (6/72). The most frequent deleted region was 9p21 (8/72). Correlating gene expression array data with aCGH identified 321 candidate oncogenes, which were overexpressed and showed frequent DNA copy number gains; and 12 candidate tumor suppressor genes which were down-regulated and showed frequent DNA copy number losses in human gastric cancers. Three networks of significantly expressed genes in gastric cancer samples were identified by ingenuity pathway analysis.

Conclusions

This study provides insight into DNA copy number variations and their contribution to altered gene expression profiles during human gastric cancer development. It provides novel candidate driver oncogenes or tumor suppressor genes for human gastric cancer, useful pathway maps for the future understanding of the molecular pathogenesis of this malignancy, and the construction of new therapeutic targets.     

Relatively Small Contribution of Methylation and Genomic Copy Number Aberration to the Aberrant Expression of Inflammation-Related Genes in HBV-Related Hepatocellular Carcinoma

BackgroundIt is well known that chronic inflammation plays a pivotal role in the development of hepatitis B virus (HBV) related hepatocellular carcinoma (HCC). However, the causes behind aberrant expression of inflammation-related genes occurred in HCC remain unclear.MethodsWe performed array-based analyses to comprehensively investigate the contributions of DNA methylation and somatic copy number aberration (SCNA) to the aberrant expression of 1,027 inflammation-related genes in 30 HCCs and paired non-tumor tissues. The results were validated in public datasets and an additional sample set of 47 paired HCCs and non-tumor tissues.ResultsWe identified 252 differentially expressed, 125 aberrantly methylated and 287 copy number changed inflammation-related genes. Despite reasonable statistical power, among them, only 11 genes and 56 genes whose aberrant expression was associated with DNA methylation or SCNA, respectively. DNA methylation and SCNA together contributed to less than 30% aberrant expression of inflammation-related genes.ConclusionThese results suggest that molecular mechanisms other than DNA methylation and SCNA might play major role in the regulation of aberrant expression of inflammation-related gene in HBV-related HCCs.     

High-resolution array copy number analyses for detection of deletion,gain, amplification and copy-neutral LOH in primary neuroblastoma tumors: Four cases of homozygous deletions of the CDKN2A gene

Background

Neuroblastoma is a very heterogeneous pediatric tumor of the sympathetic nervous system showing clinically significant patterns of genetic alterations. Favorable tumors usually have near-triploid karyotypes with few structural rearrangements. Aggressive stage 4 tumors often have near-diploid or near-tetraploid karyotypes and structural rearrangements. Whole genome approaches for analysis of genome-wide copy number have been used to analyze chromosomal abnormalities in tumor samples. We have used array-based copy number analysis using oligonucleotide single nucleotide polymorphisms (SNP) arrays to analyze the chromosomal structure of a large number of neuroblastoma tumors of different clinical and biological subsets.

Results

Ninety-two neuroblastoma tumors were analyzed with 50 K and/or 250 K SNP arrays from Affymetrix, using CNAG3.0 software. Thirty percent of the tumors harbored 1p deletion, 22% deletion of 11q, 26% had MYCN amplification and 45% 17q gain. Most of the tumors with 1p deletion were found among those with MYCN amplification. Loss of 11q was most commonly seen in tumors without MYCN amplification. In the case of MYCN amplification, two types were identified. One type displayed simple continuous amplicons; the other type harbored more complex rearrangements. MYCN was the only common gene in all cases with amplification. Complex amplification on chromosome 12 was detected in two tumors and three different overlapping regions of amplification were identified. Two regions with homozygous deletions, four cases with CDKN2A deletions in 9p and one case with deletion on 3p (the gene RBMS3) were also detected in the tumors.

Conclusion

SNP arrays provide useful tools for high-resolution characterization of significant chromosomal rearrangements in neuroblastoma tumors. The mapping arrays from Affymetrix provide both copy number and allele-specific information at a resolution of 10–12 kb. Chromosome 9p, especially the gene CDKN2A, is subject to homozygous (four cases) and heterozygous deletions (five cases) in neuroblastoma tumors.

     

SAAS-CNV: A Joint Segmentation Approach on Aggregated and Allele Specific Signals for the Identification of Somatic Copy Number Alterations with Next-Generation Sequencing Data

Cancer genomes exhibit profound somatic copy number alterations (SCNAs). Studying tumor SCNAs using massively parallel sequencing provides unprecedented resolution and meanwhile gives rise to new challenges in data analysis, complicated by tumor aneuploidy and heterogeneity as well as normal cell contamination. While the majority of read depth based methods utilize total sequencing depth alone for SCNA inference, the allele specific signals are undervalued. We proposed a joint segmentation and inference approach using both signals to meet some of the challenges. Our method consists of four major steps: 1) extracting read depth supporting reference and alternative alleles at each SNP/Indel locus and comparing the total read depth and alternative allele proportion between tumor and matched normal sample; 2) performing joint segmentation on the two signal dimensions; 3) correcting the copy number baseline from which the SCNA state is determined; 4) calling SCNA state for each segment based on both signal dimensions. The method is applicable to whole exome/genome sequencing (WES/WGS) as well as SNP array data in a tumor-control study. We applied the method to a dataset containing no SCNAs to test the specificity, created by pairing sequencing replicates of a single HapMap sample as normal/tumor pairs, as well as a large-scale WGS dataset consisting of 88 liver tumors along with adjacent normal tissues. Compared with representative methods, our method demonstrated improved accuracy, scalability to large cancer studies, capability in handling both sequencing and SNP array data, and the potential to improve the estimation of tumor ploidy and purity.     

Identification of rare genetic variants in novel loci associated with Paget’s disease of bone

In genome-wide association studies, single nucleotide polymorphisms located in five novel loci were associated with PDB. We aimed at identifying rare genetic variants of candidate genes located in these loci and search for genetic association with PDB in the French-Canadian population. Exons, promoter and exon–intron junctions from patients with familial PDB and healthy individuals were sequenced in candidate genes, located within novel loci associated with PDB in our population. Rare variant was defined by a minor allele frequency <0.05 or absent from dbSNP (NCBI). We sequenced seven genes in 1p13 locus, three genes in 7q33, three genes in 8q22, and five genes in 15q24 locus. We identified 126 rare variants in at least one patient with PDB of whom 55 were located in 1p13 locus, 32 in 7q33, 10 in 8q22 and 29 in 15q24 locus. We located 71 of these 126 rare variants in an intron, 30 in an exon and 9 in an untranslated region. 60 % of these variants were located in functionally relevant gene regions. Among the 12 missense rare variants in PDB, two (rs62620995 in TM7SF4; rs62641691 in CD276) were predicted to be damaging by in silico analysis tools. Rs62620995, which altered a conserved amino acid (p.Leu397Phe) in the TM7SF4 gene, encoding the DC-STAMP protein involved in osteoclastogenesis through RANK signaling pathway, was found to have a marginal association with PDB (p = 0.09). Rs35500845, located in the CTHRC1 gene, which encodes a regulator of collagen matrix deposition, was also associated with PDB in the French-Canadian population (p = 0.046).     

Candidate Luminal B Breast Cancer Genes Identified by Genome,Gene Expression and DNA Methylation Profiling

Breast cancers (BCs) of the luminal B subtype are estrogen receptor-positive (ER+), highly proliferative, resistant to standard therapies and have a poor prognosis. To better understand this subtype we compared DNA copy number aberrations (CNAs), DNA promoter methylation, gene expression profiles, and somatic mutations in nine selected genes, in 32 luminal B tumors with those observed in 156 BCs of the other molecular subtypes. Frequent CNAs included 8p11-p12 and 11q13.1-q13.2 amplifications, 7q11.22-q34, 8q21.12-q24.23, 12p12.3-p13.1, 12q13.11-q24.11, 14q21.1-q23.1, 17q11.1-q25.1, 20q11.23-q13.33 gains and 6q14.1-q24.2, 9p21.3-p24,3, 9q21.2, 18p11.31-p11.32 losses. A total of 237 and 101 luminal B-specific candidate oncogenes and tumor suppressor genes (TSGs) presented a deregulated expression in relation with their CNAs, including 11 genes previously reported associated with endocrine resistance. Interestingly, 88% of the potential TSGs are located within chromosome arm 6q, and seven candidate oncogenes are potential therapeutic targets. A total of 100 candidate oncogenes were validated in a public series of 5,765 BCs and the overexpression of 67 of these was associated with poor survival in luminal tumors. Twenty-four genes presented a deregulated expression in relation with a high DNA methylation level. FOXO3, PIK3CA and TP53 were the most frequent mutated genes among the nine tested. In a meta-analysis of next-generation sequencing data in 875 BCs, KCNB2 mutations were associated with luminal B cases while candidate TSGs MDN1 (6q15) and UTRN (6q24), were mutated in this subtype. In conclusion, we have reported luminal B candidate genes that may play a role in the development and/or hormone resistance of this aggressive subtype.     

Loss of retinoblastoma but not p16 function allows bypass of replicative senescence in human fibroblasts

Current models envision replicative senescence to be under dual control by the p53 and retinoblastoma (RB) tumour-suppressor pathways. The role of the p16^INK4a–RB pathway is controversial, and the function of RB in human cells has not been tested directly. We used targeted homologous recombination to knock out one copy of RB in presenescent human fibroblasts. During entry into senescence, RB^+/− cells underwent spontaneous loss of heterozygosity and the resultant RB^−/− clones bypassed senescence. The extended lifespan phase was eventually terminated by a crisis-like state. The same phenotype was documented for p21 ^CIP1/WAF1 and p53 heterozygous cells, indicating that loss of function of all three genes results in failure to establish senescence. By contrast, the abolition of p16 function by the expression of a p16-insensitive cyclin-dependent kinase 4 protein or siRNA-mediated knockdown provided only minimal lifespan extension that was terminated by senescence. We propose that p53, p21 and RB act in a linear genetic pathway to regulate cell entry into replicative senescence.     

Pan cancer patterns of allelic imbalance from chromosomal alterations in 33 tumor types

Somatic copy number alterations (SCNAs) serve as hallmarks of tumorigenesis and often result in deviations from one-to-one allelic ratios at heterozygous loci, leading to allelic imbalance (AI). The Cancer Genome Atlas (TCGA) reports SCNAs identified using a circular binary segmentation algorithm, providing segment mean copy number estimates from single-nucleotide polymorphism DNA microarray total intensities (log R ratio), but not allele-specific intensities (“B allele” frequencies) that inform of AI. Our approach provides more sensitive identification of SCNAs by modeling the “B allele” frequencies jointly, thereby bolstering the catalog of chromosomal alterations in this widely utilized resource. Here we present AI summaries for all 33 tumor sites in TCGA, including those induced by SCNAs and copy-neutral loss-of-heterozygosity (cnLOH). We identified AI in 94% of the tumors, higher than in previous reports. Recurrent events included deletions of 17p, 9q, 3p, amplifications of 8q, 1q, 7p, as well as mixed event types on 8p and 13q. We also observed both site-specific and pan-cancer (spanning 17p) cnLOH, patterns which have not been comprehensively characterized. The identification of such cnLOH events elucidates tumor suppressors and multi-hit pathways to carcinogenesis. We also contrast the landscapes inferred from AI- and total intensity-derived SCNAs and propose an automated procedure to improve and adjust SCNAs in TCGA for cases where high levels of aneuploidy obscured baseline intensity identification. Our findings support the exploration of additional methods for robust automated inference procedures and to aid empirical discoveries across TCGA.     

Molecular-genetic analysis of two cases with retinoblastoma: Benefits for disease management

Effective counselling and management of retinoblastoma families using genetic information is presently practised in many parts of the world. We studied histopathological, chromosomal and molecular-genetic data of two retinoblastoma patients from India. The two patients, one with bilateral and the other with unilateral retinoblastoma, underwent complete ophthalmic examination, cytogenetic study, retinoblastoma gene (RB1) mutational analysis andRB1 promoter region methylation screening. In the bilateral retinoblastoma patient deletion of chromosome region 13q14 in peripheral blood lymphocytes and a hemizygous novel 8-bp deletion in exon 4 ofRB1 in tumour sample were observed. In the unilaterally affected patient CGA to TGA transition protein truncation mutations were observed in exons 8 and 14 ofRB1.     

Modeling and correct the GC bias of tumor and normal WGS data for SCNA based tumor subclonal population inferring

Background

Somatic copy number alternations (SCNAs) can be utilized to infer tumor subclonal populations in whole genome seuqncing studies, where usually their read count ratios between tumor-normal paired samples serve as the inferring proxy. Existing SCNA based subclonal population inferring tools consider the GC bias of tumor and normal sample is of the same fature, and could be fully offset by read count ratio. However, we found that, the read count ratio on SCNA segments presents a Log linear biased pattern, which influence existing read count ratios based subclonal inferring tools performance. Currently no correction tools take into account the read ratio bias.

Results

We present Pre-SCNAClonal, a tool that improving tumor subclonal population inferring by correcting GC-bias at SCNAs level. Pre-SCNAClonal first corrects GC bias using Markov chain Monte Carlo probability model, then accurately locates baseline DNA segments (not containing any SCNAs) with a hierarchy clustering model. We show Pre-SCNAClonal’s superiority to exsiting GC-bias correction methods at any level of subclonal population.

Conclusions

Pre-SCNAClonal could be run independently as well as serving as pre-processing/gc-correction step in conjuntion with exsiting SCNA-based subclonal inferring tools.

     

Somatic Mutation Allelic Ratio Test Using ddPCR (SMART-ddPCR): An Accurate Method for Assessment of Preferential Allelic Imbalance in Tumor DNA

The extent to which heritable genetic variants can affect tumor development has yet to be fully elucidated. Tumor selection of single nucleotide polymorphism (SNP) risk alleles, a phenomenon called preferential allelic imbalance (PAI), has been demonstrated in some cancer types. We developed a novel application of digital PCR termed Somatic Mutation Allelic Ratio Test using Droplet Digital PCR (SMART-ddPCR) for accurate assessment of tumor PAI, and have applied this method to test the hypothesis that heritable SNPs associated with childhood acute lymphoblastic leukemia (ALL) may demonstrate tumor PAI. These SNPs are located at CDKN2A (rs3731217) and IKZF1 (rs4132601), genes frequently lost in ALL, and at CEBPE (rs2239633), ARID5B (rs7089424), PIP4K2A (rs10764338), and GATA3 (rs3824662), genes located on chromosomes gained in high-hyperdiploid ALL. We established thresholds of AI using constitutional DNA from SNP heterozygotes, and subsequently measured allelic copy number in tumor DNA from 19–142 heterozygote samples per SNP locus. We did not find significant tumor PAI at these loci, though CDKN2A and IKZF1 SNPs showed a trend towards preferential selection of the risk allele (p = 0.17 and p = 0.23, respectively). Using a genomic copy number control ddPCR assay, we investigated somatic copy number alterations (SCNA) underlying AI at CDKN2A and IKZF1, revealing a complex range of alterations including homozygous and hemizygous deletions and copy-neutral loss of heterozygosity, with varying degrees of clonality. Copy number estimates from ddPCR showed high agreement with those from multiplex ligation-dependent probe amplification (MLPA) assays. We demonstrate that SMART-ddPCR is a highly accurate method for investigation of tumor PAI and for assessment of the somatic alterations underlying AI. Furthermore, analysis of publicly available data from The Cancer Genome Atlas identified 16 recurrent SCNA loci that contain heritable cancer risk SNPs associated with a matching tumor type, and which represent candidate PAI regions warranting further investigation.     

Somatic Copy Number Alterations Associated with Japanese or Endometriosis in Ovarian Clear Cell Adenocarcinoma

When compared with other epithelial ovarian cancers, the clinical characteristics of ovarian clear cell adenocarcinoma (CCC) include 1) a higher incidence among Japanese, 2) an association with endometriosis, 3) poor prognosis in advanced stages, and 4) a higher incidence of thrombosis as a complication. We used high resolution comparative genomic hybridization (CGH) to identify somatic copy number alterations (SCNAs) associated with each of these clinical characteristics of CCC. The Human Genome CGH 244A Oligo Microarray was used to examine 144 samples obtained from 120 Japanese, 15 Korean, and nine German patients with CCC. The entire 8q chromosome (minimum corrected p-value: q = 0.0001) and chromosome 20q13.2 including the ZNF217 locus (q = 0.0078) were amplified significantly more in Japanese than in Korean or German samples. This copy number amplification of the ZNF217 gene was confirmed by quantitative real-time polymerase chain reaction (Q-PCR). ZNF217 RNA levels were also higher in Japanese tumor samples than in non-Japanese samples (P = 0.027). Moreover, endometriosis was associated with amplification of EGFR gene (q = 0.047), which was again confirmed by Q-PCR and correlated with EGFR RNA expression. However, no SCNAs were significantly associated with prognosis or thrombosis. These results indicated that there may be an association between CCC and ZNF217 amplification among Japanese patients as well as between endometriosis and EGFR gene amplifications.     

Studies on the human retinoblastoma susceptibility gene

The retinoblastoma susceptibility (RB) gene is unique among other cloned cancer genes because its causal role in a human cancer, retinoblastoma, was established by classical genetic methods before its isolation. Earlier hypotheses and experimental data suggested that inactivation of a gene in chromosome band 13q14 resulted in retinoblastoma formation. A gene in this region was identified as the RB gene on the basis of mutations found specifically in retinoblastoma tumors; however, its proposed biological activity in suppressing neoplasia has yet to be demonstrated. The RB gene product was identified as a nuclear phosphoprotein of 110 kD associated with DNA binding activity, suggesting that the RB protein may regulate other genes. Probes for the RB gene and gene product will be useful for genetic diagnosis of retinoblastoma susceptibility in affected families; for direct detection of mutant RB alleles; and, potentially, for genetic diagnosis of susceptibility to osteosarcoma and other tumors tentatively linked to RB-gene dysfunction. Continued study of the RB gene should yield further insight into mechanisms of oncogenesis, development, and gene regulation.