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.     

Whole Exome Sequencing Identifies Novel Recurrently Mutated Genes in Patients with Splenic Marginal Zone Lymphoma

The pathogenesis of splenic marginal zone lymphoma (SMZL) remains largely unknown. Recent high-throughput sequencing studies have identified recurrent mutations in key pathways, most notably NOTCH2 mutations in >25% of patients. These studies are based on small, heterogeneous discovery cohorts, and therefore only captured a fraction of the lesions present in the SMZL genome. To identify further novel pathogenic mutations within related biochemical pathways, we applied whole exome sequencing (WES) and copy number (CN) analysis to a biologically and clinically homogeneous cohort of seven SMZL patients with 7q abnormalities and IGHV1-2*04 gene usage. We identified 173 somatic non-silent variants, affecting 160 distinct genes. In additional to providing independent validation of the presence of mutation in several previously reported genes (NOTCH2, TNFAIP3, MAP3K14, MLL2 and SPEN), our study defined eight additional recurrently mutated genes in SMZL; these genes are CREBBP, CBFA2T3, AMOTL1, FAT4, FBXO11, PLA2G4D, TRRAP and USH2A. By integrating our WES and CN data we identified three mutated putative candidate genes targeted by 7q deletions (CUL1, EZH2 and FLNC), with FLNC positioned within the well-characterized 7q minimally deleted region. Taken together, this work expands the reported directory of recurrently mutated cancer genes in this disease, thereby expanding our understanding of SMZL pathogenesis. Ultimately, this work will help to establish a stratified approach to care including the possibility of targeted therapy.     

Transcriptome and DNA methylation analysis reveals molecular mechanisms underlying intrahepatic cholangiocarcinoma progression

Intrahepatic cholangiocarcinoma (iCCA) is an aggressive malignancy with increasing incidence. It has been suggested that DNA methylation drives cancer development. However, the molecular mechanisms underlying iCCA progression and the roles of DNA methylation still remain elusive. In this study, weighted correlation networks were constructed to identify gene modules and hub genes associated with the tumour stage. We identified 12 gene modules, two of which were significantly positively or negatively related to the tumour stage, respectively. Key hub genes SLC2A1, CDH3 and EFHD2 showed increased expression across the tumour stage and were correlated with poor survival, whereas decrease of FAM171A1, ONECUT1 and PHYHIPL was correlated with better survival. Pathway analysis revealed hedgehog pathway was activated in CDH3 up-regulated tumours, and chromosome separation was elevated in tumours expressing high EFHD2. JAK-STAT pathway was overrepresented in ONECUT1 down-regulated tumours, whereas Rho GTPases-formins signalling was activated in PHYHIPL down-regulated tumours. Finally, significant negative associations between expression of EFHD2, PHYHIPL and promoter DNA methylation were detected, and alterations of DNA methylation were correlated with tumour survival. In summary, we identified key genes and pathways that may participate in progression of iCCA and proposed putative roles of DNA methylation in iCCA.     

MiR-206 suppresses the deterioration of intrahepatic cholangiocarcinoma and promotes sensitivity to chemotherapy by inhibiting interactions with stromal CAFs

Background: Intrahepatic cholangiocarcinoma (iCCA) is a highly malignant subtype of cholangiocarcinoma (CCA) with poor prognosis. In iCCA, the interplay between the stroma and tumor cells results in resistance to adjuvant chemotherapy. Increasing evidence indicates that miR-206 participates in tumor progression, but its role in iCCA is still unclear. The aim of this study was to identify dysregulated miR-206 expression in iCCA and to further explore the underlying mechanism.Methods: MiR-206 expression was proven to be downregulated in iCCA tissues by qPCR, and its correlation with clinical characteristics and prognosis was investigated. iCCA-derived cancer-associated fibroblast cells (CAFs) and normal fibroblast cells (NFs) were isolated and identified. MiR-206 was knocked in or down in CAFs and CCA cells, respectively, to explore the role of miR-206, and coculture of these treated CCAs and CAFs was conducted to explore the effects of miR-206 on their mutual promoting effects. Exosomes carrying miR-206 and an orthotopic mouse model were used to determine the inhibitory effects of miR-206 on iCCA deterioration in vivo.Results: We confirmed that miR-206 is a suppressor of iCCA. Overexpressing miR-206 in CCA cells inhibited cell proliferation, migration and invasion. When cocultured with CCA cells, NFs downregulated miR-206 expression, and NFs were susceptible to transforming into CAFs. Moreover, CAFs promoted CCA cell malignant behaviors and gemcitabine resistance. Overexpressing miR-206 in CAFs or CCA cells inhibited this mutual promoting effect. Additionally, when delivered by exosomes, miR-206 suppressed tumor deterioration. And combined with gemcitabine, this treatment resulted in a longer survival time.Conclusion: Our study explained that the interaction between CCA cells and CAFs promoted iCCA deterioration. As a suppressive factor, miR-206 inhibited aggressive characteristics and gemcitabine resistance by interfering with this mutual promoting effect. This research elucidated the molecular mechanism underlying the unfavorable chemotherapeutic response of patients with iCCA, which provided a promising target for iCCA treatment.     

