MetaboAnalyst 6.0: towards a unified platform for metabolomics data processing, analysis and interpretation.

We introduce MetaboAnalyst version 6.0 as a unified platform for processing, analyzing, and interpreting data from targeted as well as untargeted metabolomics studies using liquid chromatography - mass spectrometry (LC-MS). The two main objectives in developing version 6.0 are to support tandem MS (MS2) data processing and annotation, as well as to support the analysis of data from exposomics studies and related experiments. Key features of MetaboAnalyst 6.0 include: (i) a significantly enhanced Spectra Processing module with support for MS2 data and the asari algorithm; (ii) a MS2 Peak Annotation module based on comprehensive MS2 reference databases with fragment-level annotation; (iii) a new Statistical Analysis module dedicated for handling complex study design with multiple factors or phenotypic descriptors; (iv) a Causal Analysis module for estimating metabolite - phenotype causal relations based on two-sample Mendelian randomization, and (v) a Dose-Response Analysis module for benchmark dose calculations. In addition, we have also improved MetaboAnalyst's visualization functions, updated its compound database and metabolite sets, and significantly expanded its pathway analysis support to around 130 species. MetaboAnalyst 6.0 is freely available at https://www.metaboanalyst.ca.

MicrobiomeAnalyst 2.0: comprehensive statistical, functional and integrative analysis of microbiome data.

Microbiome studies have become routine in biomedical, agricultural and environmental sciences with diverse aims, including diversity profiling, functional characterization, and translational applications. The resulting complex, often multi-omics datasets demand powerful, yet user-friendly bioinformatics tools to reveal key patterns, important biomarkers, and potential activities. Here we introduce MicrobiomeAnalyst 2.0 to support comprehensive statistics, visualization, functional interpretation, and integrative analysis of data outputs commonly generated from microbiome studies. Compared to the previous version, MicrobiomeAnalyst 2.0 features three new modules: (i) a Raw Data Processing module for amplicon data processing and taxonomy annotation that connects directly with the Marker Data Profiling module for downstream statistical analysis; (ii) a Microbiome Metabolomics Profiling module to help dissect associations between community compositions and metabolic activities through joint analysis of paired microbiome and metabolomics datasets; and (iii) a Statistical Meta-Analysis module to help identify consistent signatures by integrating datasets across multiple studies. Other important improvements include added support for multi-factor differential analysis and interactive visualizations for popular graphical outputs, updated methods for functional prediction and correlation analysis, and expanded taxon set libraries based on the latest literature. These new features are demonstrated using a multi-omics dataset from a recent type 1 diabetes study. MicrobiomeAnalyst 2.0 is freely available at microbiomeanalyst.ca.

OmicsNet 2.0: a web-based platform for multi-omics integration and network visual analytics.

Researchers are increasingly seeking to interpret molecular data within a multi-omics context to gain a more comprehensive picture of their study system. OmicsNet (www.omicsnet.ca) is a web-based tool developed to allow users to easily build, visualize, and analyze multi-omics networks to study rich relationships among lists of 'omics features of interest. Three major improvements have been introduced in OmicsNet 2.0, which include: (i) enhanced network visual analytics with eleven 2D graph layout options and a novel 3D module layout; (ii) support for three new 'omics types: single nucleotide polymorphism (SNP) list from genetic variation studies; taxon list from microbiome profiling studies, as well as liquid chromatography-mass spectrometry (LC-MS) peaks from untargeted metabolomics; and (iii) measures to improve research reproducibility by coupling R command history with the release of the companion OmicsNetR package, and generation of persistent links to share interactive network views. We performed a case study using the multi-omics data obtained from a recent large-scale investigation on inflammatory bowel disease (IBD) and demonstrated that OmicsNet was able to quickly create meaningful multi-omics context to facilitate hypothesis generation and mechanistic insights.

Integrative systems immunology uncovers molecular networks of the cell cycle that stratify COVID-19 severity.

