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.     

Discovery of a Good Responder Subtype of Esophageal Squamous Cell Carcinoma with Cytotoxic T-Lymphocyte Signatures Activated by Chemoradiotherapy

Definitive chemoradiotherapy (CRT) is a less invasive therapy for esophageal squamous cell carcinoma (ESCC). Five-year survival rate of locally advanced ESCC patients by definitive CRT were 37%. We previously reported that tumor-specific cytotoxic T-lymphocyte (CTL) activation signatures were preferentially found in long-term survivors. However, it is unknown whether the CTL activation is actually driven by CRT. We compared gene expression profiles among pre- and post-treatment biopsy specimens of 30 ESCC patients and 121 pre-treatment ESCC biopsy specimens. In the complete response (CR) cases, 999 overexpressed genes including at least 234 tumor-specific CTL-activation associated genes such as IFNG, PRF1, and GZMB, were found in post-treatment biopsy specimens. Clustering analysis using expression profiles of these 234 genes allowed us to distinguish the immune-activated cases, designating them as I-type, from other cases. However, despite the better CR rate in the I-type, overall survival was not significantly better in both these 30 cases and another 121 cases. Further comparative study identified a series of epithelial to mesenchymal transition-related genes overexpressed in the early relapse cases. Importantly, the clinical outcome of CDH2-negative cases in the I-type was significantly better than that of the CDH2-positive cases in the I-type. Furthermore, NK cells, which were activated by neutrophils-producing S100A8/S100A9, and CTLs were suggested to cooperatively enhance the effect of CRT in the CDH2-negative I-type. These results suggested that CTL gene activation may provide a prognostic advantage in ESCCs with epithelial characteristics.     

Identification of a DNA Methylome Profile of Esophageal Squamous Cell Carcinoma and Potential Plasma Epigenetic Biomarkers for Early Diagnosis

DNA methylation is a critical epigenetic mechanism involved in key cellular processes. Its deregulation has been linked to many human cancers including esophageal squamous cell carcinoma (ESCC). This study was designed to explore the whole methylation status of ESCC and to identify potential plasma biomarkers for early diagnosis. We used Infinium Methylation 450k array to analyze ESCC tissues (n = 4), paired normal surrounding tissues (n = 4) and normal mucosa from healthy individuals (n = 4), and combined these with gene expression data from the GEO database. One hundred and sixty eight genes had differentially methylated CpG sites in their promoter region and a gene expression pattern inverse to the direction of change in DNA methylation. These genes were involved in several cancer-related pathways. Three genes were validated in additional 42 ESCC tissues and paired normal surrounding tissues. The methylation frequency of EPB41L3, GPX3, and COL14A1 were higher in tumor tissues than in normal surrounding tissues (P<0.017). The higher methylation frequency of EPB41l3 was correlated with large tumor size (P = 0.044) and advanced pT tumor stage (P = 0.001). The higher methylation frequency of GPX3 and COL14A1 were correlated with advanced pN tumor stage (P = 0.001 and P<0.001). The methylation of EPB41L3, GPX3, and COL14A1 genes were only found in ESCC patients'' plasma, but not in normal individuals upon testing 42 ESCC patients and 50 healthy individuals. Diagnostic sensitivity was increased when methylation of any of the 3 genes were counted (64.3% sensitivity and 100% specificity). These differentially methylated genes in plasma may be used as biomarkers for early diagnosis of ESCC.     

Replication Study of ESCC Susceptibility Genetic Polymorphisms Locating in the ADH1B-ADH1C-ADH7 Cluster Identified by GWAS