Impact of Gene Dosage on Gene Expression,Biological Processes and Survival in Cervical Cancer: A Genome-Wide Follow-Up Study

We investigated the role of tumor copy number (CN)–altered genome (CN-AG) in the carcinogenesis of cervical cancer (CC), especially its effect on gene expression, biological processes, and patient survival. Fifty-nine human papillomavirus 16 (HPV16)-positive CCs were investigated with microarrays–31 for mapping CN-AG and 55 for global gene expression, with 27 CCs in common. Five-year survival was investigated in 55 patients. Deletions and amplifications >2.5 Mb were defined as CN alterations. The %CN-AG varied from 0 to 32.2% (mean = 8.1±8.9). Tumors were classified as low (mean = 0.5±0.6, n = 11), medium (mean = 5.4±2.4, n = 10), or high (mean = 19.2±6.6, n = 10) CN. The highest %CN-AG was found in 3q, which contributed an average of 55% of all CN alterations. Genome-wide, only 5.3% of CN-altered genes were deregulated directly by gene dosage. In contrast, the rate in fully duplicated 3q was twice as high. Amplification of 3q explained 23.2% of deregulated genes in whole tumors (r² = 0.232, p = 0.006; analysis of variance), including genes located in 3q and other chromosomes. A total of 862 genes were deregulated exclusively in high-CN tumors, but only 22.9% were CN altered. This suggests that the remaining genes are not deregulated directly by gene dosage, but by mechanisms induced in trans by CN-altered genes. Anaphase-promoting complex/cyclosome (APC/C)-dependent proteasome proteolysis, glycolysis, and apoptosis were upregulated, whereas cell adhesion and angiogenesis were downregulated exclusively in high-CN tumors. The high %CN-AG and upregulated gene expression profile of APC/C-dependent proteasome proteolysis were associated with poor patient survival (p<0.05, log-rank test). Along with glycolysis, they were linearly associated with FIGO stage (r>0.38, p<0.01, Spearman test). Therefore, inhibition of APC/C-dependent proteasome proteolysis and glycolysis could be useful for CC treatment. However, whether they are indispensable for tumor growth remains to be demonstrated.     

Adjusting Breast Cancer Patient Prognosis with Non-HER2-Gene Patterns on Chromosome 17

Background

HER2 and TOP2A gene status are assessed for diagnostic and research purposes in breast cancer with fluorescence in situ hybridization (FISH). However, FISH probes do not target only the annotated gene, while chromosome 17 (chr17) is among the most unstable chromosomes in breast cancer. Here we asked whether the status of specifically targeted genes on chr17 might help in refining prognosis of early high-risk breast cancer patients.

Methods

Copy numbers (CN) for 14 genes on chr17, 4 of which were within and 10 outside the core HER2 amplicon (HER2- and non-HER2-genes, respectively) were assessed with qPCR in 485 paraffin-embedded tumor tissue samples from breast cancer patients treated with adjuvant chemotherapy in the frame of two randomized phase III trials.

Principal Findings

HER2-genes CN strongly correlated to each other (Spearman’s rho >0.6) and were concordant with FISH HER2 status (Kappa 0.6697 for ERBB2 CN). TOP2A CN were not concordant with TOP2A FISH status (Kappa 0.1154). CN hierarchical clustering revealed distinct patterns of gains, losses and complex alterations in HER2- and non-HER2-genes associated with IHC4 breast cancer subtypes. Upon multivariate analysis, non-HER2-gene gains independently predicted for shorter disease-free survival (DFS) and overall survival (OS) in patients with triple-negative cancer, as compared to luminal and HER2-positive tumors (interaction p = 0.007 for DFS and p = 0.011 for OS). Similarly, non-HER2-gene gains were associated with worse prognosis in patients who had undergone breast-conserving surgery as compared to modified radical mastectomy (p = 0.004 for both DFS and OS). Non-HER2-gene losses were unfavorable prognosticators in patients with 1–3 metastatic nodes, as compared to those with 4 or more nodes (p = 0.017 for DFS and p = 0.001 for OS).

