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

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.

Integration of probabilistic functional networks without an external Gold Standard

BACKGROUND

Probabilistic functional integrated networks (PFINs) are designed to aid our understanding of cellular biology and can be used to generate testable hypotheses about protein function. PFINs are generally created by scoring the quality of interaction datasets against a Gold Standard dataset, usually chosen from a separate high-quality data source, prior to their integration. Use of an external Gold Standard has several drawbacks, including data redundancy, data loss and the need for identifier mapping, which can complicate the network build and impact on PFIN performance. Additionally, there typically are no Gold Standard data for non-model organisms.

RESULTS

We describe the development of an integration technique, ssNet, that scores and integrates both high-throughput and low-throughout data from a single source database in a consistent manner without the need for an external Gold Standard dataset. Using data from Saccharomyces cerevisiae we show that ssNet is easier and faster, overcoming the challenges of data redundancy, Gold Standard bias and ID mapping. In addition ssNet results in less loss of data and produces a more complete network.

CONCLUSIONS

The ssNet method allows PFINs to be built successfully from a single database, while producing comparable network performance to networks scored using an external Gold Standard source and with reduced data loss.

Predicted mouse interactome and network-based interpretation of differentially expressed genes

The house mouse or Mus musculus has become a premier mammalian model for genetic research due to its genetic and physiological similarities to humans. It brought mechanistic insights into numerous human diseases and has been routinely used to assess drug efficiency and toxicity, as well as to predict patient responses. To facilitate molecular mechanism studies in mouse, we present the Mouse Interactome Database (MID, Version 1), which includes 155,887 putative functional associations between mouse protein-coding genes inferred from functional association evidence integrated from 9 public databases. These putative functional associations are expected to cover 19.32% of all mouse protein interactions, and 26.02% of these function associations may represent protein interactions. On top of MID, we developed a gene set linkage analysis (GSLA) web tool to annotate potential functional impacts from observed differentially expressed genes. Two case studies show that the MID/GSLA system provided precise and informative annotations that other widely used gene set annotation tools, such as PANTHER and DAVID, did not. Both MID and GSLA are accessible through the website http://mouse.biomedtzc.cn.

A novel probabilistic generator for large-scale gene association networks

MOTIVATION

Gene expression data provide an opportunity for reverse-engineering gene-gene associations using network inference methods. However, it is difficult to assess the performance of these methods because the true underlying network is unknown in real data. Current benchmarks address this problem by subsampling a known regulatory network to conduct simulations. But the topology of regulatory networks can vary greatly across organisms or tissues, and reference-based generators-such as GeneNetWeaver-are not designed to capture this heterogeneity. This means, for example, benchmark results from the E. coli regulatory network will not carry over to other organisms or tissues. In contrast, probabilistic generators do not require a reference network, and they have the potential to capture a rich distribution of topologies. This makes probabilistic generators an ideal approach for obtaining a robust benchmarking of network inference methods.

RESULTS

We propose a novel probabilistic network generator that (1) provides an alternative to address the inherent limitation of reference-based generators and (2) is able to create realistic gene association networks, and (3) captures the heterogeneity found across gold-standard networks better than existing generators used in practice. Eight organism-specific and 12 human tissue-specific gold-standard association networks are considered. Several measures of global topology are used to determine the similarity of generated networks to the gold-standards. Along with demonstrating the variability of network structure across organisms and tissues, we show that the commonly used "scale-free" model is insufficient for replicating these structures.

AVAILABILITY