China was one of the countries with highest esophageal squamous cell carcinoma (ESCC) incidence and mortality worldwide. Alcohol drinking has been identified as a major environmental risk-factor related to ESCC. The alcohol dehydrogenase (ADH) family are major enzymes involved in the alcohol-metabolizing pathways, including alcohol dehydrogenase 1B (ADH1B) and ADH1C. Interestingly, ADH1B and ADH1C genes locate tandemly with ADH7 in a genomic segment as a gene cluster, and are all polymorphic. Several ESCC susceptibility single nucleotide polymorphisms (SNPs) of the ADH1B-ADH1C-ADH7 cluster have been identified previously through a genome-wide association study (GWAS). In the study, we examined the association between five ADH1B-ADH1C-ADH7 cluster SNPs (rs1042026, rs17033, rs1614972, rs1789903 and rs17028973) and risk of developing ESCC. Genotypes were determined in two independent case-control sets from two regions of China. Odds ratios (ORs) and 95% confidence intervals (CIs) were estimated by logistic regression. Our data demonstrated that these ADH1B-ADH1C-ADH7 cluster SNPs confer susceptibility to ESCC in these two case-control sets, which were consistent to results of the previous GWAS.     

Focal Chromosomal Copy Number Aberrations Identify CMTM8 and GPR177 as New Candidate Driver Genes in Osteosarcoma

Osteosarcoma is an aggressive bone tumor that preferentially develops in adolescents. The tumor is characterized by an abundance of genomic aberrations, which hampers the identification of the driver genes involved in osteosarcoma tumorigenesis. Our study aims to identify these genes by the investigation of focal copy number aberrations (CNAs, <3 Mb). For this purpose, we subjected 26 primary tumors of osteosarcoma patients to high-resolution single nucleotide polymorphism array analyses and identified 139 somatic focal CNAs. Of these, 72 had at least one gene located within or overlapping the focal CNA, with a total of 94 genes. For 84 of these genes, the expression status in 31 osteosarcoma samples was determined by expression microarray analysis. This enabled us to identify the genes of which the over- or underexpression was in more than 35% of cases in accordance to their copy number status (gain or loss). These candidate genes were subsequently validated in an independent set and furthermore corroborated as driver genes by verifying their role in other tumor types. We identified CMTM8 as a new candidate tumor suppressor gene and GPR177 as a new candidate oncogene in osteosarcoma. In osteosarcoma, CMTM8 has been shown to suppress EGFR signaling. In other tumor types, CMTM8 is known to suppress the activity of the oncogenic protein c-Met and GPR177 is known as an overexpressed upstream regulator of the Wnt-pathway. Further studies are needed to determine whether these proteins also exert the latter functions in osteosarcoma tumorigenesis.     

TWIST1, MMP-21, and HLAG-1 co-overexpression is associated with ESCC aggressiveness

Esophageal squamous cell carcinoma (ESCC) is one of the most aggressive types of cancer, requiring reliable biomarkers for prognosis and therapeutic responsiveness. TWIST1, as an important factor responsible for metastasis of several cancers, is involved in tumor invasion and metastasis through indirectly regulation of MMP-21 expression. On the other hand, NF-ĸβ which is a regulator of HLAG-1 has direct interaction with TWIST1 protein. In this retrospective study we investigated the clinical significance of TWIST1, MMP-21, and HLAG-1 expression in ESCC, and the possible correlation between these genes and progression of the disease. The gene expression analyses of TWIST1, MMP-21, and HLA-G1 were performed by relative comparative real-time polymerase chain reaction in 58 ESCCs compared with corresponding margin-normal esophageal tissues. Significant overexpression of HLAG-1, TWIST1, and MMP-21 messenger RNA was observed in 22.4%, 41.4%, and 60.3% of tumor samples, respectively. Concomitant overexpression of TWIST1/MMP-21 and TWIST1/HLAG-1 were significantly correlated to each other in various clinicopathological features, including depth of tumor invasion, stage of tumor progression, lymphatic invasion, and grade of tumor cell differentiation ( P < 0.05). The current study is the first report of coexpression of TWIST1, MMP-21, and HLAG-1 in ESCC. Such findings suggest an oncogenic role for concomitant expression of these genes in ESCC invasion and metastasis.     

Estimating genomic instability mediated by Alu retroelements in breast cancer

