Search

Differential network analysis has become an important approach in identifying driver genes in development and disease. However, most studies capture only local features of the underlying gene-regulatory network topology. These approaches are vulnerable to noise and other changes which mask driver-gene activity. Therefore, methods are urgently needed which can separate the impact of true regulatory elements from stochastic changes and downstream effects. We propose the differential network flow (DNF) method to identify key regulators of progression in development or disease. Given the network representation of consecutive biological states, DNF quantifies the essentiality of each node by differences in the distribution of network flow, which are capable of capturing comprehensive topological differences from local to global feature domains. DNF achieves more accurate driver-gene identification than other state-of-the-art methods when applied to four human datasets from The Cancer Genome Atlas and three single-cell RNA-seq datasets of murine neural and hematopoietic differentiation. Furthermore, we predict key regulators of crosstalk between separate networks underlying both neuronal differentiation and the progression of neurodegenerative disease, among which APP is predicted as a driver gene of neural stem cell differentiation. Our method is a new approach for quantifying the essentiality of genes across networks of different biological states.

Glioblastoma (GBM) is a fast-growing type of malignant primary brain tumor. To explore the mechanisms in GBM, complex biological networks are used to reveal crucial changes among different biological states, which reflect on the development of living organisms. It is critical to discover the kernel differential subgraph (KDS) that leads to drastic changes. However, identifying the KDS is similar to the Steiner Tree problem that is an NP-hard problem. In this paper, we developed a criterion to explore the KDS (CKDS), which considered the connectivity and scale of KDS, the topological difference of nodes and function relevance between genes in the KDS. The CKDS algorithm was applied to simulated datasets and three single-cell RNA sequencing (scRNA-seq) datasets including GBM, fetal human cortical neurons (FHCN) and neural differentiation. Then we performed the network topology and functional enrichment analyses on the extracted KDSs. Compared with the state-of-art methods, the CKDS algorithm outperformed on simulated datasets to discover the KDSs. In the GBM and FHCN, seventeen genes (one biomarker, nine regulatory genes, one driver genes, six therapeutic targets) and KEGG pathways in KDSs were strongly supported by literature mining that they were highly interrelated with GBM. Moreover, focused on GBM, there were fifteen genes (including ten regulatory genes, three driver genes, one biomarkers, one therapeutic target) and KEGG pathways found in the KDS of neural differentiation process from activated neural stem cells (aNSC) to neural progenitor cells (NPC), while few genes and no pathway were found in the period from NPC to astrocytes (Ast). These experiments indicated that the process from aNSC to NPC is a key differentiation period affecting the development of GBM. Therefore, the CKDS algorithm provides a unique perspective in identifying cell-type-specific genes and KDSs.

Network alignment is an important bridge to understanding human protein-protein interactions (PPIs) and functions through model organisms. However, the underlying subgraph isomorphism problem complicates and increases the time required to align protein interaction networks (PINs). Parallel computing technology is an effective solution to the challenge of aligning large-scale networks via sequential computing. In this study, the typical Hungarian-Greedy Algorithm (HGA) is used as an example for PIN alignment. The authors propose a HGA with 2-nearest neighbours (HGA-2N) and implement its graphics processing unit (GPU) acceleration. Numerical experiments demonstrate that HGA-2N can find alignments that are close to those found by HGA while dramatically reducing computing time. The GPU implementation of HGA-2N optimises the parallel pattern, computing mode and storage mode and it improves the computing time ratio between the CPU and GPU compared with HGA when large-scale networks are considered. By using HGA-2N in GPUs, conserved PPIs can be observed, and potential PPIs can be predicted. Among the predictions based on 25 common Gene Ontology terms, 42.8% can be found in the Human Protein Reference Database. Furthermore, a new method of reconstructing phylogenetic trees is introduced, which shows the same relationships among five herpes viruses that are obtained using other methods.

Understanding the molecular mechanisms underlying cell migration, which plays an important role in tumor growth and progression, is critical for the development of novel tumor therapeutics. Overexpression of transmembrane protein 30A (TMEM30A) has been shown to initiate tumor cell migration, however, the molecular mechanisms through which this takes place have not yet been reported. Thus, we propose the integration of computational and experimental approaches by first predicting potential signaling networks regulated by TMEM30A using a) computational biology methods, b) our previous mass spectrometry results of the TMEM30A complex in mouse tissue, and c) a number of migration-related genes manually collected from the literature, and subsequently performing molecular biology experiments including the in vitro scratch assay and real-time quantitative polymerase chain reaction (qPCR) to validate the reliability of the predicted network. The results verify that the genes identified in the computational signaling network are indeed regulated by TMEM30A during cell migration, indicating the effectiveness of our proposed method and shedding light on the regulatory mechanisms underlying tumor migration, which facilitates the understanding of the molecular basis of tumor invasion. [ABSTRACT FROM AUTHOR]

Drought is the one of the most important environment stresses that restricts crop yield worldwide. Cassava (Manihot esculenta Crantz) is an important food and energy crop that has many desirable traits such as drought, heat and low nutrients tolerance. However, the mechanisms underlying drought tolerance in cassava are unclear. Ethylene signaling pathway, from the upstream receptors to the downstream transcription factors, plays important roles in environmental stress responses during plant growth and development. In this study, we used bioinformatics approaches to identify and characterize candidate Manihot esculenta ethylene receptor genes and transcription factor genes. Using computational methods, we localized these genes on cassava chromosomes, constructed phylogenetic trees and identified stress-responsive cis-elements within their 5’ upstream regions. Additionally, we measured the trehalose and proline contents in cassava fresh leaves after drought, osmotic, and salt stress treatments, and then it was found that the regulation patterns of contents of proline and trehalose in response to various dehydration stresses were differential, or even the opposite, which shows that plant may take different coping strategies to deal with different stresses, when stresses come. Furthermore, expression profiles of these genes in different organs and tissues under non-stress and abiotic stress were investigated through quantitative real-time PCR (qRT-PCR) analyses in cassava. Expression profiles exhibited clear differences among different tissues under non-stress and various dehydration stress conditions. We found that the leaf and tuberous root tissues had the greatest and least responses, respectively, to drought stress through the ethylene signaling pathway in cassava. Moreover, tuber and root tissues had the greatest and least reponses to osmotic and salt stresses through ethylene signaling in cassava, respectively. These results show that these plant tissues had differential expression levels of genes involved in ethylene signaling in response to the stresses tested. Moreover, after several gene duplication events, the spatiotemporally differential expression pattern of homologous genes in response to abiotic and biotic stresses may imply their functional diversity as a mechanism for adapting to the environment. Our data provide a framework for further research on the molecular mechanisms of cassava resistance to drought stress and provide a foundation for breeding drought-resistant new cultivars. [ABSTRACT FROM AUTHOR]