DNA Copy Number Amplification Research Articles

Comprehensive single cell transcriptomics analysis of murine osteosarcoma uncovers Skp2 function in metastasis, genomic instability and immune activation and reveals additional target pathways.

Osteosarcoma (OS) is the most common primary pediatric bone malignancy. One promising new therapeutic target is SKP2, encoding a substrate recognition factor of the SCF E3 ubiquitin ligase responsible for ubiquitination and proteasome degradation of substrate p27, thus driving cellular proliferation. We have shown previously that knockout of Skp2 in an immunocompetent transgenic mouse model of OS improved survival, drove apoptosis, and induced tumor inflammation. Here, we applied single-cell RNA-sequencing (scRNA-seq) to study primary OS tumors derived from Osx-Cre driven conditional knockout of Rb1 and Trp53. We showed that murine OS models recapitulate the tumor heterogeneity and microenvironment complexity observed in patient tumors. We further compared this model with OS models with functional disruption of Skp2: one with Skp2 knockout and the other with the Skp2-p27 interaction disrupted (resulting in p27 overexpression). We found reduction of T cell exhaustion and upregulation of interferon activation, along with evidence of replicative and endoplasmic reticulum-related stress in the Skp2 disruption models, and showed that interferon induction was correlated with improved survival in OS patients. Additionally, our scRNA-seq analysis uncovered decreased activities of metastasis-related gene signatures in the Skp2-disrupted OS, which we validated by observation of a strong reduction in lung metastasis in the Skp2 knockout mice. Finally, we report several potential mechanisms of escape from targeting Skp2 in OS, including upregulation of Myc targets, DNA copy number amplification and overexpression of alternative E3 ligase genes, and potential alternative lineage activation. These mechanistic insights into OS tumor biology and Skp2 function suggest novel targets for new, synergistic therapies, while the data and our comprehensive analysis may serve as a public resource for further big data-driven OS research.

An integrative pan-cancer bioinformatics analysis of MSRB1 and its association with tumor immune microenvironment, prognosis, and immunotherapy

Methionine sulfoxide reductase B1 (MSRB1) is involved in the development and immune regulation of multiple tumors. However, the role of MSRB1 in the tumor microenvironment and its potential as a therapeutic target remain largely unknown. In this study, MSRB1 expression patterns were evaluated using pan-cancer RNA sequencing data from multiple cell lines, tissues, and single cells. The pan-cancer prognostic role of MSRB1 was assessed and the association between MSRB1 expression and certain cancer characteristics was analyzed. We showed that MSRB1 expression levels were increased in several types of cancer (P < 0.05) and in certain cell types (macrophages, dendritic cells, and malignant tumor cells). The upregulation of MSRB1 expression was due to DNA copy number amplification. Furthermore, MSRB1 was significantly associated with the activation of immune pathways (P < 0.05, NES > 0), immune cell infiltration, and expression of immune checkpoint molecules. In addition, high expression of MSRB1 was found in a series of in vivo and in vitro immunotherapy response models (P < 0.05), and showed resistance to most targeted drugs. Our results indicated that MSRB1 may regulate the tumor immune microenvironment through an immunoresponse and potentially influence cancer development. This could make it a promising predictive biomarker and therapeutic target for precise tumor immunotherapy.

Inhalable mixture of polycyclic aromatic hydrocarbons and metals, DNA oxidative stress and nasal ribosomal DNA copy number amplification: Direct and indirect effect analyses among population

The ribosomal DNA (rDNA) plays a crucial role in maintaining genome stability. So far, alterations of rDNA from airborne pollutants exposure remain unclear. Nasal epithelial cells are the earliest respiratory barrier, which has an accessible surrogate for evaluating respiratory impairment. We developed a mixture-centered biomarkers study integrated epidemiological and biological evidence among 768 subjects, a mixture of polycyclic aromatic hydrocarbons (PAHs) and metals. We identified the mixed exposure of PAHs and metals by environmental and biological monitoring, selected urinary 8-hydroxy-2′-deoxyguanosine as DNA oxidative stress marker, and measured their rDNA copy number (rDNA CN) in nasal epithelial cells. We performed linear regression, adaptive elastic net regression, BKMR, and mediation analyses to assess the direct and indirect effects. We found a 10% elevation in urinary 1-hydroxypyrene was correlated with a separate 0.31% and 0.82% amplification of nasal 5S and 45S rDNA CN, respectively (all P < 0.05). A 10% increment of urine nickel was associated with a separate 0.37% and 1.18% elevation of nasal 5S and 45S rDNA CN, respectively (all P < 0.05). BKMR results also confirmed our findings of PAHs and nickel. Our findings suggested that DNA oxidative stress might trigger rDNA instability induced by inhaled PAHs and metals.

