Integrative analysis of copy number and transcriptional expression profiles in esophageal cancer to identify a novel driver gene for therapy

From modulation of cellular plasticity to potentiation of therapeutic resistance: new and emerging roles of MYB transcription factors in human malignancies

MYB transcription factors are encoded by a large family of highly conserved genes from plants to vertebrates. There are three members of the MYB gene family in human, namely, MYB, MYBL1, and MYBL2 that encode MYB/c-MYB, MYBL1/A-MYB, and MYBL2/B-MYB, respectively. MYB was the first member to be identified as a cellular homolog of the v-myb oncogene carried by the avian myeloblastosis virus (AMV) causing leukemia in chickens. Under the normal scenario, MYB is predominantly expressed in hematopoietic tissues, colonic crypts, and neural stem cells and plays a role in maintaining the undifferentiated state of the cells. Over the years, aberrant expression of MYB genes has been reported in several malignancies and recent years have witnessed tremendous progress in understanding of their roles in processes associated with cancer development. Here, we review various MYB alterations reported in cancer along with the roles of MYB family proteins in tumor cell plasticity, therapy resistance, and other hallmarks of cancer. We also discuss studies that provide mechanistic insights into the oncogenic functions of MYB transcription factors to identify potential therapeutic vulnerabilities.

A Nutrient-Deficient Microenvironment Facilitates Ferroptosis Resistance via the FAM60A-PPAR Axis in Pancreatic Ductal Adenocarcinoma

Ferroptosis, a nonapoptotic form of cell death, is an emerging potential therapeutic target for various diseases, including cancer. However, the role of ferroptosis in pancreatic cancer remains poorly understood. Pancreatic ductal adenocarcinoma (PDAC) is characterized by a poor prognosis and chemotherapy resistance, attributed to its high Kirsten rats arcomaviral oncogene homolog mutation rate and severe nutritional deficits resulting from a dense stroma. Several studies have linked rat sarcoma (RAS) mutations to ferroptosis, suggesting that inducing ferroptosis may be an effective strategy against oncogenic RAS-bearing tumors. We investigated the role of Family With Sequence Similarity 60 Member A (FAM60A) in this study, a protein closely associated with a poor prognosis and highly expressed in PDAC and tumor tissue from KrasG12D/+;Trp53R172H/+; Pdx1-Cre mice, in regulating ferroptosis, tumor growth, and gemcitabine sensitivity in vitro and in vivo. Our results demonstrate that FAM60A regulates 3 essential metabolic enzymes, ACSL1/4 and GPX4, to protect PDAC cells from ferroptosis. Furthermore, we found that YY1 transcriptionally regulates FAM60A expression by promoting its transcription, and the Hippo-YY1 pathway is restricted in the low-amino-acid milieu in the context of nutrient deprivation, leading to downstream suppression of peroxisome proliferator-activated receptor and ACSL1/4 and activation of GPX4 pathways. Importantly, FAM60A knockdown sensitized PDAC cells to gemcitabine treatment. A new understanding of FAM60A transcriptional regulation pattern in PDAC and its dual function in ferroptosis reliever and chemotherapy resistance is provided by our study. Targeting FAM60A may therefore offer a promising therapeutic approach for PDAC by simultaneously addressing 2 major features of the disease (high RAS mutation rate and tumor microenvironment nutrient deficiency) and preventing tumor cell metabolic adaptation.

Integrative Analyses of Pyrimidine Salvage Pathway-Related Genes Revealing the Associations Between UPP1 and Tumor Microenvironment

Background

The pyrimidine salvage pathway plays a critical role in tumor progression and patient outcomes. The roles of pyrimidine salvage pathway-related genes (PSPGs) in cancer, however, are not fully understood. This study aims to depict the characteristics of PSPGs across various cancers.

Methods

An integrative pan-cancer analysis of six PSPGs (CDA, UCK1, UCK2, UCKL1, UPP1, and UPP2) was conducted using TCGA data, single-cell RNA sequencing datasets, and patient samples. Single-cell transcriptome analysis and RT-qPCR were used to validate the relation between UPP1 and cytokines. Flow cytometry was performed to validate the role of UPP1 in immune checkpoint regulation. The correlation between UPP1 and tumor associated neutrophils (TAN) were investigated and validated by single-cell transcriptome analysis and tissue microarrays (TMAs).

Results

PSPGs showed low mutation rates but significant copy number variations, particularly amplifications in UCKL1, UPP1, and UCK2 across various cancers. DNA methylation patterns varied, with notable negative correlations between methylation and gene expression in UPP1. PSPGs were broadly up-regulated in multiple cancers, with correlations to clinical staging and prognosis. Proteomic data further confirmed these findings. Functional analysis revealed PSPGs' associations with tumor proliferation, metastasis, and various signaling pathways. UPP1 showed strong correlations with the tumor microenvironment (TME), particularly with cytokines, immune checkpoints, and various immune cells. Single-cell transcriptome analysis confirmed these associations, highlighting UPP1's influence on cytokine expression and immune checkpoint regulation. In esophageal squamous cell carcinoma (ESCC), UPP1-high tumor cells were significantly associated with immunosuppressive cells in the TME. Spatial analysis using TMAs revealed that UPP1+ tumor cells were predominantly located at the invasive margin and closely associated with neutrophils, correlating with poorer patient prognosis.

