MicroRNA–mRNA interaction network using TSK-type recurrent neural fuzzy network

MicroRNAs are short non-coding RNAs that can regulate gene expression during various crucial cell processes such as differentiation, proliferation and apoptosis. Changes in expression profiles of miRNA play an important role in the development of many cancers, including CRC. Therefore, the identification of cancer related miRNAs and their target genes are important for cancer biology research. In this paper, we applied TSK-type recurrent neural fuzzy network (TRNFN) to infer miRNA–mRNA association network from paired miRNA, mRNA expression profiles of CRC patients. We demonstrated that the method we proposed achieved good performance in recovering known experimentally verified miRNA–mRNA associations. Moreover, our approach proved successful in identifying 17 validated cancer miRNAs which are directly involved in the CRC related pathways. Targeting such miRNAs may help not only to prevent the recurrence of disease but also to control the growth of advanced metastatic tumors. Our regulatory modules provide valuable insights into the pathogenesis of cancer.     

Geometric interpretation of gene coexpression network analysis

THE MERGING OF NETWORK THEORY AND MICROARRAY DATA ANALYSIS TECHNIQUES HAS SPAWNED A NEW FIELD: gene coexpression network analysis. While network methods are increasingly used in biology, the network vocabulary of computational biologists tends to be far more limited than that of, say, social network theorists. Here we review and propose several potentially useful network concepts. We take advantage of the relationship between network theory and the field of microarray data analysis to clarify the meaning of and the relationship among network concepts in gene coexpression networks. Network theory offers a wealth of intuitive concepts for describing the pairwise relationships among genes, which are depicted in cluster trees and heat maps. Conversely, microarray data analysis techniques (singular value decomposition, tests of differential expression) can also be used to address difficult problems in network theory. We describe conditions when a close relationship exists between network analysis and microarray data analysis techniques, and provide a rough dictionary for translating between the two fields. Using the angular interpretation of correlations, we provide a geometric interpretation of network theoretic concepts and derive unexpected relationships among them. We use the singular value decomposition of module expression data to characterize approximately factorizable gene coexpression networks, i.e., adjacency matrices that factor into node specific contributions. High and low level views of coexpression networks allow us to study the relationships among modules and among module genes, respectively. We characterize coexpression networks where hub genes are significant with respect to a microarray sample trait and show that the network concept of intramodular connectivity can be interpreted as a fuzzy measure of module membership. We illustrate our results using human, mouse, and yeast microarray gene expression data. The unification of coexpression network methods with traditional data mining methods can inform the application and development of systems biologic methods.     

Analysis of gene expression for studying tumor progression: the case of glucocorticoid administration

Background

Glucocorticoids are commonly used as adjuvant treatment for side-effects and have anti-proliferative activity in several tumors but, on the other hand, their proliferative effect has been reported in several studies, some of them involving the spread of cancer. We shall attempt to reconcile these incongruities from the genomic and tissue-physiology perspectives with our findings.

Methods

An accurate phenotype analysis of microarray data can help to solve multiple paradoxes derived from tumor-progression models. We have developed a new strategy to facilitate the study of interdependences among the phenotypes defined by the sample clusters obtained by common clustering methods (HC, SOTA, SOM, PAM). These interdependences are obtained by the detection of non-linear expression-relationships where each fluctuation in the relationship implies a phenotype change and each relationship typology implies a specific phenotype interdependence. As a result, multiple phenotypic changes are identified together with the genes involved in the phenotype transitions. In this way, we study the phenotypic changes from microarray data that describe common phenotypes in cancer from different tissues, and we cross our results with biomedical databases to relate the glucocorticoid activity to the phenotypic changes.

Results

11,244 significant non-linear expression relationships, classified into 11 different typologies, have been detected from the data matrix analyzed. From them, 415 non-linear expression relationships were related to glucocorticoid activity. Studying them, we have found the possible reason for opposite effects of some stressor agents like dexamethasone on tumor progression and it has been confirmed by literature. This hidden reason has resulted in being linked with the type of tumor progression of the tissues. In the first type of tumor progression found, new cells can be stressed during proliferation and stressor agents increase tumor proliferation. In the second type, cell stress and tumor proliferation are antagonists so, therefore, stressor agents stop tumor proliferation in order to stress the cells. The non-linear expression relationships among DUSP6, FERMT2, FKBP5, EGFR, NEDD4L and CITED2 genes are used to synthesize these findings.     

Empirical Bayes identification [correction of identication] of tumor progression genes from microarray data

The use of microarray data has become quite commonplace in medical and scientific experiments. We focus here on microarray data generated from cancer studies. It is potentially important for the discovery of biomarkers to identify genes whose expression levels correlate with tumor progression. In this article, we propose a simple procedure for the identification of such genes, which we term tumor progression genes. The first stage involves estimation based on the proportional odds model. At the second stage, we calculate two quantities: a q-value, and a shrinkage estimator of the test statistic is constructed to adjust for the multiple testing problem. The relationship between the proposed method with the false discovery rate is studied. The proposed methods are applied to data from a prostate cancer microarray study.     