Development of an osteosarcoma model with MYCN amplification and TP53 mutation in hiPS cell-derived neural crest cells.

Mesenchymal stem cell- or osteoblast-derived osteosarcoma is the most common malignant bone tumor. Its highly metastatic malignant phenotypes, which are often associated with a poor prognosis, have been correlated with the modulation of TP53- and cell-cycle-related pathways. MYC, which regulates the transcription of cell-cycle modulating genes, is used as a representative prognostic marker for osteosarcoma. Another member of the MYC oncoprotein family, MYCN, is highly expressed in a subset of osteosarcoma, however its roles in osteosarcoma have not been fully elucidated. Here, we attempted to create an in vitro tumorigenesis model using hiPSC-derived neural crest cells, which are precursors of mesenchymal stem cells, by overexpressing MYCN on a heterozygous TP53 hotspot mutation (c.733G>A; p.G245S) background. MYCN-expressing TP53 mutated transformed clones were isolated by soft agar colony formation, and administered subcutaneously into the periadrenal adipose tissue of immunodeficient mice, resulting in the development of chondroblastic osteosarcoma. MYCN suppression decreased the proliferation of MYCN-induced osteosarcoma cells, suggesting MYCN as a potential target for a subset of osteosarcoma treatment. Further, comprehensive analysis of gene expression and exome sequencing of MYCN-induced clones indicated osteosarcoma-specific molecular features, such as the activation of TGF-β signaling and DNA copy number amplification of GLI1. The model of MYCN-expressing chondroblastic osteosarcoma was developed from hiPSC-derived neural crest cells, providing a useful tool for the development of new tumor models using hiPSC-derived progenitor cells with gene modifications and in vitro transformation.

GAPDH: A common housekeeping gene with an oncogenic role in pan-cancer

Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) is one of the most prominent housekeeping proteins and is widely used as an internal control in some semi-quantitative assays. In addition to glycolysis, GAPDH is involved in several cancer-related biological processes and has been reported to be commonly dysregulated in multiple cancer types. Therefore, its role in the physiological process of cancer needs to be urgently elucidated. Pan-cancer analysis indicated that GAPDH is ubiquitously highly expressed in most cancer types, and that patients with a high GAPDH expression of in tumor tissues have a poor prognosis. The concordance of GAPDH expression in tumors with the infiltration of immune cells and immune checkpoints implies a certain association between GAPDH and the tumor microenvironment as well as tumor development. Gene Set Enrichment Analysis revealed that GAPDH may contribute to multiple important cancer-related pathways and biological processes. Multi-omics analysis and in vitro cell experiments revealed that GAPDH overexpression is regulated by DNA copy number amplification and promoter methylation modification. Importantly, a transcription factor, forkhead box M1 (FOXM1), which is capable of regulating GAPDH expression, was also identified and was confirmed to be an oncogene and ubiquitously highly expressed in multiple cancer types. Semi-quantitative chromatin immunoprecipitation, quantitative PCR, and dual-luciferase assays showed that FOXM1 mainly binds to the promoter region of GAPDH in two cancer cell lines. The present findings revealed the implication of GAPDH in tumor development, thus bringing attention to this important molecule and casting doubts on its role as an internal reference gene in cancer studies.

Extra Spindle Pole Bodies-Like 1 Serves as a Prognostic Biomarker and Promotes Lung Adenocarcinoma Metastasis

