Integrated analysis of gene expression and genomic aberration data in osteosarcoma (OS)

Cytogenetic analyses have revealed that complex karyotypes with numerous and highly variable genomic aberrations including single-nucleotide polymorphisms (SNPs) and copy number variants (CNVs), are observed in most of the conventional osteosarcomas (OSs). Several genome-wide studies have reported that the dysregulated expression of many genes is correlated with genomic aberrations in OS. We first compared OS gene expression in Gene Expression Omnibus (GEO) data sets and genomic aberrations in International Cancer Genome Consortium (ICGC) database to identify differentially expressed genes (DEGs) associated with SNPs or CNVs in OS. Then the function annotation of SNP- or CNV-associated DEGs was performed in terms of gene ontology analysis, pathway analysis and protein-protein interactions (PPIs). Finally, the expression of genes correlated with both SNPs and CNVs were confirmed by quantitative reverse-transcription PCR. Eight publicly available GEO data sets were obtained, and a set of 979 DEGs were identified (472 upregulated and 507 downregulated DEGs). Moreover, we obtained 1039 SNPs mapped in 938 genes, and 583 CNV sites mapped in 2915 genes. Comparing genomic aberrations and DGEs, we found 41 SNP-associated DEGs and 124 CNV-associated DEGs, in which 7 DGEs were associated with both SNPs and CNVs, including WWP1, EXT1, LDHB, C8orf59, PLEKHA5, CCT3 and VWF. The result of function annotation showed that ossification, bone development and skeletal system development were the significantly enriched terms of biological processes for DEGs. PPI network analysis showed that CCT3, COPS3 and WWP1 were the significant hub proteins. We conclude that these genes, including CCT3, COPS3 and WWP1 are candidate driver genes of importance in OS tumorigenesis.