Conclusion

Our study depicted the multi-dimensional view of PSPGs in cancer, with a particular focus on UPP1's role in the TME. Targeting UPP1 holds promise as a potential strategy for cancer therapy.

FAM60A promotes osteosarcoma development and progression

BACKGROUND

Osteosarcoma (OS) is a highly malignant primary bone tumor. Family of homology 60A (FAM60A) reportedly contributes to the malignant growth of some tumors.

METHODS

Herein we investigated the mRNA expression level of FAM60A by combining OS and non-cancer samples from public databases. Immunohistochemistry was performed to determine protein expression levels of FAM60A in patients with OS. Further, RT-qPCR and western blotting were conducted to evaluate FAM60A expression in various OS cell lines. CCK-8 assay, colony formation assay, and flow cytometry were applied to determine the function of FAM60A. Finally, functional enrichment analysis was performed based on FAM60A co-expressed genes.

RESULTS

FAM60A mRNA expression level was found to be significantly upregulated (standardized mean difference = 1.27, 95% CI [0.67-1.88]). Survival analyses suggested that higher expression of FAM60A was indicative of poor prognoses. Similarly, FAM60A protein expression level was also observed to be upregulated. Knocking down FAM60A expression inhibited OS cell proliferation, increased apoptosis, and blocked cells from entering the S phase. Besides, cell cycle was the most prominently enriched pathway, and BUB1, DTL, and EXO1 were identified as hub genes.

CONCLUSIONS

FAM60A expression was found to be markedly upregulated in OS; furthermore, FAM60A was observed to promote OS cell proliferation, inhibit apoptosis, and participate in cell cycle regulation. Besides, FAM60A may interact with hub genes to participate in the progress of OS.

Risk scoring based on DNA methylation-driven related DEGs for colorectal cancer prognosis with systematic insights

Colorectal cancer is a common malignant tumor of the digestive tract. Despite advances in diagnostic techniques and medications. Its prognosis remains challenging. DNA methylation-driven related circulating tumor cells have attracted enormous interest in diagnosing owing to their non-invasive nature and early recognition properties. However, the mechanism through which risk biomarkers act remains elusive. Here, we designed a risk model based on differentially expressed genes, DNA methylation, robust, and survival-related factors in the framework of Cox regression. The model has satisfactory performance and is independently verified by an external and isolated dataset in terms of C-index value, ROC, and tROC. The model was applied to Colorectal cancer patients who were subsequently divided into high- and low-risk groups. Functional annotations, genomic alterations, tumor immune environment, and drug sensitivity were analyzed. We observed that up-regulated genes are associated with epithelial cell differentiation and MAPK signaling pathways. The down-regulated genes are related to IL-7 signaling and apoptosis-induced DNA fragmentation. Interestingly, the immune system was inhibited in high-risk groups. High-frequency mutation genes tend to co-occur. High-risk score patients are related to copy number amplification events. To address the challenges, we suggested eleven and twenty-one drugs that are sensitive to low- and high-risk patients. Finally, an artificial neural network was provided to evaluate the immunotherapeutic efficiency. Taken together, the findings demonstrated that our risk score model is robust and reliable for evaluating the prognosis with novel diagnostic and treatment targets. It also yields benefits for the treatment and provides unique insights into developing therapeutic strategies.

Dysregulated ceRNA network modulated by copy number variation-driven lncRNAs in breast cancer: A comprehensive analysis

Breast cancer is a malignancy harmful to physical and mental health in women, with quite high mortality. Copy number variations (CNVs) are vital factors affecting the progression of breast cancer. Detecting CNVs in breast cancer to predict the prognosis of patients has become a promising approach to accurate treatment in recent years. The differential analysis was performed on CNVs of long noncoding RNAs (lncRNAs) as well as the expression of lncRNAs, microRNAs (miRNAs) and mRNAs in normal tissue and breast tumor tissue based on The Cancer Genome Atlas (TCGA) database. The CNV-driven lncRNAs were identified by the Kruskal-Wallis test. Meanwhile, a competitive endogenous RNA (ceRNA) network regulated by CNV-driven lncRNA was constructed. As the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed, the mRNAs in the dysregulated ceRNA network were mainly enriched in the biological functions and signaling pathways, including the Focal Adhesion-PI3K-Akt-mTOR-signaling pathway, the neuronal system, metapathway biotransformation Phase I and II and blood circulation, etc. The relationship between the CNVs of five lncRNAs and their gene expression in the ceRNA network was analyzed via a chi-square test, which confirmed that except for LINC00243, the expression of four lncRNAs was notably correlated with the CNVs. The survival analysis revealed that only the copy number gain of LINC00536 was evidently related to the poor prognosis of patients. The CIBERSORT algorithm showed that five lncRNAs were correlated with the abundance of immune cell infiltration and immune checkpoints. In a word, by analyzing CNV-driven lncRNAs and the ceRNA network regulated by these lncRNAs, this study explored the mechanism of breast cancer and provided novel insights into new biomarkers.