Identification of cancer genes using a statistical framework for multiexperiment analysis of nondiscretized array CGH data

Tumor formation is in part driven by DNA copy number alterations (CNAs), which can be measured using microarray-based Comparative Genomic Hybridization (aCGH). Multiexperiment analysis of aCGH data from tumors allows discovery of recurrent CNAs that are potentially causal to cancer development. Until now, multiexperiment aCGH data analysis has been dependent on discretization of measurement data to a gain, loss or no-change state. Valuable biological information is lost when a heterogeneous system such as a solid tumor is reduced to these states. We have developed a new approach which inputs nondiscretized aCGH data to identify regions that are significantly aberrant across an entire tumor set. Our method is based on kernel regression and accounts for the strength of a probe's signal, its local genomic environment and the signal distribution across multiple tumors. In an analysis of 89 human breast tumors, our method showed enrichment for known cancer genes in the detected regions and identified aberrations that are strongly associated with breast cancer subtypes and clinical parameters. Furthermore, we identified 18 recurrent aberrant regions in a new dataset of 19 p53-deficient mouse mammary tumors. These regions, combined with gene expression microarray data, point to known cancer genes and novel candidate cancer genes.     

Discriminating lysosomal membrane protein types using dynamic neural network

This work presents a dynamic artificial neural network methodology, which classifies the proteins into their classes from their sequences alone: the lysosomal membrane protein classes and the various other membranes protein classes. In this paper, neural networks-based lysosomal-associated membrane protein type prediction system is proposed. Different protein sequence representations are fused to extract the features of a protein sequence, which includes seven feature sets; amino acid (AA) composition, sequence length, hydrophobic group, electronic group, sum of hydrophobicity, R-group, and dipeptide composition. To reduce the dimensionality of the large feature vector, we applied the principal component analysis. The probabilistic neural network, generalized regression neural network, and Elman regression neural network (RNN) are used as classifiers and compared with layer recurrent network (LRN), a dynamic network. The dynamic networks have memory, i.e. its output depends not only on the input but the previous outputs also. Thus, the accuracy of LRN classifier among all other artificial neural networks comes out to be the highest. The overall accuracy of jackknife cross-validation is 93.2% for the data-set. These predicted results suggest that the method can be effectively applied to discriminate lysosomal associated membrane proteins from other membrane proteins (Type-I, Outer membrane proteins, GPI-Anchored) and Globular proteins, and it also indicates that the protein sequence representation can better reflect the core feature of membrane proteins than the classical AA composition.     

Pattern identification and classification in gene expression data using an autoassociative neural network model

The application of DNA microarray technology for analysis of gene expression creates enormous opportunities to accelerate the pace in understanding living systems and identification of target genes and pathways for drug development and therapeutic intervention. Parallel monitoring of the expression profiles of thousands of genes seems particularly promising for a deeper understanding of cancer biology and the identification of molecular signatures supporting the histological classification schemes of neoplastic specimens. However, the increasing volume of data generated by microarray experiments poses the challenge of developing equally efficient methods and analysis procedures to extract, interpret, and upgrade the information content of these databases. Herein, a computational procedure for pattern identification, feature extraction, and classification of gene expression data through the analysis of an autoassociative neural network model is described. The identified patterns and features contain critical information about gene-phenotype relationships observed during changes in cell physiology. They represent a rational and dimensionally reduced base for understanding the basic biology of the onset of diseases, defining targets of therapeutic intervention, and developing diagnostic tools for the identification and classification of pathological states. The proposed method has been tested on two different microarray datasets-Golub's analysis of acute human leukemia [Golub et al. (1999) Science 286:531-537], and the human colon adenocarcinoma study presented by Alon et al. [1999; Proc Natl Acad Sci USA 97:10101-10106]. The analysis of the neural network internal structure allows the identification of specific phenotype markers and the extraction of peculiar associations among genes and physiological states. At the same time, the neural network outputs provide assignment to multiple classes, such as different pathological conditions or tissue samples, for previously unseen instances.     

Prioritizing candidate peptides for cancer vaccines through predicting peptide presentation by HLA-I proteins

Cancer (treatment) vaccines that are made of neoantigens, or peptides unique to tumor cells due to somatic mutations, have emerged as a promising method to reinvigorate the immune response against cancer. A key step to prioritizing neoantigens for cancer vaccines is computationally predicting which neoantigens are presented on the cell surface by a human leukocyte antigen (HLA). We propose to address this challenge by training a neural network using mass spectrometry (MS) data composed of peptides presented by at least one of several HLAs of a subject. We embed the neural network within a mixture model and train the neural network by maximizing the likelihood of the mixture model. After evaluating our method using data sets where the peptide presentation status was known, we applied it to analyze somatic mutations of 60 melanoma patients and identified a group of neoantigens more immunogenic in tumor cells than in normal cells. Moreover, neoantigen burden estimated by our method was significantly associated with a measurement of the immune system activity, suggesting these neoantigens could induce an immune response.     

