A network map of apelin-mediated signaling

The apelin receptor (APLNR) is a class A (rhodopsin-like) G-protein coupled receptor with a wide distribution throughout the human body. Activation of the apelin/APLNR system regulates AMPK/PI3K/AKT/mTOR and RAF/ERK1/2 mediated signaling pathways. APLNR activation orchestrates several downstream signaling cascades, which play diverse roles in physiological effects, including effects upon vasoconstriction, heart muscle contractility, energy metabolism regulation, and fluid homeostasis angiogenesis. We consolidated a network map of the APLNR signaling map owing to its biomedical importance. The curation of literature data pertaining to the APLNR system was performed manually by the NetPath criteria. The described apelin receptor signaling map comprises 35 activation/inhibition events, 38 catalysis events, 4 molecular associations, 62 gene regulation events, 113 protein expression types, and 4 protein translocation events. The APLNR signaling pathway map data is made freely accessible through the WikiPathways Database (https://www.wikipathways.org/index.php/Pathway:WP5067).Supplementary InformationThe online version contains supplementary material available at 10.1007/s12079-021-00614-6.     

SARS-CoV-2 signaling pathway map: A functional landscape of molecular mechanisms in COVID-19

Coronavirus disease (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). COVID-19 has been declared a pandemic by WHO. The clinical manifestation and disease progression in COVID-19 patients varies from minimal symptoms to severe respiratory issues with multiple organ failure. Understanding the mechanism of SARS-CoV-2 interaction with host cells will provide key insights into the effective molecular targets for the development of novel therapeutics. Recent studies have identified virus-mediated phosphorylation or activation of some major signaling pathways, such as ERK1/2, JNK, p38, PI3K/AKT and NF-κB signaling, that potentially elicit the cytokine storm that serves as a major cause of tissue injuries. Several studies highlight the aggressive inflammatory response particularly ‘cytokine storm’ in SARS-CoV-2 patients. A depiction of host molecular dynamics triggered by SARS-CoV-2 in the form of a network of signaling molecules will be helpful for COVID-19 research. Therefore, we developed the signaling pathway map of SARS-CoV-2 infection using data mined from the recently published literature. This integrated signaling pathway map of SARS-CoV-2 consists of 326 proteins and 73 reactions. These include information pertaining to 1,629 molecular association events, 30 enzyme catalysis events, 43 activation/inhibition events, and 8,531 gene regulation events. The pathway map is publicly available through WikiPathways: https://www.wikipathways.org/index.php/Pathway:WP5115.Supplementary InformationThe online version contains supplementary material available at 10.1007/s12079-021-00632-4.     

A modular map of Bradykinin-mediated inflammatory signaling network

Bradykinin, a member of the kallikrein-kinin system (KKS), is associated with an inflammatory response pathway with diverse vascular permeability functions, including thrombosis and blood coagulation. In majority, bradykinin signals through Bradykinin Receptor B2 (B2R). B2R is a G protein-coupled receptor (GPCR) coupled to G protein family such as Gα_qs, Gα_q/Gα_11, Gα_i1, and Gβ1_γ2_. B2R stimulation leads to the activation of a signaling cascade of downstream molecules such as phospholipases, protein kinase C, Ras/Raf-1/MAPK, and PI3K/AKT and secondary messengers such as inositol-1,4,5-trisphosphate, diacylglycerol and Ca²⁺ ions. These secondary messengers modulate the production of nitric oxide or prostaglandins. Bradykinin-mediated signaling is implicated in inflammation, chronic pain, vasculopathy, neuropathy, obesity, diabetes, and cancer. Despite the biomedical importance of bradykinin, a resource of bradykinin-mediated signaling pathway is currently not available. Here, we developed a pathway resource of signaling events mediated by bradykinin. By employing data mining strategies in the published literature, we describe an integrated pathway reaction map of bradykinin consisting of 233 reactions. Bradykinin signaling pathway events included 25 enzyme catalysis reactions, 12 translocations, 83 activation/inhibition reactions, 11 molecular associations, 45 protein expression and 57 gene regulation events. The pathway map is made publicly available on the WikiPathways Database with the ID URL: https://www.wikipathways.org/index.php/Pathway:WP5132. The bradykinin-mediated signaling pathway map will facilitate the identification of novel candidates as therapeutic targets for diseases associated with dysregulated bradykinin signaling.Supplementary InformationThe online version contains supplementary material available at 10.1007/s12079-021-00652-0.     

A complete map of the Calcium/calmodulin-dependent protein kinase kinase 2 (CAMKK2) signaling pathway

Calcium/calmodulin-dependent protein kinase kinase 2 (CAMKK2) is a serine/threonine-protein kinase belonging to the Ca²⁺/calmodulin-dependent protein kinase subfamily. CAMKK2 has an autocatalytic site, which gets exposed when Ca²⁺/calmodulin (CAM) binds to it. This results in autophosphorylation and complete activation of CAMKK2. The three major known downstream targets of CAMKK2 are 5′-adenosine monophosphate (AMP)-activated protein kinase (AMPKα), calcium/calmodulin-dependent protein kinase 1 (CAMK1) and calcium/calmodulin-dependent protein kinase 4 (CAMK4). Activation of these targets by CAMKK2 is important for the maintenance of different cellular and physiological processes within the cell. CAMKK2 is found to be important in neuronal development, bone remodeling, adipogenesis, and systemic glucose homeostasis, osteoclastgensis and postnatal myogensis. CAMKK2 is reported to be involved in pathologies like Duchenne muscular dystrophy, inflammation, osteoporosis and bone remodeling and is also reported to be overexpressed in prostate cancer, hepatic cancer, ovarian and gastric cancer. CAMKK2 is involved in increased cell proliferation and migration through CAMKK2/AMPK pathway in prostate cancer and activation of AKT in ovarian cancer. Although CAMKK2 is a molecule of great importance, a public resource of the CAMKK2 signaling pathway is currently lacking. Therefore, we carried out detailed data mining and documentation of the signaling events associated with CAMKK2 from published literature and developed an integrated reaction map of CAMKK2 signaling. This resulted in the cataloging of 285 reactions belonging to the CAMKK2 signaling pathway, which includes 33 protein–protein interactions, 74 post-translational modifications, 7 protein translocation events, and 22 activation/inhibition events. Besides, 124 gene regulation events and 25 activator/inhibitors involved in CAMKK2 activation were also cataloged. The CAMKK2 signaling pathway map data is made freely accessible through WikiPathway database (https://www.wikipathways.org/index.php/Pathway:WP4874). We expect that data on a signaling map of CAMKK2 will provide the scientific community with an improved platform to facilitate further molecular as well as biomedical investigations on CAMKK2 and its utility in the development of biomarkers and therapeutic targets.Electronic supplementary materialThe online version of this article (10.1007/s12079-020-00592-1) contains supplementary material, which is available to authorized users.     

Accelerating the Original Profile Kernel

One of the most accurate multi-class protein classification systems continues to be the profile-based SVM kernel introduced by the Leslie group. Unfortunately, its CPU requirements render it too slow for practical applications of large-scale classification tasks. Here, we introduce several software improvements that enable significant acceleration. Using various non-redundant data sets, we demonstrate that our new implementation reaches a maximal speed-up as high as 14-fold for calculating the same kernel matrix. Some predictions are over 200 times faster and render the kernel as possibly the top contender in a low ratio of speed/performance. Additionally, we explain how to parallelize various computations and provide an integrative program that reduces creating a production-quality classifier to a single program call. The new implementation is available as a Debian package under a free academic license and does not depend on commercial software. For non-Debian based distributions, the source package ships with a traditional Makefile-based installer. Download and installation instructions can be found at https://rostlab.org/owiki/index.php/Fast_Profile_Kernel. Bugs and other issues may be reported at https://rostlab.org/bugzilla3/enter_bug.cgi?product=fastprofkernel.     

Using the Semantic Web for Rapid Integration of WikiPathways with Other Biological Online Data Resources

The diversity of online resources storing biological data in different formats provides a challenge for bioinformaticians to integrate and analyse their biological data. The semantic web provides a standard to facilitate knowledge integration using statements built as triples describing a relation between two objects. WikiPathways, an online collaborative pathway resource, is now available in the semantic web through a SPARQL endpoint at http://sparql.wikipathways.org. Having biological pathways in the semantic web allows rapid integration with data from other resources that contain information about elements present in pathways using SPARQL queries. In order to convert WikiPathways content into meaningful triples we developed two new vocabularies that capture the graphical representation and the pathway logic, respectively. Each gene, protein, and metabolite in a given pathway is defined with a standard set of identifiers to support linking to several other biological resources in the semantic web. WikiPathways triples were loaded into the Open PHACTS discovery platform and are available through its Web API (https://dev.openphacts.org/docs) to be used in various tools for drug development. We combined various semantic web resources with the newly converted WikiPathways content using a variety of SPARQL query types and third-party resources, such as the Open PHACTS API. The ability to use pathway information to form new links across diverse biological data highlights the utility of integrating WikiPathways in the semantic web.     

Approximate Bayesian Computation

Approximate Bayesian computation (ABC) constitutes a class of computational methods rooted in Bayesian statistics. In all model-based statistical inference, the likelihood function is of central importance, since it expresses the probability of the observed data under a particular statistical model, and thus quantifies the support data lend to particular values of parameters and to choices among different models. For simple models, an analytical formula for the likelihood function can typically be derived. However, for more complex models, an analytical formula might be elusive or the likelihood function might be computationally very costly to evaluate. ABC methods bypass the evaluation of the likelihood function. In this way, ABC methods widen the realm of models for which statistical inference can be considered. ABC methods are mathematically well-founded, but they inevitably make assumptions and approximations whose impact needs to be carefully assessed. Furthermore, the wider application domain of ABC exacerbates the challenges of parameter estimation and model selection. ABC has rapidly gained popularity over the last years and in particular for the analysis of complex problems arising in biological sciences (e.g., in population genetics, ecology, epidemiology, and systems biology).

This is a “Topic Page” article for PLOS Computational Biology.

     

Evolving Digital Ecological Networks

“It is hard to realize that the living world as we know it is just one among many possibilities” [1]. Evolving digital ecological networks are webs of interacting, self-replicating, and evolving computer programs (i.e., digital organisms) that experience the same major ecological interactions as biological organisms (e.g., competition, predation, parasitism, and mutualism). Despite being computational, these programs evolve quickly in an open-ended way, and starting from only one or two ancestral organisms, the formation of ecological networks can be observed in real-time by tracking interactions between the constantly evolving organism phenotypes. These phenotypes may be defined by combinations of logical computations (hereafter tasks) that digital organisms perform and by expressed behaviors that have evolved. The types and outcomes of interactions between phenotypes are determined by task overlap for logic-defined phenotypes and by responses to encounters in the case of behavioral phenotypes. Biologists use these evolving networks to study active and fundamental topics within evolutionary ecology (e.g., the extent to which the architecture of multispecies networks shape coevolutionary outcomes, and the processes involved).

This is a “Topic Page” article for PLOS Computational Biology.

     

Viral Phylodynamics

Viral phylodynamics is defined as the study of how epidemiological, immunological, and evolutionary processes act and potentially interact to shape viral phylogenies. Since the coining of the term in 2004, research on viral phylodynamics has focused on transmission dynamics in an effort to shed light on how these dynamics impact viral genetic variation. Transmission dynamics can be considered at the level of cells within an infected host, individual hosts within a population, or entire populations of hosts. Many viruses, especially RNA viruses, rapidly accumulate genetic variation because of short generation times and high mutation rates. Patterns of viral genetic variation are therefore heavily influenced by how quickly transmission occurs and by which entities transmit to one another. Patterns of viral genetic variation will also be affected by selection acting on viral phenotypes. Although viruses can differ with respect to many phenotypes, phylodynamic studies have to date tended to focus on a limited number of viral phenotypes. These include virulence phenotypes, phenotypes associated with viral transmissibility, cell or tissue tropism phenotypes, and antigenic phenotypes that can facilitate escape from host immunity. Due to the impact that transmission dynamics and selection can have on viral genetic variation, viral phylogenies can therefore be used to investigate important epidemiological, immunological, and evolutionary processes, such as epidemic spread [2], spatio-temporal dynamics including metapopulation dynamics [3], zoonotic transmission, tissue tropism [4], and antigenic drift [5]. The quantitative investigation of these processes through the consideration of viral phylogenies is the central aim of viral phylodynamics.

This is a “Topic Page” article for PLOS Computational Biology.

     

Flow Cytometry Bioinformatics

Flow cytometry bioinformatics is the application of bioinformatics to flow cytometry data, which involves storing, retrieving, organizing, and analyzing flow cytometry data using extensive computational resources and tools. Flow cytometry bioinformatics requires extensive use of and contributes to the development of techniques from computational statistics and machine learning. Flow cytometry and related methods allow the quantification of multiple independent biomarkers on large numbers of single cells. The rapid growth in the multidimensionality and throughput of flow cytometry data, particularly in the 2000s, has led to the creation of a variety of computational analysis methods, data standards, and public databases for the sharing of results. Computational methods exist to assist in the preprocessing of flow cytometry data, identifying cell populations within it, matching those cell populations across samples, and performing diagnosis and discovery using the results of previous steps. For preprocessing, this includes compensating for spectral overlap, transforming data onto scales conducive to visualization and analysis, assessing data for quality, and normalizing data across samples and experiments. For population identification, tools are available to aid traditional manual identification of populations in two-dimensional scatter plots (gating), to use dimensionality reduction to aid gating, and to find populations automatically in higher dimensional space in a variety of ways. It is also possible to characterize data in more comprehensive ways, such as the density-guided binary space partitioning technique known as probability binning, or by combinatorial gating. Finally, diagnosis using flow cytometry data can be aided by supervised learning techniques, and discovery of new cell types of biological importance by high-throughput statistical methods, as part of pipelines incorporating all of the aforementioned methods. Open standards, data, and software are also key parts of flow cytometry bioinformatics. Data standards include the widely adopted Flow Cytometry Standard (FCS) defining how data from cytometers should be stored, but also several new standards under development by the International Society for Advancement of Cytometry (ISAC) to aid in storing more detailed information about experimental design and analytical steps. Open data is slowly growing with the opening of the CytoBank database in 2010 and FlowRepository in 2012, both of which allow users to freely distribute their data, and the latter of which has been recommended as the preferred repository for MIFlowCyt-compliant data by ISAC. Open software is most widely available in the form of a suite of Bioconductor packages, but is also available for web execution on the GenePattern platform.

This is a “Topic Page” article for PLOS Computational Biology.

     

Genetic Networks Inducing Invasive Growth in Saccharomyces cerevisiae Identified Through Systematic Genome-Wide Overexpression

The budding yeast Saccharomyces cerevisiae can respond to nutritional and environmental stress by implementing a morphogenetic program wherein cells elongate and interconnect, forming pseudohyphal filaments. This growth transition has been studied extensively as a model signaling system with similarity to processes of hyphal development that are linked with virulence in related fungal pathogens. Classic studies have identified core pseudohyphal growth signaling modules in yeast; however, the scope of regulatory networks that control yeast filamentation is broad and incompletely defined. Here, we address the genetic basis of yeast pseudohyphal growth by implementing a systematic analysis of 4909 genes for overexpression phenotypes in a filamentous strain of S. cerevisiae. Our results identify 551 genes conferring exaggerated invasive growth upon overexpression under normal vegetative growth conditions. This cohort includes 79 genes lacking previous phenotypic characterization. Pathway enrichment analysis of the gene set identifies networks mediating mitogen-activated protein kinase (MAPK) signaling and cell cycle progression. In particular, overexpression screening suggests that nuclear export of the osmoresponsive MAPK Hog1p may enhance pseudohyphal growth. The function of nuclear Hog1p is unclear from previous studies, but our analysis using a nuclear-depleted form of Hog1p is consistent with a role for nuclear Hog1p in repressing pseudohyphal growth. Through epistasis and deletion studies, we also identified genetic relationships with the G2 cyclin Clb2p and phenotypes in filamentation induced by S-phase arrest. In sum, this work presents a unique and informative resource toward understanding the breadth of genes and pathways that collectively constitute the molecular basis of filamentation.     

Mechanisms of Ephrin Receptor Protein Kinase-Independent Signaling in Amphid Axon Guidance in Caenorhabditis elegans

Eph receptors and their ephrin ligands are key conserved regulators of axon guidance and can function in a variety of signaling modes. Here we analyze the genetic and cellular requirements for Eph signaling in a Caenorhabditis elegans axon guidance choice point, the ventral guidance of axons in the amphid commissure. The C. elegans Eph receptor EFN-1 has both kinase-dependent and kinase-independent roles in amphid ventral guidance. Of the four C. elegans ephrins, we find that only EFN-1 has a major role in amphid axon ventral guidance, and signals in both a receptor kinase-dependent and kinase-independent manner. Analysis of EFN-1 and EFN-1 expression and tissue-specific requirements is consistent with a model in which VAB-1 acts in amphid neurons, interacting with EFN-1 expressed on surrounding cells. Unexpectedly, left-hand neurons are more strongly affected than right-hand neurons by loss of Eph signaling, indicating a previously undetected left–right asymmetry in the requirement for Eph signaling. By screening candidate genes involved in Eph signaling, we find that the Eph kinase-independent pathway involves the ABL-1 nonreceptor tyrosine kinase and possibly the phosphatidylinositol 3-kinase pathway. Overexpression of ABL-1 is sufficient to rescue EFN-1 ventral guidance defects cell autonomously. Our results reveal new aspects of Eph signaling in a single axon guidance decision in vivo.     

Coordination of Recombination with Meiotic Progression in the Caenorhabditis elegans Germline by KIN-18, a TAO Kinase That Regulates the Timing of MPK-1 Signaling

Meiosis is a tightly regulated process requiring coordination of diverse events. A conserved ERK/MAPK-signaling cascade plays an essential role in the regulation of meiotic progression. The Thousand And One kinase (TAO) kinase is a MAPK kinase kinase, the meiotic role of which is unknown. We have analyzed the meiotic functions of KIN-18, the homolog of mammalian TAO kinases, in Caenorhabditis elegans. We found that KIN-18 is essential for normal meiotic progression; mutants exhibit accelerated meiotic recombination as detected both by analysis of recombination intermediates and by crossover outcome. In addition, ectopic germ-cell differentiation and enhanced levels of apoptosis were observed in kin-18 mutants. These defects correlate with ectopic activation of MPK-1 that includes premature, missing, and reoccurring MPK-1 activation. Late progression defects in kin-18 mutants are suppressed by inhibiting an upstream activator of MPK-1 signaling, KSR-2. However, the acceleration of recombination events observed in kin-18 mutants is largely MPK-1-independent. Our data suggest that KIN-18 coordinates meiotic progression by modulating the timing of MPK-1 activation and the progression of recombination events. The regulation of the timing of MPK-1 activation ensures the proper timing of apoptosis and is required for the formation of functional oocytes. Meiosis is a conserved process; thus, revealing that KIN-18 is a novel regulator of meiotic progression in C. elegans would help to elucidate TAO kinase’s role in germline development in higher eukaryotes.     

ggtreeExtra: Compact Visualization of Richly Annotated Phylogenetic Data

We present the ggtreeExtra package for visualizing heterogeneous data with a phylogenetic tree in a circular or rectangular layout (https://www.bioconductor.org/packages/ggtreeExtra). The package supports more data types and visualization methods than other tools. It supports using the grammar of graphics syntax to present data on a tree with richly annotated layers and allows evolutionary statistics inferred by commonly used software to be integrated and visualized with external data. GgtreeExtra is a universal tool for tree data visualization. It extends the applications of the phylogenetic tree in different disciplines by making more domain-specific data to be available to visualize and interpret in the evolutionary context.     

State Transitions in the TORC1 Signaling Pathway and Information Processing in Saccharomyces cerevisiae

TOR kinase complex I (TORC1) is a key regulator of cell growth and metabolism in all eukaryotes. Previous studies in yeast have shown that three GTPases—Gtr1, Gtr2, and Rho1—bind to TORC1 in nitrogen and amino acid starvation conditions to block phosphorylation of the S6 kinase Sch9 and activate protein phosphatase 2A (PP2A). This leads to downregulation of 450 Sch9-dependent protein and ribosome synthesis genes and upregulation of 100 PP2A-dependent nitrogen assimilation and amino acid synthesis genes. Here, using bandshift assays and microarray measurements, we show that the TORC1 pathway also populates three other stress/starvation states. First, in glucose starvation conditions, the AMP-activated protein kinase (AMPK/Snf1) and at least one other factor push the TORC1 pathway into an off state, in which Sch9-branch signaling and PP2A-branch signaling are both inhibited. Remarkably, the TORC1 pathway remains in the glucose starvation (PP2A inhibited) state even when cells are simultaneously starved for nitrogen and glucose. Second, in osmotic stress, the MAPK Hog1/p38 drives the TORC1 pathway into a different state, in which Sch9 signaling and PP2A-branch signaling are inhibited, but PP2A-branch signaling can still be activated by nitrogen starvation. Third, in oxidative stress and heat stress, TORC1-Sch9 signaling is blocked while weak PP2A-branch signaling occurs. Together, our data show that the TORC1 pathway acts as an information-processing hub, activating different genes in different conditions to ensure that available energy is allocated to drive growth, amino acid synthesis, or a stress response, depending on the needs of the cell.     

Multi-state Modeling of Biomolecules

Multi-state modeling of biomolecules refers to a series of techniques used to represent and compute the behavior of biological molecules or complexes that can adopt a large number of possible functional states. Biological signaling systems often rely on complexes of biological macromolecules that can undergo several functionally significant modifications that are mutually compatible. Thus, they can exist in a very large number of functionally different states. Modeling such multi-state systems poses two problems: the problem of how to describe and specify a multi-state system (the “specification problem”) and the problem of how to use a computer to simulate the progress of the system over time (the “computation problem”). To address the specification problem, modelers have in recent years moved away from explicit specification of all possible states and towards rule-based formalisms that allow for implicit model specification, including the κ-calculus [1], BioNetGen [2]–[5], the Allosteric Network Compiler [6], and others [7], [8]. To tackle the computation problem, they have turned to particle-based methods that have in many cases proved more computationally efficient than population-based methods based on ordinary differential equations, partial differential equations, or the Gillespie stochastic simulation algorithm [9], [10]. Given current computing technology, particle-based methods are sometimes the only possible option. Particle-based simulators fall into two further categories: nonspatial simulators, such as StochSim [11], DYNSTOC [12], RuleMonkey [9], [13], and the Network-Free Stochastic Simulator (NFSim) [14], and spatial simulators, including Meredys [15], SRSim [16], [17], and MCell [18]–[20]. Modelers can thus choose from a variety of tools, the best choice depending on the particular problem. Development of faster and more powerful methods is ongoing, promising the ability to simulate ever more complex signaling processes in the future.

This is a “Topic Page” article for PLOS Computational Biology.

     

ModuleRole: A Tool for Modulization,Role Determination and Visualization in Protein-Protein Interaction Networks

Rapidly increasing amounts of (physical and genetic) protein-protein interaction (PPI) data are produced by various high-throughput techniques, and interpretation of these data remains a major challenge. In order to gain insight into the organization and structure of the resultant large complex networks formed by interacting molecules, using simulated annealing, a method based on the node connectivity, we developed ModuleRole, a user-friendly web server tool which finds modules in PPI network and defines the roles for every node, and produces files for visualization in Cytoscape and Pajek. For given proteins, it analyzes the PPI network from BioGRID database, finds and visualizes the modules these proteins form, and then defines the role every node plays in this network, based on two topological parameters Participation Coefficient and Z-score. This is the first program which provides interactive and very friendly interface for biologists to find and visualize modules and roles of proteins in PPI network. It can be tested online at the website http://www.bioinfo.org/modulerole/index.php, which is free and open to all users and there is no login requirement, with demo data provided by “User Guide” in the menu Help. Non-server application of this program is considered for high-throughput data with more than 200 nodes or user’s own interaction datasets. Users are able to bookmark the web link to the result page and access at a later time. As an interactive and highly customizable application, ModuleRole requires no expert knowledge in graph theory on the user side and can be used in both Linux and Windows system, thus a very useful tool for biologist to analyze and visualize PPI networks from databases such as BioGRID.

Availability

ModuleRole is implemented in Java and C, and is freely available at http://www.bioinfo.org/modulerole/index.php. Supplementary information (user guide, demo data) is also available at this website. API for ModuleRole used for this program can be obtained upon request.     

PathVisio 3: An Extendable Pathway Analysis Toolbox

PathVisio is a commonly used pathway editor, visualization and analysis software. Biological pathways have been used by biologists for many years to describe the detailed steps in biological processes. Those powerful, visual representations help researchers to better understand, share and discuss knowledge. Since the first publication of PathVisio in 2008, the original paper was cited more than 170 times and PathVisio was used in many different biological studies. As an online editor PathVisio is also integrated in the community curated pathway database WikiPathways.Here we present the third version of PathVisio with the newest additions and improvements of the application. The core features of PathVisio are pathway drawing, advanced data visualization and pathway statistics. Additionally, PathVisio 3 introduces a new powerful extension systems that allows other developers to contribute additional functionality in form of plugins without changing the core application.PathVisio can be downloaded from http://www.pathvisio.org and in 2014 PathVisio 3 has been downloaded over 5,500 times. There are already more than 15 plugins available in the central plugin repository. PathVisio is a freely available, open-source tool published under the Apache 2.0 license (http://www.apache.org/licenses/LICENSE-2.0). It is implemented in Java and thus runs on all major operating systems. The code repository is available at http://svn.bigcat.unimaas.nl/pathvisio. The support mailing list for users is available on https://groups.google.com/forum/#!forum/wikipathways-discuss and for developers on https://groups.google.com/forum/#!forum/wikipathways-devel.

This is a PLOS Computational Biology software article.