Increased AURKA promotes cell proliferation and predicts poor prognosis in bladder cancer

Background

Bladder cancer (BC) is the most common cancer of the urinary bladder and upper tract, in which the clinical management is limited. AURKA (aurora kinase A) has been identified as an oncogene in cancer development; however, its potential role and underlying mechanisms in the progression of BC remain unknown.

Results

In this study, we evaluated Aurora kinase A (AURKA) expression in patient samples by performing gene expression profiling, and found that AURKA expression levels were significantly higher in BC tissues than in normal tissues. Increased AURKA in BC was strongly associated with stage and grade. Moreover, BC patients with elevated AURKA achieved poor overall survival rates. The experiments in vitro comprehensively validated the critical role of AURKA in promoting BC cell proliferation using the methods of gene overexpression and gene silencing. Furthermore, we proved that AURKA inhibitor MLN8237 arrested BC cell growth and induced apoptosis.

Conclusions

These findings implicate AURKA acting as an effective biomarker for BC detection and prognosis, as well as therapeutic target.

     

Identification of cancer subtypes from single-cell RNA-seq data using a consensus clustering method

Background

Human cancers are complex ecosystems composed of cells with distinct molecular signatures. Such intratumoral heterogeneity poses a major challenge to cancer diagnosis and treatment. Recent advancements of single-cell techniques such as scRNA-seq have brought unprecedented insights into cellular heterogeneity. Subsequently, a challenging computational problem is to cluster high dimensional noisy datasets with substantially fewer cells than the number of genes.

Methods

In this paper, we introduced a consensus clustering framework conCluster, for cancer subtype identification from single-cell RNA-seq data. Using an ensemble strategy, conCluster fuses multiple basic partitions to consensus clusters.

Results

Applied to real cancer scRNA-seq datasets, conCluster can more accurately detect cancer subtypes than the widely used scRNA-seq clustering methods. Further, we conducted co-expression network analysis for the identified melanoma subtypes.

Conclusions

Our analysis demonstrates that these subtypes exhibit distinct gene co-expression networks and significant gene sets with different functional enrichment.

     

Long non-coding RNA DBCCR1-003 regulate the expression of DBCCR1 via DNMT1 in bladder cancer

Background

Many long non coding RNAs have been identified as key modulators in cancer development. A lncRNA, DBCCR1-003, derived from the locus of tumor suppressor gene DBCCR1 (deleted in bladder cancer chromosome region 1), has unknown function. In the present study, we explored function and molecular mechanism of DBCCR1-003 in bladder cancer (BC) development.

Methods

We evaluated the expression levels of DBCCR1-003 in tissues and cells with western blot and quantitative real-time polymerase chain reaction. Multiple approaches including chromatin immunoprecipitation assay and RNA immunoprecipitation were used to confirm the direct binding of DBCCR1-003 to DNMT1. The recombinant vector overexpressing DBCCR1-003 was constructed. Cell proliferation assay, colony formation assay and flow cytometric analysis were employed to measure the role of DBCCR1-003 in regulation of cell proliferation, cycle and apoptosis.

Results

Firstly we detected the expression of DBCCR1-003, DBCCR1, DNMT1 (DNA methyltransferase 1) and DNA methylation in the promoter of DBCCR1. We found low expression of DBCCR1-003, same as DBCCR1, while high expression of DNMT1 and hypermethylation of DBCCR1 gene promoter in BC tissues and T24 cells line. Further studies revealed that treatment of DNMT inhibitor, 5-aza-2-deoxycytidine(DAC), or overexpression of DBCCR1-003 led to increased DBCCR1 expression by reversion of promoter hypermethylation and DNMT1 binding to DBCCR1 promoter in T24 cells. Importantly, RNA immunoprecipitation (RIP) showed that DBCCR1-003 physically associates with DNMT1. The binding of them was increased with the inhibition of DBCCR1 promoter methylation, indicating that DBCCR1-003 may bind to DNMT1 and prevent DNMT1-mediated the methylation of DBCCR1. Furthermore, overexpression of DBCCR1-003 resulted in significant inhibition of T24 cells growth through the inducing G0/G1 arrest and apoptosis.

Conclusions