Several perturbations in the number of peripheral blood leukocytes, such as neutrophilia and lymphopenia associated with Coronavirus disease 2019 (COVID-19) severity, point to systemic molecular cell cycle alterations during severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection. However, the landscape of cell cycle alterations in COVID-19 remains primarily unexplored. Here, we performed an integrative systems immunology analysis of publicly available proteome and transcriptome data to characterize global changes in the cell cycle signature of COVID-19 patients. We found significantly enriched cell cycle-associated gene co-expression modules and an interconnected network of cell cycle-associated differentially expressed proteins (DEPs) and genes (DEGs) by integrating the molecular data of 1469 individuals (981 SARS-CoV-2 infected patients and 488 controls [either healthy controls or individuals with other respiratory illnesses]). Among these DEPs and DEGs are several cyclins, cell division cycles, cyclin-dependent kinases, and mini-chromosome maintenance proteins. COVID-19 patients partially shared the expression pattern of some cell cycle-associated genes with other respiratory illnesses but exhibited some specific differential features. Notably, the cell cycle signature predominated in the patients' blood leukocytes (B, T, and natural killer cells) and was associated with COVID-19 severity and disease trajectories. These results provide a unique global understanding of distinct alterations in cell cycle-associated molecules in COVID-19 patients, suggesting new putative pathways for therapeutic intervention.

OmicsAnalyst: a comprehensive web-based platform for visual analytics of multi-omics data.

Data analysis and interpretation remain a critical bottleneck in current multi-omics studies. Here, we introduce OmicsAnalyst, a user-friendly, web-based platform that allows users to perform a wide range of well-established data-driven approaches for multi-omics integration, and visually explore their results in a clear and meaningful manner. To help navigate complex landscapes of multi-omics analysis, these approaches are organized into three visual analytics tracks: (i) the correlation network analysis track, where users choose among univariate and multivariate methods to identify important features and explore their relationships in 2D or 3D networks; (ii) the cluster heatmap analysis track, where users apply several cutting-edge multi-view clustering algorithms and explore their results via interactive heatmaps; and (iii) the dimension reduction analysis track, where users choose among several recent multivariate techniques to reveal global data structures, and explore corresponding scores, loadings and biplots in interactive 3D scatter plots. The three visual analytics tracks are equipped with comprehensive options for parameter customization, view customization and targeted analysis. OmicsAnalyst lowers the access barriers to many well-established methods for multi-omics integration via novel visual analytics. It is freely available at https://www.omicsanalyst.ca.

MetaboAnalyst 5.0: narrowing the gap between raw spectra and functional insights.

Since its first release over a decade ago, the MetaboAnalyst web-based platform has become widely used for comprehensive metabolomics data analysis and interpretation. Here we introduce MetaboAnalyst version 5.0, aiming to narrow the gap from raw data to functional insights for global metabolomics based on high-resolution mass spectrometry (HRMS). Three modules have been developed to help achieve this goal, including: (i) a LC-MS Spectra Processing module which offers an easy-to-use pipeline that can perform automated parameter optimization and resumable analysis to significantly lower the barriers to LC-MS1 spectra processing; (ii) a Functional Analysis module which expands the previous MS Peaks to Pathways module to allow users to intuitively select any peak groups of interest and evaluate their enrichment of potential functions as defined by metabolic pathways and metabolite sets; (iii) a Functional Meta-Analysis module to combine multiple global metabolomics datasets obtained under complementary conditions or from similar studies to arrive at comprehensive functional insights. There are many other new functions including weighted joint-pathway analysis, data-driven network analysis, batch effect correction, merging technical replicates, improved compound name matching, etc. The web interface, graphics and underlying codebase have also been refactored to improve performance and user experience. At the end of an analysis session, users can now easily switch to other compatible modules for a more streamlined data analysis. MetaboAnalyst 5.0 is freely available at https://www.metaboanalyst.ca.

