Urothelial cancer gene regulatory networks inferred from large-scale RNAseq, Bead and Oligo gene expression data

TRPV1 is a potential biomarker for the prediction and treatment of multiple cancers based on a pan-cancer analysis

BACKGROUND

Transient receptor potential cation channel subfamily V member 1 (TRPV1) was considered to play pivotal roles in multiple cancers; however, the expression and clinical significance of the TRPV1 remain unclear, which were explored in this study.

RESULTS

The pan-cancer analysis was performed based on 10,236 samples in 32 cancers. Differential TRPV1 expression levels were detected in 12 cancers (p < 0.05). TRPV1 demonstrated its conspicuous prognosis significance and prediction effects for some cancers (e.g., lung adenocarcinoma), indicating its potential as a valuable and novel biomarker in treating and predicting cancers. TRPV1 expression was relevant to DNA methyltransferases, mismatch repair genes, tumor mutational burden, and microsatellite instability. TRPV1 expression was associated with the immune microenvironment of some cancers, and its roles in different cancers may be mediated by affecting various immune cells. Gene set enrichment analysis discloses the significant relevance of TRPV1 expression with a series of metabolic and immunoregulatory-related pathways.

CONCLUSIONS

This study provided a comprehensive workflow of the expression, clinical significance, and underlying mechanisms of TRPV1 in pan-cancer. TRPV1 may be an underlying biomarker for predicting and treating multiple cancer.

The Performance Comparison of Gene Co-expression Networks of Breast and Prostate Cancer using Different Selection Criteria

Gene co-expression networks (GCN) present undirected relations between genes to understand molecular structures behind the diseases, including cancer. The utilization of various biological datasets and gene network inference (GNI) algorithms can reveal meaningful gene-gene interactions of GCNs. This study applies three GNI algorithms on mRNA gene expression, RNA-Seq, and miRNA-target genes datasets to infer GCNs of breast and prostate cancers. To evaluate the performance of the GCNs, we utilize overlap analysis via literature data, topological assessment, and Gene Ontology-based biological assessment. The results emphasize how the selection of biological datasets and GNI algorithms affect the performance results on different evaluation criteria. GCNs on microarray gene expression data slightly outperform in overlap analysis. Also, GCNs on RNA-Seq and gene expression datasets follow scale-free topology. The biological assessment results are close to each other on all biological datasets. C3NET algorithm-based GCNs did not contain any biological assessment modules; therefore, it is not optimal for biological assessment. GNI algorithms' selection did not change the overlap analysis and topological assessment results. Our primary objective is to compare the performance results of biological datasets and GNI algorithms based on different evaluation criteria. For this purpose, we developed the GNIAP R package that enables users to select different GNI algorithms to infer GCNs. The GNIAP R package also provides literature-based overlap analysis, and topological and biological analyses on GCNs. Users can access the GNIAP R package via https://github.com/ozgurcingiz/GNIAP .

RNAseq expression patterns of canine invasive urothelial carcinoma reveal two distinct tumor clusters and shared regions of dysregulation with human bladder tumors

BACKGROUND

Invasive urothelial carcinoma (iUC) is highly similar between dogs and humans in terms of pathologic presentation, molecular subtypes, response to treatment and age at onset. Thus, the dog is an established and relevant model for testing and development of targeted drugs benefiting both canine and human patients. We sought to identify gene expression patterns associated with two primary types of canine iUC tumors: those that express a common somatic mutation in the BRAF gene, and those that do not.

METHODS

We performed RNAseq on tumor and normal tissues from pet dogs. Analysis of differential expression and clustering, and positional and individual expression was used to develop gene set enrichment profiles distinguishing iUC tumors with and without BRAFV595E mutations, as well as genomic regions harboring excessive numbers of dysregulated genes.

RESULTS