Taken together, these findings demonstrated that a novel tumor suppressor DBCCR1-003 regulates the expression of DBCCR1 via binding to DNMT1 and preventing DNMT1-mediated the methylation of DBCCR1 in BC. LncRNA DBCCR1-003 may serve as a novel biomarker and therapeutic target for BC in future cancer clinic.

     

Analysis of metabolites and metabolic pathways in breast cancer in a Korean prospective cohort: the Korean Cancer Prevention Study-II

Introduction

Since blood is in contact with all tissues in the body and is considered to dynamically reflect the body’s pathophysiological status, serum metabolomics changes are important and have diagnostic value in early cancer detection.

Objectives

In this prospective study, we investigated the application of metabolomics to differentiate subjects with incident breast cancer (BC) from subjects who remained free of cancer during a mean follow-up period of 7 years with the aim of identifying valuable biomarkers for BC.

Methods

Baseline serum samples from 84 female subjects with incident BC (BC group) and 88 cancer-free female subjects (control group) were used. Metabolic alterations associated with BC were investigated via metabolomics analysis of the baseline serum samples using ultra-performance liquid chromatography-linear-trap quadrupole-Orbitrap mass spectrometry.

Results

A total of 57 metabolites were identified through the metabolic analysis. Among them, 20 metabolite levels were significantly higher and 22 metabolite levels were significantly lower in the BC group than in the control group at baseline. Ten metabolic pathways, including amino acid metabolism, arachidonic acid (AA) metabolism, fatty acid metabolism, linoleic acid metabolism, and retinol metabolism, showed significant differences between the BC group and the control group. Logistic regression revealed that the incidence of BC was affected by leucine, AA, prostaglandin (PG)J₂, PGE₂, and γ-linolenic acid (GLA).

Conclusions

This prospective study showed the clinical relevance of dysregulation of various metabolisms on the incidence of BC. Additionally, leucine, AA, PGJ₂, PGE₂, and GLA were identified as independent variables affecting the incidence of BC.

     

Statistical analysis of fractionation resistance by functional category and expression

Background

The current literature establishes the importance of gene functional category and expression in promoting or suppressing duplicate gene loss after whole genome doubling in plants, a process known as fractionation. Inspired by studies that have reported gene expression to be the dominating factor in preventing duplicate gene loss, we analyzed the relative effect of functional category and expression.

Methods

We use multivariate methods to study data sets on gene retention, function and expression in rosids and asterids to estimate effects and assess their interaction.

Results

Our results suggest that the effect on duplicate gene retention fractionation by functional category and expression are independent and have no statistical interaction.

Conclusion

In plants, functional category is the more dominant factor in explaining duplicate gene loss.

     

Optimizing gene set annotations combining GO structure and gene expression data

Background

With the rapid accumulation of genomic data, it has become a challenge issue to annotate and interpret these data. As a representative, Gene set enrichment analysis has been widely used to interpret large molecular datasets generated by biological experiments. The result of gene set enrichment analysis heavily relies on the quality and integrity of gene set annotations. Although several methods were developed to annotate gene sets, there is still a lack of high quality annotation methods. Here, we propose a novel method to improve the annotation accuracy through combining the GO structure and gene expression data.

Results

We propose a novel approach for optimizing gene set annotations to get more accurate annotation results. The proposed method filters the inconsistent annotations using GO structure information and probabilistic gene set clusters calculated by a range of cluster sizes over multiple bootstrap resampled datasets. The proposed method is employed to analyze p53 cell lines, colon cancer and breast cancer gene expression data. The experimental results show that the proposed method can filter a number of annotations unrelated to experimental data and increase gene set enrichment power and decrease the inconsistent of annotations.

Conclusions

A novel gene set annotation optimization approach is proposed to improve the quality of gene annotations. Experimental results indicate that the proposed method effectively improves gene set annotation quality based on the GO structure and gene expression data.

     

SOX2 as a novel contributor of oxidative metabolism in melanoma cells

Background

