COVID-GWAB: A Web-Based Prediction of COVID-19 Host Genes via Network Boosting of Genome-Wide Association Data.

Host genetics affect both the susceptibility and response to viral infection. Searching for host genes that contribute to COVID-19, the Host Genetics Initiative (HGI) was formed to investigate the genetic factors involved in COVID-19 via genome-wide association studies (GWAS). The GWAS suffer from limited statistical power and in general, only a few genes can pass the conventional significance thresholds. This statistical limitation may be overcome by boosting weak association signals through integrating independent functional information such as molecular interactions. Additionally, the boosted results can be evaluated by various independent data for further connections to COVID-19. We present COVID-GWAB, a web-based tool to boost original GWAS signals from COVID-19 patients by taking the signals of the interactome neighbors. COVID-GWAB takes summary statistics from the COVID-19 HGI or user input data and reprioritizes candidate host genes for COVID-19 using HumanNet, a co-functional human gene network. The current version of COVID-GWAB provides the pre-processed data of releases 5, 6, and 7 of the HGI. Additionally, COVID-GWAB provides web interfaces for a summary of augmented GWAS signals, prediction evaluations by appearance frequency in COVID-19 literature, single-cell transcriptome data, and associated pathways. The web server also enables browsing the candidate gene networks.

Comparison of Postoperative Gastrointestinal Motility of Sugammadex and Neostigmine in Patients Undergoing Robotic Thyroidectomy: A Retrospective Study.

Postoperative bowel dysfunction poses difficulty to patients during their recovery from surgery, and reversal agents may affect bowel function. This study aimed to investigate and compare the effects of sugammadex and a neostigmine/glycopyrrolate combination on postoperative bowel movement in patients undergoing robotic thyroidectomy. The electronic medical records of 122 patients, who underwent robotic thyroidectomy between March 2018 and December 2020, were retrospectively reviewed. Demographic, clinical, and laboratory findings and the first gas-passing time after surgery were assessed. The number of patients with a first gas emission time over 24 h was significantly higher in the neostigmine group than in the sugammadex group (p = 0.008). Multivariate logistic regression analysis indicated that sugammadex was a prognostic factor for the first gas-passing time within 24 h (odds ratio = 4.60, 95% confidence interval 1.47-14.36, p = 0.005). Although postoperative bowel motility, based on the first gas emission time, was comparable, the number of patients with a first gas emission time within 24 h was significantly higher in the sugammadex group than in the neostigmine group. This shows that the use of sugammadex did not affect the delayed recovery of postoperative bowel motility after robotic thyroidectomy.

Effects of Remimazolam vs. Sevoflurane Anesthesia on Intraoperative Hemodynamics in Patients with Gastric Cancer Undergoing Robotic Gastrectomy: A Propensity Score-Matched Analysis.

Remimazolam has been suggested to improve the maintenance of hemodynamic stability when compared with other agents used for general anesthesia. This study aimed to compare the effects of remimazolam and sevoflurane anesthesia on hemodynamic stability in patients undergoing robotic gastrectomy. We retrospectively reviewed the electronic medical records of 199 patients who underwent robotic gastrectomy with sevoflurane (n = 135) or remimazolam (n = 64) anesthesia from January to November 2021. Propensity scores were used for 1:1 matching between the groups. The primary outcome was the difference in use of intraoperative vasopressors between groups. Secondary outcomes included differences in incidence and dose of vasopressors, as well as intraoperative hemodynamic variables, between groups. Remimazolam anesthesia was associated with a significantly less frequent use of ephedrine (odds ratio (OR): 0.13; 95% confidence interval (CI): 0.05−0.38, p < 0.001), phenylephrine (OR: 0.12; 95% CI: 0.04−0.40, p < 0.001), and any vasopressor (OR: 0.06; 95% CI: 0.02−0.25, p < 0.001) compared with sevoflurane anesthesia. Remimazolam anesthesia enables better maintenance of hemodynamic stability than sevoflurane anesthesia. Thus, remimazolam anesthesia may be beneficial for patients who are expected to experience hypotension due to the combined effects of CO2 pneumoperitoneum and the head-up position utilized during robotic gastrectomy.

HumanNet v3: an improved database of human gene networks for disease research.