Alu-PCR is a relatively simple technique that can be used to investigate genomic instability in cancer. This technique allows identification of the loss, gain or amplification of gene sequences based on the analysis of segments between two Alu elements coupled with quantitative and qualitative analyses of the profiles obtained from tumor samples, surgical margins and blood. In this work, we used Alu-PCR to identify gene alterations in ten patients with invasive ductal breast cancer. Several deletions and insertions were identified, indicating genomic instability in the tumor and adjacent normal tissue. Although not associated with specific genes, the alterations, which involved chromosomal bands 1p36.23, 1q41, 11q14.3, 13q14.2, occurred in areas of well-known genomic instability in breast and other types of cancer. These results indicate the potential usefulness of Alu-PCR in identifying altered gene sequences in breast cancer. However, caution is required in its application since the Alu primer can produce non-specific amplification.     

Peroxiredoxin 1 is involved in disassembly of flagella and cilia

Cilia/flagella are evolutionarily conserved cellular organelles. In this study, we demonstrated that Dunaliella salina Peroxiredoxin 1 (DsPrdx1) localized to the flagella and basal bodies, and was involved in flagellar disassembly. The link between DsPrdx1 and flagella of Dunaliella salina (D. salina) encouraged us to explore the function of its human homologue, Homo sapiens Peroxiredoxin 1 (HsPrdx1) in development and physiology. Our results showed that HsPrdx1 was overexpressed, and cilia were lost in esophageal squamous cell carcinoma (ESCC) cells compared with the non-cancerous esophageal epithelial cells Het-1A. Furthermore, when HsPrdx1 was knocked down by short hairpin RNA (shRNA) lentivirus in ESCC cells, the phenotype of cilia lost can be reversed, and the expression levels of tumor suppressor genes LKB1 and p-AMPK were increased, and the activity of the oncogene Aurora A was inhibited compared with those in cells transfected with scrambe-shRNA lentivirus. These findings firstly showed that Prdx1 is involved in disassembly of flagella and cilia, and suggested that the abnormal expression of the cilia-related gene including Prdx1 may affect both ciliogenesis and cancernogenesis.     

Cadherin-6 is a putative tumor suppressor and target of epigenetically dysregulated miR-429 in cholangiocarcinoma

Cholangiocarcinoma (CC) is a rare malignancy of the extrahepatic or intrahepatic biliary tract with an outstanding poor prognosis. Non-surgical therapeutic regimens result in minimally improved survival of CC patients. Global genomic analyses identified a few recurrently mutated genes, some of them in genes involved in epigenetic patterning. In a previous study, we demonstrated global DNA methylation changes in CC, indicating major contribution of epigenetic alterations to cholangiocarcinogenesis. Here, we aimed at the identification and characterization of CC-related, differentially methylated regions (DMRs) in potential microRNA promoters and of genes targeted by identified microRNAs. Twenty-seven hypermethylated and 13 hypomethylated potential promoter regions of microRNAs, known to be associated with cancer-related pathways like Wnt, ErbB, and PI3K-Akt signaling, were identified. Selected DMRs were confirmed in 2 independent patient cohorts. Inverse correlation between promoter methylation and expression suggested miR-129-2 and members of the miR-200 family (miR-200a, miR-200b, and miR-429) as novel tumor suppressors and oncomiRs, respectively, in CC. Tumor suppressor genes deleted in liver cancer 1 (DLC1), F-box/WD-repeat-containing protein 7 (FBXW7), and cadherin-6 (CDH6) were identified as presumed targets in CC. Tissue microarrays of a representative and well-characterized cohort of biliary tract cancers (n=212) displayed stepwise downregulation of CDH6 and association with poor patient outcome. Ectopic expression of CDH6 on the other hand, delayed growth in the CC cell lines EGI-1 and TFK-1, together suggesting a tumor suppressive function of CDH6. Our work represents a valuable repository for the study of epigenetically altered miRNAs in cholangiocarcinogenesis and novel putative, CC-related tumor suppressive miRNAs and oncomiRs.     