Identification of a ceRNA Network Driven by Copy Number Variations in Esophageal Cancer

BACKGROUND

Copy number variation (CNV) is associated with progression of esophageal cancer (EC), a common gastrointestinal neoplasm.

METHODS

Using sequencing data, CNV data, and clinical data of EC transcriptome samples obtained from public databases, we performed differential expression analysis on sequencing data. Differentially expressed CNV-driven lncRNAs were screened using the chi-square test, and CNV-driven lncRNA-associated miRNAs and mRNAs were predicted. Cytoscape software was then used to construct ceRNA networks. Gene ontology and Kyoto Encyclopedia of Genes and Genomes analyses were performed to investigate biological functions of mRNAs in the ceRNA network. Survival curves were plotted to explore correlations between lncRNAs in the ceRNA network and overall survival of CNV patients. Multiple databases were used to predict lncRNAs-related drugs.

RESULTS

A dysregulated lncRNA-associated ceRNA network driven by CNV in EC, including 11 lncRNAs, 11 miRNAs and 159 mRNAs, was constructed. Downstream enrichment of mRNAs was related to biological processes such as extracellular matrix organization, indicating that these mRNAs mainly participate in intercellular exchange between tumor cells. Additionally, expression of all lncRNAs in the ceRNA network, except LINC00950, LINC01270 and MIR181A1HG, was correlated with patients' CNV. In addition, none of the 11 lncRNAs was significantly correlated with overall survival of CNV patients. CH5424802 and PD-033299CNV mainly affected the RTK signaling pathway and the cell cycle of tumor cells via RP11-180N14.1 and RP11-273 G15.2 in the ceRNA network.

CONCLUSIONS

This study identified 11 CNV-driven lncRNAs that might affect EC development, 2 of which have promising effects if applied to drug treatment. These findings might assist in identifying novel treatments for EC.

Chromosomal Instability, Selection and Competition: Factors That Shape the Level of Karyotype Intra-Tumor Heterogeneity

Simple Summary

Each cancer consists of billions of cells. These cells are far from identical; hence, the population of cells that constitute a tumor is heterogeneous. A salient property that varies between cells in a tumor is their karyotype, the number and configuration of the chromosomes. The level of karyotype heterogeneity can be used to predict the survival of a patient. In this review, we describe the processes that shape the level of karyotype heterogeneity in a cancer.

Abstract

Intra-tumor heterogeneity (ITH) is a pan-cancer predictor of survival, with high ITH being correlated to a dismal prognosis. The level of ITH is, hence, a clinically relevant characteristic of a malignancy. ITH of karyotypes is driven by chromosomal instability (CIN). However, not all new karyotypes generated by CIN are viable or competitive, which limits the amount of ITH. Here, we review the cellular processes and ecological properties that determine karyotype ITH. We propose a framework to understand karyotype ITH, in which cells with new karyotypes emerge through CIN, are selected by cell intrinsic and cell extrinsic selective pressures, and propagate through a cancer in competition with other malignant cells. We further discuss how CIN modulates the cell phenotype and immune microenvironment, and the implications this has for the subsequent selection of karyotypes. Together, we aim to provide a comprehensive overview of the biological processes that shape the level of karyotype heterogeneity.

Restraint of FAM60A has a cancer-inhibiting role in pancreatic carcinoma via the effects on the Akt/GSK-3β/β-catenin signaling pathway

Family with sequence similarity 60A (FAM60A) has been reported as a new cancer-related protein that affects the malignant progression of some cancers. However, whether FAM60A plays a part in pancreatic carcinoma is undetermined. This work was designed to examine the impact of FAM60A in pancreatic carcinoma. Abundant expression of FAM60A was observed in the primary tumor tissue of pancreatic carcinoma. Moreover, a high FAM60A level was related to a poor overall survival in pancreatic carcinoma patients. Malignant behaviors of pancreatic carcinoma cells, such as proliferation and invasiveness, were markedly affected by FAM60A depletion. In addition, FAM60A depletion enhanced the drug sensitivity of pancreatic carcinoma cells to gemcitabine. Further study revealed that FAM60A depletion impaired the activities of Akt and β-catenin. Inhibiting the activity of Akt abolished FAM60A-mediated β-catenin activation. Re-expression of β-catenin partially diminished the FAM60A-depletion-mediated cancer suppressive effect in pancreatic carcinoma cells. In vivo experiments demonstrated that FAM60A depletion prohibited the xenograft formation of pancreatic carcinoma cells, with concurrent reductions of Akt and β-catenin activities. Collectively, our findings indicate that FAM60A exerts a cancer-promoting role in pancreatic carcinoma through affection of the Akt/β-catenin pathway. This work indicates that FAM60A acts as a tumor promoter in pancreatic carcinoma and can be utilized as a potential target for anti-pancreatic carcinoma therapy development.

Identification and validation of the N6-methyladenosine RNA methylation regulator ZC3H13 as a novel prognostic marker and potential target for hepatocellular carcinoma