Extra spindle pole bodies-like 1 (ESPL1), a cysteine endopeptidase, plays a vital role in chromosome inheritance. However, the association of ESPL1 with prognosis and immune infiltration in lung adenocarcinoma (LUAD) has not yet been explored. Here, we analyzed the expression level, prognostic values, diagnostic value, and immune infiltration level in LUAD using various databases. Immunohistochemistry (IHC) and quantitative real-time PCR (qRT-PCR) assays were used to detect the expression of ESPL1 in LUAD tissues and cell lines. In this study, we found that ESPL1 was upregulated in LUAD and a higher expression of ESPL1 was correlated with unfavorable prognosis in LUAD. Meanwhile, Cox hazard regression analysis results suggested that ESPL1 may be an independent prognostic factor for LUAD. Moreover, we demonstrated that ESPL1 expression was significantly correlated with immune infiltration of Th2 and dendritic cells in LUAD. We also confirmed that DNA copy number amplification and DNA hypo-methylation were positively correlated with ESPL1 expression in LUAD. Additionally, DNA copy number amplification was significantly associated with adverse clinical outcomes in LUAD. Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene set enrichment analysis (GSEA) confirmed that ESPL1 was mainly involved in the DNA replication and glycolysis signaling pathway. Finally, we revealed that ESPL1 was highly expressed in LUAD tissues and cell lines. Knockdown of ESPL1 significantly inhibited cell migration and the invasion abilities of LUAD. Our study comprehensively confirmed that ESPL1 expression may serve as a novel prognostic biomarker for both the clinical outcome and immune cell infiltration in LUAD.

LINC00467 Is Upregulated by DNA Copy Number Amplification and Hypomethylation and Shows ceRNA Potential in Lung Adenocarcinoma.

Lung adenocarcinoma (LUAD) is the most common histological lung cancer, and it is the leading cause of cancer-related deaths worldwide. Long noncoding RNAs (lncRNAs) have been implicated in tumorigenesis. LINC00467 is a novel lncRNA that is abnormally expressed in several cancer types including LUAD. However, its function and regulatory mechanism in LUAD progression remain unclear. In this study, based on The Cancer Genome Atlas data mining, we demonstrated that DNA copy number amplification and hypomethylation was positively correlated with LINC00467 expression in LUAD. In addition, DNA copy number amplification was significantly associated with distant metastasis, immune infiltration and poor survival. Microarray analysis demonstrated that LINC00467 knockdown in the LUAD A549 cell line led to a distinct microRNA expression profile that impacted various target genes involved in multiple biological processes. This finding suggests that LINC00467 may regulate LUAD progression by functioning as a competing endogenous RNA (ceRNA). Finally, we constructed a ceRNA network that included two microRNAs (hsa-miR-1225-5p, hsa-miR-575) and five mRNAs (BARX2, BCL9, KCNK1, KIAA1324, TMEM182) specific to LINC00467 in LUAD. Subsequent Kaplan-Meier survival analysis in both The Cancer Genome Atlas and Gene Expression Omnibus databases revealed that two genes, BARX2 and BCL9, were potential prognostic biomarkers for LUAD patients. In conclusion, our data provide possible mechanisms underlying the abnormal upregulation of LINC00467 as well as a comprehensive view of the LINC00467-mediated ceRNA network in LUAD, thereby highlighting its potential role in diagnosis and therapy.

Overexpression of Replication-Dependent Histone Signifies a Subset of Dedifferentiated Liposarcoma with Increased Aggressiveness.

Simple SummaryAlthough the genetic event that initiates liposarcoma is relatively well studied, understanding how the tumor progresses and displays differential severities will shed new insights into the diagnosis and prognosis determination of patients with liposarcoma. We analyzed the genome and transcriptome of liposarcoma samples and found that the overexpression of the histone gene cluster is an important determinant for dedifferentiated liposarcoma in predicting genome instability and cell proliferation. The overexpression of histone genes was caused by increased HMGA2 with amplification of chromosome 12q13-14. This dedifferentiated liposarcoma-specific genetic signature may serve as a biomarker for prognosis prediction.Liposarcoma (LPS) is an adult soft tissue malignancy that arises from fat tissue, where well-differentiated (WD) and dedifferentiated (DD) forms are the most common. DDLPS represents the progression of WDLPS into a more aggressive high-grade and metastatic form. Although a few DNA copy-number amplifications are known to be specifically found in WD- or DDLPS, systematic genetic differences that signify subtype determination between WDLPS and DDLPS remain unclear. Here, we profiled the genome and transcriptome of 38 LPS tumors to uncover the genetic signatures of subtype differences. Replication-dependent histone (RD-HIST) mRNAs were highly elevated and their regulation was disrupted in a subset of DDLPS, increasing cellular histone molecule levels, as measured using RNA-seq (the averaged fold change of 53 RD-HIST genes between the DD and WD samples was 10.9) and immunohistochemistry. The change was not observed in normal tissues. Integrated whole-exome sequencing, RNA-seq, and methylation analyses revealed that the overexpressed HMGA2 (the fold change between DD and WD samples was 7.3) was responsible for the increased RD-HIST level, leading to aberrant cell proliferation. Therefore, HMGA2-mediated elevation of RD-HISTs were crucial events in determining the aggressiveness of DDLPS, which may serve as a biomarker for prognosis prediction for liposarcoma patients.

