Current topics in signal transduction in bacteria

Among the signal transfer systems in bacteria two types predominate: two-component regulatory systems and quorum sensing systems. Both types of system can mediate signal transfer across the bacterial cell envelope; however, the signalling molecule typically is not taken up into the cells in the former type of system, whereas it usually is in the latter. The Two-component systems include the recently described (eukaryotic) phosphorelay systems; quorum sensing systems can be based upon autoinducers of the N-acylated homoserine lactones, and on autoinducers of a peptidic nature. A single bacterial cell contains many signalling modules that primarily operate in parallel. This may give rise to neural-network behaviour. Recently, however, for both types of basic signal transfer modules, it has been demonstrated that they also can be organised in series (i.e. in a hierarchical order). Besides their hierarchical position in the signal transduction network of the cell, the spatial distribution of individual signalling modules may also be an important factor in their efficiency in signal transfer. Many challenges lie hidden in future work to understand these signal transfer processes in more detail. These are discussed here, with emphasis on the mutual interactions between different signal transfer processes. Successful contributions to this work will require rigorous mathematical modelling of the performance of signal transduction components, and -networks, as well as studies on light-sensing signal transduction systems, because of the unsurpassed time resolution obtainable in those latter systems, the opportunity to apply repeated reproducible stimuli, etc. The increased understanding of bacterial behaviour that already has resulted – and may further result – from these studies, can be used to fine-tune the beneficial activities of bacteria and/or more efficiently inhibit their deleterious ones.     

LncRNA and mRNA integration network reconstruction reveals novel key regulators in esophageal squamous-cell carcinoma

Many experimental and computational studies have identified key protein coding genes in initiation and progression of esophageal squamous cell carcinoma (ESCC). However, the number of researches that tried to reveal the role of long non-coding RNAs (lncRNAs) in ESCC has been limited. LncRNAs are one of the important regulators of cancers which are transcribed dominantly in the genome and in various conditions. The main goal of this study was to use a systems biology approach to predict novel lncRNAs as well as protein coding genes associated with ESCC and assess their prognostic values. By using microarray expression data for mRNAs and lncRNAs from a large number of ESCC patients, we utilized “Weighted Gene Co-expression Network Analysis” (WGCNA) method to make a big coding-non-coding gene co-expression network, and discovered important functional modules. Gene set enrichment and pathway analysis revealed major biological processes and pathways involved in these modules. After selecting some protein coding genes involved in biological processes and pathways related to cancer, we used “LncTar”, a computational tool to predict potential interactions between these genes and lncRNAs. By combining interaction results with Pearson correlations, we introduced some novel lncRNAs with putative key regulatory roles in the network. Survival analysis with Kaplan-Meier estimator and Log-rank test statistic confirmed that most of the introduced genes are associated with poor prognosis in ESCC. Overall, our study reveals novel protein coding genes and lncRNAs associated with ESCC, along with their predicted interactions. Based on the promising results of survival analysis, these genes can be used as good estimators of patients' survival, or even can be analyzed further as new potential signatures or targets for the therapy of ESCC disease.     

Deciphering Ascorbic Acid Regulatory Pathways in Ripening Tomato Fruit Using a Weighted Gene Correlation Network Analysis Approach

Genotype is generally determined by the co-expression of diverse genes and multiple regulatory pathways in plants. Gene co-expression analysis combining with physiological trait data provides very important information about the gene function and regulatory mechanism. L-Ascorbic acid （AsA）, which is an essential nutrient component for human health and plant metabolism, plays key roles in diverse biological processes such as cell cycle, cell expansion, stress resistance, hormone synthesis, and signaling. Here, we applied a weighted gene correlation network analysis approach based on gene expression values and AsA content data in ripening tomato （Solanum lycopersicum L.） fruit with different AsA content levels, which leads to identification of AsA relevant modules and vital genes in AsA regulatory pathways. Twenty- four modules were compartmentalized according to gene expression profiling. Among these modules, one negatively related module containing genes involved in redox processes and one positively related module enriched with genes involved in AsA biosynthetic and recycling pathways were further analyzed. The present work herein indicates that redox pathways as well as hormone-signal pathways are closely correlated with AsA accumulation in ripening tomato fruit, and allowed us to prioritize candidate genes for follow-up studies to dissect this interplay at the biochemical and molecular level.     

Towards mechanistic models of plant organ growth