N6-methyladenosine (m6A) RNA methylation has been implicated in various malignancies. This study aimed to identify prognostic signature based on m6A methylation regulators for hepatocellular carcinoma (HCC) and provide candidate targets for HCC treatment. In this study, the expression levels, prognostic values, correlation with tumor grades and genetic variations of m6A-related genes in HCC were evaluated using bioinformatics analyses. Interestingly, the results show that methyltransferase zinc finger CCCH-type containing 13 (ZC3H13) was expressed at a significantly low level in HCC. Survival outcome analysis suggested that significant correlations existed between ZC3H13 downregulation and poor overall survival (OS) and poor recurrence-free survival (RFS) in HCC patients. Therefore, ZC3H13 was chosen for further experimental validation. The expression of ZC3H13 in HCC cell lines was investigated by western blotting. Knockdown of ZC3H13 significantly enhanced the migration and invasion of HCC cells, as demonstrated by wound healing and transwell assays. Moreover, upregulating ZC3H13 repressed the growth of xenograft tumors in vivo. Functional and pathway enrichment analyses indicated that ZC3H13 might be involved in transcriptional dysregulation or the JAK-STAT signaling pathway in cancer. Additionally, ZC3H13 expression was significantly correlated with lymphocytes and immunomodulators. Therefore, ZC3H13 is a promising candidate as a novel biomarker and therapeutic target for HCC.

CNViz: An R/Shiny Application for Interactive Copy Number Variant Visualization in Cancer

Copy number variants (CNVs) comprise a class of mutation which includes deletion, duplication, or amplification events that range in size from smaller than a single-gene or exon, to the size of a full chromosome. These changes can affect gene expression levels and are thus implicated in disease, including cancer. Although a variety of tools and methodologies exist to detect CNVs using data from massively parallel sequencing (also referred to as next-generation sequencing), it can be difficult to appreciate the copy number profile in a list format or as a static image. CNViz is a freely accessible R/Bioconductor package that launches an interactive R/Shiny visualization tool to facilitate review of copy number data. As inputs, it requires genomic locations and corresponding copy number ratios for probe, gene, and/or segment-level data. If supplied, loss of heterozygosity (LOH), focal variant data [single nucleotide variants (SNVs) and small insertions and deletions (indels)], and metadata (e.g., specimen purity and ploidy) can also be incorporated into the visualization. The CNViz R/Bioconductor package is an easy-to-use tool built with the intent of encouraging visualization and exploration of copy number variation. CNViz can be used in a clinical setting as well as for research to study patterns in human cancers more broadly. The intuitive interface allows users to visualize the copy number profile of a specimen, dynamically change resolution to explore gene and probe-level copy number changes, and simultaneously integrate LOH, SNV, and indel findings. CNViz is available for download as an R package via Bioconductor. An example of the application is available at rebeccagreenblatt.shinyapps.io/cnviz_example.

Trafficking Protein TMED3 Promotes Esophageal Squamous Cell Carcinoma

BACKGROUND

The molecular mechanisms of esophageal squamous cell carcinoma (ESCC) remain poorly understood. Transmembrane emp24 trafficking protein 3 (TMED3) acts as an oncogene or tumor suppressor gene in different cancers. Our study was to explore the clinicopathological significance and functional roles of TMED3 in ESCC.

METHODS

Immunohistochemistry, qPCR, and western blotting were used to analyze the expression of TMED3 in ESCC tissues and cells. Statistical analysis was performed to analyze the relationship between TMED3 expression and tumor characteristics in patients with ESCC. The role of TMED3 in vitro and in vivo was investigated by performing functional verification experiments and using a xenograft mouse model. Proteins that are functionally related to TMED3 were recognized by Affymetrix microarray and Ingenuity Pathway Analysis (IPA). Functional verification experiments were performed to analyze the role of FAM60A (a protein functionally related to TMED3) in vitro.

RESULTS

We confirmed the overexpression of TMED3 was correlated with poor prognosis in ESCC patients. When TMED3 was knocked down, ESCC cell proliferation, migration, and invasion were inhibited whereas cell apoptosis was promoted in vitro, and tumorigenicity was inhibited in vivo. We further revealed significant changes in gene expression profile in TMED3 knockdown cells. Among these differentially expressed genes, FAM60A was overexpressed in ESCC tissues. Furthermore, knocking down FAM60A significantly weakened the proliferation ability of ESCC cells and reversed TMED3's tumorigenicity of ESCC cells.

CONCLUSION

Our study revealed an oncogenic role of TMED3 in ESCC.

A single-nucleotide-polymorphism in the 5′-flanking region of MSX1 gene as a predictive marker candidate for platinum-based therapy of esophageal carcinoma

Background:

Platinum derivatives are important treatment options for patients with esophageal carcinoma (EC), and a predictive marker for platinum-based therapy is needed for precision medicine.

Patients and methods:

This study contained two cohorts consisting of EC patients treated using platinum-based chemoradiation therapy (CRT) as the first-line and another external cohort of nationwide clinicogenomic data from the BioBank Japan (BBJ).

Results:

Genome-wide association study (GWAS) of therapeutic outcomes, refractory disease or not, following platinum-based CRT as first-line in 94 patients in the first cohort suggested the association of 89 SNPs using p < 0.0001. The top 10 SNPs selected from each chromosomal region by odds ratio were evaluated for progression-free survival (PFS) and overall survival (OS) hazard ratios in the first cohort, resulting in four candidates (p < 0.0025). The four selected candidates were re-evaluated in another cohort of 24 EC patients, which included patients prospectively enrolled in this study to fulfill the sample size statistically suggested by the results of the first cohort, and of the four, only rs3815544 was replicated (p < 0.0125). Furthermore, this candidate genotype of rs3815544 proceeded to the re-evaluation study in an external cohort consisting of EC patients treated with platinum derivatives and/or by radiation therapy as the first-line treatment in BBJ, which confirmed that the alternative allele (G) of rs3815544 was statistically associated with non-response (SD or PD) to platinum-based therapy in EC patients (odds ratio = 1.801, p = 0.048). The methylation QTL database as well as online clinicogenomic databases suggested that the region including rs3815544 may regulate MSX1 expression through CpG methylation, and this down-regulation was statistically associated with poor prognosis after platinum-based therapies for EC.

Conclusion:

rs3815544 is a novel candidate predictive marker for platinum-based EC therapy.

Overexpression of Epithelial Splicing Regulatory Protein 1 in Metastatic Lesions of Serous Ovarian Carcinoma Correlates with Poor Patient Prognosis

Background: Epithelial splicing regulatory proteins (ESRPs) can regulate alternative splicing of RNA and play roles in tumorigenesis and development of various malignancies. In this study, bioinformatic analyses and immunohistochemistry (IHC) were used to investigate the function of ESPRs in serous ovarian carcinoma (SOC) oncogenesis and metastasis. Materials and Methods: The mRNA levels of ESRPs were analyzed by Oncomine and gene expression profiling interactive analysis (GEPIA). Prognostic values of ESRPs were analyzed by GEPIA and the UALCAN website. Genetic variations of ESRPs were analyzed by cBioPortal. ESRP1 was selected for further research. The relationship between ESRP1 and immunoregulatory molecules was studied by using the TISIDB database. ESRP1 protein expression in OC was investigated via IHC assays. Results: ESRP1 and ESRP2 mRNA were significantly upregulated in SOC (p < 0.05). The prognostic value of ESRP1 mRNA in SOC was inconsistent, and ESRP2 mRNA level did not relate to prognosis for OC patients. The IHC results showed higher ESRP1 expression in OC tissues than in normal ovarian tissues (p = 0.002), and ESRP1 expression in metastatic lesions of OC patients was higher than in paired primary OC tissues (p = 0.035). The ESRP1 expression was related to FIGO stage, differentiation, and peritoneal metastasis (p = 0.016; 0.031; 0.038, respectively). The ESRP1 switch (the differential expression of ESRP1 between metastatic and primary tumor of ovarian carcinoma) was significantly associated with E-cadherin expression in metastatic OC tumors (p = 0.012). The ESRP1 expression in both metastasis and ESRP1 switch significantly correlated with poor prognosis of OC patients (p = 0.045; 0.038, respectively), and ESRP1 switch and FIGO stage were independent risk factors for OC patient prognosis (p = 0.033; 0.009, respectively). Conclusions: The ESRP1 may promote OC metastasis by promoting OC cell colonization via the mesenchymal-epithelial transition (MET) process. The ESRP1 expression in metastatic lesions of OC patients may be a biomarker for predicting prognosis and a potential therapeutic target in OC.

Loss of FAM60A attenuates cell proliferation in glioma via suppression of PI3K/Akt/mTOR signaling pathways

BACKGROUND

Glioma is a common malignant tumor of the central nervous system with a high incidence and mortality. Family with sequence similarity 60 member A (FAM60A) is a new subunit of the Sin3 deacetylase complex. The clinical significance and biologic role of FAM60A in glioma remain unclear.

METHODS

The expression of FAM60A in normal glial cells, glioma cells, and five-paired gliomas, and adjacent noncancerous tissues was quantified using real-time polymerase chain reaction (PCR) and western blotting. FAM60A protein expression in 179 archived, paraffin-embedded glioma samples was analyzed using immunohistochemistry. The roles of FAM60A in glioma cell proliferation and tumorigenicity were explored in vitro and in vivo. The underlying molecular mechanisms were elucidated using Western blot assay. Serum exosomal FAM60A levels of glioma patients were detected using electron microscopy, western blot, and real-time PCR.

RESULTS

FAM60A expression was significantly up-regulated in glioma tissues and cell lines and positively associated with a worse outcome in glioma. Knockdown of FAM60A could inhibit glioma cell proliferation and tumorigenicity in vitro and in vivo. Besides, FAM60A expression was detectable in extracted serum exosomes with a higher expression in the glioma cancer group than in the normal group.

CONCLUSIONS

Loss of FAM60A attenuates cell proliferation in glioma by suppressing PI3K/Akt/mTOR signaling pathways. Therefore, FAM60A may act as a prognostic biomarker and therapeutic target for glioma.

A five-gene signature for predicting overall survival of esophagus adenocarcinoma