Conclusions

TOP2A FISH and qPCR may not identify the same pathology on chr17q. Non-HER2 chr17 CN patterns may further predict outcome in breast cancer patients with known favorable and unfavorable prognosis.     

Molecular characterization of acquired tolerance of tumor cells to picropodophyllin (PPP)

Background

Picropodophyllin (PPP) is a promising novel anti-neoplastic agent that efficiently kills tumor cells in vitro and causes tumor regression and increased survival in vivo. We have previously reported that PPP treatment induced moderate tolerance in two out of 10 cell lines only, and here report the acquired genomic and expression alterations associated with PPP selection over 1.5 years of treatment.

Methodology/Principal Findings

Copy number alterations monitored using metaphase and array-based comparative genomic hybridization analyses revealed largely overlapping alterations in parental and maximally tolerant cells. Gain/ amplification of the MYC and PVT1 loci in 8q24.21 were verified on the chromosome level. Abnormalities observed in connection to PPP treatment included regular gains and losses, as well as homozygous losses in 10q24.1-q24.2 and 12p12.3-p13.2 in one of the lines and amplification at 5q11.2 in the other. Abnormalities observed in both tolerant derivatives include amplification/gain of 5q11.2, gain of 11q12.1-q14.3 and gain of 13q33.3-qter. Using Nexus software analysis we combined the array-CGH data with data from gene expression profilings and identified genes that were altered in both inputs. A subset of genes identified as downregulated (ALDH1A3, ANXA1, TLR4 and RAB5A) or upregulated (COX6A1, NFIX, ME1, MAPK and TAP2) were validated by siRNA in the tolerant or parental cells to alter sensitivity to PPP and confirmed to alter sensitivity to PPP in further cell lines.

Conclusions

Long-term PPP selection lead to altered gene expression in PPP tolerant cells with increase as well as decrease of genes involved in cell death such as PTEN and BCL2. In addition, acquired genomic copy number alterations were observed that were often reflected by altered mRNA expression levels for genes in the same regions.     

Oncogene Panel Sequencing Analysis Identifies Candidate Actionable Genes in Advanced Well-Differentiated Gastroenteropancreatic Neuroendocrine Tumors

Objective: To report the rate of candidate actionable somatic mutations in patients with locally advanced and metastatic gastro-enteropancreatic (GEP) neuroendocrine tumors (NET) and of other genetic alterations that may be associated with tumorigenesis.Methods: A phase II mutation targeted therapy trial was conducted in patients with advanced well-differentiated G1/G2 GEP-NET. Mutations found in the mTOR pathway-associated genes led to treatment with the mTOR inhibitor everolimus, and were defined as actionable. Tumor deoxyribonucleic acid (DNA) from GEP-NET were sequenced and compared with germline DNA, using the OncoVAR-NET assay, designed for hybrid capture sequencing of 500 tumor suppressor genes and oncogenes. Somatic variants were called and copy-number (CN) variant analysis was performed.Results: Thirty patients (14 small-intestine, 8 pancreatic, 3 unknown primary NET, and 5 of other primary sites) harbored 37 lesions (4 patients had DNA of multiple lesions sequenced). Only 2 patients with sporadic NET (n = 26) had an actionable mutation leading to treatment with everolimus. Driver somatic mutations were detected in 18 of 30 patients (21/37 lesions sequenced). In the remaining samples without a driver mutation, CN alterations were found in 11/16 tumors (10/12 patients), including CN loss of chromosome (Chr) 18 (P<.05), CN gain of Chr 5, and loss of Chr 13. CN losses in Chr 18 were more common in patients without driver mutations detected. Pronounced genetic heterogeneity was detected in patients with multiple lesions sequenced.Conclusion: Genome-wide DNA sequencing may identify candidate actionable genes and lead to the identification of novel target genes for advanced well-differentiated GEP-NET.Abbreviations: Chr = chromosome; CN = copy number; DNA = deoxyribonucleic acid; FDA = Food and Drug Administration; GEP = gastro-enteropancreatic; MEN-1 = multiple endocrine neoplasia syndrome type 1; mTOR = mammalian target of rapamycin; NET = neuroendocrine tumor; PFS = progression-free survival; PNET = pancreatic neuroendocrine tumors; SINET = small-intestine neuroendocrine tumor     