Mutation and Copy Number Alterations Analysis of KIF23 in Glioma.

In glioma, kinesin family member 23 (KIF23) is up-regulated and plays a vital role in oncogenesis. However, the mechanism underlying KIF23 overexpression in malignant glioma remains to be elucidated. This study aims to find potential causes of KIF23 high expression at genome level. To clarify this issue, we obtained point mutation and copy number alterations (CNAs) of KIF23 in 319 gliomas using whole-exome sequencing. Only two glioma samples with missense mutations in KIF23 coding region were identified, while 7 patients were detected with amplification of KIF23. Additional analysis showed that KIF23 amplification was significantly associated with higher expression of KIF23. Gene ontology analysis indicated that higher copy number of KIF23 was associated TNF-α signaling pathway and mitotic cell circle checkpoint, which probably caused by subsequent upregulated expression of KIF23. Moreover, pan-cancer analysis showed that gaining of copy number was significantly associated with higher expression of KIF23, consolidating our findings in glioma. Thus, it was deduced that elevated KIF23 expression in glioma tended to be caused by DNA copy number amplification, instead of mutation.

Comprehensive analysis on the expression levels and prognostic values of LOX family genes in kidney renal clear cell carcinoma.

BackgroundsKidney renal clear cell carcinoma (KIRC) is a major pathological type of renal cell carcinoma (RCC), and the prognosis of advanced KIRC patients is often unsatisfactory. Some lysine oxidase (LOX) family genes have been proven to be upregulated in some malignancies and play pivotal roles in the carcinogenesis. However, their roles in KIRC remain unclear.Materials and MethodsHere, we used some online databases (eg, ONCOMINE, GEPIA, UALCAN, c‐BioPortal, Human Protein Altas) to comprehensively explored the expression levels and the prognostic values of LOX family genes in KIRC using bioinformatic methods.ResultsThe results revealed that lysyl oxidase (LOX) and lysyl oxidase‐like 2 (LOXL2) were significantly overexpressed in KIRC at the level of mRNA expression, protein expression, and RCC cell lines. Further analysis demonstrated that higher mRNA expression of LOX and LOXL2 were significantly correlated with poor survival, tumor grade, individual cancer stages, and nodal metastasis status. DNA copy number amplifications and mRNA upregulation, DNA deep deletion, and mRNA upregulation were the main genetic mutations of LOX and LOXL2, respectively. Prognostic analysis showed that the altered group had significantly poorer overall survival (OS) compared to the unaltered group (p = .0387). Co‐expression analysis showed CP, PLOD2, and COL5A1 were significantly correlated with LOX, and COL1A2 was positively correlated with LOXL2. Further analysis confirmed that these co‐expressed genes were significantly upregulated and predicted unfavorable prognosis in KIRC.ConclusionMulti‐level analysis demonstrated that LOX and LOXL2 were significantly upregulated and predicted poor survival in KIRC, which may apply as promising biomarkers for diagnosis and therapy of KIRC in the future.

Comprehensive analysis of the value of RAB family genes in prognosis of breast invasive carcinoma.