Esophageal adenocarcinoma (EAC) is common and aggressive with increasing trend of incidence. Urgent need for an effective signature to assess EAC prognosis and facilitate tailored treatment is required.Differentially expressed mRNAs (DEMs) were identified by analyzing EAC tissues and adjacent normal samples from The Cancer Genome Atlas (TCGA). Then univariate regression analyses were performed to confirm prognostic DEMs. We used least absolute shrinkage and selection operator (LASSO) to build a prognostic mRNA signature whose performance was assessed by Kaplan-Meier curve, receiver operating characteristic (ROC). GSE72874 were used as an external test set. The performances of the signature were also validated in internal TCGA and external test sets. Gene set enrichment analysis (GSEA) and tumor immunity analysis were performed to decipher the biological mechanisms of the signature.A 5-mRNA signature consisted of SLC26A9, SINHCAF, MICB, KRT19, and MT1X was developed to predict prognosis of EAC. The 5-mRNA signature was promising as a biomarker for predicting 3-year survival rate of EAC in the internal test set, the entire TCGA set, and the external test set with areas under the curve (AUC) = 0.849, 0.924, and 0.747, respectively. Patients were divided into low- and high-risk groups based on risk scores of the signature. The high-risk group was mainly associated with cancer-related pathways and low levels of B cell infiltration.The 5-mRNA prognostic signature we identified can reliably predict prognosis and facilitate individualized treatment decisions for EAC patients.

Identification of DNA-Repair-Related Five-Gene Signature to Predict Prognosis in Patients with Esophageal Cancer

Esophageal cancer (ESCA) is a leading cause of cancer-related mortality, with poor prognosis worldwide. DNA damage repair is one of the hallmarks of cancer. Loss of genomic integrity owing to inactivation of DNA repair genes can increase the risk of cancer progression and lead to poor prognosis. We aimed to identify a novel gene signature related to DNA repair to predict the prognosis of ESCA patients. Based on gene expression profiles of ESCA patients from The Cancer Genome Atlas and gene set enrichment analysis, 102 genes related to DNA repair were identified as candidates. After stepwise Cox regression analysis, we established a five-gene prognostic model comprising DGCR8, POM121, TAF9, UPF3B, and BCAP31. Kaplan-Meier survival analysis confirmed a strong correlation between the prognostic model and survival. Moreover, we verified the clinical value of the prognostic signature under the influence of different clinical parameters. We found that small-molecule drugs (trametinib, selumetinib, and refametinib) could help to improve patient survival. In summary, our study provides a novel and promising prognostic signature based on DNA-repair-related genes to predict survival of patients with ESCA. Systematic data mining provides a theoretical basis for further exploring the molecular pathogenesis of ESCA and identifying therapeutic targets.

Pan-cancer driver copy number alterations identified by joint expression/CNA data analysis

Analysis of large gene expression datasets from biopsies of cancer patients can identify co-expression signatures representing particular biomolecular events in cancer. Some of these signatures involve genomically co-localized genes resulting from the presence of copy number alterations (CNAs), for which analysis of the expression of the underlying genes provides valuable information about their combined role as oncogenes or tumor suppressor genes. Here we focus on the discovery and interpretation of such signatures that are present in multiple cancer types due to driver amplifications and deletions in particular regions of the genome after doing a comprehensive analysis combining both gene expression and CNA data from The Cancer Genome Atlas.

Future directions in advanced penile cancer - mechanisms of carcinogenesis and a search for targeted therapy

Penile squamous cell carcinoma (SCC) is a rare and aggressive urological malignancy. Advanced penile SCC requires multimodal management, including surgery and systemic therapy. Given its rarity, there have been few substantial advances in our understanding of the molecular and genomic drivers of penile SCC, especially for patients with relapsed or advanced disease. In this review, we discuss the molecular and genomic landscape of penile SCC, clinical trials in progress and implications for novel therapeutic targets. Future work should focus on preclinical models to provide a platform for investigation and validation of new molecular pathways for testing of therapeutics.

Systematic Establishment of Robustness and Standards in Patient-Derived Xenograft Experiments and Analysis

Patient-derived xenografts (PDX) are tumor-in-mouse models for cancer. PDX collections, such as the NCI PDXNet, are powerful resources for preclinical therapeutic testing. However, variations in experimental and analysis procedures have limited interpretability. To determine the robustness of PDX studies, the PDXNet tested temozolomide drug response for three prevalidated PDX models (sensitive, resistant, and intermediate) across four blinded PDX Development and Trial Centers using independently selected standard operating procedures. Each PDTC was able to correctly identify the sensitive, resistant, and intermediate models, and statistical evaluations were concordant across all groups. We also developed and benchmarked optimized PDX informatics pipelines, and these yielded robust assessments across xenograft biological replicates. These studies show that PDX drug responses and sequence results are reproducible across diverse experimental protocols. In addition, we share the range of experimental procedures that maintained robustness, as well as standardized cloud-based workflows for PDX exome-sequencing and RNA-sequencing analyses and for evaluating growth. SIGNIFICANCE: The PDXNet Consortium shows that PDX drug responses and sequencing results are reproducible across diverse experimental protocols, establishing the potential for multisite preclinical studies to translate into clinical trials.

FAM60A, increased by Helicobacter pylori, promotes proliferation and suppresses apoptosis of gastric cancer cells by targeting the PI3K/AKT pathway

