AutophagyNet: high-resolution data source for the analysis of autophagy and its regulation.

Macroautophagy/autophagy is a highly-conserved catabolic procss eliminating dysfunctional cellular components and invading pathogens. Autophagy malfunction contributes to disorders such as cancer, neurodegenerative and inflammatory diseases. Understanding autophagy regulation in health and disease has been the focus of the last decades. We previously provided an integrated database for autophagy research, the Autophagy Regulatory Network (ARN). For the last eight years, this resource has been used by thousands of users. Here, we present a new and upgraded resource, AutophagyNet. It builds on the previous database but contains major improvements to address user feedback and novel needs due to the advancement in omics data availability. AutophagyNet contains updated interaction curation and integration of over 280,000 experimentally verified interactions between core autophagy proteins and their protein, transcriptional and post-transcriptional regulators as well as their potential upstream pathway connections. AutophagyNet provides annotations for each core protein about their role: 1) in different types of autophagy (mitophagy, xenophagy, etc.); 2) in distinct stages of autophagy (initiation, expansion, termination, etc.); 3) with subcellular and tissue-specific localization. These annotations can be used to filter the dataset, providing customizable download options tailored to the user's needs. The resource is available in various file formats (e.g. CSV, BioPAX and PSI-MI), and data can be analyzed and visualized directly in Cytoscape. The multi-layered regulation of autophagy can be analyzed by combining AutophagyNet with tissue- or cell type-specific (multi-)omics datasets (e.g. transcriptomic or proteomic data). The resource is publicly accessible at http://autophagynet.org.Abbreviations: ARN: Autophagy Regulatory Network; ATG: autophagy related; BCR: B cell receptor pathway; BECN1: beclin 1; GABARAP: GABA type A receptor-associated protein; IIP: innate immune pathway; LIR: LC3-interacting region; lncRNA: long non-coding RNA; MAP1LC3B: microtubule associated protein 1 light chain 3 beta; miRNA: microRNA; NHR: nuclear hormone receptor; PTM: post-translational modification; RTK: receptor tyrosine kinase; TCR: T cell receptor; TLR: toll like receptor.

Immune checkpoint inhibitor-induced colitis is mediated by polyfunctional lymphocytes and is dependent on an IL23/IFNγ axis.

Immune checkpoint inhibitors (CPIs) are a relatively newly licenced cancer treatment, which make a once previously untreatable disease now amenable to a potential cure. Combination regimens of anti-CTLA4 and anti-PD-1 show enhanced efficacy but are prone to off-target immune-mediated tissue injury, particularly at the barrier surfaces. To probe the impact of immune checkpoints on intestinal homoeostasis, mice are challenged with anti-CTLA4 and anti-PD-1 immunotherapy and manipulation of the intestinal microbiota. The immune profile of the colon of these mice with CPI-colitis is analysed using bulk RNA sequencing, single-cell RNA sequencing and flow cytometry. CPI-colitis in mice is dependent on the composition of the intestinal microbiota and by the induction of lymphocytes expressing interferon-γ (IFNγ), cytotoxicity molecules and other pro-inflammatory cytokines/chemokines. This pre-clinical model of CPI-colitis could be attenuated following blockade of the IL23/IFNγ axis. Therapeutic targeting of IFNγ-producing lymphocytes or regulatory networks, may hold the key to reversing CPI-colitis.

The gut microbiota and metabolome are associated with diminished COVID-19 vaccine-induced antibody responses in immunosuppressed inflammatory bowel disease patients.

