Use of a combined antibacterial synergy approach and the ANNOgesic tool to identify novel targets within the gene networks of multidrug-resistant Klebsiella pneumoniae.

Since the 1980s, the development of new drug classes for the treatment of multidrug-resistant Klebsiella pneumoniae has become limited, highlighting the urgent need for novel antibiotics. To address this challenge, this study aimed to explore the synergistic interactions between chemical compounds and representative antibiotics, such as carbapenem and colistin. The primary objective of this study was not only to mitigate the adverse impact of multidrug-resistant K. pneumoniae on public health but also to establish a sustainable balance among humans, animals, and the environment. Phenotypical measurements were conducted using the broth microdilution technique to determine the drug sensitivity of bacterial strains. Additionally, a genotypical approach was employed, involving traditional RNA sequencing analysis to identify differentially expressed genes and the computational ANNOgesic tool to detect noncoding RNAs. This study revealed the existence of various pathways and regulatory RNA elements that form a functional network. These pathways, characterized by the expression of specific genes, contribute to the combined treatment effect and bacterial survival strategies. The connections between pathways are facilitated by regulatory RNA elements that respond to environmental changes. These findings suggest an adaptive response of bacteria to harsh environmental conditions.IMPORTANCENoncoding RNAs were identified as key players in post-transcriptional regulation. Moreover, this study predicted the presence of novel small regulatory RNAs that interact with target genes, as well as the involvement of riboswitches and RNA thermometers in conjunction with associated genes. These findings will contribute to the discovery of potential antimicrobial therapeutic candidates. Overall, this study offers valuable insights into the synergistic effects of chemical compounds and antibiotics, highlighting the role of regulatory RNA elements in bacterial response, and survival strategies. The identification of novel noncoding RNAs and their interactions with target genes, riboswitches, and RNA thermometers holds promise for the development of antimicrobial therapies.

WNT5A drives interleukin-6-dependent epithelial-mesenchymal transition via the JAK/STAT pathway in keloid pathogenesis.

Background: Keloid scarring is a fibroproliferative disease caused by aberrant genetic activation with an unclear underlying mechanism. Genetic predisposition, aberrant cellular responses to environmental factors, increased inflammatory cytokines and epithelial-mesenchymal transition (EMT) phenomena are known as major contributors. In this study, we aimed to identify the molecular drivers that initiate keloid pathogenesis. Methods: Bulk tissue RNA sequencing analyses of keloid and normal tissues along with ex vivo and in vitro tests were performed to identify the contributing genes to keloid pathogenesis. An animal model of inflammatory keloid scarring was reproduced by replication of a skin fibrosis model with intradermal bleomycin injection in C57BL/6 mice. Results: Gene set enrichment analysis revealed upregulation of Wnt family member 5A (WNT5A) expression and genes associated with EMT in keloid tissues. Consistently, human keloid tissues and the bleomycin-induced skin fibrosis animal model showed significantly increased expression of WNT5A and EMT markers. Increased activation of the interleukin (IL)-6/Janus kinase (JAK)/signal transducer and activator of transcription (STAT) pathway and subsequent elevation of EMT markers was also observed in keratinocytes co-cultured with WNT5A-activated fibroblasts or keloid fibroblasts. Furthermore, WNT5A silencing and the blockage of IL-6 secretion via neutralizing IL-6 antibody reversed hyperactivation of the STAT pathway and EMT markers in keratinocytes. Lastly, STAT3 silencing significantly reduced the EMT-like phenotypes in both keratinocytes and IL-6-stimulated keratinocytes. Conclusions: Intercellular communication via the WNT5A and STAT pathways possibly underlies a partial mechanism of EMT-like phenomena in keloid pathogenesis. IL-6 secreted from WNT5A-activated fibroblasts or keloid fibroblasts activates the JAK/STAT signaling pathway in adjacent keratinocytes which in turn express EMT markers. A better understanding of keloid development and the role of WNT5A in EMT will promote the development of next-generation targeted treatments for keloid scars.