We identified two expression clusters that are defined by the presence/absence of a BRAFV595E (BRAFV600E in humans) somatic mutation. BRAFV595E tumors shared significantly more dysregulated genes than BRAF wild-type tumors, and vice versa, with 398 genes differentiating the two clusters. Key genes fall into clades of limited function: tissue development, cell cycle regulation, immune response, and membrane transport. The genomic site with highest number of dysregulated genes overall lies in a locus corresponding to human chromosome 8q24, a region frequently amplified in human urothelial cancers.

CONCLUSIONS

These data identify critical sets of genes that are differently regulated in association with an activating mutation in the MAPK/ERK pathway in canine iUC tumors. The experiments also highlight the value of the canine system in identifying expression patterns associated with a common, shared cancer.

Two-tier combinatorial structure to integrate various gene co-expression networks of prostate cancer

Advances in DNA sequencing technologies enable researchers to integrate various biological datasets in order to reveal hidden relations at the molecular level. In this study, we present a two-tiered combinatorial structure (TTCS) to integrate gene co-expression networks (GCNs) that are inferred from microarray gene expression, RNA-Seq and miRNA-target gene data. In the initial phase of TTCS, we derive GCNs by using gene network inference (GNI) algorithms for each dataset. In the first and second integration phases, we use straightforward methods: intersection, union and simple majority voting to combine GCNs. We use overlap, topological and biological analyses in performance evaluation and investigate the integration effects of GCNs separately for all phases. Our results prove that the first integration phase has limited contribution on performance. However, combining the biological datasets in the second phase significantly enhances the overlap and topological performance analyses.

Identification of Key Genes and Pathways Associated with RUNX1 Mutations in Acute Myeloid Leukemia Using Bioinformatics Analysis

BACKGROUND RUNXl plays a key regulatory role in the process of hematopoiesis and is a common target for multiple chromosomal translocations in human acute leukemia. Mutations of RUNX1 gene can lead to acute leukemia and affect the prognosis of AML patients. We aimed to identify pivotal genes and pathways involved in RUNX1-mutated patients of with acute myeloid leukemia (AML) and to explore possible molecular markers for novel therapeutic targets of the disease. MATERIAL AND METHODS The RNA sequencing datasets of 151 cases of AML were obtained from the Cancer Genome Atlas database. Differentially expressed genes (DEGs) were identified using edgeR of the R platform. PPI (protein-protein interaction) network clustering modules were analyzed with ClusterONE, and the KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analyses for modules were performed. RESULTS A total of 379 genes were identified as DEGs. The KEGG enrichment analysis of DEGs showed significantly enriched pathways in cancer, extracellular matrix (ECM)-receptor interaction pathway, and cyclic adenosine monophosphate (cAMP) signaling pathway. The top 10 genes ranked by degree were PRKACG, ANKRD7, RNFL7, ROPN11, TEX14, PRMT8, OTOA, CFAP99, NRXN1, and DMRT1, which were identified as hub genes from the protein-protein interaction network (PPI). Statistical analysis revealed that RUNX1-mutated patients with AML had a shorter median survival time (MST) with poor clinical outcome and an increased risk of death when compared with those without RUNX1 mutations. CONCLUSIONS DEGs and pathways identified in the present study will help understand the molecular mechanisms underlying RUNX1 mutations in AML and develop effective therapeutic strategies for RUNX1-mutation AML.

Identification of Key Genes and Pathways Associated with RUNX1 Mutations in Acute Myeloid Leukemia Using Bioinformatics Analysis

