Comparison of RNA-Seq and microarray in the prediction of protein expression and survival prediction

Gene expression profiling using RNA-sequencing (RNA-seq) and microarray technologies is widely used in cancer research to identify biomarkers for clinical endpoint prediction. We compared the performance of these two methods in predicting protein expression and clinical endpoints using The Cancer Genome Atlas (TCGA) datasets of lung cancer, colorectal cancer, renal cancer, breast cancer, endometrial cancer, and ovarian cancer. We calculated the correlation coefficients between gene expression measured by RNA-seq or microarray and protein expression measured by reverse phase protein array (RPPA). In addition, after selecting the top 103 survival-related genes, we compared the random forest survival prediction model performance across test platforms and cancer types. Both RNA-seq and microarray data were retrieved from TCGA dataset. Most genes showed similar correlation coefficients between RNA-seq and microarray, but 16 genes exhibited significant differences between the two methods. The BAX gene was recurrently found in colorectal cancer, renal cancer, and ovarian cancer, and the PIK3CA gene belonged to renal cancer and breast cancer. Furthermore, the survival prediction model using microarray was better than the RNA-seq model in colorectal cancer, renal cancer, and lung cancer, but the RNA-seq model was better in ovarian and endometrial cancer. Our results showed good correlation between mRNA levels and protein measured by RPPA. While RNA-seq and microarray performance were similar, some genes showed differences, and further clinical significance should be evaluated. Additionally, our survival prediction model results were controversial

Networks of intergenic long-range enhancers and snpRNAs drive castration-resistant phenotype of prostate cancer and contribute to pathogenesis of multiple common human disorders

Biological and mechanistic relevance of intergenic disease-associated genetic loci (IDAGL) containing highly statistically significant disease-linked SNPs remains unknown. Here we present the experimental and clinical evidence revealing important role of IDAGL in human diseases. Targeted RT-PCR screen coupled with sequencing of purified PCR products detects widespread transcription at multiple intergenic disease-associated genomic loci (IDAGL) and identifies 96 small non-coding trans-regulatory RNAs of ~ 100-300 nt in length containing SNPs associated with 21 common human disorders (snpRNAs). Functionality of snpRNAs is supported by multiple independent lines of experimental evidence demonstrating their cell-type-specific expression and evolutionary conservation of sequences, genomic coordinates, and biological effects. Analysis of chromatin state signatures, expression profiling experiments using microarray and Q-PCR technologies, and luciferase reporter assays indicate that many IDAGL are Polycomb-regulated long-range enhancers. Expression of snpRNAs in human and mouse cells markedly affects cellular behavior and induces allele-specific clinically-relevant phenotypic changes: NLRP1-locus snpRNAs exert regulatory effects on monocyte/macrophage trans-differentiation, induce prostate cancer (PC) susceptibility snpRNAs, and transform low-malignancy hormone-dependent human PC cells into highly malignant androgen-independent PC. Q-PCR analysis and luciferase reporter assays demonstrate that snpRNA sequences represent allele-specific “decoy” targets of microRNAs which function as SNP-allele-specific modifiers of microRNA expression and activity. We demonstrate that trans-acting RNA molecules facilitating androgen depletion-independent growth (ADIG) in vitro and castration-resistant (CR) phenotype in vivo of PC contain intergenic 8q24-locus SNP variants which were recently linked with increased risk of developing PC. Expression level of 8q24-locus PC susceptibility snpRNAs is regulated by NLRP1-locus snpRNAs, which are transcribed from the intergenic long-range enhancer sequence located in 17p13 region at ~ 30 kb distance from the NLRP1 gene. Q-PCR analysis of clinical PC samples reveals markedly increased snpRNA expression levels in tumor tissues compared to the adjacent normal prostate [122-fold and 45-fold in Gleason 7 tumors (p = 0.03); 370-fold and 127-fold in Gleason 8 tumors (p = 0.0001); for NLRP1-locus and 8q24-locus SnpRNAs, respectively]. Highly concordant expression profiles of the NLRP1-locus snpRNAs and 8q24 CR-locus snpRNAs (r = 0.896; p < 0.0001) in clinical PC samples and experimental evidence of trans-regulatory effects of NLRP1-locus snpRNAs on expression of 8q24-locus SnpRNAs indicate that ADIG and CR phenotype of human PC cells can be triggered by RNA molecules transcribed from the NLRP1-locus intergenic enhancer and down-stream activation of the 8q24-locus snpRNAs. Our results define the intergenic NLRP1 and 8q24 regions as regulatory loci of ADIG and CR phenotype of human PC, reveal previously unknown molecular links between the innate immunity/inflammasome system and development of hormone-independent PC, and identify novel diagnostic and therapeutic targets exploration of which should be highly beneficial for clinical management of PC

Genomic and proteomic profiling for cancer diagnosis in dogs

Global gene expression, whereby tumours are classified according to similar gene expression patterns or ‘signatures’ regardless of cell morphology or tissue characteristics, is being increasingly used in both the human and veterinary fields to assist in cancer diagnosis and prognosis. Many studies on canine tumours have focussed on RNA expression using techniques such as microarrays or next generation sequencing. However, proteomic studies combining two-dimensional polyacrylamide gel electrophoresis or two-dimensional differential gel electrophoresis with mass spectrometry have also provided a wealth of data on gene expression in tumour tissues. In addition, proteomics has been instrumental in the search for tumour biomarkers in blood and other body fluids