Purpose: Several RAB family genes have been studied extensively and proven to play pivotal roles in the occurrence and development of certain cancers. Here, we explored commonly expressed RAB family genes in humans and their prognostic significance using bioinformatics, and then identified potential biomarkers of breast invasive carcinoma (BRCA). Materials and methods: The prognostic values (overall survival) of RAB family genes in BRCA were obtained using Gene Expression Profiling Interactive Analysis (GEPIA). The expression patterns of RAB family genes and their relationships with clinicopathological parameters in BRCA were measured using the ONCOMINE and UALCAN databases, respectively. Genetic mutations and survival analysis were investigated using the cBio Cancer Genomics Portal (c-BioPortal). Interacting genes of potential biomarkers were identified using STRING, and functional enrichment analyses were performed using FunRich v3.1.3. Results: In total, 64 RAB genes were identified and analyzed in our study. Results showed that RAB1B, RAB2A, and RAB18 were up-regulated and significantly associated with poor overall survival in BRCA. Furthermore, their higher expression was positively correlated with clinicopathological parameters (e.g. cancer stage and nodal metastasis status). DNA copy number amplifications and mRNA up-regulation were the main genetic mutations, and the altered group showed significantly poorer overall survival compared with the unaltered group. Functional enrichment analysis of RAB1B, RAB2A, and RAB18 indicated they were closely involved in GTPase activity. Conclusions: RAB1B, RAB2A, and RAB18 were up-regulated and significantly correlated with poor prognosis in BRCA. Thus, they could be applied as novel biomarkers of BRCA in future studies.

Abstract P4-10-26: Protein-protein interaction of HER2 predicts the prognosis of hormone receptor-positive HER2-negative breast cancer

Abstract Introduction: Patients with Hormone receptor (HR)-positive / HER2-negtive breast cancer have not received targeted therapy for Her2 based on criteria of the protein expression or DNA copy number amplification. However, the clinical meaning of the level of Her2 signaling activity has not been fully elucidated yet. We recently developed a novel technique of single-molecule protein-protein interaction (PPI) profiling to investigate the importance of HER2 signaling activity in this major subgroup of breast cancer. Methods: This study is a study with the biospecimen and registered medical record from the prospective collected cohort of Samsung Medical Center, and all of these biospecimens were provided by Samsung Medical Center BioBank. First, We selected 138 of 8463 patients who were diagnosed with invasive breast cancer by curative surgery between January 2008 and December 2013. Specimens were selected based on modified Allred Immunohistochemistry (IHC) scores; 6 for hormone receptors (ER/PR+) and 1 for Her2. Then, we applied the single-molecule co-immunoprecipitation, a novel technology determining protein-protein interactions (PPIs), to each selected biospecimen. PPI and expression level for HER2 protein of each tumor extract were analyzed with HER2-PPI-Dx kit. All the measurement were performed using a single-molecule fluorescence imaging equipment. Based on HER2 PPI and expression level analysis, we developed an algorithm based on the 72 patients’ samples which can predict the prognosis of these patients. This algorithm was validated with the samples from the other 66 patients. Results: In total, 138 samples from each patients (72 in training set and 66 in validation set) were analyzed. The median follow-up time was 68.0 (range: 7-123) months. Single-molecule measurement for HER2 level showed that clinical specimens showed significant differences in HER2-PPI level even though they had the same IHC score (65-fold from 117 patients with ‘HER2 IHC score: 0’ and 33-fold from 21 patients with ‘HER2-IHC score: 1’, respectively). For the training set, by combining HER2 PPI and HER2 protein expression data, we could select 12 patients (17%) as a high level HER2-PPI group. This high level HER2-PPI group showed the significantly poorer disease free survival (DFS) rate compared to that of the low level group (42% for high-level vs. 83% for low-level, p=0.001). For the validation set with the same cutoff, 17 patients (26%) with identified as high level HER2-PPI group. We still found a significantly poorer DFS rate in this group compared to that of the low level group (82% vs. 88%, p=0.033). For total 138 patients set, we found a significantly poorer DFS rate in this group compared to that of the low level group (p=0.014). Conclusion: Our results showed that HER2 PPI may have an important clinical meaning to identify a high risk of recurrence in the patients with strong HR-positive and HER2-negative breast cancers defined by conventional method. For those patients with high level HER2-PPI, HER2-targeted therapies may be beneficial despite of their HER2-negativity. This will be proved for the future. Citation Format: Hee Jun Choi, Hong Won Lee, Sangmin Kim, Tae-Young Yoon, Hyunwoo Kim, Seul Lee, Seok Jin Nam, Seok Won Kim, Sooyoun Cho, Yeon Hee Park, Jeong Eon Lee. Protein-protein interaction of HER2 predicts the prognosis of hormone receptor-positive HER2-negative breast cancer [abstract]. In: Proceedings of the 2019 San Antonio Breast Cancer Symposium; 2019 Dec 10-14; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(4 Suppl):Abstract nr P4-10-26.