miRNet 2.0: network-based visual analytics for miRNA functional analysis and systems biology.

miRNet is an easy-to-use, web-based platform designed to help elucidate microRNA (miRNA) functions by integrating users' data with existing knowledge via network-based visual analytics. Since its first release in 2016, miRNet has been accessed by >20 000 researchers worldwide, with ∼100 users on a daily basis. While version 1.0 was focused primarily on miRNA-target gene interactions, it has become clear that in order to obtain a global view of miRNA functions, it is necessary to bring other important players into the context during analysis. Driven by this concept, in miRNet version 2.0, we have (i) added support for transcription factors (TFs) and single nucleotide polymorphisms (SNPs) that affect miRNAs, miRNA-binding sites or target genes, whilst also greatly increased (>5-fold) the underlying knowledgebases of miRNAs, ncRNAs and disease associations; (ii) implemented new functions to allow creation and visual exploration of multipartite networks, with enhanced support for in situ functional analysis and (iii) revamped the web interface, optimized the workflow, and introduced microservices and web application programming interface (API) to sustain high-performance, real-time data analysis. The underlying R package is also released in tandem with version 2.0 to allow more flexible data analysis for R programmers. The miRNet 2.0 website is freely available at https://www.mirnet.ca.

NetworkAnalyst 3.0: a visual analytics platform for comprehensive gene expression profiling and meta-analysis.

The growing application of gene expression profiling demands powerful yet user-friendly bioinformatics tools to support systems-level data understanding. NetworkAnalyst was first released in 2014 to address the key need for interpreting gene expression data within the context of protein-protein interaction (PPI) networks. It was soon updated for gene expression meta-analysis with improved workflow and performance. Over the years, NetworkAnalyst has been continuously updated based on community feedback and technology progresses. Users can now perform gene expression profiling for 17 different species. In addition to generic PPI networks, users can now create cell-type or tissue specific PPI networks, gene regulatory networks, gene co-expression networks as well as networks for toxicogenomics and pharmacogenomics studies. The resulting networks can be customized and explored in 2D, 3D as well as Virtual Reality (VR) space. For meta-analysis, users can now visually compare multiple gene lists through interactive heatmaps, enrichment networks, Venn diagrams or chord diagrams. In addition, users have the option to create their own data analysis projects, which can be saved and resumed at a later time. These new features are released together as NetworkAnalyst 3.0, freely available at https://www.networkanalyst.ca.

OmicsNet: a web-based tool for creation and visual analysis of biological networks in 3D space.

Biological networks play increasingly important roles in omics data integration and systems biology. Over the past decade, many excellent tools have been developed to support creation, analysis and visualization of biological networks. However, important limitations remain: most tools are standalone programs, the majority of them focus on protein-protein interaction (PPI) or metabolic networks, and visualizations often suffer from 'hairball' effects when networks become large. To help address these limitations, we developed OmicsNet - a novel web-based tool that allows users to easily create different types of molecular interaction networks and visually explore them in a three-dimensional (3D) space. Users can upload one or multiple lists of molecules of interest (genes/proteins, microRNAs, transcription factors or metabolites) to create and merge different types of biological networks. The 3D network visualization system was implemented using the powerful Web Graphics Library (WebGL) technology that works natively in most major browsers. OmicsNet supports force-directed layout, multi-layered perspective layout, as well as spherical layout to help visualize and navigate complex networks. A rich set of functions have been implemented to allow users to perform coloring, shading, topology analysis, and enrichment analysis. OmicsNet is freely available at http://www.omicsnet.ca.

MADS-box transcription factor Mcm1 controls cell cycle, fungal development, cell integrity and virulence in the filamentous insect pathogenic fungus Beauveria bassiana.