Targeted Knockdown of IQGAP1 Inhibits the Progression of Esophageal Squamous Cell Carcinoma In Vitro and In Vivo

IQGAP1 is a scaffolding protein that can regulate several distinct signaling pathways. The accumulating evidence has demonstrated that IQGAP1 plays an important role in tumorigenesis and tumor progression. However, the function of IQGAP1 in esophageal squamous cell carcinoma (ESCC) has not been thoroughly investigated. In the present study, we showed that IQGAP1 was overexpressed in ESCC tumor tissues, and its overexpression was correlated with the invasion depth of ESCC. Importantly, by using RNA interference (RNAi) technology we successfully silenced IQGAP1 gene in two ESCC cell lines, EC9706 and KYSE150, and for the first time found that suppressing IQGAP1 expression not only obviously reduced the tumor cell growth, migration and invasion in vitro but also markedly inhibited the tumor growth, invasion, lymph node and lung metastasis in xenograft mice. Furthermore, Knockdown of IQGAP1 expression in ESCC cell lines led to a reversion of epithelial to mesenchymal transition (EMT) progress. These results suggest that IQGAP1 plays crucial roles in regulating ESCC occurrence and progression. IQGAP1 silencing may therefore develop into a promising novel anticancer therapy.     

Long intergenic non-protein coding RNA 00858 participates in the occurrence and development of esophageal squamous cell carcinoma through the activation of the FTO-m6A-MYC axis by recruiting ZNF184

ObjectivesWe aimed at probing impact of LINC00858 on esophageal squamous cell carcinoma (ESCC) progression via ZNF184-FTO-m⁶A-MYC axis.MethodsExpression of related genes (LINC00858, ZNF184, FTO, and MYC) was detected in ESCC tissues or cells and their relationships were assessed. After expression alterations in ESCC cells, cell proliferation, invasion, migration, and apoptosis were detected. Tumor formation in nude mice was conducted.ResultsLINC00858, ZNF184, FTO, and MYC were overexpressed in ESCC tissues and cells. LINC00858 enhanced ZNF184 expression to upregulate FTO, which augmented MYC expression. LINC00858 knockdown diminished ESCC cell proliferative, migratory, and invasive properties while elevating apoptosis, which was negated by FTO overexpression. FTO knockdown exerted similar functions of LINC00858 knockdown on ESCC cell movements, which was annulled by MYC upregulation. Silencing LINC00858 repressed tumor growth and related gene expression in nude mice.ConclusionsLINC00858 modulated MYC m⁶A modification via FTO by recruiting ZNF184, thus facilitating ESCC progression.     

Epigenetic inactivation of SPINT2 is associated with tumor suppressive function in esophageal squamous cell carcinoma

Hepatocyte growth factor activator inhibitor type 2 (SPINT2), a Kunitz-type serine proteinase inhibitor, has been identified as a putative tumor suppressor gene silenced by promoter methylation. We aimed to investigate whether SPINT2 might act as an esophageal squamous cell carcinoma (ESCC) tumor suppressor gene. Four ESCC cell lines, Fifty-two ESCC tissues and twenty-nine neighboring non-cancerous tissues were included in this study. The expression of SPINT2 was monitored by real time PCR. Bisulfite genomic sequencing and methylation-specific PCR were used to analyze methylation status. The effect of SPINT2 on cell proliferation and apoptosis in EC109 and EC9706 cells was observed by CCK-8 assay and flow cytometric analysis. We found that silencing of SPINT2 was associated with promoter methylation in ESCC cell lines. The densely methylated SPINT2 promoter region was confirmed by bisulfite genomic sequencing. Ectopic expression of SPINT2 inhibited cell proliferation through inducing cell apoptosis in vitro. Furthermore, methylation-specific PCR analysis revealed that SPINT2 promoter methylation was prominent in carcinoma tissues (52.08%) compared with neighboring non-cancerous tissues (22.58%). Kaplan–Meier analysis showed that patients with SPINT2 hypermethylation had shorter survival time. The tumor suppressor gene of SPINT2 is commonly silenced by promoter hypermethylation in human ESCC and SPINT2 hypermethylation is correlated with poor overall survival, implicating SPINT2 is an underlying prognostic marker for human ESCC.     