BACKGROUND RUNXl plays a key regulatory role in the process of hematopoiesis and is a common target for multiple chromosomal translocations in human acute leukemia. Mutations of RUNX1 gene can lead to acute leukemia and affect the prognosis of AML patients. We aimed to identify pivotal genes and pathways involved in RUNX1-mutated patients of with acute myeloid leukemia (AML) and to explore possible molecular markers for novel therapeutic targets of the disease. MATERIAL AND METHODS The RNA sequencing datasets of 151 cases of AML were obtained from the Cancer Genome Atlas database. Differentially expressed genes (DEGs) were identified using edgeR of the R platform. PPI (protein-protein interaction) network clustering modules were analyzed with ClusterONE, and the KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analyses for modules were performed. RESULTS A total of 379 genes were identified as DEGs. The KEGG enrichment analysis of DEGs showed significantly enriched pathways in cancer, extracellular matrix (ECM)-receptor interaction pathway, and cyclic adenosine monophosphate (cAMP) signaling pathway. The top 10 genes ranked by degree were PRKACG, ANKRD7, RNFL7, ROPN11, TEX14, PRMT8, OTOA, CFAP99, NRXN1, and DMRT1, which were identified as hub genes from the protein-protein interaction network (PPI). Statistical analysis revealed that RUNX1-mutated patients with AML had a shorter median survival time (MST) with poor clinical outcome and an increased risk of death when compared with those without RUNX1 mutations. CONCLUSIONS DEGs and pathways identified in the present study will help understand the molecular mechanisms underlying RUNX1 mutations in AML and develop effective therapeutic strategies for RUNX1-mutation AML.

Analytical Strategy to Prioritize Alzheimer's Disease Candidate Genes in Gene Regulatory Networks Using Public Expression Data

Alzheimer’s disease (AD) progressively destroys cognitive abilities in the aging population with tremendous effects on memory. Despite recent progress in understanding the underlying mechanisms, high drug attrition rates have put a question mark behind our knowledge about its etiology. Re-evaluation of past studies could help us to elucidate molecular-level details of this disease. Several methods to infer such networks exist, but most of them do not elaborate on context specificity and completeness of the generated networks, missing out on lesser-known candidates. In this study, we present a novel strategy that corroborates common mechanistic patterns across large scale AD gene expression studies and further prioritizes potential biomarker candidates. To infer gene regulatory networks (GRNs), we applied an optimized version of the BC3Net algorithm, named BC3Net10, capable of deriving robust and coherent patterns. In principle, this approach initially leverages the power of literature knowledge to extract AD specific genes for generating viable networks. Our findings suggest that AD GRNs show significant enrichment for key signaling mechanisms involved in neurotransmission. Among the prioritized genes, well-known AD genes were prominent in synaptic transmission, implicated in cognitive deficits. Moreover, less intensive studied AD candidates (STX2, HLA-F, HLA-C, RAB11FIP4, ARAP3, AP2A2, ATP2B4, ITPR2, and ATP2A3) are also involved in neurotransmission, providing new insights into the underlying mechanism. To our knowledge, this is the first study to generate knowledge-instructed GRNs that demonstrates an effective way of combining literature-based knowledge and data-driven analysis to identify lesser known candidates embedded in stable and robust functional patterns across disparate datasets.

A comparative analytical assay of gene regulatory networks inferred using microarray and RNA-seq datasets

A Gene Regulatory Network (GRN) is a collection of interactions between molecular regulators and their targets in cells governing gene expression level. Omics data explosion generated from high-throughput genomic assays such as microarray and RNA-Seq technologies and the emergence of a number of pre-processing methods demands suitable guidelines to determine the impact of transcript data platforms and normalization procedures on describing associations in GRNs. In this study exploiting publically available microarray and RNA-Seq datasets and a gold standard of transcriptional interactions in Arabidopsis, we performed a comparison between six GRNs derived by RNA-Seq and microarray data and different normalization procedures. As a result we observed that compared algorithms were highly data-specific and Networks reconstructed by RNA-Seq data revealed a considerable accuracy against corresponding networks captured by microarrays. Topological analysis showed that GRNs inferred from two platforms were similar in several of topological features although we observed more connectivity in RNA-Seq derived genes network. Taken together transcriptional regulatory networks obtained by Robust Multiarray Averaging (RMA) and Variance-Stabilizing Transformed (VST) normalized data demonstrated predicting higher rate of true edges over the rest of methods used in this comparison.