Deregulated metabolism is a hallmark of cancer and recent evidence underlines that targeting tumor energetics may improve therapy response and patient outcome. Despite the general attitude of cancer cells to exploit the glycolytic pathway even in the presence of oxygen (aerobic glycolysis or “Warburg effect”), tumor metabolism is extremely plastic, and such ability to switch from glycolysis to oxidative phosphorylation (OxPhos) allows cancer cells to survive under hostile microenvironments. Recently, OxPhos has been related with malignant progression, chemo-resistance and metastasis. OxPhos is induced under extracellular acidosis, a well-known characteristic of most solid tumors, included melanoma.

Methods

To evaluate whether SOX2 modulation is correlated with metabolic changes under standard or acidic conditions, SOX2 was silenced and overexpressed in several melanoma cell lines. To demonstrate that SOX2 directly represses HIF1A expression we used chromatin immunoprecipitation (ChIP) and luciferase assay.

Results

In A375-M6 melanoma cells, extracellular acidosis increases SOX2 expression, that sustains the oxidative cancer metabolism exploited under acidic conditions. By studying non-acidic SSM2c and 501-Mel melanoma cells (high- and very low-SOX2 expressing cells, respectively), we confirmed the metabolic role of SOX2, attributing SOX2-driven OxPhos reprogramming to HIF1α pathway disruption.

Conclusions

SOX2 contributes to the acquisition of an aggressive oxidative tumor phenotype, endowed with enhanced drug resistance and metastatic ability.

     

Chromosome 9p copy number gains involving PD-L1 are associated with a specific proliferation and immune-modulating gene expression program active across major cancer types

Background

Inhibition of the PD-L1/PD-1 immune checkpoint axis represents one of the most promising approaches of immunotherapy for various cancer types. However, immune checkpoint inhibition is successful only in subpopulations of patients emphasizing the need for powerful biomarkers that adequately reflect the complex interaction between the tumor and the immune system. Recently, recurrent copy number gains (CNG) in chromosome 9p involving PD-L1 were detected in many cancer types including lung cancer, melanoma, bladder cancer, head and neck cancer, cervical cancer, soft tissue sarcoma, prostate cancer, gastric cancer, ovarian cancer, and triple-negative breast cancer.

Methods

Here, we applied functional genomics to analyze global mRNA expression changes associated with chromosome 9p gains. Using the TCGA data set, we identified a list of 75 genes that were strongly up-regulated in tumors with chromosome 9p gains across many cancer types.

Results

As expected, the gene set was enriched for chromosome 9p and in particular chromosome 9p24 (36 genes and 23 genes). Furthermore, we found enrichment of two expression programs derived from genes within and beyond 9p: one implicated in cell cycle regulation (22 genes) and the other implicated in modulation of the immune system (16 genes). Among these were specific cytokines and chemokines, e.g. CCL4, CCL8, CXCL10, CXCL11, other immunoregulatory genes such as IFN-G and IDO1 as well as highly expressed proliferation-related kinases and genes including PLK1, TTK, MELK and CDC20 that represent potential drug targets.

Conclusions

Collectively, these data shed light on mechanisms of immune escape and stimulation of proliferation in cancer with PD-L1 CNG and highlight additional vulnerabilities that may be therapeutically exploitable.

     

Analysis of PD-1 related immune transcriptional profile in different cancer types

Background

Programmed cell death 1 (PD-1) functions as an immune checkpoint in the process of anti-tumor immune response. The PD-1 blockade is now becoming a fundamental part in cancer immunotherapy. So it’s essential to elicit the PD-1 related immune process in different types of cancer.

Methods

The Cancer Genome Atlas was used to collect the RNA-seq data of 33 cancer types. The microenvironment cell populations-counter was used to analyze the immune cell infiltrates. KEGG and GO analysis were performed to investigate PD-1 associated biological process. Kaplan–Meier survival curves and Cox’s proportional hazards model were performed for prognostic value analysis.

Results

We demonstrated that PD-1 expression varied in different cancer types. The uveal melanoma had a low PD-1 expression and poor infiltrated with immune cells. But it showed the strong correlation of PD-1 with the most types of immune cells. The PD-1 demonstrated a robust relationship with other immunomodulators and showed its involvement in critical functions correlated with anti-tumor immune pathways. Survival analysis indicated the PD-1 expression suggested different prognosis in different cancer types.

Conclusions

Our investigations promote a better understanding of the PD-1 blockade and provide PD-1 related personized combined immunotherapy for different types of cancer patients.

     