Modelling and simulation are increasingly used as tools in the study of plant growth and developmental processes. By formulating experimentally obtained knowledge as a system of interacting mathematical equations, it becomes feasible for biologists to gain a mechanistic understanding of the complex behaviour of biological systems. In this review, the modelling tools that are currently available and the progress that has been made to model plant development, based on experimental knowledge, are described. In terms of implementation, it is argued that, for the modelling of plant organ growth, the cellular level should form the cornerstone. It integrates the output of molecular regulatory networks to two processes, cell division and cell expansion, that drive growth and development of the organ. In turn, these cellular processes are controlled at the molecular level by hormone signalling. Therefore, combining a cellular modelling framework with regulatory modules for the regulation of cell division, expansion, and hormone signalling could form the basis of a functional organ growth simulation model. The current state of progress towards this aim is that the regulation of the cell cycle and hormone transport have been modelled extensively and these modules could be integrated. However, much less progress has been made on the modelling of cell expansion, which urgently needs to be addressed. A limitation of the current generation models is that they are largely qualitative. The possibilities to characterize existing and future models more quantitatively will be discussed. Together with experimental methods to measure crucial model parameters, these modelling techniques provide a basis to develop a Systems Biology approach to gain a fundamental insight into the relationship between gene function and whole organ behaviour.     

Reconstructing past ecological networks: the reconfiguration of seed-dispersal interactions after megafaunal extinction

The late Quaternary megafaunal extinction impacted ecological communities worldwide, and affected key ecological processes such as seed dispersal. The traits of several species of large-seeded plants are thought to have evolved in response to interactions with extinct megafauna, but how these extinctions affected the organization of interactions in seed-dispersal systems is poorly understood. Here, we combined ecological and paleontological data and network analyses to investigate how the structure of a species-rich seed-dispersal network could have changed from the Pleistocene to the present and examine the possible consequences of such changes. Our results indicate that the seed-dispersal network was organized into modules across the different time periods but has been reconfigured in different ways over time. The episode of megafaunal extinction and the arrival of humans changed how seed dispersers were distributed among network modules. However, the recent introduction of livestock into the seed-dispersal system partially restored the original network organization by strengthening the modular configuration. Moreover, after megafaunal extinctions, introduced species and some smaller native mammals became key components for the structure of the seed-dispersal network. We hypothesize that such changes in network structure affected both animal and plant assemblages, potentially contributing to the shaping of modern ecological communities. The ongoing extinction of key large vertebrates will lead to a variety of context-dependent rearranged ecological networks, most certainly affecting ecological and evolutionary processes.     

Protein complexes are central in the yeast genetic landscape

If perturbing two genes together has a stronger or weaker effect than expected, they are said to genetically interact. Genetic interactions are important because they help map gene function, and functionally related genes have similar genetic interaction patterns. Mapping quantitative (positive and negative) genetic interactions on a global scale has recently become possible. This data clearly shows groups of genes connected by predominantly positive or negative interactions, termed monochromatic groups. These groups often correspond to functional modules, like biological processes or complexes, or connections between modules. However it is not yet known how these patterns globally relate to known functional modules. Here we systematically study the monochromatic nature of known biological processes using the largest quantitative genetic interaction data set available, which includes fitness measurements for ～5.4 million gene pairs in the yeast Saccharomyces cerevisiae. We find that only 10% of biological processes, as defined by Gene Ontology annotations, and less than 1% of inter-process connections are monochromatic. Further, we show that protein complexes are responsible for a surprisingly large fraction of these patterns. This suggests that complexes play a central role in shaping the monochromatic landscape of biological processes. Altogether this work shows that both positive and negative monochromatic patterns are found in known biological processes and in their connections and that protein complexes play an important role in these patterns. The monochromatic processes, complexes and connections we find chart a hierarchical and modular map of sensitive and redundant biological systems in the yeast cell that will be useful for gene function prediction and comparison across phenotypes and organisms. Furthermore the analysis methods we develop are applicable to other species for which genetic interactions will progressively become more available.     

A literature-based similarity metric for biological processes

Background  

Recent analyses in systems biology pursue the discovery of functional modules within the cell. Recognition of such modules requires the integrative analysis of genome-wide experimental data together with available functional schemes. In this line, methods to bridge the gap between the abstract definitions of cellular processes in current schemes and the interlinked nature of biological networks are required.     

Key genes and functional coexpression modules involved in the pathogenesis of systemic lupus erythematosus