GWAS Meets Microarray: Are the Results of Genome-Wide Association Studies and Gene-Expression Profiling Consistent? Prostate Cancer as an Example

Background

Genome-wide association studies (GWASs) and global profiling of gene expression (microarrays) are two major technological breakthroughs that allow hypothesis-free identification of candidate genes associated with tumorigenesis. It is not obvious whether there is a consistency between the candidate genes identified by GWAS (GWAS genes) and those identified by profiling gene expression (microarray genes).

Methodology/Principal Findings

We used the Cancer Genetic Markers Susceptibility database to retrieve single nucleotide polymorphisms from candidate genes for prostate cancer. In addition, we conducted a large meta-analysis of gene expression data in normal prostate and prostate tumor tissue. We identified 13,905 genes that were interrogated by both GWASs and microarrays. On the basis of P values from GWASs, we selected 1,649 most significantly associated genes for functional annotation by the Database for Annotation, Visualization and Integrated Discovery. We also conducted functional annotation analysis using same number of the top genes identified in the meta-analysis of the gene expression data. We found that genes involved in cell adhesion were overrepresented among both the GWAS and microarray genes.

Conclusions/Significance

We conclude that the results of these analyses suggest that combining GWAS and microarray data would be a more effective approach than analyzing individual datasets and can help to refine the identification of candidate genes and functions associated with tumor development.     

Discovering biological progression underlying microarray samples

In biological systems that undergo processes such as differentiation, a clear concept of progression exists. We present a novel computational approach, called Sample Progression Discovery (SPD), to discover patterns of biological progression underlying microarray gene expression data. SPD assumes that individual samples of a microarray dataset are related by an unknown biological process (i.e., differentiation, development, cell cycle, disease progression), and that each sample represents one unknown point along the progression of that process. SPD aims to organize the samples in a manner that reveals the underlying progression and to simultaneously identify subsets of genes that are responsible for that progression. We demonstrate the performance of SPD on a variety of microarray datasets that were generated by sampling a biological process at different points along its progression, without providing SPD any information of the underlying process. When applied to a cell cycle time series microarray dataset, SPD was not provided any prior knowledge of samples' time order or of which genes are cell-cycle regulated, yet SPD recovered the correct time order and identified many genes that have been associated with the cell cycle. When applied to B-cell differentiation data, SPD recovered the correct order of stages of normal B-cell differentiation and the linkage between preB-ALL tumor cells with their cell origin preB. When applied to mouse embryonic stem cell differentiation data, SPD uncovered a landscape of ESC differentiation into various lineages and genes that represent both generic and lineage specific processes. When applied to a prostate cancer microarray dataset, SPD identified gene modules that reflect a progression consistent with disease stages. SPD may be best viewed as a novel tool for synthesizing biological hypotheses because it provides a likely biological progression underlying a microarray dataset and, perhaps more importantly, the candidate genes that regulate that progression.     

circRNA-associated ceRNA network construction reveals the circRNAs involved in the progression and prognosis of breast cancer

Recently, increasing evidences show that circular RNAs (circRNAs) are important regulators of various diseases, especially cancer. However, the regulatory role and the potential mechanism of action of circRNAs in breast cancer remain largely unknown. In this study, weighted gene co-expression network analysis was conducted with the differentially expressed miRNAs and mRNAs in breast cancer from The Cancer Genome Atlas database to identify the key modules associated with the carcinogenesis of breast cancer. In the significant turquoise and brown modules, 22 miRNAs and 1877 mRNAs were identified, respectively. Then, We compared and predicted the target genes and performed survival analysis to identify the miRNAs and mRNAs related to the prognosis of breast cancer. A circRNA-related competitive endogenous RNA network was identified by database co-screening, and deleted in liver cancer 1 (DLC1) was identified as a key gene. Finally, to assess how genes in key modules and key genes contribute to the development of breast cancer, relevant pathway information was obtained through DAVID and Gene Set Enrichment Analysis. These data demonstrated that three circRNAs (hsa-circ-0083373, hsa-circ-0083374, and hsa-circ-0083375) that regulate DLC1 expression via hsa-mir-511 and are involved in the pathogenesis and development of breast cancer.     

A regulatory gene network related to the porcine umami taste receptor (TAS1R1/TAS1R3)

Taste perception plays an important role in the mediation of food choices in mammals. The first porcine taste receptor genes identified, sequenced and characterized, TAS1R1 and TAS1R3, were related to the dimeric receptor for umami taste. However, little is known about their regulatory network. The objective of this study was to unfold the genetic network involved in porcine umami taste perception. We performed a meta‐analysis of 20 gene expression studies spanning 480 porcine microarray chips and screened 328 taste‐related genes by selective mining steps among the available 12 320 genes. A porcine umami taste‐specific regulatory network was constructed based on the normalized coexpression data of the 328 genes across 27 tissues. From the network, we revealed the ‘taste module’ and identified a coexpression cluster for the umami taste according to the first connector with the TAS1R1/TAS1R3 genes. Our findings identify several taste‐related regulatory genes and extend previous genetic background of porcine umami taste.