Helicobacter pylori (H. pylori) infection can promote the development of gastric cancer (GC); however, the underlying mechanism is not clear. FAM60A has been found showing high levels in some cancer cells, including lung cancer (A549), and pancreatic cancer (Capan-2) cell lines. Data in oncomine showed that FAM60A overexpression was an critical prognostic factor in GC. In this study, we showed that knockdown of FAM60A could revert the increase of proliferation and the decrease of apoptosis caused by H.pylori infection in HGC-27 and AGS cells. Conversely, FAM60A upregulation promoted proliferation and inhibited apoptosis in HGC-27 and AGS cells. We also found that the PI3K/AKT pathway inhibitor LY294002 could revert the changes caused by FAM60A upregulation in HGC-27 and AGS cells. Thus, our study provides evidence that FAM60A act as a carcinogen and suggests that H. pylori-induced upregulation of FAM60A may contribute to the development of gastric cancer.

Copy number variation is highly correlated with differential gene expression: a pan-cancer study

BACKGROUND

Cancer is a heterogeneous disease with many genetic variations. Lines of evidence have shown copy number variations (CNVs) of certain genes are involved in development and progression of many cancers through the alterations of their gene expression levels on individual or several cancer types. However, it is not quite clear whether the correlation will be a general phenomenon across multiple cancer types.

METHODS

In this study we applied a bioinformatics approach integrating CNV and differential gene expression mathematically across 1025 cell lines and 9159 patient samples to detect their potential relationship.

RESULTS

Our results showed there is a close correlation between CNV and differential gene expression and the copy number displayed a positive linear influence on gene expression for the majority of genes, indicating that genetic variation generated a direct effect on gene transcriptional level. Another independent dataset is utilized to revalidate the relationship between copy number and expression level. Further analysis show genes with general positive linear influence on gene expression are clustered in certain disease-related pathways, which suggests the involvement of CNV in pathophysiology of diseases.

CONCLUSIONS

This study shows the close correlation between CNV and differential gene expression revealing the qualitative relationship between genetic variation and its downstream effect, especially for oncogenes and tumor suppressor genes. It is of a critical importance to elucidate the relationship between copy number variation and gene expression for prevention, diagnosis and treatment of cancer.

Annexin A8 can serve as potential prognostic biomarker and therapeutic target for ovarian cancer: based on the comprehensive analysis of Annexins

BACKGROUND

Annexins are involved in vesicle trafficking, cell proliferation and apoptosis, but their functional mechanisms in ovarian cancer remain unclear. In this study, we analyzed Annexins in ovarian cancer using different databases and selected Annexin A8 (ANXA8), which showed the greatest prognostic value, for subsequent validation in immunohistochemical (IHC) assays.

METHODS

The mRNA expression levels, genetic variations, prognostic values and gene-gene interaction network of Annexins in ovarian cancer were analyzed using the Oncomine, Gene Expression Profiling Interactive Analysis (GEPIA), cBioPortal, Kaplan-Meier plotter and GeneMANIA database. ANXA8 was selected for analyzing the biological functions and pathways of its co-expressed genes, and its correlation with immune system responses via the Database for Annotation, Visualization, and Integrated Discovery (DAVID) and the TISIDB database, respectively. We validated the expression of ANXA8 in ovarian cancer via IHC assays and analyzed its correlation with clinicopathological parameters and prognosis.

RESULTS

ANXA2/3/8/11 mRNA expression levels were significantly upregulated in ovarian cancer, and ANXA5/6/7 mRNA expression levels were significantly downregulated. Prognostic analysis suggested that significant correlations occurred between ANXA2/4/8/9 mRNA upregulation and poor overall survival, and between ANXA8/9/11 mRNA upregulation and poor progression-free survival in patients with ovarian serous tumors. Taken together, results suggested that ANXA8 was most closely associated with ovarian cancer tumorigenesis and progression. Further analyses indicated that ANXA8 may be involved in cell migration, cell adhesion, and vasculature development, as well as in the regulation of PI3K-Akt, focal adhesion, and proteoglycans. Additionally, ANXA8 expression was significantly correlated with lymphocytes and immunomodulators. The IHC results showed that ANXA8 expression was higher in the malignant tumor group than in the borderline and benign tumor groups and normal ovary group, and high ANXA8 expression was an independent risk factor for survival and prognosis of ovarian cancer patients (P = 0.013).

CONCLUSIONS

Members of the Annexin family display varying degrees of abnormal expressions in ovarian cancer. ANXA8 was significantly highly expressed in ovarian cancer, and high ANXA8 expression was significantly correlated with poor prognosis. Therefore, ANXA8 is a high candidate as a novel biomarker and therapeutic target for ovarian cancer.

Genomic data analysis workflows for tumors from patient-derived xenografts (PDXs): challenges and guidelines

BACKGROUND