We performed a systematic review of genome‐wide gene expression datasets to identify key genes and functional modules involved in the pathogenesis of systemic lupus erythematosus (SLE) at a systems level. Genome‐wide gene expression datasets involving SLE patients were searched in Gene Expression Omnibus and ArrayExpress databases. Robust rank aggregation (RRA) analysis was used to integrate those public datasets and identify key genes associated with SLE. The weighted gene coexpression network analysis (WGCNA) was adapted to identify functional modules involved in SLE pathogenesis, and the gene ontology enrichment analysis was utilized to explore their functions. The aberrant expressions of several randomly selected key genes were further validated in SLE patients through quantitative real‐time polymerase chain reaction. Fifteen genome‐wide gene expression datasets were finally included, which involved a total of 1,778 SLE patients and 408 healthy controls. A large number of significantly upregulated or downregulated genes were identified through RRA analysis, and some of those genes were novel SLE gene signatures and their molecular roles in etiology of SLE remained vague. WGCNA further successfully identified six main functional modules involved in the pathogenesis of SLE. The most important functional module involved in SLE included 182 genes and mainly enriched in biological processes, including defense response to virus, interferon signaling pathway, and cytokine‐mediated signaling pathway. This study identifies a number of key genes and functional coexpression modules involved in SLE, which provides deepening insights into the molecular mechanism of SLE at a systems level and also provides some promising therapeutic targets.     

Characterizing the role benthos plays in large coastal seas and estuaries: A modular approach

Ecologists studying coastal and estuarine benthic communities have long taken a macroecological view, by relating benthic community patterns to environmental factors across several spatial scales. Although many general ecological patterns have been established, often a significant amount of the spatial and temporal variation in soft-sediment communities within and among systems remains unexplained. Here we propose a framework that may aid in unraveling the complex influence of environmental factors associated with the different components of coastal systems (i.e. the terrestrial and benthic landscapes, and the hydrological seascape) on benthic communities, and use this information to assess the role played by benthos in coastal ecosystems. A primary component of the approach is the recognition of system modules (e.g. marshes, dendritic systems, tidal rivers, enclosed basins, open bays, lagoons). The modules may differentially interact with key forcing functions (e.g. temperature, salinity, currents) that influence system processes and in turn benthic responses and functions. Modules may also constrain benthic characteristics and related processes within certain ecological boundaries and help explain their overall spatio-temporal variation. We present an example of how benthic community characteristics are related to the modular structure of 14 coastal seas and estuaries, and show that benthic functional group composition is significantly related to the modular structure of these systems. We also propose a framework for exploring the role of benthic communities in coastal systems using this modular approach and offer predictions of how benthic communities may vary depending on the modular composition and characteristics of a coastal system.     

Modularity in the evolution of yeast protein interaction network

Protein interaction networks are known to exhibit remarkable structures: scale-free and small-world and modular structures. To explain the evolutionary processes of protein interaction networks possessing scale-free and small-world structures, preferential attachment and duplication-divergence models have been proposed as mathematical models. Protein interaction networks are also known to exhibit another remarkable structural characteristic, modular structure. How the protein interaction networks became to exhibit modularity in their evolution? Here, we propose a hypothesis of modularity in the evolution of yeast protein interaction network based on molecular evolutionary evidence. We assigned yeast proteins into six evolutionary ages by constructing a phylogenetic profile. We found that all the almost half of hub proteins are evolutionarily new. Examining the evolutionary processes of protein complexes, functional modules and topological modules, we also found that member proteins of these modules tend to appear in one or two evolutionary ages. Moreover, proteins in protein complexes and topological modules show significantly low evolutionary rates than those not in these modules. Our results suggest a hypothesis of modularity in the evolution of yeast protein interaction network as systems evolution.     

Building the cellular puzzle: control in multi-level reaction networks

Quantitative conceptual tools dealing with control and regulation of cellular processes have been mostly developed for and applied to the pathways of intermediary metabolism. Yet, cellular processes are organized in different levels, metabolism forming the lowest level in a cascade of processes. Well-known examples are the DNA-mRNA-enzyme-metabolism cascade and the signal transduction cascades consisting of covalent modification cycles. The reaction network that constitutes each level can be viewed as a "module" in which reactions are linked by mass transfer. Although in principle all of these cellular modules are ultimately linked by mass transfer, in practice they can often be regarded as "isolated" from each other in terms of mass transfer. Here modules can interact with each other only by means of regulatory or catalytic effects-a chemical species in one module may affect the rate of a reaction in another module by binding to an enzyme or transport system or by acting as a catalyst. This paper seeks to answer two questions about the control and regulation of such multi-level reaction networks: (i) How can the control properties of the system as a whole be expressed in terms of the control properties of individual modules and the effects between modules? (ii) How do the control properties of a module in its isolated state change when it is embedded in the whole system through its connections with the other modules? In order to answer these questions a quantitative theoretical framework is developed and applied to systems containing two, three or four fully interacting modules; it is shown how it can be extended in principle to n modules. This newly developed theory therefore makes it possible to quantitatively dissect intermodular, internal and external regulation in multi-level systems.     