RNA sequencing and bioinformatics analysis revealed PACSIN3 as a potential novel biomarker for platinum resistance in epithelial ovarian cancer.

BACKGROUND: Failure to respond to treatment in epithelial ovarian cancer can often be attributed to platinum-based chemotherapy resistance. However, the possible mechanisms or candidate biomarkers associated with platinum resistance are yet to be elucidated, even though many researchers have performed related studies. METHODS: We performed RNA sequencing of clinical specimens obtained from patients with platinum-sensitive or resistant epithelial ovarian cancer (EOC). Furthermore, various bioinformatics approaches, including spatial analysis of functional enrichment, were used to identify key regulators and associated underlying mechanisms of platinum resistance in EOC. RESULTS: Through RNA-sequencing, we identified 263 differentially expressed genes (98 upregulated and 165 downregulated) and subjected them to Gene Oncology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses, which were characterized to the traditional platinum-resistant characteristics. Subsequently, the gene interaction network and module analysis by spatial analysis of functional enrichment software demonstrated protein kinase C and casein kinase substrate in neurons 3 (PACSIN3) as the only upregulated hub gene, and neurotensin (NTS) and KIAA0319 as downregulated hub genes in platinum-resistant EOC. We selected PACSIN3 for further analysis because it has not been studied in relation to response to platinum-based chemotherapy. PACSIN3 was significantly upregulated in ovarian cancer cells compared to immortalized human ovarian surface epithelial cells. In addition, cisplatin-induced apoptosis was measured in PACSIN3 knockout OVCA433 and BRCA-mutated EOC cell line, SNU251, by a fluorescence-activated cell sorting-based Annexin-V/propium iodide double staining assay, which revealed a significant increase in apoptosis. CONCLUSIONS: Taken together, the present study presents PACSIN3 as a promising predictive biomarker associated with platinum resistance, especially in BRCA-mutated epithelial ovarian cancers.

Relationship between Egg Consumption and Body Composition as Well as Serum Cholesterol Level: Korea National Health and Nutrition Examination Survey 2008-2011.

We analyzed the relationship between egg consumption, body composition, and serum cholesterol levels. We obtained data on egg consumption by using a food frequency questionnaire (FFQ) (13,132 adults) and the 24-h dietary recall (24HR) (13,366 adults) from the fourth and fifth Korea National Health and Nutrition Examination Surveys (2008-2011). In men, consuming 2-3 eggs/week was associated with higher fat mass (FM), percentage body fat (PBF), and fat-to-muscle ratio (FtoM), compared to consuming <1 egg/week. In women, consuming 1-6 eggs/week was associated with higher low-density lipoprotein cholesterol, consuming 2-6 eggs/week was associated with higher total cholesterol, and consuming 4-6 eggs/week was associated with higher FM and high-density lipoprotein cholesterol, compared to consuming <1 egg/week. There was no relationship between egg consumption and the prevalence of dyslipidemia, and there was no relationship between egg consumption, body composition, and serum cholesterol levels according to the 24HR. However, there was some association with other cardiovascular diseases and consumption of certain amounts of eggs. Egg consumption investigated by FFQ was associated with body composition and serum cholesterol levels. However, the egg consumption investigated by the 24HR resulted in no health benefit or harm with respect to body composition and cholesterol.

Construction of Functional Protein Networks Using Domain Profile Associations.

Proteins are major functional molecules that physically and functionally interact to carry out cellular processes. The physical interactions are generally mediated by domain-level interactions. Thus, novel protein-protein interactions can be predicted using various computational methods based on domain-domain interactions, using resolved structures of protein complexes. Functional protein interactions can be inferred based on shared domains between proteins, since proteins involved in the same biological processes tend to harbor common domains. We recently developed a method of inferring functional interactions between proteins using associations between their domain compositions, which can be represented as domain profiles. Since the method requires only protein domain annotations, it can be easily applied to any species with a sequenced genome. Here, we describe in detail the method of generating domain profiles for proteins and measuring the association between them to infer functional interactions between proteins. We also demonstrate that domain profile association can be used to successfully construct a large-scale functional network of human proteins.