GTF2IRD1 on chromosome 7 is a novel oncogene regulating the tumor-suppressor gene TGFβR2 in colorectal cancer.

Chromosome 7q (Ch.7q) is clonally amplified in colorectal cancer (CRC). We aimed to identify oncogenes on Ch.7q that are overexpressed through DNA copy number amplification and determine the biological and clinical significance of these oncogenes in CRC. We identified general transcription factor 2I repeat domain‐containing protein 1 (GTF2IRD1) as a potential oncogene using a CRC dataset from The Cancer Genome Atlas with a bioinformatics approach. We measured the expression of GTF2IRD1 in 98 patients with CRC using immunohistochemistry and RT‐quantitative PCR (RT‐qPCR). The biological effects of GTF2IRD1 expression were explored by gene set enrichment analysis (GSEA). Next, we undertook in vitro cell proliferation and cell cycle assays using siGTF2IRD1‐transfected CRC cells. We further investigated the oncogenic mechanisms through which GTF2IRD1 promoted CRC progression. Finally, we assessed the clinical significance of GTF2IRD1 expression by RT‐qPCR. GTF2IRD1 was overexpressed in tumor cells and liver metastatic lesions. The GSEA revealed a positive correlation between GTF2IRD1 expression and cell cycle progression‐related genes. GTF2IRD1 knockdown inhibited cell proliferation and induced cell cycle arrest in Smad4‐mutated CRC. GTF2IRD1 downregulated the expression of the gene encoding transforming growth factor β receptor 2 (TGFβR2), a tumor‐suppressor gene in Smad4‐mutated CRC. On multivariate analysis, high GTF2IRD1 expression was an independent poor prognostic factor. Clinicopathological analysis showed that GTF2IRD1 expression was positively correlated with liver metastasis. In conclusion, GTF2IRD1 promoted CRC progression by downregulating TGFβR2 and could be a prognostic biomarker on Ch.7q in CRC. GTF2IRD1 could also be a novel oncogene in CRC.

DNA copy number gain-mediated lncRNA LINC01061 upregulation predicts poor prognosis and promotes papillary thyroid cancer progression

Several DNA copy number amplifications (CNAs) have been reported in papillary thyroid cancer (PTC). However, the functional role of CNAs in PTC remains very unclear. And whether there is a correlation between long noncoding RNA (lncRNA) and CNA requires to be explored. Here, we identified a novel lncRNA LINC01061. The genomic copy number of LINC01061 is amplified, which leads to its elevated expression level in PTC tissues. Moreover, increased level of LINC01061 was correlated with aggressive clinicopathological characteristics. Functional study indicated that LINC01061 silence significantly inhibited the proliferation, cell-cycle and invasion of PTC cells in vitro and tumor growth in vivo. Mechanistically, we showed that LINC01061 interacted with miR-4316 to promote E2F6 expression. The expression of miR-4316 was downregulated in PTC tissues while that of E2F6 was upregulated. Through rescue assay, we demonstrated that LINC01061 promoted PTC cell proliferation, cell-cycle progression and invasion by regulating miR-4316/E2F6 signaling pathway. In conclusion, our research indicated that LINC01061 might be a target for PTC therapy.

Ribosomal DNA copy number amplification and loss in human cancers is linked to tumor genetic context, nucleolus activity, and proliferation.