Adipose tissue gene expression analysis reveals changes in inflammatory, mitochondrial respiratory and lipid metabolic pathways in obese insulin-resistant subjects

Background

To elucidate gene expression associated with copy number changes, we performed a genome-wide copy number and expression microarray analysis of 25 pairs of gastric tissues.

Methods

We applied laser capture microdissection (LCM) to obtain samples for microarray experiments and profiled DNA copy number and gene expression using 244K CGH Microarray and Human Exon 1.0 ST Microarray.

Results

Obviously, gain at 8q was detected at the highest frequency (70%) and 20q at the second (63%). We also identified molecular genetic divergences for different TNM-stages or histological subtypes of gastric cancers. Interestingly, the C20orf11 amplification and gain at 20q13.33 almost separated moderately differentiated (MD) gastric cancers from poorly differentiated (PD) type. A set of 163 genes showing the correlations between gene copy number and expression was selected and the identified genes were able to discriminate matched adjacent noncancerous samples from gastric cancer samples in an unsupervised two-way hierarchical clustering. Quantitative RT-PCR analysis for 4 genes (C20orf11, XPO5, PUF60, and PLOD3) of the 163 genes validated the microarray results. Notably, some candidate genes (MCM4 and YWHAZ) and its adjacent genes such as PRKDC, UBE2V2, ANKRD46, ZNF706, and GRHL2, were concordantly deregulated by genomic aberrations.

Conclusions

Taken together, our results reveal diverse chromosomal region alterations for different TNM-stages or histological subtypes of gastric cancers, which is helpful in researching clinicopathological classification, and highlight several interesting genes as potential biomarkers for gastric cancer.     

Array-comparative genomic hybridization reveals loss of SOCS6 is associated with poor prognosis in primary lung squamous cell carcinoma

Background

Primary tumor recurrence commonly occurs after surgical resection of lung squamous cell carcinoma (SCC). Little is known about the genes driving SCC recurrence.

Methods

We used array comparative genomic hybridization (aCGH) to identify genes affected by copy number alterations that may be involved in SCC recurrence. Training and test sets of resected primary lung SCC were assembled. aCGH was used to determine genomic copy number in a training set of 62 primary lung SCCs (28 with recurrence and 34 with no evidence of recurrence) and the altered copy number of candidate genes was confirmed by quantitative PCR (qPCR). An independent test set of 72 primary lung SCCs (20 with recurrence and 52 with no evidence of recurrence) was used for biological validation. mRNA expression of candidate genes was studied using qRT-PCR. Candidate gene promoter methylation was evaluated using methylation microarrays and Sequenom EpiTYPER analysis.

Results

18q22.3 loss was identified by aCGH as being significantly associated with recurrence (p = 0.038). Seven genes within 18q22.3 had aCGH copy number loss associated with recurrence but only SOCS6 copy number was both technically replicated by qPCR and biologically validated in the test set. SOCS6 copy number loss correlated with reduced mRNA expression in the study samples and in the samples with copy number loss, there was a trend for increased methylation, albeit non-significant. Overall survival was significantly poorer in patients with SOCS6 loss compared to patients without SOCS6 loss in both the training (30 vs. 43 months, p = 0.023) and test set (27 vs. 43 months, p = 0.010).

Conclusion

Reduced copy number and mRNA expression of SOCS6 are associated with disease recurrence in primary lung SCC and may be useful prognostic biomarkers.     

A Meta-Analysis of Retinoblastoma Copy Numbers Refines the List of Possible Driver Genes Involved in Tumor Progression