Integrative pathway-based survival prediction utilizing the interaction between gene expression and DNA methylation in breast cancer

Background

Integrative analysis on multi-omics data has gained much attention recently. To investigate the interactive effect of gene expression and DNA methylation on cancer, we propose a directed random walk-based approach on an integrated gene-gene graph that is guided by pathway information.

Methods

Our approach first extracts a single pathway profile matrix out of the gene expression and DNA methylation data by performing the random walk over the integrated graph. We then apply a denoising autoencoder to the pathway profile to further identify important pathway features and genes. The extracted features are validated in the survival prediction task for breast cancer patients.

Results

The results show that the proposed method substantially improves the survival prediction performance compared to that of other pathway-based prediction methods, revealing that the combined effect of gene expression and methylation data is well reflected in the integrated gene-gene graph combined with pathway information. Furthermore, we show that our joint analysis on the methylation features and gene expression profile identifies cancer-specific pathways with genes related to breast cancer.

Conclusions

In this study, we proposed a DRW-based method on an integrated gene-gene graph with expression and methylation profiles in order to utilize the interactions between them. The results showed that the constructed integrated gene-gene graph can successfully reflect the combined effect of methylation features on gene expression profiles. We also found that the selected features by DA can effectively extract topologically important pathways and genes specifically related to breast cancer.

     

Integrated analysis of microRNA-target interactions with clinical outcomes for cancers

Background

Clinical statement alone is not enough to predict the progression of disease. Instead, the gene expression profiles have been widely used to forecast clinical outcomes. Many genes related to survival have been identified, and recently miRNA expression signatures predicting patient survival have been also investigated for several cancers. However, miRNAs and their target genes associated with clinical outcomes have remained largely unexplored.

Methods

Here, we demonstrate a survival analysis based on the regulatory relationships of miRNAs and their target genes. The patient survivals for the two major cancers, ovarian cancer and glioblastoma multiforme (GBM), are investigated through the integrated analysis of miRNA-mRNA interaction pairs.

Results

We found that there is a larger survival difference between two patient groups with an inversely correlated expression profile of miRNA and mRNA. It supports the idea that signatures of miRNAs and their targets related to cancer progression can be detected via this approach.

Conclusions

This integrated analysis can help to discover coordinated expression signatures of miRNAs and their target mRNAs that can be employed for therapeutics in human cancers.

     

Volatile metabolomic signature of bladder cancer cell lines based on gas chromatography–mass spectrometry

Introduction

Recent studies provide a convincing support that the presence of cancer cells in the body leads to the alteration of volatile organic compounds (VOCs) emanating from biological samples, particularly of those closely related with tumoral tissues. Thus, a great interest emerged for the study of cancer volatilome and subsequent attempts to confirm VOCs as potential diagnostic biomarkers.

Objectives

The aim of this study was to determine the volatile metabolomic signature of bladder cancer (BC) cell lines and provide an in vitro proof-of-principle that VOCs emanated into the extracellular medium may discriminate BC cells from normal bladder epithelial cells.

Methods

VOCs in the culture media of three BC cell lines (Scaber, J82, 5637) and one normal bladder cell line (SV-HUC-1) were extracted by headspace-solid phase microextraction and analysed by gas chromatography-mass spectrometry (HS-SPME/GC–MS). Two different pH (pH 2 and 7) were used for VOCs extraction to infer the best pH to be used in in vitro metabolomic studies.

Results

Multivariate analysis revealed a panel of volatile metabolites that discriminated cancerous from normal bladder cells, at both pHs, although a higher number of discriminative VOCs was obtained at neutral pH. Most of the altered metabolites were ketones and alkanes, which were generally increased in BC compared to normal cells, and alcohols, which were significantly decreased in BC cells. Among them, three metabolites, namely 2-pentadecanone, dodecanal and γ-dodecalactone (the latter only tentatively identified), stood out as particularly important metabolites and promising volatile biomarkers for BC detection. Furthermore, our results also showed the potential of VOCs in discriminating BC cell lines according to tumour grade and histological subtype.

Conclusions