Ribosomal RNAs (rRNAs) are transcribed from two multicopy DNA arrays: the 5S ribosomal DNA (rDNA) array residing in a single human autosome and the 45S rDNA array residing in five human autosomes. The arrays are among the most variable segments of the genome, exhibit concerted copy number variation (cCNV), encode essential components of the ribosome, and modulate global gene expression. Here we combined whole genome data from >700 tumors and paired normal tissues to provide a portrait of rDNA variation in human tissues and cancers of diverse mutational signatures, including stomach and lung adenocarcinomas, ovarian cancers, and others of the TCGA panel. We show that cancers undergo coupled 5S rDNA array expansion and 45S rDNA loss that is accompanied by increased estimates of proliferation rate and nucleolar activity. These somatic changes in rDNA CN occur in a background of over 10-fold naturally occurring rDNA CN variation across individuals and cCNV of 5S-45S arrays in some but not all tissues. Analysis of genetic context revealed associations between cancer rDNA CN amplification or loss and the presence of specific somatic alterations, including somatic SNPs and copy number gain/losses in protein coding genes across the cancer genome. For instance, somatic inactivation of the tumor suppressor gene TP53 emerged with a strong association with coupled 5S expansion / 45S loss in several cancers. Our results uncover frequent and contrasting changes in the 5S and 45S rDNA along rapidly proliferating cell lineages with high nucleolar activity. We suggest that 5S rDNA amplification facilitates increased proliferation, nucleolar activity, and ribosomal synthesis in cancer, whereas 45S rDNA loss emerges as a byproduct of transcription-replication conflict in rapidly replicating tumor cells. The observations raise the prospects of using the rDNA arrays as re-emerging targets for the design of novel strategies in cancer therapy.

Multi-omics of 34 colorectal cancer cell lines - a resource for biomedical studies

BackgroundColorectal cancer (CRC) cell lines are widely used pre-clinical model systems. Comprehensive insights into their molecular characteristics may improve model selection for biomedical studies.MethodsWe have performed DNA, RNA and protein profiling of 34 cell lines, including (i) targeted deep sequencing (n = 612 genes) to detect single nucleotide variants and insertions/deletions; (ii) high resolution DNA copy number profiling; (iii) gene expression profiling at exon resolution; (iv) small RNA expression profiling by deep sequencing; and (v) protein expression analysis (n = 297 proteins) by reverse phase protein microarrays.ResultsThe cell lines were stratified according to the key molecular subtypes of CRC and data were integrated at two or more levels by computational analyses. We confirm that the frequencies and patterns of DNA aberrations are associated with genomic instability phenotypes and that the cell lines recapitulate the genomic profiles of primary carcinomas. Intrinsic expression subgroups are distinct from genomic subtypes, but consistent at the gene-, microRNA- and protein-level and dominated by two distinct clusters; colon-like cell lines characterized by expression of gastro-intestinal differentiation markers and undifferentiated cell lines showing upregulation of epithelial-mesenchymal transition and TGFβ signatures. This sample split was concordant with the gene expression-based consensus molecular subtypes of primary tumors. Approximately ¼ of the genes had consistent regulation at the DNA copy number and gene expression level, while expression of gene-protein pairs in general was strongly correlated. Consistent high-level DNA copy number amplification and outlier gene- and protein- expression was found for several oncogenes in individual cell lines, including MYC and ERBB2.ConclusionsThis study expands the view of CRC cell lines as accurate molecular models of primary carcinomas, and we present integrated multi-level molecular data of 34 widely used cell lines in easily accessible formats, providing a resource for preclinical studies in CRC.

Phosphoserine Phosphatase Is a Novel Prognostic Biomarker on Chromosome 7 in Colorectal Cancer.

Amplification of chromosome 7p (Ch.7p) is common in colorectal cancer (CRC). The aim of this study was to identify potential driver genes on Ch.7p that are overexpressed due to DNA copy number amplification and determine their clinical significance in CRC. We identified phosphoserine phosphatase (PSPH) as a potential driver gene using a CRC dataset from The Cancer Genome Atlas (TCGA) using a bioinformatics approach. The expression of PSPH in 124 primary CRCs was examined by quantitative reverse transcription polymerase chain reaction (PCR) and immunohistochemistry. The biological effect of PSPH expression was explored by Gene Set Enrichment Analysis (GSEA) using the TCGA dataset. PSPH was overexpressed in tumor tissues and PSPH positively correlated with depth of invasion and distant metastasis. On multivariate analysis, high PSPH expression was an independent poor prognostic factor. These results were supported by GSEA. PSPH could be a novel prognostic biomarker with malignant potential on Ch.7p in CRC.

Detection of wild-type EGFR amplification and EGFRvIII mutation in CSF-derived extracellular vesicles of glioblastoma patients.