MADS-box transcription factor Mcm1 plays crucial roles in regulating mating processes and pathogenesis in some fungi. However, its roles are varied in fungal species, and its function remains unclear in entomopathogenic fungi. Here, Mcm1 orthologue, Bbmcm1, was characterized in a filamentous entomopathogenic fungus, Beauveria bassiana. Disruption of Bbmcm1 resulted in a distinct reduction in growth with abnormal conidiogenesis, and a significant decrease in conidial viability with abnormal germination. ΔBbmcm1 displayed impaired cell integrity, with distorted cell wall structure and altered cell wall component. Abnormal cell cycle was detected in ΔBbmcm1 with longer G2 /M phase but shorter G1 /G0 and S phases in unicellular blastospores, and sparser septa in multicellular hyphae, which might be responsible for defects in development and differentiation due to the regulation of cell cycle-involved genes, as well as the corresponding cellular events-associated genes. Dramatically decreased virulence was examined in ΔBbmcm1, with impaired ability to escape haemocyte encapsulation, which was consistent with markedly reduced cuticle-degrading enzyme production by repressing their coding genes, and downregulated fungal effector protein-coding genes, suggesting a novel role of Mcm1 in interaction with host insect. These data indicate that Mcm1 is a crucial regulator of development, cell integrity, cell cycle and virulence in insect fungal pathogens.

Evaluation of ligand-based NMR screening methods to characterize small molecule binding to HIV-1 glycoprotein-41.

Small molecule inhibitors of glycoprotein-41 (gp41) are able to prevent HIV infection by binding to a hydrophobic pocket (HP) contained within the gp41 ectodomain, and preventing progression of fusion. There is little structural information on gp41-ligand complexes, owing to hydrophobicity of the ligands, occlusion of the HP in folded gp41 ectodomain, and failure to form crystals of complexes. Here we used an engineered gp41 ectodomain protein containing an exposed HP and a small molecule designed to bind with weak affinity to the HP. We evaluated NMR methods, including WaterLOGSY, Saturation Transfer Difference spectroscopy (STD-NMR) and 1H relaxation rate difference spectroscopy with and without target irradiation (DIRECTION) for their ability to probe complex formation and structure. WaterLOGSY was the most sensitive technique for monitoring formation of the complex. STD-NMR and DIRECTION experiments gave similar pharmacophore mapping profiles, although the low dynamic range of the DIRECTION experiment limited its discrimination and sensitivity. A unique binding pose was identified from the STD data and provided clues for future optimization. Advantages and disadvantages of the techniques are discussed. This is the first example of the use of STD for structural analysis of a gp41-small molecule complex.

Synthesis and evaluation of the diarylthiourea analogs as novel anti-cancer agents.

Ten p-nitrodiarylthiourea analogs were designed, synthesized and evaluated in breast (MCF-7, T-47D, MDA-MB-453) and prostate (DU-145, PC-3, LNCaP) cancer cell lines for their anticancer activities. The majority of the compounds were able to inhibit the growth of these six cancer cell lines at low micromolar concentrations. Compound 7 was found to be the most potent anticancer agent in this series with GI50 values of 3.16µM for MCF-7, 2.53µM for T-47D, 4.77µM for MDA-MB-453 breast cancer lines and 3.54µM for LNCaP prostate cancer cell line. These GI50 values were comparable to the parent compound, SHetA2.

Web-based multi-omics integration using the Analyst software suite.

The growing number of multi-omics studies demands clear conceptual workflows coupled with easy-to-use software tools to facilitate data analysis and interpretation. This protocol covers three key components involved in multi-omics analysis, including single-omics data analysis, knowledge-driven integration using biological networks and data-driven integration through joint dimensionality reduction. Using the dataset from a recent multi-omics study of human pancreatic islet tissue and plasma samples, the first section introduces how to perform transcriptomics/proteomics data analysis using ExpressAnalyst and lipidomics data analysis using MetaboAnalyst. On the basis of significant features detected in these workflows, the second section demonstrates how to perform knowledge-driven integration using OmicsNet. The last section illustrates how to perform data-driven integration from the normalized omics data and metadata using OmicsAnalyst. The complete protocol can be executed in ~2 h. Compared with other available options for multi-omics integration, the Analyst software suite described in this protocol enables researchers to perform a wide range of omics data analysis tasks via a user-friendly web interface.