Patient-derived xenograft (PDX) models are in vivo models of human cancer that have been used for translational cancer research and therapy selection for individual patients. The Jackson Laboratory (JAX) PDX resource comprises 455 models originating from 34 different primary sites (as of 05/08/2019). The models undergo rigorous quality control and are genomically characterized to identify somatic mutations, copy number alterations, and transcriptional profiles. Bioinformatics workflows for analyzing genomic data obtained from human tumors engrafted in a mouse host (i.e., Patient-Derived Xenografts; PDXs) must address challenges such as discriminating between mouse and human sequence reads and accurately identifying somatic mutations and copy number alterations when paired non-tumor DNA from the patient is not available for comparison.

RESULTS

We report here data analysis workflows and guidelines that address these challenges and achieve reliable identification of somatic mutations, copy number alterations, and transcriptomic profiles of tumors from PDX models that lack genomic data from paired non-tumor tissue for comparison. Our workflows incorporate commonly used software and public databases but are tailored to address the specific challenges of PDX genomics data analysis through parameter tuning and customized data filters and result in improved accuracy for the detection of somatic alterations in PDX models. We also report a gene expression-based classifier that can identify EBV-transformed tumors. We validated our analytical approaches using data simulations and demonstrated the overall concordance of the genomic properties of xenograft tumors with data from primary human tumors in The Cancer Genome Atlas (TCGA).

CONCLUSIONS

The analysis workflows that we have developed to accurately predict somatic profiles of tumors from PDX models that lack normal tissue for comparison enable the identification of the key oncogenic genomic and expression signatures to support model selection and/or biomarker development in therapeutic studies. A reference implementation of our analysis recommendations is available at https://github.com/TheJacksonLaboratory/PDX-Analysis-Workflows .

SINHCAF/FAM60A and SIN3A specifically repress HIF-2α expression

The SIN3A-HDAC (histone deacetylase) complex is a master transcriptional repressor, required for development but often deregulated in disease. Here, we report that the recently identified new component of this complex, SINHCAF (SIN3A and HDAC-associated factor)/FAM60A (family of homology 60A), links the SIN3A-HDAC co-repressor complex function to the hypoxia response. We show that SINHCAF specifically represses HIF-2α mRNA and protein expression, via its interaction with the transcription factor SP1 (specificity protein 1) and recruitment of HDAC1 to the HIF-2α promoter. SINHCAF control over HIF-2α results in functional cellular changes in in vitro angiogenesis and viability. Our analysis reveals an unexpected link between SINHCAF and the regulation of the hypoxia response.

A five-gene based risk score with high prognostic value in colorectal cancer

Colorectal cancer (CRC) is one of the most frequently occurring malignancies worldwide. The outcomes of patients with similar clinical symptoms or at similar pathological stages remain unpredictable. This inherent clinical diversity is most likely due to the genetic heterogeneity. The present study aimed to create a predicting tool to evaluate patient survival based on genetic profile. Firstly, three Gene Expression Omnibus (GEO) datasets (GSE9348, GSE44076 and GSE44861) were utilized to identify and validate differentially expressed genes (DEGs) in CRC. The GSE14333 dataset containing survival information was then introduced in order to screen and verify prognosis-associated genes. Of the 66 DEGs, the present study screened out 46 biomarkers closely associated to patient overall survival. By Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis, it was demonstrated that these genes participated in multiple biological processes which were highly associated with cancer proliferation, drug-resistance and metastasis, thus further affecting patient survival. The five most important genes, MET proto-oncogene, receptor tyrosine kinase, carboxypeptidase M, serine hydroxymethyltransferase 2, guanylate cyclase activator 2B and sodium voltage-gated channel a subunit 9 were selected by a random survival forests algorithm, and were further made up to a linear risk score formula by multivariable cox regression. Finally, the present study tested and verified this risk score within three independent GEO datasets (GSE14333, GSE17536 and GSE29621), and observed that patients with a high risk score had a lower overall survival (P<0.05). Furthermore, this risk score was the most significant compared with other predicting factors including age and American Joint Committee on Cancer stage, in the model, and was able to predict patient survival independently and directly. The findings suggest that this survival associated DEGs-based risk score is a powerful and accurate prognostic tool and is promisingly implemented in a clinical setting.

Microarray analyses reveal genes related to progression and prognosis of esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma is a high morbidity and mortality cancer in China. Here are few biomarkers and therapeutic targets. Our study was aimed to identify candidate genes correlated to ESCC. Oncomine, The Cancer Genome Atlas, Gene Expression Omnibus were retrieved for eligible ESCC data. Deregulated genes were identified by meta-analysis and validated by an independent dataset. Survival analyses and bioinformatics analyses were used to explore potential mechanisms. Copy number variant analyses identified upstream mechanisms of candidate genes. In our study, top 200 up/down-regulated genes were identified across two microarrays. A total of 139 different expression genes were validated in GSE53625. Survival analysis found that nine genes were closely related to prognosis. Furthermore, Gene Ontology analyses and Kyoto Encyclopedia of Genes and Genomes analyses showed that different expression genes were mainly enriched in cell division, cell cycle and cell-cell adhesion pathways. Copy number variant analyses indicated that overexpression of ECT2 and other five genes were correlated with copy number amplification. The current study demonstrated that ECT2 and other eight candidate genes were correlated to progression and prognosis of esophageal squamous cell carcinoma, which might provide novel insights to the mechanisms.