Randomness and preserved patterns in cancer network

Architecture and topologies of gene regulatory networks associated with breast cancer, adjacent normal, and normal tissues

Most cancer studies employ adjacent normal tissues to tumors (ANTs) as controls, which are not completely normal and represent a pre-cancerous state. However, the regulatory landscape of ANTs compared to tumor and non-tumor-bearing normal tissues is largely unexplored. Among cancers, breast cancer is the most commonly diagnosed cancer and a leading cause of death in women worldwide, with a lack of sufficient treatment regimens for various reasons. Hence, we aimed to gain deeper insights into normal, pre-cancerous, and cancerous regulatory systems of breast tissues towards identifying ANT and subtype-specific candidate genes. For this, we constructed and analyzed eight gene regulatory networks (GRNs), including five subtypes (viz., Basal, Her2, Luminal A, Luminal B, and Normal-Like), one ANT, and two normal tissue networks. Whereas several topological properties of these GRNs enabled us to identify tumor-related features of ANT, escape velocity centrality (EVC+) identified 24 functionally significant common genes, including well-known genes such as E2F1, FOXA1, JUN, BRCA1, GATA3, ERBB2, and ERBB3 across all six tissues including subtypes and ANT. Similarly, the EVC+ also helped us to identify tissue-specific key genes (Basal: 18, Her2: 6, Luminal A: 5, Luminal B: 5, Normal-Like: 2, and ANT: 7). Additionally, differentially correlated switching gene pairs along with functional, pathway, and disease annotations highlighted the cancer-associated role of these genes. In a nutshell, the present study revealed ANT and subtype-specific regulatory features and key candidate genes, which can be explored further using in vitro and in vivo experiments for better and effective disease management at an early stage.

Random matrix theory tools for the predictive analysis of functional magnetic resonance imaging examinations

Purpose

Random matrix theory (RMT) is an increasingly useful tool for understanding large, complex systems. Prior studies have examined functional magnetic resonance imaging (fMRI) scans using tools from RMT, with some success. However, RMT computations are highly sensitive to a number of analytic choices, and the robustness of findings involving RMT remains in question. We systematically investigate the usefulness of RMT on a wide variety of fMRI datasets using a rigorous predictive framework.

Approach

We develop open-source software to efficiently compute RMT features from fMRI images and examine the cross-validated predictive potential of eigenvalue and RMT-based features ("eigenfeatures") with classic machine-learning classifiers. We systematically vary pre-processing extent, normalization procedures, RMT unfolding procedures, and feature selection and compare the impact of these analytic choices on the distributions of cross-validated prediction performance for each combination of dataset binary classification task, classifier, and feature. To deal with class imbalance, we use the area under the receiver operating characteristic curve (AUROC) as the main performance metric.

Results

Across all classification tasks and analytic choices, we find RMT- and eigenvalue-based "eigenfeatures" to have predictive utility more often than not (82.4% of median AUROCs > 0.5 ; median AUROC range across classification tasks 0.47 to 0.64). Simple baseline reductions on source timeseries, by contrast, were less useful (58.8% of median AUROCs > 0.5 , median AUROC range across classification tasks 0.42 to 0.62). Additionally, eigenfeature AUROC distributions were overall more right-tailed than baseline features, suggesting greater predictive potential. However, performance distributions were wide and often significantly affected by analytic choices.

Conclusions

Eigenfeatures clearly have potential for understanding fMRI functional connectivity in a wide variety of scenarios. The utility of these features is strongly dependent on analytic decisions, suggesting caution when interpreting past and future studies applying RMT to fMRI. However, our study demonstrates that the inclusion of RMT statistics in fMRI investigations could improve prediction performances across a wide variety of phenomena.

Mapping drug-target interactions and synergy in multi-molecular therapeutics for pressure-overload cardiac hypertrophy

Advancements in systems biology have resulted in the development of network pharmacology, leading to a paradigm shift from "one-target, one-drug" to "target-network, multi-component therapeutics". We employ a chimeric approach involving in-vivo assays, gene expression analysis, cheminformatics, and network biology to deduce the regulatory actions of a multi-constituent Ayurvedic concoction, Amalaki Rasayana (AR) in animal models for its effect in pressure-overload cardiac hypertrophy. The proteomics analysis of in-vivo assays for Aorta Constricted and Biologically Aged rat models identify proteins expressed under each condition. Network analysis mapping protein-protein interactions and synergistic actions of AR using multi-component networks reveal drug targets such as ACADM, COX4I1, COX6B1, HBB, MYH14, and SLC25A4, as potential pharmacological co-targets for cardiac hypertrophy. Further, five out of eighteen AR constituents potentially target these proteins. We propose a distinct prospective strategy for the discovery of network pharmacological therapies and repositioning of existing drug molecules for treating pressure-overload cardiac hypertrophy.

Analysis of human mitochondrial genome co-occurrence networks of Asian population at varying altitudes