ERCC1 Single Nucleotide Polymorphism C8092A,but Not Its Expression Is Associated with Survival of Esophageal Squamous Cell Carcinoma Patients from Fujian Province,China

Esophageal carcinoma is one of the world''s deadliest cancers. Esophageal squamous cell carcinoma (ESCC) is more frequent than adenocarcenoma (AC) in China. Platinum-based chemotherapy with surgical resection is a common treatment approach for ESCC; however, the treatment response is uncertain. Evidence suggests polymorphisms in genes encoding excision repair cross-complementing group 1 (ERCC1), a protein involved in nuclear excision repair (NER), may help predict response to cisplatin and other platinum-based chemotherapeutics. Multiple ERCC1 single nucleotide polymorphisms (SNPs) have been associated with platinum chemotherapy response. Two common SNPs occur at the C8092A and C118T loci. Our study aimed to determine if 1) an association exists between ERCC1 tumor expression and patient survival, 2) whether adjuvant therapy influence on survival is related to histological ERCC1 presence in tumor cell nuclei, and 3) whether other clinicopathological characteristics in a cohort of patients following surgery for various stages of ESCC are associated with tumor ERCC1 expression. One hundred eight patients were included in the study, and tumor biopsy was collected for genotyping and immunohistochemical analysis of ERCC1. Sixty-seven patients (62%) received no adjuvant therapy, and the rest had either platinum-based chemotherapy (28.5%), radiotherapy (6.5%) or both treatments (2.8%). Log-rank analysis revealed no significant connection between tumor ERCC1 expression (P = 0.12) or adjuvant therapy (P = 0.56) on patient survival. Also, non-parametric Mann-Whitney analysis showed no significant link between tumor size or nodus tumor formation and ERCC1 presence in patients in the study. Interestingly, C8092A SNP showed significant association with patient survival (P = 0.01), with patients homozygous for the mutant allele showing the most significantly reduced survival (P = 0.04) compared to those homozygous for the dominant allele (CC). Our results provide novel insight into the genotypic variation of patients from Quanzhou, Fujian province China.     

Clinical Relevance of Genomic Changes in Recurrent Pediatric Solid Tumors

PURPOSE: Relapsed/refractory pediatric cancers show poor prognosis; however, their genomic patterns remain unknown. To investigate the genetic mechanisms of tumor relapse and therapy resistance, we characterized genomic alterations in diagnostic and relapsed lesions in patients with relapsed/refractory pediatric solid tumors using targeted deep sequencing. PATIENTS AND METHODS: A targeted sequencing panel covering the exons of 381 cancer genes was used to characterize 19 paired diagnostic and relapsed samples from patients with relapsed/refractory pediatric solid tumors. RESULTS: The mean coverage for all samples was 930.6× (SD = 213.8). Among the 381 genes, 173 single nucleotide variations (SNVs)/insertion-deletions (InDels), 100 copy number alterations, and 1 structural variation were detected. A total of 72.6% of SNVs in primary tumors were also found in recurrent lesions, and 27.2% of SNVs in recurrent tumors had newly occurred. Among SNVs/InDels detected only in recurrent lesions, 71% had a low variant allele fraction (<10%). Patients were classified into three categories based on the mutation patterns after cancer treatment. A significant association between the major mutation patterns and clinical outcome was observed. Patients whose relapsed tumor had fewer mutations than the diagnostic sample tended to be older, had longer progression-free survival, and achieved complete remission after relapse. Contrastingly, patients whose genetic profile only had concordant mutations without any change had the worst outcome. CONCLUSIONS: We characterized genomic changes in recurrent pediatric solid tumors. These findings could help to understand the biology of relapsed childhood cancer and to develop personalized treatment based on their genetic profile.