RNAs within extracellular vesicles (EVs) have potential as diagnostic biomarkers for patients with cancer and are identified in a variety of biofluids. Glioblastomas (GBMs) release EVs containing RNA into cerebrospinal fluid (CSF). Here we describe a multi-institutional study of RNA extracted from CSF-derived EVs of GBM patients to detect the presence of tumor-associated amplifications and mutations in epidermal growth factor receptor (EGFR). CSF and matching tumor tissue were obtained from patients undergoing resection of GBMs. We determined wild-type (wt)EGFR DNA copy number amplification, as well as wtEGFR and EGFR variant (v)III RNA expression in tumor samples. We also characterized wtEGFR and EGFRvIII RNA expression in CSF-derived EVs. EGFRvIII-positive tumors had significantly greater wtEGFR DNA amplification (P = 0.02) and RNA expression (P = 0.03), and EGFRvIII-positive CSF-derived EVs had significantly more wtEGFR RNA expression (P = 0.004). EGFRvIII was detected in CSF-derived EVs for 14 of the 23 EGFRvIII tissue-positive GBM patients. Conversely, only one of the 48 EGFRvIII tissue-negative patients had the EGFRvIII mutation detected in their CSF-derived EVs. These results yield a sensitivity of 61% and a specificity of 98% for the utility of CSF-derived EVs to detect an EGFRvIII-positive GBM. Our results demonstrate CSF-derived EVs contain RNA signatures reflective of the underlying molecular genetic status of GBMs in terms of wtEGFR expression and EGFRvIII status. The high specificity of the CSF-derived EV diagnostic test gives us an accurate determination of positive EGFRvIII tumor status and is essentially a less invasive "liquid biopsy" that might direct mutation-specific therapies for GBMs.

MA17.10 YES1 Kinase is a New Therapeutic Target in Non-small Cell Lung Cancer

Next-generation sequencing techniques have allowed the discovery of driver mutations in non-small cell lung cancer (NSCLC) that can be translated into advances in cancer diagnosis and treatment. However, specific oncogenic alterations are still unknown in a high proportion of NSCLC patients that therefore cannot benefit from targeted therapies. The challenge is to identify new genetic alterations that allow the use of molecular-targeted therapies. In previous studies from our group (Aramburu et al. BMC Genomics 2015), the analysis of tumor molecular profiles from patients with NSCLC allowed us to identify the DNA copy number amplification of YES1 kinase (v-YES-1 Yamaguchi sarcoma viral oncogene homolog 1) as a prognostic marker in lung cancer. YES1 kinase is member of the Src family of non-receptor protein tyrosine kinases that are involved in the regulation of cell growth, apoptosis, cell-cell adhesion, cytoskeleton remodeling, and differentiation. The aim of this project is to evaluate if YES1 is a driver gene in NSCLC, and if targeting its activation may be a potential new therapeutic strategy. We first evaluated the prognostic role of YES1 protein expression in two independent series of 76 and 234 NSCLC patients, respectively. In both series, the multivariate analysis revealed that high YES1 expression is an independent poor prognostic factor for overall survival (CUN series HR: 3.416 [0.933-12.508]; MD Anderson series HR: 1.570 [1.032-2.391]). We next evaluated the effect of YES1 knockdown in 5 NSCLC cell lines with YES1 amplification and overexpression, and in 3 cell lines without YES1 amplification and with low protein expression. YES1 downregulation by two specific siRNAs decreased proliferation and cell survival only in those cells overexpressing YES1. Congruently, YES1 inhibition led to apoptosis only in those cells. Consistent with these results, constitutive overexpression of YES1 in cells with low YES1 expression significantly enhanced cell proliferation. We next evaluated the effect of the multi-target Src kinase inhibitor dasatinib on the proliferation of NSCLC cell lines with high (8 cell lines) or low (4 cell lines) YES1 expression. Dasatinib dramatically inhibited proliferation in high YES1-expressing cell lines, whereas low YES1 cell lines were more resistant to dasatinib treatment (GI50s were four orders of magnitude higher in resistant cells). In conclusion, our results indicate that YES1 is a promising therapeutic target in NSCLC. Furthermore, amplification and high expression of YES1 may define a subset of patients who may potentially benefit from dasatinib treatment.

Cancer therapeutic target genes identified on chromosome 20q

Cancer therapeutic target genes identified on chromosome 20q