BackgroundWhile RB1 loss initiates retinoblastoma development, additional somatic copy number alterations (SCNAs) can drive tumor progression. Although SCNAs have been identified with good concordance between studies at a cytoband resolution, accurate identification of single genes for all recurrent SCNAs is still challenging. This study presents a comprehensive meta-analysis of genome-wide SCNAs integrated with gene expression profiling data, narrowing down the list of plausible retinoblastoma driver genes.MethodsWe performed SCNA profiling of 45 primary retinoblastoma samples and eight retinoblastoma cell lines by high-resolution microarrays. We combined our data with genomic, clinical and histopathological data of ten published genome-wide SCNA studies, which strongly enhanced the power of our analyses (N = 310).ResultsComprehensive recurrence analysis of SCNAs in all studies integrated with gene expression data allowed us to reduce candidate gene lists for 1q, 2p, 6p, 7q and 13q to a limited gene set. Besides the well-established driver genes RB1 (13q-loss) and MYCN (2p-gain) we identified CRB1 and NEK7 (1q-gain), SOX4 (6p-gain) and NUP205 (7q-gain) as novel retinoblastoma driver candidates. Depending on the sample subset and algorithms used, alternative candidates were identified including MIR181 (1q-gain) and DEK (6p gain). Remarkably, our study showed that copy number gains rarely exceeded change of one copy, even in pure tumor samples with 100% homozygosity at the RB1 locus (N = 34), which is indicative for intra-tumor heterogeneity. In addition, profound between-tumor variability was observed that was associated with age at diagnosis and differentiation grades.InterpretationSince focal alterations at commonly altered chromosome regions were rare except for 2p24.3 (MYCN), further functional validation of the oncogenic potential of the described candidate genes is now required. For further investigations, our study provides a refined and revised set of candidate retinoblastoma driver genes.     

Increased MET Gene Copy Number but Not mRNA Level Predicts Postoperative Recurrence in Patients with Non–Small Cell Lung Cancer

The aim of the present study was to investigate the relationship of MET copy number (CN) and MET mRNA expression to other molecular alterations, clinicopathologic characteristics, and survival of patients with resected non–small cell lung cancer. One hundred fifty-one paired surgical samples of tumor and tumor-distant normal lung tissues were analyzed by comparative quantitative polymerase chain reaction (PCR) methods with commercially available assays and the CopyCaller software v. 1.0 for post-PCR data processing (downloadable from www.appliedbiosystems.com). MET copy gain (set as more than 3.0 copies per cell) was found in 18.5% of the samples and occurred more frequently in the adenocarcinomas (ADCs) with an increased epidermal growth factor receptor (EGFR) or human epidermal growth factor receptor 2 (HER2) CN (P = .001 and .030 for EGFR and HER2, respectively) and in the ADCs with EGFR activating mutations (P = .051) but did not correlate with KRAS dosage or mutational status. MET mRNA level was 1.76-fold higher [95% confidence interval (CI), 1.29-2.40] in the tumor compared to unaffected lung tissue and associated significantly with MET CN (beta coefficient, 1.51; 95% CI, 1.22-1.87; P < .001). In the multivariable analysis, patients diagnosed with ADC with increased MET CN had a significantly higher risk of disease recurrence (hazard ratio, 1.76; 95% CI, 1.20-2.57; P = .004). An increased MET CN in combination with histologic type appears to be a prognostic factor in patients with ADC after a curative surgery.     

A 1.1 Mb deletion in distal 13q deletion syndrome region with congenital heart defect and postaxial polydactyly: Additional support for a CHD locus at distal 13q34 region

13q deletion syndrome is a rare genetic disorder, especially for group 3 deletion (13q33–q34 deletion). Previously we described a patient with congenital heart defect and mental retardation and proposed that a distal 6 Mb region might contain the causative gene of congenital heart defect. Here we present a new patient with congenital heart defects (CHD), hand and foot anomalies and mild mental retardation. We identified a 1.1 Mb deletion at chromosome 13q34 with high resolution SNP-array BeadChips (HumanOmni1-Quad, Illumina, USA). This chromosome region contains ten annotated genes, including GRK1, TFDP1, RASA3 and GAS6. To our knowledge, this represents the smallest 13q34 deletion identified to date. Our study provides additional support that distal 13q34 deletion region might contain key gene(s) responsible for cardiac development.     

A transcriptomic computational analysis of mastic oil-treated Lewis lung carcinomas reveals molecular mechanisms targeting tumor cell growth and survival

Background

DNA copy number alterations are frequently observed in ovarian cancer, but it remains a challenge to identify the most relevant alterations and the specific causal genes in those regions.

Methods

We obtained high-resolution 500K SNP array data for 52 ovarian tumors and identified the most statistically significant minimal genomic regions with the most prevalent and highest-level copy number alterations (recurrent CNAs). Within a region of recurrent CNA, comparison of expression levels in tumors with a given CNA to tumors lacking that CNA and to whole normal ovary samples was used to select genes with CNA-specific expression patterns. A public expression array data set of laser capture micro-dissected (LCM) non-malignant fallopian tube epithelia and LCM ovarian serous adenocarcinoma was used to evaluate the effect of cell-type mixture biases.