mGWAS-Explorer 2.0: Causal Analysis and Interpretation of Metabolite-Phenotype Associations.

Metabolomics-based genome-wide association studies (mGWAS) are key to understanding the genetic regulations of metabolites in complex phenotypes. We previously developed mGWAS-Explorer 1.0 to link single-nucleotide polymorphisms (SNPs), metabolites, genes and phenotypes for hypothesis generation. It has become clear that identifying potential causal relationships between metabolites and phenotypes, as well as providing deep functional insights, are crucial for further downstream applications. Here, we introduce mGWAS-Explorer 2.0 to support the causal analysis between >4000 metabolites and various phenotypes. The results can be interpreted within the context of semantic triples and molecular quantitative trait loci (QTL) data. The underlying R package is released for reproducible analysis. Using two case studies, we demonstrate that mGWAS-Explorer 2.0 is able to detect potential causal relationships between arachidonic acid and Crohn's disease, as well as between glycine and coronary heart disease.

Using ExpressAnalyst for Comprehensive Gene Expression Analysis in Model and Non-Model Organisms.

ExpressAnalyst is a web-based platform that enables intuitive, end-to-end transcriptomics and proteomics data analysis. Users can start from FASTQ files, gene/protein abundance tables, or gene/protein lists. ExpressAnalyst will perform read quantification, gene expression table processing and normalization, differential expression analysis, or meta-analysis with complex study designs. The results are presented via various interactive visualizations such as volcano plots, heatmaps, networks, and ridgeline charts, with built-in functional enrichment analysis to allow flexible data exploration and understanding. ExpressAnalyst currently contains built-in support for 29 common organisms. For non-model organisms without good reference genomes, it can perform comprehensive transcriptome profiling directly from RNA-seq reads. These common tasks are covered in 11 Basic Protocols. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: RNA-seq count table uploading, processing, and normalization Basic Protocol 2: Differential expression analysis with linear models Basic Protocol 3: Functional analysis with volcano plot, enrichment network, and ridgeline visualization Basic Protocol 4: Hierarchical clustering analysis of transcriptomics data using interactive heatmaps Basic Protocol 5: Cross-species gene expression analysis based on ortholog mapping results Basic Protocol 6: Proteomics and microarray data processing and normalization Basic Protocol 7: Preparing multiple gene expression tables for meta-analysis Basic Protocol 8: Statistical and functional meta-analysis of gene expression data Basic Protocol 9: Functional analysis of transcriptomics signatures Basic Protocol 10: Dose-response and time-series data analysis Basic Protocol 11: RNA-seq reads processing and quantification with and without reference transcriptomes.

ExpressAnalyst: A unified platform for RNA-sequencing analysis in non-model species.

The increasing application of RNA sequencing to study non-model species demands easy-to-use and efficient bioinformatics tools to help researchers quickly uncover biological and functional insights. We developed ExpressAnalyst ( www.expressanalyst.ca ), a web-based platform for processing, analyzing, and interpreting RNA-sequencing data from any eukaryotic species. ExpressAnalyst contains a series of modules that cover from processing and annotation of FASTQ files to statistical and functional analysis of count tables or gene lists. All modules are integrated with EcoOmicsDB, an ortholog database that enables comprehensive analysis for species without a reference transcriptome. By coupling ultra-fast read mapping algorithms with high-resolution ortholog databases through a user-friendly web interface, ExpressAnalyst allows researchers to obtain global expression profiles and gene-level insights from raw RNA-sequencing reads within 24 h. Here, we present ExpressAnalyst and demonstrate its utility with a case study of RNA-sequencing data from multiple non-model salamander species, including two that do not have a reference transcriptome.