Pathway-specific protein domains are predictive for human diseases.

Protein domains are basic functional units of proteins. Many protein domains are pervasive among diverse biological processes, yet some are associated with specific pathways. Human complex diseases are generally viewed as pathway-level disorders. Therefore, we hypothesized that pathway-specific domains could be highly informative for human diseases. To test the hypothesis, we developed a network-based scoring scheme to quantify specificity of domain-pathway associations. We first generated domain profiles for human proteins, then constructed a co-pathway protein network based on the associations between domain profiles. Based on the score, we classified human protein domains into pathway-specific domains (PSDs) and non-specific domains (NSDs). We found that PSDs contained more pathogenic variants than NSDs. PSDs were also enriched for disease-associated mutations that disrupt protein-protein interactions (PPIs) and tend to have a moderate number of domain interactions. These results suggest that mutations in PSDs are likely to disrupt within-pathway PPIs, resulting in functional failure of pathways. Finally, we demonstrated the prediction capacity of PSDs for disease-associated genes with experimental validations in zebrafish. Taken together, the network-based quantitative method of modeling domain-pathway associations presented herein suggested underlying mechanisms of how protein domains associated with specific pathways influence mutational impacts on diseases via perturbations in within-pathway PPIs, and provided a novel genomic feature for interpreting genetic variants to facilitate the discovery of human disease genes.

Evaluation of an Automated Screening Assay, Compared to Indirect Immunofluorescence, an Extractable Nuclear Antigen Assay, and a Line Immunoassay in a Large Cohort of Asian Patients with Antinuclear Antibody-Associated Rheumatoid Diseases: A Multicenter Retrospective Study.

We assessed the diagnostic utility of the connective tissue disease (CTD) screen as an automated screening test, in comparison with the indirect immunofluorescence (IIF), EliA extractable nuclear antigen (ENA), and line immunoassay (LIA) for patients with antinuclear antibody- (ANA-) associated rheumatoid disease (AARD). A total of 1115 serum samples from two university hospitals were assayed using these four autoantibody-based methods. The AARD group consisted of patients with systemic lupus erythematosus (SLE), systemic sclerosis (SSc), Sjögren's syndrome (SS), and mixed connective tissue disease (MCTD). The qualitative results of all four autoantibody assays showed a significant association with AARDs, compared to controls (P < 0.0001 for all). The areas under the receiver operating characteristic curves (ROC-AUCs) of the CTD screen for differentiating total AARDs, SLE, SSc, SS, and MCTD from controls were 0.89, 0.93, 0.73, 0.93, and 0.95, respectively. The ROC-AUCs of combination testing with LIA were slightly higher in patients with AARDs (0.92) than those of CTD screen alone. Multivariate analysis indicated that all four autoantibody assays could independently predict AARDs. CTD screening alone and in combination with IIF, EliA ENA, and LIA are potentially valuable diagnostic approaches for predicting AARDs. Combining CTD screen with LIA might be effective for AARD patients.

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.

From sequencing data to gene functions: co-functional network approaches.

Advanced high-throughput sequencing technology accumulated massive amount of genomics and transcriptomics data in the public databases. Due to the high technical accessibility, DNA and RNA sequencing have huge potential for the study of gene functions in most species including animals and crops. A proven analytic platform to convert sequencing data to gene functional information is co-functional network. Because all genes exert their functions through interactions with others, network analysis is a legitimate way to study gene functions. The workflow of network-based functional study is composed of three steps: (i) inferencing co-functional links, (ii) evaluating and integrating the links into genome-scale networks, and (iii) generating functional hypotheses from the networks. Co-functional links can be inferred from DNA sequencing data by using phylogenetic profiling, gene neighborhood, domain profiling, associalogs, and co-expression analysis from RNA sequencing data. The inferred links are then evaluated and integrated into a genome-scale network with aid from gold-standard co-functional links. Functional hypotheses can be generated from the network based on (i) network connectivity, (ii) network propagation, and (iii) subnetwork analysis. The functional analysis pipeline described here requires only sequencing data which can be readily available for most species by next-generation sequencing technology. Therefore, co-functional networks will greatly potentiate the use of the sequencing data for the study of genetics in any cellular organism.