Systematic identification of common functional modules related to heart failure with different etiologies

The development of heart failure (HF) is a complex process that can be initiated by multiple etiologies. Identifying common functional modules associated with HF is a challenging task. Here, we developed a systems method to identify these common functional modules by integrating multiple expression profiles, protein interactions from four species, gene function annotations, and text information. We identified 1439 consistently differentially expressed genes (CDEGs) across HF with different etiologies by applying three meta-analysis methods to multiple HF-related expression profiles. Using a weighted human interaction network constructed by combining interaction data from multiple species, we extracted 60 candidate CDEG modules. We further evaluated the functional relevance of each module by using expression, interaction network, functional annotations, and text information together. Finally, five functional modules with significant biological relevance were identified. We found that almost half of the genes in these modules are hubs in the weighted network, and that these modules can accurately classify HF patients from healthy subjects. We also identified many significantly enriched biological processes that contribute to the pathophysiology of HF, including two new ones, RNA splicing and vesicle-mediated protein transport. In summary, we proposed a novel framework to analyze common functional modules related to HF with different etiologies. Our findings provide important insights into the complex mechanism of HF. Further biological experimentations should be required to validate these novel biological processes.     

Network cluster analysis of protein–protein interaction network identified biomarker for early onset colorectal cancer

Colorectal cancer (CRC) is a major cause of morbidity and mortality throughout the world. However, the genetic alterations and molecular mechanism of the early onset CRCs are not fully investigated. The present study aimed to characterize early onset CRC by analyzing its gene expression compared with normal controls and to identify network-based biomarkers of early onset CRC. The gene expression profiles of early onset CRC were downloaded from Gene Expression Omnibus and the differentially expressed genes (DEGs) in CRC patients were identified. Then, a protein–protein interaction (PPI) network was constructed and the clusters in PPI were analyzed by ClusterONE. Furthermore, the gene ontology functional analysis and pathway enrichment analysis were conducted to the modules in PPI network. A systems biology approach integrating microarray data and PPI was further applied to construct a PPI network in CRC. Total 631 DEGs were identified from the early onset CRC compared to healthy controls. These genes were found to be involved in several biological processes, including cell communication, cell proliferation, cell shape and apoptosis. Five functional modules which may play important roles in the initiation of early onset CRC were identified from the PPI network. Functional annotation revealed that these five modules were involved in the pathways of signal transduction, carcinogenesis and metastasis. The hub nodes of these five modules, CDC42, TEX11, QKI, CAV1 and FN1, may serve as the biomarkers of early onset CRC and have the potential to be targets for therapeutic intervention. However, further investigations are still needed to confirm our findings.     

Identifying responsive modules by mathematical programming: an application to budding yeast cell cycle

High-throughput biological data offer an unprecedented opportunity to fully characterize biological processes. However, how to extract meaningful biological information from these datasets is a significant challenge. Recently, pathway-based analysis has gained much progress in identifying biomarkers for some phenotypes. Nevertheless, these so-called pathway-based methods are mainly individual-gene-based or molecule-complex-based analyses. In this paper, we developed a novel module-based method to reveal causal or dependent relations between network modules and biological phenotypes by integrating both gene expression data and protein-protein interaction network. Specifically, we first formulated the identification problem of the responsive modules underlying biological phenotypes as a mathematical programming model by exploiting phenotype difference, which can also be viewed as a multi-classification problem. Then, we applied it to study cell-cycle process of budding yeast from microarray data based on our biological experiments, and identified important phenotype- and transition-based responsive modules for different stages of cell-cycle process. The resulting responsive modules provide new insight into the regulation mechanisms of cell-cycle process from a network viewpoint. Moreover, the identification of transition modules provides a new way to study dynamical processes at a functional module level. In particular, we found that the dysfunction of a well-known module and two new modules may directly result in cell cycle arresting at S phase. In addition to our biological experiments, the identified responsive modules were also validated by two independent datasets on budding yeast cell cycle.     

Modular organization of protein interaction networks