Networks have been established as an extremely powerful framework to understand and predict the behavior of many large-scale complex systems. We studied network motifs, the basic structural elements of networks, to describe the possible role of co-occurrence of genomic variations behind high altitude adaptation in the Asian human population. Mitochondrial DNA (mtDNA) variations have been acclaimed as one of the key players in understanding the biological mechanisms behind adaptation to extreme conditions. To explore the cumulative effects of variations in the mitochondrial genome with the variation in the altitude, we investigated human mt-DNA sequences from the NCBI database at different altitudes under the co-occurrence motifs framework. Analysis of the co-occurrence motifs using similarity clustering revealed a clear distinction between lower and higher altitude regions. In addition, the previously known high altitude markers 3394 and 7697 (which are definitive sites of haplogroup M9a1a1c1b) were found to co-occur within their own gene complexes indicating the impact of intra-genic constraint on co-evolution of nucleotides. Furthermore, an ancestral 'RSRS50' variant 10,398 was found to co-occur only at higher altitudes supporting the fact that a separate route of colonization at these altitudes might have taken place. Overall, our analysis revealed the presence of co-occurrence interactions specific to high altitude at a whole mitochondrial genome level. This study, combined with the classical haplogroups analysis is useful in understanding the role of co-occurrence of mitochondrial variations in high altitude adaptation.

The Human Body as a Super Network: Digital Methods to Analyze the Propagation of Aging

Biological aging is a complex process involving multiple biological processes. These can be understood theoretically though considering them as individual networks-e.g., epigenetic networks, cell-cell networks (such as astroglial networks), and population genetics. Mathematical modeling allows the combination of such networks so that they may be studied in unison, to better understand how the so-called "seven pillars of aging" combine and to generate hypothesis for treating aging as a condition at relatively early biological ages. In this review, we consider how recent progression in mathematical modeling can be utilized to investigate aging, particularly in, but not exclusive to, the context of degenerative neuronal disease. We also consider how the latest techniques for generating biomarker models for disease prediction, such as longitudinal analysis and parenclitic analysis can be applied to as both biomarker platforms for aging, as well as to better understand the inescapable condition. This review is written by a highly diverse and multi-disciplinary team of scientists from across the globe and calls for greater collaboration between diverse fields of research.

Network spectra for drug-target identification in complex diseases: new guns against old foes

The fundamental understanding of altered complex molecular interactions in a diseased condition is the key to its cure. The overall functioning of these molecules is kind of jugglers play in the cell orchestra and to anticipate these relationships among the molecules is one of the greatest challenges in modern biology and medicine. Network science turned out to be providing a successful and simple platform to understand complex interactions among healthy and diseased tissues. Furthermore, much information about the structure and dynamics of a network is concealed in the eigenvalues of its adjacency matrix. In this review, we illustrate rapid advancements in the field of network science in combination with spectral graph theory that enables us to uncover the complexities of various diseases. Interpretations laid by network science approach have solicited insights into molecular relationships and have reported novel drug targets and biomarkers in various complex diseases.

Spectral properties of complex networks

This review presents an account of the major works done on spectra of adjacency matrices drawn on networks and the basic understanding attained so far. We have divided the review under three sections: (a) extremal eigenvalues, (b) bulk part of the spectrum, and (c) degenerate eigenvalues, based on the intrinsic properties of eigenvalues and the phenomena they capture. We have reviewed the works done for spectra of various popular model networks, such as the Erdős-Rényi random networks, scale-free networks, 1-d lattice, small-world networks, and various different real-world networks. Additionally, potential applications of spectral properties for natural processes have been reviewed.

Random Matrix Analysis for Gene Interaction Networks in Cancer Cells

Investigations of topological uniqueness of gene interaction networks in cancer cells are essential for understanding the disease. Although cancer is considered to originate from the topological alteration of a huge molecular interaction network in cellular systems, the theoretical study to investigate such complex networks is still insufficient. It is necessary to predict the behavior of a huge complex interaction network from the behavior of a finite size network. Based on the random matrix theory, we study the distribution of the nearest neighbor level spacings P(s) of interaction matrices of gene networks in human cancer cells. The interaction matrices are computed using the Cancer Network Galaxy (TCNG) database which is a repository of gene interactions inferred by a Bayesian network model. 256 NCBI GEO entries regarding gene expressions in human cancer cells have been used for the inference. We observe the Wigner distribution of P(s) when the gene networks are dense networks that have more than ~38,000 edges. In the opposite case, when the networks have smaller numbers of edges, the distribution P(s) becomes the Poisson distribution. We investigate relevance of P(s) both to the sparseness of the networks and to edge frequency factor which is the reliance (likelihood) of the inferred gene interactions.

Dysregulated genes targeted by microRNAs and metabolic pathways in bladder cancer revealed by bioinformatics methods

