Deletion of PD-1 destabilizes the lineage identity and metabolic fitness of tumor-infiltrating regulatory T cells.

Regulatory T (Treg) cells have an immunosuppressive function and highly express the immune checkpoint receptor PD-1 in the tumor microenvironment; however, the function of PD-1 in tumor-infiltrating (TI) Treg cells remains controversial. Here, we showed that conditional deletion of PD-1 in Treg cells delayed tumor progression. In Pdcd1fl/flFoxp3eGFP-Cre-ERT2(+/-) mice, in which both PD-1-expressing and PD-1-deficient Treg cells coexisted in the same tissue environment, conditional deletion of PD-1 in Treg cells resulted in impairment of the proliferative and suppressive capacity of TI Treg cells. PD-1 antibody therapy reduced the TI Treg cell numbers, but did not directly restore the cytokine production of TI CD8+ T cells in TC-1 lung cancer. Single-cell analysis indicated that PD-1 signaling promoted lipid metabolism, proliferation and suppressive pathways in TI Treg cells. These results suggest that PD-1 ablation or inhibition can enhance antitumor immunity by weakening Treg cell lineage stability and metabolic fitness in the tumor microenvironment.

A metabolic switch orchestrated by IL-18 and the cyclic dinucleotide cGAMP programs intestinal tolerance.

Tissues are exposed to diverse inflammatory challenges that shape future inflammatory responses. While cellular metabolism regulates immune function, how metabolism programs and stabilizes immune states within tissues and tunes susceptibility to inflammation is poorly understood. Here, we describe an innate immune metabolic switch that programs long-term intestinal tolerance. Intestinal interleukin-18 (IL-18) stimulation elicited tolerogenic macrophages by preventing their proinflammatory glycolytic polarization via metabolic reprogramming to fatty acid oxidation (FAO). FAO reprogramming was triggered by IL-18 activation of SLC12A3 (NCC), leading to sodium influx, release of mitochondrial DNA, and activation of stimulator of interferon genes (STING). FAO was maintained in macrophages by a bistable switch that encoded memory of IL-18 stimulation and by intercellular positive feedback that sustained the production of macrophage-derived 2'3'-cyclic GMP-AMP (cGAMP) and epithelial-derived IL-18. Thus, a tissue-reinforced metabolic switch encodes durable immune tolerance in the gut and may enable reconstructing compromised immune tolerance in chronic inflammation.

scHumanNet: a single-cell network analysis platform for the study of cell-type specificity of disease genes.