This generator is implemented in the R package "SeqNet" and is available on CRAN (https://cran.r-project.org/web/packages/SeqNet/index.html).

Metabolic remodelling during early mouse embryo development

During early mammalian embryogenesis, changes in cell growth and proliferation depend on strict genetic and metabolic instructions. However, our understanding of metabolic reprogramming and its influence on epigenetic regulation in early embryo development remains elusive. Here we show a comprehensive metabolomics profiling of key stages in mouse early development and the two-cell and blastocyst embryos, and we reconstructed the metabolic landscape through the transition from totipotency to pluripotency. Our integrated metabolomics and transcriptomics analysis shows that while two-cell embryos favour methionine, polyamine and glutathione metabolism and stay in a more reductive state, blastocyst embryos have higher metabolites related to the mitochondrial tricarboxylic acid cycle, and present a more oxidative state. Moreover, we identify a reciprocal relationship between α-ketoglutarate (α-KG) and the competitive inhibitor of α-KG-dependent dioxygenases, L-2-hydroxyglutarate (L-2-HG), where two-cell embryos inherited from oocytes and one-cell zygotes display higher L-2-HG, whereas blastocysts show higher α-KG. Lastly, increasing 2-HG availability impedes erasure of global histone methylation markers after fertilization. Together, our data demonstrate dynamic and interconnected metabolic, transcriptional and epigenetic network remodelling during early mouse embryo development.

Biomolecule and Bioentity Interaction Databases in Systems Biology: A Comprehensive Review

Technological advances in high-throughput techniques have resulted in tremendous growth of complex biological datasets providing evidence regarding various biomolecular interactions. To cope with this data flood, computational approaches, web services, and databases have been implemented to deal with issues such as data integration, visualization, exploration, organization, scalability, and complexity. Nevertheless, as the number of such sets increases, it is becoming more and more difficult for an end user to know what the scope and focus of each repository is and how redundant the information between them is. Several repositories have a more general scope, while others focus on specialized aspects, such as specific organisms or biological systems. Unfortunately, many of these databases are self-contained or poorly documented and maintained. For a clearer view, in this article we provide a comprehensive categorization, comparison and evaluation of such repositories for different bioentity interaction types. We discuss most of the publicly available services based on their content, sources of information, data representation methods, user-friendliness, scope and interconnectivity, and we comment on their strengths and weaknesses. We aim for this review to reach a broad readership varying from biomedical beginners to experts and serve as a reference article in the field of Network Biology.

Systems biology analysis identifies TNFRSF9 as a functional marker of tumor-infiltrating regulatory T-cell enabling clinical outcome prediction in lung cancer

Regulatory T cells (Tregs) are enriched in the tumor microenvironment and play key roles in immune evasion of cancer cells. Cell surface markers specific for tumor-infiltrating Tregs (TI-Tregs) can be effectively targeted to enhance antitumor immunity and used for stratification of immunotherapy outcomes. Here, we present a systems biology approach to identify functional cell surface markers for TI-Tregs. We selected differentially expressed genes for surface proteins of TI-Tregs and compared these with other CD4⁺ T cells using bulk RNA-sequencing data from murine lung cancer models. Thereafter, we filtered for human orthologues with conserved expression in TI-Tregs using single-cell transcriptome data from patients with non-small cell lung cancer (NSCLC). To evaluate the functional importance of expression-based markers of TI-Tregs, we utilized network-based measure of context-associated centrality in a Treg-specific coregulatory network. We identified TNFRSF9 (also known as 4-1BB or CD137), a previously reported target for enhancing antitumor immunity, among the final candidates for TI-Treg markers with high functional importance score. We found that the low TNFRSF9 expression level in Tregs was associated with enhanced overall survival rate and response to anti-PD-1 immunotherapy in patients with NSCLC, proposing that TNFRSF9 promotes immune suppressive activity of Tregs in tumor. Collectively, these results demonstrated that integrative transcriptome and network analysis can facilitate the discovery of functional markers of tumor-specific immune cells to develop novel therapeutic targets and biomarkers for boosting cancer immunotherapy.

CoCoCoNet: conserved and comparative co-expression across a diverse set of species

Co-expression analysis has provided insight into gene function in organisms from Arabidopsis to zebrafish. Comparison across species has the potential to enrich these results, for example by prioritizing among candidate human disease genes based on their network properties or by finding alternative model systems where their co-expression is conserved. Here, we present CoCoCoNet as a tool for identifying conserved gene modules and comparing co-expression networks. CoCoCoNet is a resource for both data and methods, providing gold standard networks and sophisticated tools for on-the-fly comparative analyses across 14 species. We show how CoCoCoNet can be used in two use cases. In the first, we demonstrate deep conservation of a nucleolus gene module across very divergent organisms, and in the second, we show how the heterogeneity of autism mechanisms in humans can be broken down by functional groups and translated to model organisms. CoCoCoNet is free to use and available to all at https://milton.cshl.edu/CoCoCoNet, with data and R scripts available at ftp://milton.cshl.edu/data.

Computational resources associating diseases with genotypes, phenotypes and exposures

The causes of a disease and its therapies are not only related to genotypes, but also associated with other factors, including phenotypes, environmental exposures, drugs and chemical molecules. Distinguishing disease-related factors from many neutral factors is critical as well as difficult. Over the past two decades, bioinformaticians have developed many computational resources to integrate the omics data and discover associations among these factors. However, researchers and clinicians are experiencing difficulties in choosing appropriate resources from hundreds of relevant databases and software tools. Here, in order to assist the researchers and clinicians, we systematically review the public computational resources of human diseases related to genotypes, phenotypes, environment factors, drugs and chemical exposures. We briefly describe the development history of these computational resources, followed by the details of the relevant databases and software tools. We finally conclude with a discussion of current challenges and future opportunities as well as prospects on this topic.

Genetic cooperativity in multi-layer networks implicates cell survival and senescence in the striatum of Huntington's disease mice synchronous to symptoms

Motivation

Huntington’s disease (HD) may evolve through gene deregulation. However, the impact of gene deregulation on the dynamics of genetic cooperativity in HD remains poorly understood. Here, we built a multi-layer network model of temporal dynamics of genetic cooperativity in the brain of HD knock-in mice (allelic series of Hdh mice). To enhance biological precision and gene prioritization, we integrated three complementary families of source networks, all inferred from the same RNA-seq time series data in Hdh mice, into weighted-edge networks where an edge recapitulates path-length variation across source-networks and age-points.

Results

Weighted edge networks identify two consecutive waves of tight genetic cooperativity enriched in deregulated genes (critical phases), pre-symptomatically in the cortex, implicating neurotransmission, and symptomatically in the striatum, implicating cell survival (e.g. Hipk4) intertwined with cell proliferation (e.g. Scn4b) and cellular senescence (e.g. Cdkn2a products) responses. Top striatal weighted edges are enriched in modulators of defective behavior in invertebrate models of HD pathogenesis, validating their relevance to neuronal dysfunction in vivo. Collectively, these findings reveal highly dynamic temporal features of genetic cooperativity in the brain of Hdh mice where a 2-step logic highlights the importance of cellular maintenance and senescence in the striatum of symptomatic mice, providing highly prioritized targets.

Availability and implementation

Weighted edge network analysis (WENA) data and source codes for performing spectral decomposition of the signal (SDS) and WENA analysis, both written using Python, are available at http://www.broca.inserm.fr/HD-WENA/.

Supplementary information

Supplementary data are available at Bioinformatics online.

CORUM: the comprehensive resource of mammalian protein complexes-2019

CORUM is a database that provides a manually curated repository of experimentally characterized protein complexes from mammalian organisms, mainly human (67%), mouse (15%) and rat (10%). Given the vital functions of these macromolecular machines, their identification and functional characterization is foundational to our understanding of normal and disease biology. The new CORUM 3.0 release encompasses 4274 protein complexes offering the largest and most comprehensive publicly available dataset of mammalian protein complexes. The CORUM dataset is built from 4473 different genes, representing 22% of the protein coding genes in humans. Protein complexes are described by a protein complex name, subunit composition, cellular functions as well as the literature references. Information about stoichiometry of subunits depends on availability of experimental data. Recent developments include a graphical tool displaying known interactions between subunits. This allows the prediction of structural interconnections within protein complexes of unknown structure. In addition, we present a set of 58 protein complexes with alternatively spliced subunits. Those were found to affect cellular functions such as regulation of apoptotic activity, protein complex assembly or define cellular localization. CORUM is freely accessible at http://mips.helmholtz-muenchen.de/corum/.

BraInMap Elucidates the Macromolecular Connectivity Landscape of Mammalian Brain

Connectivity webs mediate the unique biology of the mammalian brain. Yet, while cell circuit maps are increasingly available, knowledge of their underlying molecular networks remains limited. Here, we applied multi-dimensional biochemical fractionation with mass spectrometry and machine learning to survey endogenous macromolecules across the adult mouse brain. We defined a global "interactome" comprising over one thousand multi-protein complexes. These include hundreds of brain-selective assemblies that have distinct physical and functional attributes, show regional and cell-type specificity, and have links to core neurological processes and disorders. Using reciprocal pull-downs and a transgenic model, we validated a putative 28-member RNA-binding protein complex associated with amyotrophic lateral sclerosis, suggesting a coordinated function in alternative splicing in disease progression. This brain interaction map (BraInMap) resource facilitates mechanistic exploration of the unique molecular machinery driving core cellular processes of the central nervous system. It is publicly available and can be explored here https://www.bu.edu/dbin/cnsb/mousebrain/.

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.

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.

netGO: R-Shiny package for network-integrated pathway enrichment analysis

Summary

We present an R-Shiny package, netGO, for novel network-integrated pathway enrichment analysis. The conventional Fisher’s exact test (FET) considers the extent of overlap between target genes and pathway gene-sets, while recent network-based analysis tools consider only network interactions between the two. netGO implements an intuitive framework to integrate both the overlap and networks into a single score, and adaptively resamples genes based on network degrees to assess the pathway enrichment. In benchmark tests for gene expression and genome-wide association study (GWAS) data, netGO captured the relevant gene-sets better than existing tools, especially when analyzing a small number of genes. Specifically, netGO provides user-interactive visualization of the target genes, enriched gene-set and their network interactions for both netGO and FET results for further analysis. For this visualization, we also developed a standalone R-Shiny package shinyCyJS to connect R-shiny and the JavaScript version of cytoscape.

Availability and implementation

netGO R-Shiny package is freely available from github, https://github.com/unistbig/netGO.

Supplementary information

Supplementary data are available at Bioinformatics online.

Epithelial RABGEF1 deficiency promotes intestinal inflammation by dysregulating intrinsic MYD88-dependent innate signaling

Intestinal epithelial cells (IECs) contribute to the regulation of intestinal homeostasis and inflammation through their interactions with the environment and host immune responses. Yet our understanding of IEC-intrinsic regulatory pathways remains incomplete. Here, we identify the guanine nucleotide exchange factor RABGEF1 as a regulator of intestinal homeostasis and innate pathways dependent on IECs. Mice with IEC-specific Rabgef1 deletion (called Rabgef1^IEC-KO mice) developed a delayed spontaneous colitis associated with the local upregulation of IEC chemokine expression. In mouse models of colitis based on Interleukin-10 deficiency or dextran sodium sulfate (DSS) exposure, we found that IEC-intrinsic RABGEF1 deficiency exacerbated development of intestinal pathology and dysregulated IEC innate pathways and chemokine expression. Mechanistically, we showed that RABGEF1 deficiency in mouse IECs in vitro was associated with an impairment of early endocytic events, an increased activation of the p38 mitogen-activated protein kinase (MAPK)-dependent pathway, and increased chemokine secretion. Moreover, we provided evidence that the development of spontaneous colitis was dependent on microbiota-derived signals and intrinsic MYD88-dependent pathways in vivo. Our study identifies mouse RABGEF1 as an important regulator of intestinal inflammation, MYD88-dependent IEC-intrinsic signaling, and chemokine production. This suggests that RABGEF1-dependent pathways represent interesting therapeutic targets for inflammatory conditions in the gut.

NGSEA: Network-Based Gene Set Enrichment Analysis for Interpreting Gene Expression Phenotypes with Functional Gene Sets

Gene set enrichment analysis (GSEA) is a popular tool to identify underlying biological processes in clinical samples using their gene expression phenotypes. GSEA measures the enrichment of annotated gene sets that represent biological processes for differentially expressed genes (DEGs) in clinical samples. GSEA may be suboptimal for functional gene sets; however, because DEGs from the expression dataset may not be functional genes per se but dysregulated genes perturbed by bona fide functional genes. To overcome this shortcoming, we developed network-based GSEA (NGSEA), which measures the enrichment score of functional gene sets using the expression difference of not only individual genes but also their neighbors in the functional network. We found that NGSEA outperformed GSEA in identifying pathway gene sets for matched gene expression phenotypes. We also observed that NGSEA substantially improved the ability to retrieve known anti-cancer drugs from patient-derived gene expression data using drug-target gene sets compared with another method, Connectivity Map. We also repurposed FDA-approved drugs using NGSEA and experimentally validated budesonide as a chemical with anti-cancer effects for colorectal cancer. We, therefore, expect that NGSEA will facilitate both pathway interpretation of gene expression phenotypes and anti-cancer drug repositioning. NGSEA is freely available at www.inetbio.org/ngsea.

UNIPred-Web: a web tool for the integration and visualization of biomolecular networks for protein function prediction

Background

One of the main issues in the automated protein function prediction (AFP) problem is the integration of multiple networked data sources. The UNIPred algorithm was thereby proposed to efficiently integrate —in a function-specific fashion— the protein networks by taking into account the imbalance that characterizes protein annotations, and to subsequently predict novel hypotheses about unannotated proteins. UNIPred is publicly available as R code, which might result of limited usage for non-expert users. Moreover, its application requires efforts in the acquisition and preparation of the networks to be integrated. Finally, the UNIPred source code does not handle the visualization of the resulting consensus network, whereas suitable views of the network topology are necessary to explore and interpret existing protein relationships.

Results

We address the aforementioned issues by proposing UNIPred-Web, a user-friendly Web tool for the application of the UNIPred algorithm to a variety of biomolecular networks, already supplied by the system, and for the visualization and exploration of protein networks. We support different organisms and different types of networks —e.g., co-expression, shared domains and physical interaction networks. Users are supported in the different phases of the process, ranging from the selection of the networks and the protein function to be predicted, to the navigation of the integrated network. The system also supports the upload of user-defined protein networks. The vertex-centric and the highly interactive approach of UNIPred-Web allow a narrow exploration of specific proteins, and an interactive analysis of large sub-networks with only a few mouse clicks.

Conclusions

UNIPred-Web offers a practical and intuitive (visual) guidance to biologists interested in gaining insights into protein biomolecular functions. UNIPred-Web provides facilities for the integration of networks, and supplies a framework for the imbalance-aware protein network integration of nine organisms, the prediction of thousands of GO protein functions, and a easy-to-use graphical interface for the visual analysis, navigation and interpretation of the integrated networks and of the functional predictions.

Large-scale neuroanatomical study uncovers 198 gene associations in mouse brain morphogenesis

Brain morphogenesis is an important process contributing to higher-order cognition, however our knowledge about its biological basis is largely incomplete. Here we analyze 118 neuroanatomical parameters in 1,566 mutant mouse lines and identify 198 genes whose disruptions yield NeuroAnatomical Phenotypes (NAPs), mostly affecting structures implicated in brain connectivity. Groups of functionally similar NAP genes participate in pathways involving the cytoskeleton, the cell cycle and the synapse, display distinct fetal and postnatal brain expression dynamics and importantly, their disruption can yield convergent phenotypic patterns. 17% of human unique orthologues of mouse NAP genes are known loci for cognitive dysfunction. The remaining 83% constitute a vast pool of genes newly implicated in brain architecture, providing the largest study of mouse NAP genes and pathways. This offers a complementary resource to human genetic studies and predict that many more genes could be involved in mammalian brain morphogenesis.

Brain morphogenesis is an important process contributing to higher-order cognition, however our knowledge about its biological basis is largely incomplete. Here, authors analyzed 118 neuroanatomical parameters in 1,566 mutant mouse lines to identify 198 genes whose disruptions yield neuroanatomical phenotypes

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.

CNS cell type-specific gene profiling of P301S tau transgenic mice identifies genes dysregulated by progressive tau accumulation

The microtubule-associated protein tau undergoes aberrant modification resulting in insoluble brain deposits in various neurodegenerative diseases, including frontotemporal dementia (FTD), progressive supranuclear palsy, and corticobasal degeneration. Tau aggregates can form in different cell types of the central nervous system (CNS) but are most prevalent in neurons. We have previously recapitulated aspects of human FTD in mouse models by overexpressing mutant human tau in CNS neurons, including a P301S tau variant in TAU58/2 mice, characterized by early-onset and progressive behavioral deficits and FTD-like neuropathology. The molecular mechanisms underlying the functional deficits of TAU58/2 mice remain mostly elusive. Here, we employed functional genomics (i.e. RNAseq) to determine differentially expressed genes in young and aged TAU58/2 mice to identify alterations in cellular processes that may contribute to neuropathy. We identified genes in cortical brain samples differentially regulated between young and old TAU58/2 mice relative to nontransgenic littermates and by comparative analysis with a dataset of CNS cell type-specific genes expressed in nontransgenic mice. Most differentially-regulated genes had known or putative roles in neurons and included presynaptic and excitatory genes. Specifically, we observed changes in presynaptic factors, glutamatergic signaling, and protein scaffolding. Moreover, in the aged mice, expression levels of several genes whose expression was annotated to occur in other brain cell types were altered. Immunoblotting and immunostaining of brain samples from the TAU58/2 mice confirmed altered expression and localization of identified and network-linked proteins. Our results have revealed genes dysregulated by progressive tau accumulation in an FTD mouse model.

Integrative analysis of multiple genomic variables using a hierarchical Bayesian model

Motivation

Genes showing congruent differences in several genomic variables between two biological conditions are crucial to unravel causalities behind phenotypes of interest. Detecting such genes is important in biomedical research, e.g. when identifying genes responsible for cancer development. Small sample sizes common in next-generation sequencing studies are a key challenge, and there are still only very few statistical methods to analyze more than two genomic variables in an integrative, model-based way. Here, we present a novel bioinformatics approach to detect congruent differences between two biological conditions in a larger number of different measurements such as various epigenetic marks or mRNA transcript levels.

Results

We propose a coefficient quantifying the degree to which genes present consistent alterations in multiple (more than two) genomic variables when comparing samples presenting a condition of interest (e.g. cancer) to a reference group. A hierarchical Bayesian model is employed to assess uncertainty on a gene level, incorporating information on functional relationships between genes. We demonstrate the approach on different data sets containing RNA-seq gene transcripton and up to four ChIP-seq histone modification measurements. Both the coefficient-based ranking and the inference based on the model lead to a plausible prioritizing of candidate genes when analyzing multiple genomic variables.

Availability and implementation

BUGS code in the Supplement.

Contact

m.schaefer@uni-duesseldorf.de.

Supplementary information

Supplementary data are available at Bioinformatics online.

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.

Naming CRISPR alleles: endonuclease-mediated mutation nomenclature across species

The widespread use of CRISPR/Cas and other targeted endonuclease technologies in many species has led to an explosion in the generation of new mutations and alleles. The ability to generate many different mutations from the same target sequence either by homology-directed repair with a donor sequence or non-homologous end joining-induced insertions and deletions necessitates a means for representing these mutations in literature and databases. Standardized nomenclature can be used to generate unambiguous, concise, and specific symbols to represent mutations and alleles. The research communities of a variety of species using CRISPR/Cas and other endonuclease-mediated mutation technologies have developed different approaches to naming and identifying such alleles and mutations. While some organism-specific research communities have developed allele nomenclature that incorporates the method of generation within the official allele or mutant symbol, others use metadata tags that include method of generation or mutagen. Organism-specific research community databases together with organism-specific nomenclature committees are leading the way in providing standardized nomenclature and metadata to facilitate the integration of data from alleles and mutations generated using CRISPR/Cas and other targeted endonucleases.

Network-Based Gene Function Prediction in Mouse and Other Model Vertebrates Using MouseNet Server

The mouse, Mus musculus, is a popular model organism for the study of human genes involved in development, immunology, and disease phenotypes. Despite recent revolutions in gene-knockout technologies in mouse, identification of candidate genes for functions of interest can further accelerate the discovery of novel gene functions. The collaborative nature of genetic functions allows for the inference of gene functions based on the principle of guilt-by-association. Genome-scale co-functional networks could therefore provide functional predictions for genes via network analysis. We recently constructed such a network for mouse (MouseNet), which interconnects over 88% of protein-coding genes with 788,080 functional relationships. The companion web server ( www.inetbio.org/mousenet ) enables researchers with no bioinformatics expertise to generate predictions that facilitate discovery of novel gene functions. In this chapter, we present the theoretical framework for MouseNet, as well as step-by-step instructions and technical tips for functional prediction of genes and pathways in mouse and other model vertebrates.

KERIS: kaleidoscope of gene responses to inflammation between species

A cornerstone of modern biomedical research is the use of animal models to study disease mechanisms and to develop new therapeutic approaches. In order to help the research community to better explore the similarities and differences of genomic response between human inflammatory diseases and murine models, we developed KERIS: kaleidoscope of gene responses to inflammation between species (available at http://www.igenomed.org/keris/). As of June 2016, KERIS includes comparisons of the genomic response of six human inflammatory diseases (burns, trauma, infection, sepsis, endotoxin and acute respiratory distress syndrome) and matched mouse models, using 2257 curated samples from the Inflammation and the Host Response to Injury Glue Grant studies and other representative studies in Gene Expression Omnibus. A researcher can browse, query, visualize and compare the response patterns of genes, pathways and functional modules across different diseases and corresponding murine models. The database is expected to help biologists choosing models when studying the mechanisms of particular genes and pathways in a disease and prioritizing the translation of findings from disease models into clinical studies.

PoplarGene: poplar gene network and resource for mining functional information for genes from woody plants

Poplar is not only an important resource for the production of paper, timber and other wood-based products, but it has also emerged as an ideal model system for studying woody plants. To better understand the biological processes underlying various traits in poplar, e.g., wood development, a comprehensive functional gene interaction network is highly needed. Here, we constructed a genome-wide functional gene network for poplar (covering ~70% of the 41,335 poplar genes) and created the network web service PoplarGene, offering comprehensive functional interactions and extensive poplar gene functional annotations. PoplarGene incorporates two network-based gene prioritization algorithms, neighborhood-based prioritization and context-based prioritization, which can be used to perform gene prioritization in a complementary manner. Furthermore, the co-functional information in PoplarGene can be applied to other woody plant proteomes with high efficiency via orthology transfer. In addition to poplar gene sequences, the webserver also accepts Arabidopsis reference gene as input to guide the search for novel candidate functional genes in PoplarGene. We believe that PoplarGene (http://bioinformatics.caf.ac.cn/PoplarGene and http://124.127.201.25/PoplarGene) will greatly benefit the research community, facilitating studies of poplar and other woody plants.