We demonstrate that a GC–MS metabolomics-based approach for analysis of VOCs is a valuable strategy for identifying new and specific biomarkers that may improve BC diagnosis. Future studies should entail the validation of volatile signature found for BC cell lines in biofluids from BC patients.

     

Characterizing biomarkers in osteosarcoma metastasis based on an ego-network

Objectives

To characterize biomarkers that underlie osteosarcoma (OS) metastasis based on an ego-network.

Results

From the microarray data, we obtained 13,326 genes. By combining PPI data and microarray data, 10,520 shared genes were found and constructed into ego-networks. 17 significant ego-networks were identified with p < 0.05. In the pathway enrichment analysis, seven ego-networks were identified with the most significant pathway.

Conclusions

These significant ego-modules were potential biomarkers that reveal the potential mechanisms in OS metastasis, which may contribute to understanding cancer prognoses and providing new perspectives in the treatment of cancer.

     

Study of the tumor microenvironment during breast cancer progression

Background

Different cells and mediators in the tumor microenvironment play important roles in the progression of breast cancer. The aim of this study was to determine the composition of the microenvironment during tumor progression in order to discover new related biomarkers and potentials for targeted therapy.

Methods

In this study, breast cancer biopsies from four different stages, and control breast biopsies were collected. Then, the mRNA expression of several markers related to different CD4⁺ T cell subsets including regulatory T cells (Treg), T helper (Th) type 1, 2 and 17 were determined. In addition, we investigated the expression of two inflammatory cytokines (TNF-α and IL-6) and inflammatory mediators including FASL, IDO, SOCS1, VEGF, and CCR7.

Results

The results showed that the expression of Th1 and Th17 genes was decreased in tumor tissues compared to control tissues. In addition, we found that the gene expression related to these two cell subsets decreased during cancer progression. Moreover, the expression level of TNF-α increased with tumor progression.

Conclusion

We conclude that the expression of genes related to immune response and inflammation is different between tumor tissues and control tissues. In addition, this difference was perpetuated through the different stages of cancer.

     

Weighted set enrichment of gene expression data

Background

Sets of genes that are known to be associated with each other can be used to interpret microarray data. This gene set approach to microarray data analysis can illustrate patterns of gene expression which may be more informative than analyzing the expression of individual genes. Various statistical approaches exist for the analysis of gene sets. There are three main classes of these methods: over-representation analysis, functional class scoring, and pathway topology based methods.

Methods

We propose weighted hypergeometric and weighted chi-squared methods in order to assign a rank to the degree to which each gene participates in the enrichment. Each gene is assigned a weight determined by the absolute value of its log fold change, which is then raised to a certain power. The power value can be adjusted as needed. Datasets from the Gene Expression Omnibus are used to test the method. The significantly enriched pathways are validated through searching the literature in order to determine their relevance to the dataset.

Results

Although these methods detect fewer significantly enriched pathways, they can potentially produce more relevant results. Furthermore, we compare the results of different enrichment methods on a set of microarray studies all containing data from various rodent neuropathic pain models.

Discussion

Our method is able to produce more consistent results than other methods when evaluated on similar datasets. It can also potentially detect relevant pathways that are not identified by the standard methods. However, the lack of biological ground truth makes validating the method difficult.

     

Stress amplifies sex differences in primate prefrontal profiles of gene expression

Background

Stress is a recognized risk factor for mood and anxiety disorders that occur more often in women than men. Prefrontal brain regions mediate stress coping, cognitive control, and emotion. Here, we investigate sex differences and stress effects on prefrontal cortical profiles of gene expression in squirrel monkey adults.

Methods

Dorsolateral, ventrolateral, and ventromedial prefrontal cortical regions from 18 females and 12 males were collected after stress or no-stress treatment conditions. Gene expression profiles were acquired using HumanHT-12v4.0 Expression BeadChip arrays adapted for squirrel monkeys.

Results

Extensive variation between prefrontal cortical regions was discerned in the expression of numerous autosomal and sex chromosome genes. Robust sex differences were also identified across prefrontal cortical regions in the expression of mostly autosomal genes. Genes with increased expression in females compared to males were overrepresented in mitogen-activated protein kinase and neurotrophin signaling pathways. Many fewer genes with increased expression in males compared to females were discerned, and no molecular pathways were identified. Effect sizes for sex differences were greater in stress compared to no-stress conditions for ventromedial and ventrolateral prefrontal cortical regions but not dorsolateral prefrontal cortex.