A major challenge in single-cell biology is identifying cell-type-specific gene functions, which may substantially improve precision medicine. Differential expression analysis of genes is a popular, yet insufficient approach, and complementary methods that associate function with cell type are required. Here, we describe scHumanNet (https://github.com/netbiolab/scHumanNet), a single-cell network analysis platform for resolving cellular heterogeneity across gene functions in humans. Based on cell-type-specific gene networks (CGNs) constructed under the guidance of the HumanNet reference interactome, scHumanNet displayed higher functional relevance to the cellular context than CGNs built by other methods on single-cell transcriptome data. Cellular deconvolution of gene signatures based on network compactness across cell types revealed breast cancer prognostic markers associated with T cells. scHumanNet could also prioritize genes associated with particular cell types using CGN centrality and identified the differential hubness of CGNs between disease and healthy conditions. We demonstrated the usefulness of scHumanNet by uncovering T-cell-specific functional effects of GITR, a prognostic gene for breast cancer, and functional defects in autism spectrum disorder genes specific for inhibitory neurons. These results suggest that scHumanNet will advance our understanding of cell-type specificity across human disease genes.

Transcriptomic and metabolomic analysis unveils nanoplastic-induced gut barrier dysfunction via STAT1/6 and ERK pathways.

The widespread prevalence of micro and nanoplastics in the environment raises concerns about their potential impact on human health. Recent evidence demonstrates the presence of nanoplastics in human blood and tissues following ingestion and inhalation, yet the specific risks and mechanisms of nanoplastic toxicity remain inadequately understood. In this study, we aimed to explore the molecular mechanisms underlying the toxicity of nanoplastics at both systemic and molecular levels by analyzing the transcriptomic/metabolomic responses and signaling pathways in the intestines of mice after oral administration of nanoplastics. Transcriptome analysis in nanoplastic-administered mice revealed a notable upregulation of genes involved in pro-inflammatory immune responses. In addition, nanoplastics substantially reduced the expression of tight junction proteins, including occludin, zonula occluden-1, and tricellulin, which are crucial for maintaining gut barrier integrity and function. Importantly, nanoplastic administration increased gut permeability and exacerbated dextran sulfate sodium-induced colitis. Further investigation into the underlying molecular mechanisms highlighted significant activation of signaling transsducer and activator of transcription (STAT)1 and STAT6 by nanoplastic administration, which was in line with the elevation of interferon and JAK-STAT pathway signatures identified through transcriptome enrichment analysis. Additionally, the consumption of nanoplastics specifically induced nuclear factor kappa-B (NF-κB) and extracellular signal-regulated kinase (ERK)1/2 signaling pathways in the intestines. Collectively, this study identifies molecular mechanisms contributing to adverse effects mediated by nanoplastics in the intestine, providing novel insights into the pathophysiological consequences of nanoplastic exposure.

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.

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.

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.

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.

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.

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.

SoyNet: a database of co-functional networks for soybean Glycine max.

Soybean (Glycine max) is a legume crop with substantial economic value, providing a source of oil and protein for humans and livestock. More than 50% of edible oils consumed globally are derived from this crop. Soybean plants are also important for soil fertility, as they fix atmospheric nitrogen by symbiosis with microorganisms. The latest soybean genome annotation (version 2.0) lists 56 044 coding genes, yet their functional contributions to crop traits remain mostly unknown. Co-functional networks have proven useful for identifying genes that are involved in a particular pathway or phenotype with various network algorithms. Here, we present SoyNet (available at www.inetbio.org/soynet), a database of co-functional networks for G. max and a companion web server for network-based functional predictions. SoyNet maps 1 940 284 co-functional links between 40 812 soybean genes (72.8% of the coding genome), which were inferred from 21 distinct types of genomics data including 734 microarrays and 290 RNA-seq samples from soybean. SoyNet provides a new route to functional investigation of the soybean genome, elucidating genes and pathways of agricultural importance.

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.

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.

Network perturbation by recurrent regulatory variants in cancer.

Cancer driving genes have been identified as recurrently affected by variants that alter protein-coding sequences. However, a majority of cancer variants arise in noncoding regions, and some of them are thought to play a critical role through transcriptional perturbation. Here we identified putative transcriptional driver genes based on combinatorial variant recurrence in cis-regulatory regions. The identified genes showed high connectivity in the cancer type-specific transcription regulatory network, with high outdegree and many downstream genes, highlighting their causative role during tumorigenesis. In the protein interactome, the identified transcriptional drivers were not as highly connected as coding driver genes but appeared to form a network module centered on the coding drivers. The coding and regulatory variants associated via these interactions between the coding and transcriptional drivers showed exclusive and complementary occurrence patterns across tumor samples. Transcriptional cancer drivers may act through an extensive perturbation of the regulatory network and by altering protein network modules through interactions with coding driver genes.

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.

MouseNet v2: a database of gene networks for studying the laboratory mouse and eight other model vertebrates.

Laboratory mouse, Mus musculus, is one of the most important animal tools in biomedical research. Functional characterization of the mouse genes, hence, has been a long-standing goal in mammalian and human genetics. Although large-scale knockout phenotyping is under progress by international collaborative efforts, a large portion of mouse genome is still poorly characterized for cellular functions and associations with disease phenotypes. A genome-scale functional network of mouse genes, MouseNet, was previously developed in context of MouseFunc competition, which allowed only limited input data for network inferences. Here, we present an improved mouse co-functional network, MouseNet v2 (available at http://www.inetbio.org/mousenet), which covers 17 714 genes (>88% of coding genome) with 788 080 links, along with a companion web server for network-assisted functional hypothesis generation. The network database has been substantially improved by large expansion of genomics data. For example, MouseNet v2 database contains 183 co-expression networks inferred from 8154 public microarray samples. We demonstrated that MouseNet v2 is predictive for mammalian phenotypes as well as human diseases, which suggests its usefulness in discovery of novel disease genes and dissection of disease pathways. Furthermore, MouseNet v2 database provides functional networks for eight other vertebrate models used in various research fields.

An integrated systems biology approach identifies positive cofactor 4 as a factor that increases reprogramming efficiency.

Spermatogonial stem cells (SSCs) can spontaneously dedifferentiate into embryonic stem cell (ESC)-like cells, which are designated as multipotent SSCs (mSSCs), without ectopic expression of reprogramming factors. Interestingly, SSCs express key pluripotency genes such as Oct4, Sox2, Klf4 and Myc. Therefore, molecular dissection of mSSC reprogramming may provide clues about novel endogenous reprogramming or pluripotency regulatory factors. Our comparative transcriptome analysis of mSSCs and induced pluripotent stem cells (iPSCs) suggests that they have similar pluripotency states but are reprogrammed via different transcriptional pathways. We identified 53 genes as putative pluripotency regulatory factors using an integrated systems biology approach. We demonstrated a selected candidate, Positive cofactor 4 (Pc4), can enhance the efficiency of somatic cell reprogramming by promoting and maintaining transcriptional activity of the key reprograming factors. These results suggest that Pc4 has an important role in inducing spontaneous somatic cell reprogramming via up-regulation of key pluripotency genes.

Generation of Integration-free Induced Neural Stem Cells from Mouse Fibroblasts.

The viral vector-mediated overexpression of the defined transcription factors, Brn4/Pou3f4, Sox2, Klf4, and c-Myc (BSKM), could induce the direct conversion of somatic fibroblasts into induced neural stem cells (iNSCs). However, viral vectors may be randomly integrated into the host genome thereby increasing the risk for undesired genotoxicity, mutagenesis, and tumor formation. Here we describe the generation of integration-free iNSCs from mouse fibroblasts by non-viral episomal vectors containing BSKM. The episomal vector-derived iNSCs (e-iNSCs) closely resemble control NSCs, and iNSCs generated by retrovirus (r-iNSCs) in morphology, gene expression profile, epigenetic status, and self-renewal capacity. The e-iNSCs are functionally mature, as they could differentiate into all the neuronal cell types both in vitro and in vivo Our study provides a novel concept for generating functional iNSCs using a non-viral, non-integrating, plasmid-based system that could facilitate their biomedical applicability.

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.