Network medicine has proven useful for dissecting genetic organization of complex human diseases. We have previously published HumanNet, an integrated network of human genes for disease studies. Since the release of the last version of HumanNet, many large-scale protein-protein interaction datasets have accumulated in public depositories. Additionally, the numbers of research papers and functional annotations for gene-phenotype associations have increased significantly. Therefore, updating HumanNet is a timely task for further improvement of network-based research into diseases. Here, we present HumanNet v3 (https://www.inetbio.org/humannet/, covering 99.8% of human protein coding genes) constructed by means of the expanded data with improved network inference algorithms. HumanNet v3 supports a three-tier model: HumanNet-PI (a protein-protein physical interaction network), HumanNet-FN (a functional gene network), and HumanNet-XC (a functional network extended by co-citation). Users can select a suitable tier of HumanNet for their study purpose. We showed that on disease gene predictions, HumanNet v3 outperforms both the previous HumanNet version and other integrated human gene networks. Furthermore, we demonstrated that HumanNet provides a feasible approach for selecting host genes likely to be associated with COVID-19.

Effect of Dexmedetomidine on Biochemical Recurrence in Patients after Robot-Assisted Laparoscopic Radical Prostatectomy: A Retrospective Study.

The usage of dexmedetomidine during cancer surgery in current clinical practice is debatable, largely owing to the differing reports of its efficacy based on cancer type. This study aimed to investigate the effects of dexmedetomidine on biochemical recurrence (BCR) and radiographic progression in patients with prostate cancer, who have undergone robot-assisted laparoscopic radical prostatectomy (RALP). Using follow-up data from two prospective randomized controlled studies, BCR and radiographic progression were compared between individuals who received dexmedetomidine (n = 58) and those who received saline (n = 56). Patients with complete follow-up records between July 2013 and June 2019 were enrolled in this study. There were no significant between-group differences in the number of patients who developed BCR and those who showed positive radiographic progression. Based on the Cox regression analysis, age (p = 0.015), Gleason score ≥ 8 (p < 0.001), and pathological tumor stage 3a and 3b (both p < 0.001) were shown to be significant predictors of post-RALP BCR. However, there was no impact on the dexmedetomidine or control groups. Low-dose administration of dexmedetomidine at a rate of 0.3-0.4 µg/kg/h did not significantly affect BCR incidence following RALP. In addition, no beneficial effect was noted on radiographic progression.

Human reference gut microbiome catalog including newly assembled genomes from under-represented Asian metagenomes.

BACKGROUND: Metagenome sampling bias for geographical location and lifestyle is partially responsible for the incomplete catalog of reference genomes of gut microbial species. Thus, genome assembly from currently under-represented populations may effectively expand the reference gut microbiome and improve taxonomic and functional profiling. METHODS: We assembled genomes using public whole-metagenomic shotgun sequencing (WMS) data for 110 and 645 fecal samples from India and Japan, respectively. In addition, we assembled genomes from newly generated WMS data for 90 fecal samples collected from Korea. Expecting genome assembly for low-abundance species may require a much deeper sequencing than that usually employed, so we performed ultra-deep WMS (> 30 Gbp or > 100 million read pairs) for the fecal samples from Korea. We consequently assembled 29,082 prokaryotic genomes from 845 fecal metagenomes for the three under-represented Asian countries and combined them with the Unified Human Gastrointestinal Genome (UHGG) to generate an expanded catalog, the Human Reference Gut Microbiome (HRGM). RESULTS: HRGM contains 232,098 non-redundant genomes for 5414 representative prokaryotic species including 780 that are novel, > 103 million unique proteins, and > 274 million single-nucleotide variants. This is an over 10% increase from the UHGG. The new 780 species were enriched for the Bacteroidaceae family, including species associated with high-fiber and seaweed-rich diets. Single-nucleotide variant density was positively associated with the speciation rate of gut commensals. We found that ultra-deep sequencing facilitated the assembly of genomes for low-abundance taxa, and deep sequencing (e.g., > 20 million read pairs) may be needed for the profiling of low-abundance taxa. Importantly, the HRGM significantly improved the taxonomic and functional classification of sequencing reads from fecal samples. Finally, analysis of human self-antigen homologs on the HRGM species genomes suggested that bacterial taxa with high cross-reactivity potential may contribute more to the pathogenesis of gut microbiome-associated diseases than those with low cross-reactivity potential by promoting inflammatory condition. CONCLUSIONS: By including gut metagenomes from previously under-represented Asian countries, Korea, India, and Japan, we developed a substantially expanded microbiome catalog, HRGM. Information of the microbial genomes and coding genes is publicly available ( www.mbiomenet.org/HRGM/ ). HRGM will facilitate the identification and functional analysis of disease-associated gut microbiota.

BarleyNet: A Network-Based Functional Omics Analysis Server for Cultivated Barley, Hordeum vulgare L.

Cultivated barley (Hordeum vulgare L.) is one of the most produced cereal crops worldwide after maize, bread wheat, and rice. Barley is an important crop species not only as a food source, but also in plant genetics because it harbors numerous stress response alleles in its genome that can be exploited for crop engineering. However, the functional annotation of its genome is relatively poor compared with other major crops. Moreover, bioinformatics tools for system-wide analyses of omics data from barley are not yet available. We have thus developed BarleyNet, a co-functional network of 26,145 barley genes, along with a web server for network-based predictions (http://www.inetbio.org/barleynet). We demonstrated that BarleyNet's prediction of biological processes is more accurate than that of an existing barley gene network. We implemented three complementary network-based algorithms for prioritizing genes or functional concepts to study genetic components of complex traits such as environmental stress responses: (i) a pathway-centric search for candidate genes of pathways or complex traits; (ii) a gene-centric search to infer novel functional concepts for genes; and (iii) a context-centric search for novel genes associated with stress response. We demonstrated the usefulness of these network analysis tools in the study of stress response using proteomics and transcriptomics data from barley leaves and roots upon drought or heat stresses. These results suggest that BarleyNet will facilitate our understanding of the underlying genetic components of complex traits in barley.

Understanding the immunopathogenesis of autoimmune diseases by animal studies using gene modulation: A comprehensive review.

Autoimmune diseases are clinical syndromes that result from pathogenic inflammatory responses driven by inadequate immune activation by T- and B-cells. Although the exact mechanisms of autoimmune diseases are still elusive, genetic factors also play an important role in the pathogenesis. Recently, with the advancement of understanding of the immunological and molecular basis of autoimmune diseases, gene modulation has become a potential approach for the tailored treatment of autoimmune disorders. Gene modulation can be applied to regulate the levels of interleukins (IL), tumor necrosis factor (TNF), cytotoxic T-lymphocyte-associated antigen 4 (CTLA-4), interferon-γ and other inflammatory cytokines by inhibiting these cytokine expressions using short interfering ribonucleic acid (siRNA) or by inhibiting cytokine signaling using small molecules. In addition, gene modulation delivering anti-inflammatory cytokines or cytokine antagonists showed effectiveness in regulating autoimmunity. In this review, we summarize the potential target genes for gene or immunomodulation in autoimmune diseases including rheumatoid arthritis (RA), systemic lupus erythematosus (SLE), inflammatory bowel diseases (IBD) and multiple sclerosis (MS). This article will give a new perspective on understanding immunopathogenesis of autoimmune diseases not only in animals but also in human. Emerging approaches to investigate cytokine regulation through gene modulation may be a potential approach for the tailored immunomodulation of some autoimmune diseases near in the future.

VirtualCytometry: a webserver for evaluating immune cell differentiation using single-cell RNA sequencing data.

MOTIVATION: The immune system has diverse types of cells that are differentiated or activated via various signaling pathways and transcriptional regulation upon challenging conditions. Immunophenotyping by flow and mass cytometry are the major approaches for identifying key signaling molecules and transcription factors directing the transition between the functional states of immune cells. However, few proteins can be evaluated by flow cytometry in a single experiment, preventing researchers from obtaining a comprehensive picture of the molecular programs involved in immune cell differentiation. Recent advances in single-cell RNA sequencing (scRNA-seq) have enabled unbiased genome-wide quantification of gene expression in individual cells on a large scale, providing a new and versatile analytical pipeline for studying immune cell differentiation. RESULTS: We present VirtualCytometry, a web-based computational pipeline for evaluating immune cell differentiation by exploiting cell-to-cell variation in gene expression with scRNA-seq data. Differentiating cells often show a continuous spectrum of cellular states rather than distinct populations. VirtualCytometry enables the identification of cellular subsets for different functional states of differentiation based on the expression of marker genes. Case studies have highlighted the usefulness of this subset analysis strategy for discovering signaling molecules and transcription factors for human T-cell exhaustion, a state of T-cell dysfunction, in tumor and mouse dendritic cells activated by pathogens. With more than 226 scRNA-seq datasets precompiled from public repositories covering diverse mouse and human immune cell types in normal and disease tissues, VirtualCytometry is a useful resource for the molecular dissection of immune cell differentiation. AVAILABILITY AND IMPLEMENTATION: www.grnpedia.org/cytometry.

BiomeNet: a database for construction and analysis of functional interaction networks for any species with a sequenced genome.

MOTIVATION: Owing to advanced DNA sequencing and genome assembly technology, the number of species with sequenced genomes is rapidly increasing. The aim of the recently launched Earth BioGenome Project is to sequence genomes of all eukaryotic species on Earth over the next 10 years, making it feasible to obtain genomic blueprints of the majority of animal and plant species by this time. Genetic models of the sequenced species will later be subject to functional annotation, and a comprehensive molecular network should facilitate functional analysis of individual genes and pathways. However, network databases are lagging behind genome sequencing projects as even the largest network database provides gene networks for less than 10% of sequenced eukaryotic genomes, and the knowledge gap between genomes and interactomes continues to widen. RESULTS: We present BiomeNet, a database of 95 scored networks comprising over 8 million co-functional links, which can build and analyze gene networks for any species with the sequenced genome. BiomeNet transfers functional interactions between orthologous proteins from source networks to the target species within minutes and automatically constructs gene networks with the quality comparable to that of existing networks. BiomeNet enables assembly of the first-in-species gene networks not available through other databases, which are highly predictive of diverse biological processes and can also provide network analysis by extracting subnetworks for individual biological processes and network-based gene prioritizations. These data indicate that BiomeNet could enhance the benefits of decoding the genomes of various species, thus improving our understanding of the Earth' biodiversity. AVAILABILITY AND IMPLEMENTATION: The BiomeNet is freely available at http://kobic.re.kr/biomenet/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Network Integrative Genomic and Transcriptomic Analysis of Carbapenem-Resistant Klebsiella pneumoniae Strains Identifies Genes for Antibiotic Resistance and Virulence.

Global increases in the use of carbapenems have resulted in several strains of Gram-negative bacteria acquiring carbapenem resistance, thereby limiting treatment options. Klebsiella pneumoniae is a common carbapenem-resistant pathogenic bacterium that is widely studied to identify novel antibiotic resistance mechanisms and drug targets. Antibiotic-resistant clinical isolates generally harbor many genetic alterations, and the identification of responsible mutations would provide insights into the molecular mechanisms of antibiotic resistance. We propose a method to prioritize mutated genes responsible for antibiotic resistance on the basis of expression changes in their local subnetworks, hypothesizing that mutated genes that show significant expression changes among the corresponding functionally associated genes are more likely to be involved in the carbapenem resistance. For network-based gene prioritization, we developed KlebNet (www.inetbio.org/klebnet), a genome-scale cofunctional network of K. pneumoniae genes. Using KlebNet, we reconstructed the functional modules for carbapenem resistance and virulence and identified the functional association between antibiotic resistance and virulence. Using complementation assays with the top candidate genes, we were able to validate a novel gene that negatively regulated carbapenem resistance and four novel genes that positively regulated virulence in Galleria mellonella larvae. Therefore, our study demonstrated the feasibility of network-based identification of genes required for antibiotic resistance and virulence of human-pathogenic bacteria.IMPORTANCE Klebsiella pneumoniae is a major bacterial pathogen that causes pneumonia and urinary tract infections in human. K. pneumoniae infections are treated with carbapenem, but carbapenem-resistant K. pneumoniae has been spreading worldwide. We are able to identify antimicrobial-resistant genes among mutated genes of the antibiotic-resistant clinical isolates. However, they usually harbor many mutated genes, including those that cause weak or neutral functional effects. Therefore, we need to prioritize the mutated genes to identify the more likely candidates for the follow-up functional analysis. For this study, we present a functional network of K. pneumoniae genes and propose a network-based method of prioritizing the mutated genes of the resistant clinical isolates. We also reconstructed the network-based functional modules for carbapenem resistance and virulence and retrieved the functional association between antibiotic resistance and virulence. This study demonstrated the feasibility of network-based analysis of clinical genomics data for the study of K. pneumoniae infection.

MaizeNet: a co-functional network for network-assisted systems genetics in Zea mays.

Maize (Zea mays) has multiple uses in human food, animal fodder, starch and sweetener production and as a biofuel, and is accordingly the most extensively cultivated cereal worldwide. To enhance maize production, genetic factors underlying important agricultural traits, including stress tolerance and flowering, have been explored through forward and reverse genetics approaches. Co-functional gene networks are systems biology resources useful in identifying trait-associated genes in plants by prioritizing candidate genes. Here, we present MaizeNet (http://www.inetbio.org/maizenet/), a genome-scale co-functional network of Z. mays genes, and a companion web server for network-assisted systems genetics. We describe the validation of MaizeNet network quality and its ability to functionally predict molecular pathways and complex traits in maize. Furthermore, we demonstrate that MaizeNet-based prioritization of candidate genes can facilitate the identification of cell wall biosynthesis genes and detect network communities associated with flowering-time candidate genes derived from genome-wide association studies. The demonstrated gene prioritization and subnetwork analysis can be conducted by simply submitting maize gene models based on the commonly used B73 RefGen_v3 and the latest B73 RefGen_v4 reference genomes on the MaizeNet web server. MaizeNet-based network-assisted systems genetics will substantially accelerate the discovery of trait-associated genes for crop improvement.

A Genome-Scale Co-Functional Network of Xanthomonas Genes Can Accurately Reconstruct Regulatory Circuits Controlled by Two-Component Signaling Systems.

Bacterial species in the genus Xanthomonas infect virtually all crop plants. Although many genes involved in Xanthomonas virulence have been identified through molecular and cellular studies, the elucidation of virulence-associated regulatory circuits is still far from complete. Functional gene networks have proven useful in generating hypotheses for genetic factors of biological processes in various species. Here, we present a genome-scale co-functional network of Xanthomonas oryze pv. oryzae (Xoo) genes, XooNet (www.inetbio.org/xoonet/), constructed by integrating heterogeneous types of genomics data derived from Xoo and other bacterial species. XooNet contains 106,000 functional links, which cover approximately 83% of the coding genome. XooNet is highly predictive for diverse biological processes in Xoo and can accurately reconstruct cellular pathways regulated by two-component signaling transduction systems (TCS). XooNet will be a useful in silico research platform for genetic dissection of virulence pathways in Xoo.

HumanNet v2: human gene networks for disease research.

Human gene networks have proven useful in many aspects of disease research, with numerous network-based strategies developed for generating hypotheses about gene-disease-drug associations. The ability to predict and organize genes most relevant to a specific disease has proven especially important. We previously developed a human functional gene network, HumanNet, by integrating diverse types of omics data using Bayesian statistics framework and demonstrated its ability to retrieve disease genes. Here, we present HumanNet v2 (http://www.inetbio.org/humannet), a database of human gene networks, which was updated by incorporating new data types, extending data sources and improving network inference algorithms. HumanNet now comprises a hierarchy of human gene networks, allowing for more flexible incorporation of network information into studies. HumanNet performs well in ranking disease-linked gene sets with minimal literature-dependent biases. We observe that incorporating model organisms' protein-protein interactions does not markedly improve disease gene predictions, suggesting that many of the disease gene associations are now captured directly in human-derived datasets. With an improved interactive user interface for disease network analysis, we expect HumanNet will be a useful resource for network medicine.

TRRUST v2: an expanded reference database of human and mouse transcriptional regulatory interactions.

Transcription factors (TFs) are major trans-acting factors in transcriptional regulation. Therefore, elucidating TF-target interactions is a key step toward understanding the regulatory circuitry underlying complex traits such as human diseases. We previously published a reference TF-target interaction database for humans-TRRUST (Transcriptional Regulatory Relationships Unraveled by Sentence-based Text mining)-which was constructed using sentence-based text mining, followed by manual curation. Here, we present TRRUST v2 (www.grnpedia.org/trrust) with a significant improvement from the previous version, including a significantly increased size of the database consisting of 8444 regulatory interactions for 800 TFs in humans. More importantly, TRRUST v2 also contains a database for TF-target interactions in mice, including 6552 TF-target interactions for 828 mouse TFs. TRRUST v2 is also substantially more comprehensive and less biased than other TF-target interaction databases. We also improved the web interface, which now enables prioritization of key TFs for a physiological condition depicted by a set of user-input transcriptional responsive genes. With the significant expansion in the database size and inclusion of the new web tool for TF prioritization, we believe that TRRUST v2 will be a versatile database for the study of the transcriptional regulation involved in human diseases.

GWAB: a web server for the network-based boosting of human genome-wide association data.

During the last decade, genome-wide association studies (GWAS) have represented a major approach to dissect complex human genetic diseases. Due in part to limited statistical power, most studies identify only small numbers of candidate genes that pass the conventional significance thresholds (e.g. P ≤ 5 × 10-8). This limitation can be partly overcome by increasing the sample size, but this comes at a higher cost. Alternatively, weak association signals can be boosted by incorporating independent data. Previously, we demonstrated the feasibility of boosting GWAS disease associations using gene networks. Here, we present a web server, GWAB (www.inetbio.org/gwab), for the network-based boosting of human GWAS data. Using GWAS summary statistics (P-values) for SNPs along with reference genes for a disease of interest, GWAB reprioritizes candidate disease genes by integrating the GWAS and network data. We found that GWAB could more effectively retrieve disease-associated reference genes than GWAS could alone. As an example, we describe GWAB-boosted candidate genes for coronary artery disease and supporting data in the literature. These results highlight the inherent value in sub-threshold GWAS associations, which are often not publicly released. GWAB offers a feasible general approach to boost such associations for human disease genetics.

COEXPEDIA: exploring biomedical hypotheses via co-expressions associated with medical subject headings (MeSH).

The use of high-throughput array and sequencing technologies has produced unprecedented amounts of gene expression data in central public depositories, including the Gene Expression Omnibus (GEO). The immense amount of expression data in GEO provides both vast research opportunities and data analysis challenges. Co-expression analysis of high-dimensional expression data has proven effective for the study of gene functions, and several co-expression databases have been developed. Here, we present a new co-expression database, COEXPEDIA (www.coexpedia.org), which is distinctive from other co-expression databases in three aspects: (i) it contains only co-functional co-expressions that passed a rigorous statistical assessment for functional association, (ii) the co-expressions were inferred from individual studies, each of which was designed to investigate gene functions with respect to a particular biomedical context such as a disease and (iii) the co-expressions are associated with medical subject headings (MeSH) that provide biomedical information for anatomical, disease, and chemical relevance. COEXPEDIA currently contains approximately eight million co-expressions inferred from 384 and 248 GEO series for humans and mice, respectively. We describe how these MeSH-associated co-expressions enable the identification of diseases and drugs previously unknown to be related to a gene or a gene group of interest.

Function-driven discovery of disease genes in zebrafish using an integrated genomics big data resource.

Whole exome sequencing (WES) accelerates disease gene discovery using rare genetic variants, but further statistical and functional evidence is required to avoid false-discovery. To complement variant-driven disease gene discovery, here we present function-driven disease gene discovery in zebrafish (Danio rerio), a promising human disease model owing to its high anatomical and genomic similarity to humans. To facilitate zebrafish-based function-driven disease gene discovery, we developed a genome-scale co-functional network of zebrafish genes, DanioNet (www.inetbio.org/danionet), which was constructed by Bayesian integration of genomics big data. Rigorous statistical assessment confirmed the high prediction capacity of DanioNet for a wide variety of human diseases. To demonstrate the feasibility of the function-driven disease gene discovery using DanioNet, we predicted genes for ciliopathies and performed experimental validation for eight candidate genes. We also validated the existence of heterozygous rare variants in the candidate genes of individuals with ciliopathies yet not in controls derived from the UK10K consortium, suggesting that these variants are potentially involved in enhancing the risk of ciliopathies. These results showed that an integrated genomics big data for a model animal of diseases can expand our opportunity for harnessing WES data in disease gene discovery.

Network-assisted investigation of virulence and antibiotic-resistance systems in Pseudomonas aeruginosa.

Pseudomonas aeruginosa is a Gram-negative bacterium of clinical significance. Although the genome of PAO1, a prototype strain of P. aeruginosa, has been extensively studied, approximately one-third of the functional genome remains unknown. With the emergence of antibiotic-resistant strains of P. aeruginosa, there is an urgent need to develop novel antibiotic and anti-virulence strategies, which may be facilitated by an approach that explores P. aeruginosa gene function in systems-level models. Here, we present a genome-wide functional network of P. aeruginosa genes, PseudomonasNet, which covers 98% of the coding genome, and a companion web server to generate functional hypotheses using various network-search algorithms. We demonstrate that PseudomonasNet-assisted predictions can effectively identify novel genes involved in virulence and antibiotic resistance. Moreover, an antibiotic-resistance network based on PseudomonasNet reveals that P. aeruginosa has common modular genetic organisations that confer increased or decreased resistance to diverse antibiotics, which accounts for the pervasiveness of cross-resistance across multiple drugs. The same network also suggests that P. aeruginosa has developed mechanism of trade-off in resistance across drugs by altering genetic interactions. Taken together, these results clearly demonstrate the usefulness of a genome-scale functional network to investigate pathogenic systems in P. aeruginosa.