Remodeling bacterial polysaccharides by metabolic pathway engineering.

Introducing structural modifications into biomolecules represents a powerful approach to dissect their functions and roles in biological processes. Bacterial polysaccharides, despite their rich structural information and essential roles in bacterium-host interactions and bacterial virulence, have largely been unexplored for in vivo structural modifications. In this study, we demonstrate the incorporation of a panel of monosaccharide analogs into bacterial polysaccharides in a highly homogenous manner via metabolic engineering of a promiscuous sugar nucleotide biosynthetic pathway. In addition, the bioorthorgonal functional groups metabolically incorporated were exploited for cell surface labeling using in vitro selective chemical ligation reactions. In summary, our study presents a general, facile and effective approach for in vivo generation of novel tailor-made bacterial polysaccharides.

mGWAS-Explorer: Linking SNPs, Genes, Metabolites, and Diseases for Functional Insights.

Tens of thousands of single-nucleotide polymorphisms (SNPs) have been identified to be significantly associated with metabolite abundance in over 65 genome-wide association studies with metabolomics (mGWAS) to date. Obtaining mechanistic or functional insights from these associations for translational applications has become a key research area in the mGWAS community. Here, we introduce mGWAS-Explorer, a user-friendly web-based platform to help connect SNPs, metabolites, genes, and their known disease associations via powerful network visual analytics. The application of the mGWAS-Explorer was demonstrated using a COVID-19 and a type 2 diabetes case studies.

Using MetaboAnalyst 5.0 for LC-HRMS spectra processing, multi-omics integration and covariate adjustment of global metabolomics data.

Liquid chromatography coupled with high-resolution mass spectrometry (LC-HRMS) has become a workhorse in global metabolomics studies with growing applications across biomedical and environmental sciences. However, outstanding bioinformatics challenges in terms of data processing, statistical analysis and functional interpretation remain critical barriers to the wider adoption of this technology. To help the user community overcome these barriers, we have made major updates to the well-established MetaboAnalyst platform ( www.metaboanalyst.ca ). This protocol extends the previous 2011 Nature Protocol by providing stepwise instructions on how to use MetaboAnalyst 5.0 to: optimize parameters for LC-HRMS spectra processing; obtain functional insights from peak list data; integrate metabolomics data with transcriptomics data or combine multiple metabolomics datasets; conduct exploratory statistical analysis with complex metadata. Parameter optimization may take ~2 h to complete depending on the server load, and the remaining three stages may be executed in ~60 min.

Targeting a Conserved Lysine in the Hydrophobic Pocket of HIV-1 gp41 Improves Small Molecule Antiviral Activity.

Human Immunodeficiency virus (HIV-1) fusion is mediated by glycoprotein-41, a protein that has not been widely exploited as a drug target. Small molecules directed at the gp41 ectodomain have proved to be poorly drug-like, having moderate efficacy, high hydrophobicity and/or high molecular weight. We recently investigated conversion of a fairly potent hydrophobic inhibitor into a covalent binder, by modifying it to react with a lysine residue on the protein. We demonstrated a 10-fold improvement in antiviral efficacy. Here, we continue this study, utilizing instead molecules with better inherent drug-like properties. Molecules possessing low to no antiviral activity as equilibrium binders were converted into µM inhibitors upon addition of an electrophilic warhead in the form of a sulfotetrafluorophenyl (STP) activated ester. We confirmed specificity for gp41 and for entry. The small size of the inhibitors described here offers an opportunity to expand their reach into neighboring pockets while retaining drug-likeness. STP esterification of equilibrium binders is a promising avenue to explore for inhibiting HIV-1 entry. Many gp41 targeting molecules studied over the years possess carboxylic acid groups which can be easily converted into the corresponding STP ester. It may be worth the effort to evaluate a library of such inhibitors as a way forward to small molecule inhibition of fusion of HIV and possibly other enveloped viruses.