Conclusions

Stress amplifies sex differences in gene expression profiles for prefrontal cortical regions involved in stress coping and emotion regulation. Results suggest molecular targets for new treatments of stress disorders in human mental health.

     

Comparisons of gene coexpression network modules in breast cancer and ovarian cancer

Background

Breast cancer and ovarian cancer are hormone driven and are known to have some predisposition genes in common such as the two well known cancer genes BRCA1 and BRCA2. The objective of this study is to compare the coexpression network modules of both cancers, so as to infer the potential cancer-related modules.

Methods

We applied the eigen-decomposition to the matrix that integrates the gene coexpression networks of both breast cancer and ovarian cancer. With hierarchical clustering of the related eigenvectors, we obtained the network modules of both cancers simultaneously. Enrichment analysis on Gene Ontology (GO), KEGG pathway, Disease Ontology (DO), and Gene Set Enrichment Analysis (GSEA) in the identified modules was performed.

Results

We identified 43 modules that are enriched by at least one of the four types of enrichments. 31, 25, and 18 modules are enriched by GO terms, KEGG pathways, and DO terms, respectively. The structure of 29 modules in both cancers is significantly different with p-values less than 0.05, of which 25 modules have larger densities in ovarian cancer. One module was found to be significantly enriched by the terms related to breast cancer from GO, KEGG and DO enrichment. One module was found to be significantly enriched by ovarian cancer related terms.

Conclusion

Breast cancer and ovarian cancer share some common properties on the module level. Integration of both cancers helps identifying the potential cancer associated modules.

     

Bone morphogenetic protein 2 (BMP2) induces growth suppression and enhances chemosensitivity of human colon cancer cells

Background

Molecular profiling of colorectal cancer (CRC) based on global gene expression has revealed multiple dysregulated signalling pathways associated with drug resistance and poor prognosis. However, the role of BMP2 signaling in CRC is not fully characterised.

Methods

Bioinformatics data analysis were conducted on the GSE21510 dataset. Leniviral technology was utilized to stably express BMP2 in the HCT116 CRC model. Gene expression profiling was conducted using Agilent microarray platform while data normalization and bioinformatics were conducted using GeneSpring software. Changes in gene expression were assessed using qRT-PCR. AlamarBlue assay was used to assess cell viability in vitro. In vivo experiments were conducted using SCID mice.

Results

Our data revealed frequent downregulation of BMP2 in primary CRC tissues. Additionally, interrogation of publically available gene expression datasets revealed significant downregulation of BMP2 in metastatic recurrent compared to non-metastatic cancer (p = 0.02). Global gene expression analysis in CRC cells over-expressing BMP2 revealed multiple dysregulated pathways mostly affecting cell cycle and DNA damage response. Concordantly, lentiviral-mediated re-expression of BMP2 inhibited HCT116 CRC growth, sphere formation, clonogenic potential, cell migration, and sensitized CRC cells to 5-fluorouracil (5-FU) in vitro. Additionally, BMP2 inhibited CRC tumor formation in SCID mice.

Conclusions

Our data revealed an inhibitory role for BMP2 in CRC, suggesting that restoration of BMP2 expression could be a potential therapeutic strategy for CRC.

     

Metabolic response of porcine colon explants to in vitro infection by <Emphasis Type="Italic">Brachyspira hyodysenteriae</Emphasis>: a leap into disease pathophysiology

Introduction

Swine dysentery caused by Brachyspira hyodysenteriae is a production limiting disease in pig farming. Currently antimicrobial therapy is the only treatment and control method available.

Objective

The aim of this study was to characterize the metabolic response of porcine colon explants to infection by B. hyodysenteriae.

Methods

Porcine colon explants exposed to B. hyodysenteriae were analyzed for histopathological, metabolic and pro-inflammatory gene expression changes.

Results

Significant epithelial necrosis, increased levels of l-citrulline and IL-1α were observed on explants infected with B. hyodysenteriae.

Conclusions

The spirochete induces necrosis in vitro likely through an inflammatory process mediated by IL-1α and NO.