MOTIVATION: Accumulating evidence suggests that biological systems are composed of interacting, separable, functional modules. Identifying these modules is essential to understand the organization of biological systems. RESULT: In this paper, we present a framework to identify modules within biological networks. In this approach, the concept of degree is extended from the single vertex to the sub-graph, and a formal definition of module in a network is used. A new agglomerative algorithm was developed to identify modules from the network by combining the new module definition with the relative edge order generated by the Girvan-Newman (G-N) algorithm. A JAVA program, MoNet, was developed to implement the algorithm. Applying MoNet to the yeast core protein interaction network from the database of interacting proteins (DIP) identified 86 simple modules with sizes larger than three proteins. The modules obtained are significantly enriched in proteins with related biological process Gene Ontology terms. A comparison between the MoNet modules and modules defined by Radicchi et al. (2004) indicates that MoNet modules show stronger co-clustering of related genes and are more robust to ties in betweenness values. Further, the MoNet output retains the adjacent relationships between modules and allows the construction of an interaction web of modules providing insight regarding the relationships between different functional modules. Thus, MoNet provides an objective approach to understand the organization and interactions of biological processes in cellular systems. AVAILABILITY: MoNet is available upon request from the authors.     

PRO40 Is a Scaffold Protein of the Cell Wall Integrity Pathway,Linking the MAP Kinase Module to the Upstream Activator Protein Kinase C

Mitogen-activated protein kinase (MAPK) pathways are crucial signaling instruments in eukaryotes. Most ascomycetes possess three MAPK modules that are involved in key developmental processes like sexual propagation or pathogenesis. However, the regulation of these modules by adapters or scaffolds is largely unknown. Here, we studied the function of the cell wall integrity (CWI) MAPK module in the model fungus Sordaria macrospora. Using a forward genetic approach, we found that sterile mutant pro30 has a mutated mik1 gene that encodes the MAPK kinase kinase (MAPKKK) of the proposed CWI pathway. We generated single deletion mutants lacking MAPKKK MIK1, MAPK kinase (MAPKK) MEK1, or MAPK MAK1 and found them all to be sterile, cell fusion-deficient and highly impaired in vegetative growth and cell wall stress response. By searching for MEK1 interaction partners via tandem affinity purification and mass spectrometry, we identified previously characterized developmental protein PRO40 as a MEK1 interaction partner. Although fungal PRO40 homologs have been implicated in diverse developmental processes, their molecular function is currently unknown. Extensive affinity purification, mass spectrometry, and yeast two-hybrid experiments showed that PRO40 is able to bind MIK1, MEK1, and the upstream activator protein kinase C (PKC1). We further found that the PRO40 N-terminal disordered region and the central region encompassing a WW interaction domain are sufficient to govern interaction with MEK1. Most importantly, time- and stress-dependent phosphorylation studies showed that PRO40 is required for MAK1 activity. The sum of our results implies that PRO40 is a scaffold protein for the CWI pathway, linking the MAPK module to the upstream activator PKC1. Our data provide important insights into the mechanistic role of a protein that has been implicated in sexual and asexual development, cell fusion, symbiosis, and pathogenicity in different fungal systems.     

细胞黏附分子互作网络的构建与分析

借助网络分析可对基因调控、蛋白质互作和信号转导等细胞活动进行全局和局部性质分析.以细胞黏附的蛋白质相互作用为对象,通过数据挖掘和可视化软件构建了整合蛋白介导的黏附分子互作网络,该分子互作网络由156种蛋白质通过690种相互作用相连,其平均节点度为8.66、平均聚集系数为0.24,平均路径长度为2.6.黏附分子互作网络中包含数个功能模块,这些模块涉及网络内部多种分子相互作用的启动与停止,并进一步影响细胞的黏附、迁移和骨架组织.对黏附分子网络进行模体筛选和比较,发现一些数量相对较少、以三元复合物为主要结构的关键模体,同时对各网络模块和模体对细胞黏附的调控作用进行了探讨.     

The Modular Organization of Protein Interactions in Escherichia coli

Escherichia coli serves as an excellent model for the study of fundamental cellular processes such as metabolism, signalling and gene expression. Understanding the function and organization of proteins within these processes is an important step towards a ‘systems’ view of E. coli. Integrating experimental and computational interaction data, we present a reliable network of 3,989 functional interactions between 1,941 E. coli proteins (∼45% of its proteome). These were combined with a recently generated set of 3,888 high-quality physical interactions between 918 proteins and clustered to reveal 316 discrete modules. In addition to known protein complexes (e.g., RNA and DNA polymerases), we identified modules that represent biochemical pathways (e.g., nitrate regulation and cell wall biosynthesis) as well as batteries of functionally and evolutionarily related processes. To aid the interpretation of modular relationships, several case examples are presented, including both well characterized and novel biochemical systems. Together these data provide a global view of the modular organization of the E. coli proteome and yield unique insights into structural and evolutionary relationships in bacterial networks.