Results

Fourteen recurrent deletions were detected on chromosomes 4, 6, 9, 12, 13, 15, 16, 17, 18, 22 and most prevalently on X and 8. Copy number and expression data suggest several apoptosis mediators as candidate drivers of the 8p deletions. Sixteen recurrent gains were identified on chromosomes 1, 2, 3, 5, 8, 10, 12, 15, 17, 19, and 20, with the most prevalent gains localized to 8q and 3q. Within the 8q amplicon, PVT1, but not MYC, was strongly over-expressed relative to tumors lacking this CNA and showed over-expression relative to normal ovary. Likewise, the cell polarity regulators PRKCI and ECT2 were identified as putative drivers of two distinct amplicons on 3q. Co-occurrence analyses suggested potential synergistic or antagonistic relationships between recurrent CNAs. Genes within regions of recurrent CNA showed an enrichment of Cancer Census genes, particularly when filtered for CNA-specific expression.

Conclusion

These analyses provide detailed views of ovarian cancer genomic changes and highlight the benefits of using multiple reference sample types for the evaluation of CNA-specific expression changes.     

Chromosome instability region analysis and identification of the driver genes of the epithelial ovarian cancer cell lines A2780 and SKOV3

Epithelial ovarian cancer (EOC) is one of the most prevalent gynaecological cancers worldwide. The molecular mechanisms of serous ovarian cancer (SOC) remain unclear and not well understood. SOC cases are primarily diagnosed at the late stage, resulting in a poor prognosis. Advances in molecular biology techniques allow us to obtain a better understanding of precise molecular mechanisms and to identify the chromosome instability region and key driver genes in the carcinogenesis and progression of SOC. Whole-exome sequencing was performed on the normal ovarian cell line IOSE80 and the EOC cell lines SKOV3 and A2780. The single-nucleotide variation burden, distribution, frequency and signature followed the known ovarian mutation profiles, without chromosomal bias. Recurrently mutated ovarian cancer driver genes, including LRP1B, KMT2A, ARID1A, KMT2C and ATRX were also found in two cell lines. The genome distribution of copy number alterations was found by copy number variation (CNV) analysis, including amplification of 17q12 and 4p16.1 and deletion of 10q23.33. The CNVs of MED1, GRB7 and MIEN1 located at 17q12 were found to be correlated with the overall survival of SOC patients (MED1: p = 0.028, GRB7: p = 0.0048, MIEN1: p = 0.0051), and the expression of the three driver genes in the ovarian cell line IOSE80 and EOC cell lines SKOV3 and A2780 was confirmed by western blot and cell immunohistochemistry.     

Focal DNA Copy Number Changes in Neuroblastoma Target MYCN Regulated Genes

Neuroblastoma is an embryonic tumor arising from immature sympathetic nervous system cells. Recurrent genomic alterations include MYCN and ALK amplification as well as recurrent patterns of gains and losses of whole or large partial chromosome segments. A recent whole genome sequencing effort yielded no frequently recurring mutations in genes other than those affecting ALK. However, the study further stresses the importance of DNA copy number alterations in this disease, in particular for genes implicated in neuritogenesis. Here we provide additional evidence for the importance of focal DNA copy number gains and losses, which are predominantly observed in MYCN amplified tumors. A focal 5 kb gain encompassing the MYCN regulated miR-17∼92 cluster as sole gene was detected in a neuroblastoma cell line and further analyses of the array CGH data set demonstrated enrichment for other MYCN target genes in focal gains and amplifications. Next we applied an integrated genomics analysis to prioritize MYCN down regulated genes mediated by MYCN driven miRNAs within regions of focal heterozygous or homozygous deletion. We identified RGS5, a negative regulator of G-protein signaling implicated in vascular normalization, invasion and metastasis, targeted by a focal homozygous deletion, as a new MYCN target gene, down regulated through MYCN activated miRNAs. In addition, we expand the miR-17∼92 regulatory network controlling TGFß signaling in neuroblastoma with the ring finger protein 11 encoding gene RNF11, which was previously shown to be targeted by the miR-17∼92 member miR-19b. Taken together, our data indicate that focal DNA copy number imbalances in neuroblastoma (1) target genes that are implicated in MYCN signaling, possibly selected to reinforce MYCN oncogene addiction and (2) serve as a resource for identifying new molecular targets for treatment.