MUFFINN: cancer gene discovery via network analysis of somatic mutation data.

A major challenge for distinguishing cancer-causing driver mutations from inconsequential passenger mutations is the long-tail of infrequently mutated genes in cancer genomes. Here, we present and evaluate a method for prioritizing cancer genes accounting not only for mutations in individual genes but also in their neighbors in functional networks, MUFFINN (MUtations For Functional Impact on Network Neighbors). This pathway-centric method shows high sensitivity compared with gene-centric analyses of mutation data. Notably, only a marginal decrease in performance is observed when using 10 % of TCGA patient samples, suggesting the method may potentiate cancer genome projects with small patient populations.

Weighted mutual information analysis substantially improves domain-based functional network models.

MOTIVATION: Functional protein-protein interaction (PPI) networks elucidate molecular pathways underlying complex phenotypes, including those of human diseases. Extrapolation of domain-domain interactions (DDIs) from known PPIs is a major domain-based method for inferring functional PPI networks. However, the protein domain is a functional unit of the protein. Therefore, we should be able to effectively infer functional interactions between proteins based on the co-occurrence of domains. RESULTS: Here, we present a method for inferring accurate functional PPIs based on the similarity of domain composition between proteins by weighted mutual information (MI) that assigned different weights to the domains based on their genome-wide frequencies. Weighted MI outperforms other domain-based network inference methods and is highly predictive for pathways as well as phenotypes. A genome-scale human functional network determined by our method reveals numerous communities that are significantly associated with known pathways and diseases. Domain-based functional networks may, therefore, have potential applications in mapping domain-to-pathway or domain-to-phenotype associations. AVAILABILITY AND IMPLEMENTATION: Source code for calculating weighted mutual information based on the domain profile matrix is available from www.netbiolab.org/w/WMI CONTACT: Insuklee@yonsei.ac.kr SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

TRRUST: a reference database of human transcriptional regulatory interactions.