BACKGROUND: Patients with inflammatory bowel disease (IBD) treated with anti-TNF therapy exhibit attenuated humoral immune responses to vaccination against SARS-CoV-2. The gut microbiota and its functional metabolic output, which are perturbed in IBD, play an important role in shaping host immune responses. We explored whether the gut microbiota and metabolome could explain variation in anti-SARS-CoV-2 vaccination responses in immunosuppressed IBD patients. METHODS: Faecal and serum samples were prospectively collected from infliximab-treated patients with IBD in the CLARITY-IBD study undergoing vaccination against SARS-CoV-2. Antibody responses were measured following two doses of either ChAdOx1 nCoV-19 or BNT162b2 vaccine. Patients were classified as having responses above or below the geometric mean of the wider CLARITY-IBD cohort. 16S rRNA gene amplicon sequencing, nuclear magnetic resonance (NMR) spectroscopy and bile acid profiling with ultra-high-performance liquid chromatography mass spectrometry (UHPLC-MS) were performed on faecal samples. Univariate, multivariable and correlation analyses were performed to determine gut microbial and metabolomic predictors of response to vaccination. FINDINGS: Forty-three infliximab-treated patients with IBD were recruited (30 Crohn's disease, 12 ulcerative colitis, 1 IBD-unclassified; 26 with concomitant thiopurine therapy). Eight patients had evidence of prior SARS-CoV-2 infection. Seventeen patients (39.5%) had a serological response below the geometric mean. Gut microbiota diversity was lower in below average responders (p = 0.037). Bilophila abundance was associated with better serological response, while Streptococcus was associated with poorer response. The faecal metabolome was distinct between above and below average responders (OPLS-DA R2X 0.25, R2Y 0.26, Q2 0.15; CV-ANOVA p = 0.038). Trimethylamine, isobutyrate and omega-muricholic acid were associated with better response, while succinate, phenylalanine, taurolithocholate and taurodeoxycholate were associated with poorer response. INTERPRETATION: Our data suggest that there is an association between the gut microbiota and variable serological response to vaccination against SARS-CoV-2 in immunocompromised patients. Microbial metabolites including trimethylamine may be important in mitigating anti-TNF-induced attenuation of the immune response. FUNDING: JLA is the recipient of an NIHR Academic Clinical Lectureship (CL-2019-21-502), funded by Imperial College London and The Joyce and Norman Freed Charitable Trust. BHM is the recipient of an NIHR Academic Clinical Lectureship (CL-2019-21-002). The Division of Digestive Diseases at Imperial College London receives financial and infrastructure support from the NIHR Imperial Biomedical Research Centre (BRC) based at Imperial College Healthcare NHS Trust and Imperial College London. Metabolomics studies were performed at the MRC-NIHR National Phenome Centre at Imperial College London; this work was supported by the Medical Research Council (MRC), the National Institute of Health Research (NIHR) (grant number MC_PC_12025) and infrastructure support was provided by the NIHR Imperial Biomedical Research Centre (BRC). The NIHR Exeter Clinical Research Facility is a partnership between the University of Exeter Medical School College of Medicine and Health, and Royal Devon and Exeter NHS Foundation Trust. This project is supported by the National Institute for Health Research (NIHR) Exeter Clinical Research Facility. The views expressed are those of the authors and not necessarily those of the NIHR or the UK Department of Health and Social Care.

Multilayered Networks of SalmoNet2 Enable Strain Comparisons of the Salmonella Genus on a Molecular Level.

Serovars of the genus Salmonella primarily evolved as gastrointestinal pathogens in a wide range of hosts. Some serotypes later evolved further, adopting a more invasive lifestyle in a narrower host range associated with systemic infections. A system-level knowledge of these pathogens could identify the complex adaptations associated with the evolution of serovars with distinct pathogenicity, host range, and risk to human health. This promises to aid the design of interventions and serve as a knowledge base in the Salmonella research community. Here, we present SalmoNet2, a major update to SalmoNet1, the first multilayered interaction resource for Salmonella strains, containing protein-protein, transcriptional regulatory, and enzyme-enzyme interactions. The new version extends the number of Salmonella networks from 11 to 20. We now include a strain from the second species in the Salmonella genus, a strain from the Salmonella enterica subspecies arizonae and additional strains of importance from the subspecies enterica, including S. Typhimurium strain D23580, an epidemic multidrug-resistant strain associated with invasive nontyphoidal salmonellosis (iNTS). The database now uses strain specific metabolic models instead of a generalized model to highlight differences between strains. The update has increased the coverage of high-quality protein-protein interactions, and enhanced interoperability with other computational resources by adopting standardized formats. The resource website has been updated with tutorials to help researchers analyze their Salmonella data using molecular interaction networks from SalmoNet2. SalmoNet2 is accessible at http://salmonet.org/. IMPORTANCE Multilayered network databases collate interaction information from multiple sources, and are powerful both as a knowledge base and subject of analysis. Here, we present SalmoNet2, an integrated network resource containing protein-protein, transcriptional regulatory, and metabolic interactions for 20 Salmonella strains. Key improvements to the update include expanding the number of strains, strain-specific metabolic networks, an increase in high-quality protein-protein interactions, community standard computational formats to help interoperability, and online tutorials to help users analyze their data using SalmoNet2.

Mapping the epithelial-immune cell interactome upon infection in the gut and the upper airways.

Increasing evidence points towards the key role of the epithelium in the systemic and over-activated immune response to viral infection, including SARS-CoV-2 infection. Yet, how viral infection alters epithelial-immune cell interactions regulating inflammatory responses, is not well known. Available experimental approaches are insufficient to properly analyse this complex system, and computational predictions and targeted data integration are needed as an alternative approach. In this work, we propose an integrated computational biology framework that models how infection alters intracellular signalling of epithelial cells and how this change impacts the systemic immune response through modified interactions between epithelial cells and local immune cell populations. As a proof-of-concept, we focused on the role of intestinal and upper-airway epithelial infection. To characterise the modified epithelial-immune interactome, we integrated intra- and intercellular networks with single-cell RNA-seq data from SARS-CoV-2 infected human ileal and colonic organoids as well as from infected airway ciliated epithelial cells. This integrated methodology has proven useful to point out specific epithelial-immune interactions driving inflammation during disease response, and propose relevant molecular targets to guide focused experimental analysis.

Analysing miRNA-Target Gene Networks in Inflammatory Bowel Disease and Other Complex Diseases Using Transcriptomic Data.

Patients with inflammatory bowel disease (IBD) are known to have perturbations in microRNA (miRNA) levels as well as altered miRNA regulation. Although experimental methods have provided initial insights into the functional consequences that may arise due to these changes, researchers are increasingly utilising novel bioinformatics approaches to further dissect the role of miRNAs in IBD. The recent exponential increase in transcriptomics datasets provides an excellent opportunity to further explore the role of miRNAs in IBD pathogenesis. To effectively understand miRNA-target gene interactions from gene expression data, multiple database resources are required, which have become available in recent years. In this technical note, we provide a step-by-step protocol for utilising these state-of-the-art resources, as well as systems biology approaches to understand the role of miRNAs in complex disease pathogenesis. We demonstrate through a case study example how to combine the resulting miRNA-target gene networks with transcriptomics data to find potential disease-specific miRNA regulators and miRNA-target genes in Crohn's disease. This approach could help to identify miRNAs that may have important disease-modifying effects in IBD and other complex disorders, and facilitate the discovery of novel therapeutic targets.

SignaLink3: a multi-layered resource to uncover tissue-specific signaling networks.

Signaling networks represent the molecular mechanisms controlling a cell's response to various internal or external stimuli. Most currently available signaling databases contain only a part of the complex network of intertwining pathways, leaving out key interactions or processes. Hence, we have developed SignaLink3 (http://signalink.org/), a value-added knowledge-base that provides manually curated data on signaling pathways and integrated data from several types of databases (interaction, regulation, localisation, disease, etc.) for humans, and three major animal model organisms. SignaLink3 contains over 400 000 newly added human protein-protein interactions resulting in a total of 700 000 interactions for Homo sapiens, making it one of the largest integrated signaling network resources. Next to H. sapiens, SignaLink3 is the only current signaling network resource to provide regulatory information for the model species Caenorhabditis elegans and Danio rerio, and the largest resource for Drosophila melanogaster. Compared to previous versions, we have integrated gene expression data as well as subcellular localization of the interactors, therefore uniquely allowing tissue-, or compartment-specific pathway interaction analysis to create more accurate models. Data is freely available for download in widely used formats, including CSV, PSI-MI TAB or SQL.

CytokineLink: A Cytokine Communication Map to Analyse Immune Responses-Case Studies in Inflammatory Bowel Disease and COVID-19.

Intercellular communication mediated by cytokines is critical to the development of immune responses, particularly in the context of infectious and inflammatory diseases. By releasing these small molecular weight peptides, the source cells can influence numerous intracellular processes in the target cells, including the secretion of other cytokines downstream. However, there are no readily available bioinformatic resources that can model cytokine-cytokine interactions. In this effort, we built a communication map between major tissues and blood cells that reveals how cytokine-mediated intercellular networks form during homeostatic conditions. We collated the most prevalent cytokines from the literature and assigned the proteins and their corresponding receptors to source tissue and blood cell types based on enriched consensus RNA-Seq data from the Human Protein Atlas database. To assign more confidence to the interactions, we integrated the literature information on cell-cytokine interactions from two systems of immunology databases, immuneXpresso and ImmunoGlobe. From the collated information, we defined two metanetworks: a cell-cell communication network connected by cytokines; and a cytokine-cytokine interaction network depicting the potential ways in which cytokines can affect the activity of each other. Using expression data from disease states, we then applied this resource to reveal perturbations in cytokine-mediated intercellular signalling in inflammatory and infectious diseases (ulcerative colitis and COVID-19, respectively). For ulcerative colitis, with CytokineLink, we demonstrated a significant rewiring of cytokine-mediated intercellular communication between non-inflamed and inflamed colonic tissues. For COVID-19, we were able to identify cell types and cytokine interactions following SARS-CoV-2 infection, highlighting important cytokine interactions that might contribute to severe illness in a subgroup of patients. Such findings have the potential to inform the development of novel, cytokine-targeted therapeutic strategies. CytokineLink is freely available for the scientific community through the NDEx platform and the project github repository.

SARS-CoV-2 Causes a Different Cytokine Response Compared to Other Cytokine Storm-Causing Respiratory Viruses in Severely Ill Patients.

Hyper-induction of pro-inflammatory cytokines, also known as a cytokine storm or cytokine release syndrome (CRS), is one of the key aspects of the currently ongoing SARS-CoV-2 pandemic. This process occurs when a large number of innate and adaptive immune cells activate and start producing pro-inflammatory cytokines, establishing an exacerbated feedback loop of inflammation. It is one of the factors contributing to the mortality observed with coronavirus 2019 (COVID-19) for a subgroup of patients. CRS is not unique to the SARS-CoV-2 infection; it was prevalent in most of the major human coronavirus and influenza A subtype outbreaks of the past two decades (H5N1, SARS-CoV, MERS-CoV, and H7N9). With a comprehensive literature search, we collected changing the cytokine levels from patients upon infection with the viral pathogens mentioned above. We analyzed published patient data to highlight the conserved and unique cytokine responses caused by these viruses. Our curation indicates that the cytokine response induced by SARS-CoV-2 is different compared to other CRS-causing respiratory viruses, as SARS-CoV-2 does not always induce specific cytokines like other coronaviruses or influenza do, such as IL-2, IL-10, IL-4, or IL-5. Comparing the collated cytokine responses caused by the analyzed viruses highlights a SARS-CoV-2-specific dysregulation of the type-I interferon (IFN) response and its downstream cytokine signatures. The map of responses gathered in this study could help specialists identify interventions that alleviate CRS in different diseases and evaluate whether they could be used in the COVID-19 cases.

Integrated intra- and intercellular signaling knowledge for multicellular omics analysis.

Molecular knowledge of biological processes is a cornerstone in omics data analysis. Applied to single-cell data, such analyses provide mechanistic insights into individual cells and their interactions. However, knowledge of intercellular communication is scarce, scattered across resources, and not linked to intracellular processes. To address this gap, we combined over 100 resources covering interactions and roles of proteins in inter- and intracellular signaling, as well as transcriptional and post-transcriptional regulation. We added protein complex information and annotations on function, localization, and role in diseases for each protein. The resource is available for human, and via homology translation for mouse and rat. The data are accessible via OmniPath's web service (https://omnipathdb.org/), a Cytoscape plug-in, and packages in R/Bioconductor and Python, providing access options for computational and experimental scientists. We created workflows with tutorials to facilitate the analysis of cell-cell interactions and affected downstream intracellular signaling processes. OmniPath provides a single access point to knowledge spanning intra- and intercellular processes for data analysis, as we demonstrate in applications studying SARS-CoV-2 infection and ulcerative colitis.

ViralLink: An integrated workflow to investigate the effect of SARS-CoV-2 on intracellular signalling and regulatory pathways.

The SARS-CoV-2 pandemic of 2020 has mobilised scientists around the globe to research all aspects of the coronavirus virus and its infection. For fruitful and rapid investigation of viral pathomechanisms, a collaborative and interdisciplinary approach is required. Therefore, we have developed ViralLink: a systems biology workflow which reconstructs and analyses networks representing the effect of viruses on intracellular signalling. These networks trace the flow of signal from intracellular viral proteins through their human binding proteins and downstream signalling pathways, ending with transcription factors regulating genes differentially expressed upon viral exposure. In this way, the workflow provides a mechanistic insight from previously identified knowledge of virally infected cells. By default, the workflow is set up to analyse the intracellular effects of SARS-CoV-2, requiring only transcriptomics counts data as input from the user: thus, encouraging and enabling rapid multidisciplinary research. However, the wide-ranging applicability and modularity of the workflow facilitates customisation of viral context, a priori interactions and analysis methods. Through a case study of SARS-CoV-2 infected bronchial/tracheal epithelial cells, we evidence the functionality of the workflow and its ability to identify key pathways and proteins in the cellular response to infection. The application of ViralLink to different viral infections in a context specific manner using different available transcriptomics datasets will uncover key mechanisms in viral pathogenesis.

Evolution of regulatory networks associated with traits under selection in cichlids.

BACKGROUND: Seminal studies of vertebrate protein evolution speculated that gene regulatory changes can drive anatomical innovations. However, very little is known about gene regulatory network (GRN) evolution associated with phenotypic effect across ecologically diverse species. Here we use a novel approach for comparative GRN analysis in vertebrate species to study GRN evolution in representative species of the most striking examples of adaptive radiations, the East African cichlids. We previously demonstrated how the explosive phenotypic diversification of East African cichlids can be attributed to diverse molecular mechanisms, including accelerated regulatory sequence evolution and gene expression divergence. RESULTS: To investigate these mechanisms across species at a genome-wide scale, we develop a novel computational pipeline that predicts regulators for co-extant and ancestral co-expression modules along a phylogeny, and candidate regulatory regions associated with traits under selection in cichlids. As a case study, we apply our approach to a well-studied adaptive trait-the visual system-for which we report striking cases of network rewiring for visual opsin genes, identify discrete regulatory variants, and investigate their association with cichlid visual system evolution. In regulatory regions of visual opsin genes, in vitro assays confirm that transcription factor binding site mutations disrupt regulatory edges across species and segregate according to lake species phylogeny and ecology, suggesting GRN rewiring in radiating cichlids. CONCLUSIONS: Our approach reveals numerous novel potential candidate regulators and regulatory regions across cichlid genomes, including some novel and some previously reported associations to known adaptive evolutionary traits.

Sherlock: an open-source data platform to store, analyze and integrate Big Data for computational biologists.

In the era of Big Data, data collection underpins biological research more than ever before. In many cases, this can be as time-consuming as the analysis itself. It requires downloading multiple public databases with various data structures, and in general, spending days preparing the data before answering any biological questions. Here, we introduce Sherlock, an open-source, cloud-based big data platform ( https://earlham-sherlock.github.io/) to solve this problem. Sherlock provides a gap-filling way for computational biologists to store, convert, query, share and generate biology data while ultimately streamlining bioinformatics data management. The Sherlock platform offers a simple interface to leverage big data technologies, such as Docker and PrestoDB. Sherlock is designed to enable users to analyze, process, query and extract information from extremely complex and large data sets. Furthermore, Sherlock can handle different structured data (interaction, localization, or genomic sequence) from several sources and convert them to a common optimized storage format, for example, the Optimized Row Columnar (ORC). This format facilitates Sherlock's ability to quickly and efficiently execute distributed analytical queries on extremely large data files and share datasets between teams. The Sherlock platform is freely available on GitHub, and contains specific loader scripts for structured data sources of genomics, interaction and expression databases. With these loader scripts, users can easily and quickly create and work with specific file formats, such as JavaScript Object Notation (JSON) or ORC. For computational biology and large-scale bioinformatics projects, Sherlock provides an open-source platform empowering data management, analytics, integration and collaboration through modern big data technologies.

ULK1 and ULK2 are less redundant than previously thought: computational analysis uncovers distinct regulation and functions of these autophagy induction proteins.

Macroautophagy, the degradation of cytoplasmic content by lysosomal fusion, is an evolutionary conserved process promoting homeostasis and intracellular defence. Macroautophagy is initiated primarily by a complex containing ULK1 or ULK2 (two paralogs of the yeast Atg1 protein). To understand the differences between ULK1 and ULK2, we compared the human ULK1 and ULK2 proteins and their regulation. Despite the similarity in their enzymatic domain, we found that ULK1 and ULK2 have major differences in their autophagy-related interactors and their post-translational and transcriptional regulators. We identified 18 ULK1-specific and 7 ULK2-specific protein motifs serving as different interaction interfaces. We found that interactors of ULK1 and ULK2 all have different tissue-specific expressions partially contributing to diverse and ULK-specific interaction networks in various tissues. We identified three ULK1-specific and one ULK2-specific transcription factor binding sites, and eight sites shared by the regulatory region of both genes. Importantly, we found that both their post-translational and transcriptional regulators are involved in distinct biological processes-suggesting separate functions for ULK1 and ULK2. Unravelling differences between ULK1 and ULK2 could lead to a better understanding of how ULK-type specific dysregulation affects autophagy and other cellular processes that have been implicated in diseases such as inflammatory bowel disease and cancer.

Network Biology Approaches to Identify Molecular and Systems-Level Differences Between Salmonella Pathovars.

The field of systems biology endeavors to map, study, and simulate cellular systems and their underlying mechanisms. The internal mechanisms of biological systems can be represented with networks comprising nodes and edges. Nodes denote the constituents of the biological system whereas edges indicate the relationships among them. Likewise, every layer of cellular organization can be represented by networks. Multilayered networks capture interactions between various network types, such as transcriptional regulatory networks, protein-protein interaction networks, and metabolic networks from the same biological system. This property makes multilayered networks representative of the system while its internal mechanisms are investigated. However, there are not many multilayered networks containing integrated data for nonmodel organisms including the bacterial pathogens Salmonella. Here, we outline the steps to create such an integrated network database, through the example of SalmoNet, the first integrated multilayered data resource for multiple strains belonging to distinct Salmonella serovars.

SignaLink: Multilayered Regulatory Networks.

Biological networks are graphs used to represent the inner workings of a biological system. Networks describe the relationships of the elements of biological systems using edges and nodes. However, the resulting representation of the system can sometimes be too simplistic to usefully model reality. By combining several different interaction types within one larger multilayered biological network, tools such as SignaLink provide a more nuanced view than those relying on single-layer networks (where edges only describe one kind of interaction). Multilayered networks display connections between multiple networks (i.e., protein-protein interactions and their transcriptional and posttranscriptional regulators), each one of them describing a specific set of connections. Multilayered networks also allow us to depict cross talk between cellular systems, which is a more realistic way of describing molecular interactions. They can be used to collate networks from different sources into one multilayered structure, which makes them useful as an analytic tool as well.

SalmoNet, an integrated network of ten Salmonella enterica strains reveals common and distinct pathways to host adaptation.

Salmonella enterica is a prominent bacterial pathogen with implications on human and animal health. Salmonella serovars could be classified as gastro-intestinal or extra-intestinal. Genome-wide comparisons revealed that extra-intestinal strains are closer relatives of gastro-intestinal strains than to each other indicating a parallel evolution of this trait. Given the complexity of the differences, a systems-level comparison could reveal key mechanisms enabling extra-intestinal serovars to cause systemic infections. Accordingly, in this work, we introduce a unique resource, SalmoNet, which combines manual curation, high-throughput data and computational predictions to provide an integrated network for Salmonella at the metabolic, transcriptional regulatory and protein-protein interaction levels. SalmoNet provides the networks separately for five gastro-intestinal and five extra-intestinal strains. As a multi-layered, multi-strain database containing experimental data, SalmoNet is the first dedicated network resource for Salmonella. It comprehensively contains interactions between proteins encoded in Salmonella pathogenicity islands, as well as regulatory mechanisms of metabolic processes with the option to zoom-in and analyze the interactions at specific loci in more detail. Application of SalmoNet is not limited to strain comparisons as it also provides a Salmonella resource for biochemical network modeling, host-pathogen interaction studies, drug discovery, experimental validation of novel interactions, uncovering new pathological mechanisms from emergent properties and epidemiological studies. SalmoNet is available at http://salmonet.org.