The present study aimed to identify bladder cancer-associated microRNAs (miRNAs) and target genes, and further analyze the potential molecular mechanisms involved in bladder cancer. The mRNA and miRNA expression profiling dataset GSE40355 was downloaded from the Gene Expression Omnibus database. The Limma package in R was used to identify differential expression levels. The Human microRNA Disease Database was used to identify bladder cancer-associated miRNAs and Target prediction programs were used to screen for miRNA target genes. Enrichment analysis was performed to identify biological functions. The Database for Annotation, Visualization and Integration Discovery was used to perform OMIM_DISEASE analysis, and then protein-protein interaction (PPI) analysis was performed to identify hubs with biological essentiality. ClusterONE plugins in cytoscape were used to screen modules and the InterPro database was used to perform protein domain enrichment analysis. A group of 573 disease dysregulated genes were identified in the present study. Enrichment analysis indicated that the muscle organ development and vascular smooth muscle contraction pathways were significantly enriched in terms of disease dysregulated genes. miRNAs targets (frizzled class receptor 8, EYA transcriptional coactivator and phosphatase 4, sacsin molecular chaperone, calcium voltage-gated channel auxiliary subunit β2, peptidase inhibitor 15 and catenin α2) were mostly associated with bladder cancer. PPI analysis revealed that calmodulin 1 (CALM1), Jun proto-oncogene, AP-1 transcription factor subunit (JUN) and insulin like growth factor 1 (IGF1) were the important hub nodes. Additionally, protein domain enrichment analysis indicated that the serine/threonine protein kinase active site was enriched in module 1 extracted from the PPI network. Overall, the results suggested that the IGF signaling pathway and RAS/MEK/extracellular signal-regulated kinase transduction signaling may exert vital molecular mechanisms in bladder cancer, and that CALM1, JUN and IGF1 may be used as novel potential therapeutic targets.

Understanding cancer complexome using networks, spectral graph theory and multilayer framework

Cancer complexome comprises a heterogeneous and multifactorial milieu that varies in cytology, physiology, signaling mechanisms and response to therapy. The combined framework of network theory and spectral graph theory along with the multilayer analysis provides a comprehensive approach to analyze the proteomic data of seven different cancers, namely, breast, oral, ovarian, cervical, lung, colon and prostate. Our analysis demonstrates that the protein-protein interaction networks of the normal and the cancerous tissues associated with the seven cancers have overall similar structural and spectral properties. However, few of these properties implicate unsystematic changes from the normal to the disease networks depicting difference in the interactions and highlighting changes in the complexity of different cancers. Importantly, analysis of common proteins of all the cancer networks reveals few proteins namely the sensors, which not only occupy significant position in all the layers but also have direct involvement in causing cancer. The prediction and analysis of miRNAs targeting these sensor proteins hint towards the possible role of these proteins in tumorigenesis. This novel approach helps in understanding cancer at the fundamental level and provides a clue to develop promising and nascent concept of single drug therapy for multiple diseases as well as personalized medicine.

Graph-theoretical model of global human interactome reveals enhanced long-range communicability in cancer networks

Malignant transformation is known to involve substantial rearrangement of the molecular genetic landscape of the cell. A common approach to analysis of these alterations is a reductionist one and consists of finding a compact set of differentially expressed genes or associated signaling pathways. However, due to intrinsic tumor heterogeneity and tissue specificity, biomarkers defined by a small number of genes/pathways exhibit substantial variability. As an alternative to compact differential signatures, global features of genetic cell machinery are conceivable. Global network descriptors suggested in previous works are, however, known to potentially be biased by overrepresentation of interactions between frequently studied genes-proteins. Here, we construct a cellular network of 74538 directional and differential gene expression weighted protein-protein and gene regulatory interactions, and perform graph-theoretical analysis of global human interactome using a novel, degree-independent feature—the normalized total communicability (NTC). We apply this framework to assess differences in total information flow between different cancer (BRCA/COAD/GBM) and non-cancer interactomes. Our experimental results reveal that different cancer interactomes are characterized by significant enhancement of long-range NTC, which arises from circulation of information flow within robustly organized gene subnetworks. Although enhancement of NTC emerges in different cancer types from different genomic profiles, we identified a subset of 90 common genes that are related to elevated NTC in all studied tumors. Our ontological analysis shows that these genes are associated with enhanced cell division, DNA replication, stress response, and other cellular functions and processes typically upregulated in cancer. We conclude that enhancement of long-range NTC manifested in the correlated activity of genes whose tight coordination is required for survival and proliferation of all tumor cells, and, thus, can be seen as a graph-theoretical equivalent to some hallmarks of cancer. The computational framework for differential network analysis presented herein is of potential interest for a wide range of network perturbation problems given by single or multiple gene-protein activation-inhibition.

Interplay of degree correlations and cluster synchronization

We study the evolution of coupled chaotic dynamics on networks and investigate the role of degree-degree correlation in the networks' cluster synchronizability. We find that an increase in the disassortativity can lead to an increase or a decrease in the cluster synchronizability depending on the degree distribution and average connectivity of the network. Networks with heterogeneous degree distribution exhibit significant changes in cluster synchronizability as well as in the phenomena behind cluster synchronization as compared to those of homogeneous networks. Interestingly, cluster synchronizability of a network may be very different from global synchronizability due to the presence of the driven phenomenon behind the cluster formation. Furthermore, we show how degeneracy at the zero eigenvalues provides an understanding of the occurrence of the driven phenomenon behind the synchronization in disassortative networks. The results demonstrate the importance of degree-degree correlations in determining cluster synchronization behavior of complex networks and hence have potential applications in understanding and predicting dynamical behavior of complex systems ranging from brain to social systems.