Comprehensive profiling of zebrafish hepatic proximal promoter CpG island methylation and its modification during chemical carcinogenesis

Background\ud DNA methylation is an epigenetic mechanism associated with regulation of gene expression and it is modulated during chemical carcinogenesis. The zebrafish is increasingly employed as a human disease model; however there is a lack of information on DNA methylation in zebrafish and during fish tumorigenesis. \ud \ud Results\ud A novel CpG island tiling array containing 44,000 probes, in combination with immunoprecipitation of methylated DNA, was used to achieve the first comprehensive methylation profiling of normal adult zebrafish liver. DNA methylation alterations were detected in zebrafish liver tumors induced by the environmental carcinogen 7, 12-dimethylbenz(a)anthracene. Genes significantly hypomethylated in tumors were associated particularly with proliferation, glycolysis, transcription, cell cycle, apoptosis, growth and metastasis. Hypermethylated genes included those associated with anti-angiogenesis and cellular adhesion. Of 49 genes that were altered in expression within tumors, and which also had appropriate CpG islands and were co-represented on the tiling array, approximately 45% showed significant changes in both gene expression and methylation. \ud \ud Conclusion\ud The functional pathways containing differentially methylated genes in zebrafish hepatocellular carcinoma have also been reported to be aberrantly methylated during tumorigenesis in humans. These findings increase the confidence in the use of zebrafish as a model for human cancer in addition to providing the first comprehensive mapping of DNA methylation in the normal adult zebrafish liver. \ud \u

Comparability of Microarray Data between Amplified and Non Amplified RNA in Colorectal Carcinoma

Microarray analysis reaches increasing popularity during the investigation of prognostic gene clusters in oncology. The standardisation of technical procedures will be essential to compare various datasets produced by different research groups. In several projects the amount of available tissue is limited. In such cases the preamplification of RNA might be necessary prior to microarray hybridisation. To evaluate the comparability of microarray results generated either by amplified or non amplified RNA we isolated RNA from colorectal cancer samples (stage UICC IV) following tumour tissue enrichment by macroscopic manual dissection (CMD). One part of the RNA was directly labelled and hybridised to GeneChips (HG-U133A, Affymetrix), the other part of the RNA was amplified according to the ?Eberwine? protocol and was then hybridised to the microarrays. During unsupervised hierarchical clustering the samples were divided in groups regarding the RNA pre-treatment and 5.726 differentially expressed genes were identified. Using independent microarray data of 31 amplified vs. 24 non amplified RNA samples from colon carcinomas (stage UICC III) in a set of 50 predictive genes we validated the amplification bias. In conclusion microarray data resulting from different pre-processing regarding RNA pre-amplification can not be compared within one analysis

MMpred: functional miRNA – mRNA interaction analyses by miRNA expression prediction

Background: MicroRNA (miRNA) directed gene repression is an important mechanism of posttranscriptional regulation. Comprehensive analyses of how microRNA influence biological processes requires paired miRNA-mRNA expression datasets. However, a review of both GEO and ArrayExpress repositories revealed few such datasets, which was in stark contrast to the large number of messenger RNA (mRNA) only datasets. It is of interest that numerous primary miRNAs (precursors of microRNA) are known to be co-expressed with coding genes (host genes). Results: We developed a miRNA-mRNA interaction analyses pipeline. The proposed solution is based on two miRNA expression prediction methods – a scaling function and a linear model. Additionally, miRNA-mRNA anticorrelation analyses are used to determine the most probable miRNA gene targets (i.e. the differentially expressed genes under the influence of up- or down-regulated microRNA). Both the consistency and accuracy of the prediction method is ensured by the application of stringent statistical methods. Finally, the predicted targets are subjected to functional enrichment analyses including GO, KEGG and DO, to better understand the predicted interactions. Conclusions: The MMpred pipeline requires only mRNA expression data as input and is independent of third party miRNA target prediction methods. The method passed extensive numerical validation based on the binding energy between the mature miRNA and 3’ UTR region of the target gene. We report that MMpred is capable of generating results similar to that obtained using paired datasets. For the reported test cases we generated consistent output and predicted biological relationships that will help formulate further testable hypotheses

Elephant Search with Deep Learning for Microarray Data Analysis

Even though there is a plethora of research in Microarray gene expression data analysis, still, it poses challenges for researchers to effectively and efficiently analyze the large yet complex expression of genes. The feature (gene) selection method is of paramount importance for understanding the differences in biological and non-biological variation between samples. In order to address this problem, a novel elephant search (ES) based optimization is proposed to select best gene expressions from the large volume of microarray data. Further, a promising machine learning method is envisioned to leverage such high dimensional and complex microarray dataset for extracting hidden patterns inside to make a meaningful prediction and most accurate classification. In particular, stochastic gradient descent based Deep learning (DL) with softmax activation function is then used on the reduced features (genes) for better classification of different samples according to their gene expression levels. The experiments are carried out on nine most popular Cancer microarray gene selection datasets, obtained from UCI machine learning repository. The empirical results obtained by the proposed elephant search based deep learning (ESDL) approach are compared with most recent published article for its suitability in future Bioinformatics research.Comment: 12 pages, 5 Tabl