The reconstruction of transcriptional regulatory networks (TRNs) is a long-standing challenge in human genetics. Numerous computational methods have been developed to infer regulatory interactions between human transcriptional factors (TFs) and target genes from high-throughput data, and their performance evaluation requires gold-standard interactions. Here we present a database of literature-curated human TF-target interactions, TRRUST (transcriptional regulatory relationships unravelled by sentence-based text-mining, http://www.grnpedia.org/trrust), which currently contains 8,015 interactions between 748 TF genes and 1,975 non-TF genes. A sentence-based text-mining approach was employed for efficient manual curation of regulatory interactions from approximately 20 million Medline abstracts. To the best of our knowledge, TRRUST is the largest publicly available database of literature-curated human TF-target interactions to date. TRRUST also has several useful features: i) information about the mode-of-regulation; ii) tests for target modularity of a query TF; iii) tests for TF cooperativity of a query target; iv) inferences about cooperating TFs of a query TF; and v) prioritizing associated pathways and diseases with a query TF. We observed high enrichment of TF-target pairs in TRRUST for top-scored interactions inferred from high-throughput data, which suggests that TRRUST provides a reliable benchmark for the computational reconstruction of human TRNs.

Pathway-Dependent Effectiveness of Network Algorithms for Gene Prioritization.

A network-based approach has proven useful for the identification of novel genes associated with complex phenotypes, including human diseases. Because network-based gene prioritization algorithms are based on propagating information of known phenotype-associated genes through networks, the pathway structure of each phenotype might significantly affect the effectiveness of algorithms. We systematically compared two popular network algorithms with distinct mechanisms--direct neighborhood which propagates information to only direct network neighbors, and network diffusion which diffuses information throughout the entire network--in prioritization of genes for worm and human phenotypes. Previous studies reported that network diffusion generally outperforms direct neighborhood for human diseases. Although prioritization power is generally measured for all ranked genes, only the top candidates are significant for subsequent functional analysis. We found that high prioritizing power of a network algorithm for all genes cannot guarantee successful prioritization of top ranked candidates for a given phenotype. Indeed, the majority of the phenotypes that were more efficiently prioritized by network diffusion showed higher prioritizing power for top candidates by direct neighborhood. We also found that connectivity among pathway genes for each phenotype largely determines which network algorithm is more effective, suggesting that the network algorithm used for each phenotype should be chosen with consideration of pathway gene connectivity.

FlyNet: a versatile network prioritization server for the Drosophila community.

Drosophila melanogaster (fruit fly) has been a popular model organism in animal genetics due to the high accessibility of reverse-genetics tools. In addition, the close relationship between the Drosophila and human genomes rationalizes the use of Drosophila as an invertebrate model for human neurobiology and disease research. A platform technology for predicting candidate genes or functions would further enhance the usefulness of this long-established model organism for gene-to-phenotype mapping. Recently, the power of network prioritization for gene-to-phenotype mapping has been demonstrated in many organisms. Here we present a network prioritization server dedicated to Drosophila that covers ∼95% of the coding genome. This server, dubbed FlyNet, has several distinctive features, including (i) prioritization for both genes and functions; (ii) two complementary network algorithms: direct neighborhood and network diffusion; (iii) spatiotemporal-specific networks as an additional prioritization strategy for traits associated with a specific developmental stage or tissue and (iv) prioritization for human disease genes. FlyNet is expected to serve as a versatile hypothesis-generation platform for genes and functions in the study of basic animal genetics, developmental biology and human disease. FlyNet is available for free at http://www.inetbio.org/flynet.

RiceNet v2: an improved network prioritization server for rice genes.

Rice is the most important staple food crop and a model grass for studies of bioenergy crops. We previously published a genome-scale functional network server called RiceNet, constructed by integrating diverse genomics data and demonstrated the use of the network in genetic dissection of rice biotic stress responses and its usefulness for other grass species. Since the initial construction of the network, there has been a significant increase in the amount of publicly available rice genomics data. Here, we present an updated network prioritization server for Oryza sativa ssp. japonica, RiceNet v2 (http://www.inetbio.org/ricenet), which provides a network of 25 765 genes (70.1% of the coding genome) and 1 775 000 co-functional links. Ricenet v2 also provides two complementary methods for network prioritization based on: (i) network direct neighborhood and (ii) context-associated hubs. RiceNet v2 can use genes of the related subspecies O. sativa ssp. indica and the reference plant Arabidopsis for versatility in generating hypotheses. We demonstrate that RiceNet v2 effectively identifies candidate genes involved in rice root/shoot development and defense responses, demonstrating its usefulness for the grass research community.

Network-assisted genetic dissection of pathogenicity and drug resistance in the opportunistic human pathogenic fungus Cryptococcus neoformans.

Cryptococcus neoformans is an opportunistic human pathogenic fungus that causes meningoencephalitis. Due to the increasing global risk of cryptococcosis and the emergence of drug-resistant strains, the development of predictive genetics platforms for the rapid identification of novel genes governing pathogenicity and drug resistance of C. neoformans is imperative. The analysis of functional genomics data and genome-scale mutant libraries may facilitate the genetic dissection of such complex phenotypes but with limited efficiency. Here, we present a genome-scale co-functional network for C. neoformans, CryptoNet, which covers ~81% of the coding genome and provides an efficient intermediary between functional genomics data and reverse-genetics resources for the genetic dissection of C. neoformans phenotypes. CryptoNet is the first genome-scale co-functional network for any fungal pathogen. CryptoNet effectively identified novel genes for pathogenicity and drug resistance using guilt-by-association and context-associated hub algorithms. CryptoNet is also the first genome-scale co-functional network for fungi in the basidiomycota phylum, as Saccharomyces cerevisiae belongs to the ascomycota phylum. CryptoNet may therefore provide insights into pathway evolution between two distinct phyla of the fungal kingdom. The CryptoNet web server (www.inetbio.org/cryptonet) is a public resource that provides an interactive environment of network-assisted predictive genetics for C. neoformans.

EcoliNet: a database of cofunctional gene network for Escherichia coli.

During the past several decades, Escherichia coli has been a treasure chest for molecular biology. The molecular mechanisms of many fundamental cellular processes have been discovered through research on this bacterium. Although much basic research now focuses on more complex model organisms, E. coli still remains important in metabolic engineering and synthetic biology. Despite its long history as a subject of molecular investigation, more than one-third of the E. coli genome has no pathway annotation supported by either experimental evidence or manual curation. Recently, a network-assisted genetics approach to the efficient identification of novel gene functions has increased in popularity. To accelerate the speed of pathway annotation for the remaining uncharacterized part of the E. coli genome, we have constructed a database of cofunctional gene network with near-complete genome coverage of the organism, dubbed EcoliNet. We find that EcoliNet is highly predictive for diverse bacterial phenotypes, including antibiotic response, indicating that it will be useful in prioritizing novel candidate genes for a wide spectrum of bacterial phenotypes. We have implemented a web server where biologists can easily run network algorithms over EcoliNet to predict novel genes involved in a pathway or novel functions for a gene. All integrated cofunctional associations can be downloaded, enabling orthology-based reconstruction of gene networks for other bacterial species as well. Database URL: http://www.inetbio.org/ecolinet.

AraNet v2: an improved database of co-functional gene networks for the study of Arabidopsis thaliana and 27 other nonmodel plant species.

Arabidopsis thaliana is a reference plant that has been studied intensively for several decades. Recent advances in high-throughput experimental technology have enabled the generation of an unprecedented amount of data from A. thaliana, which has facilitated data-driven approaches to unravel the genetic organization of plant phenotypes. We previously published a description of a genome-scale functional gene network for A. thaliana, AraNet, which was constructed by integrating multiple co-functional gene networks inferred from diverse data types, and we demonstrated the predictive power of this network for complex phenotypes. More recently, we have observed significant growth in the availability of omics data for A. thaliana as well as improvements in data analysis methods that we anticipate will further enhance the integrated database of co-functional networks. Here, we present an updated co-functional gene network for A. thaliana, AraNet v2 (available at http://www.inetbio.org/aranet), which covers approximately 84% of the coding genome. We demonstrate significant improvements in both genome coverage and accuracy. To enhance the usability of the network, we implemented an AraNet v2 web server, which generates functional predictions for A. thaliana and 27 nonmodel plant species using an orthology-based projection of nonmodel plant genes on the A. thaliana gene network.

YeastNet v3: a public database of data-specific and integrated functional gene networks for Saccharomyces cerevisiae.

Saccharomyces cerevisiae, i.e. baker's yeast, is a widely studied model organism in eukaryote genetics because of its simple protocols for genetic manipulation and phenotype profiling. The high abundance of publicly available data that has been generated through diverse 'omics' approaches has led to the use of yeast for many systems biology studies, including large-scale gene network modeling to better understand the molecular basis of the cellular phenotype. We have previously developed a genome-scale gene network for yeast, YeastNet v2, which has been used for various genetics and systems biology studies. Here, we present an updated version, YeastNet v3 (available at http://www.inetbio.org/yeastnet/), that significantly improves the prediction of gene-phenotype associations. The extended genome in YeastNet v3 covers up to 5818 genes (∼99% of the coding genome) wired by 362 512 functional links. YeastNet v3 provides a new web interface to run the tools for network-guided hypothesis generations. YeastNet v3 also provides edge information for all data-specific networks (∼2 million functional links) as well as the integrated networks. Therefore, users can construct alternative versions of the integrated network by applying their own data integration algorithm to the same data-specific links.