PathwayPCA: an R/Bioconductor Package for Pathway Based Integrative Analysis of Multi-Omics Data

Integrated Systems Analysis Deciphers Transcriptome and Glycoproteome Links in Alzheimer's Disease

Glycosylation is increasingly recognized as a potential therapeutic target in Alzheimer's disease. In recent years, evidence of Alzheimer's disease-specific glycoproteins has been established. However, the mechanisms underlying their dysregulation, including tissue- and cell-type specificity, are not fully understood. We aimed to explore the upstream regulators of aberrant glycosylation by integrating multiple data sources using a glycogenomics approach. We identified dysregulation of the glycosyltransferase PLOD3 in oligodendrocytes as an upstream regulator of cerebral vessels and found that it is involved in COL4A5 synthesis, which is strongly correlated with amyloid fiber formation. Furthermore, COL4A5 has been suggested to interact with astrocytes via extracellular matrix receptors as a ligand. This study suggests directions for new therapeutic strategies for Alzheimer's disease targeting glycosyltransferases.

Molecular Insights into Transcranial Direct Current Stimulation Effects: Metabolomics and Transcriptomics Analyses

INTRODUCTION

Transcranial direct current stimulation (tDCS) is an evolving non-invasive neurostimulation technique. Despite multiple studies, its underlying molecular mechanisms are still unclear. Several previous human studies of the effect of tDCS suggest that it generates metabolic effects. The induction of metabolic effects by tDCS could provide an explanation for how it generates its long-term beneficial clinical outcome.

AIM

Given these hints of tDCS metabolic effects, we aimed to delineate the metabolic pathways involved in its mode of action.

METHODS

To accomplish this, we utilized a broad analytical approach of co-analyzing metabolomics and transcriptomic data generated from anodal tDCS in rat models. Since no metabolomic dataset was available, we performed a tDCS experiment of bilateral anodal stimulation of 200 µA for 20 min and for 5 consecutive days, followed by harvesting the brain tissue below the stimulating electrode and generating a metabolomics dataset using LC-MS/MS. The analysis of the transcriptomic dataset was based on a publicly available dataset.

RESULTS

Our analyses revealed that tDCS alters the metabolic profile of brain tissue, affecting bioenergetic-related pathways, such as glycolysis and mitochondrial functioning. In addition, we found changes in calcium-related signaling.

CONCLUSIONS

We conclude that tDCS affects metabolism by modulating energy production-related processes. Given our findings concerning calcium-related signaling, we suggest that the immediate effects of tDCS on calcium dynamics drive modifications in distinct metabolic pathways. A thorough understanding of the underlying molecular mechanisms of tDCS has the potential to revolutionize its applicability, enabling the generation of personalized medicine in the field of neurostimulation and thus contributing to its optimization.

CLIN_SKAT: an R package to conduct association analysis using functionally relevant variants

Background

Availability of next generation sequencing data, allows low-frequency and rare variants to be studied through strategies other than the commonly used genome-wide association studies (GWAS). Rare variants are important keys towards explaining the heritability for complex diseases that remains to be explained by common variants due to their low effect sizes. However, analysis strategies struggle to keep up with the huge amount of data at disposal therefore creating a bottleneck. This study describes CLIN_SKAT, an R package, that provides users with an easily implemented analysis pipeline with the goal of (i) extracting clinically relevant variants (both rare and common), followed by (ii) gene-based association analysis by grouping the selected variants.

Results

CLIN_SKAT offers four simple functions that can be used to obtain clinically relevant variants, map them to genes or gene sets, calculate weights from global healthy populations and conduct weighted case–control analysis. CLIN_SKAT introduces improvements by adding certain pre-analysis steps and customizable features to make the SKAT results clinically more meaningful. Moreover, it offers several plot functions that can be availed towards obtaining visualizations for interpretation of the analyses results. CLIN_SKAT is available on Windows/Linux/MacOS and is operative for R version 4.0.4 or later. It can be freely downloaded from https://github.com/ShihChingYu/CLIN_SKAT, installed through devtools::install_github("ShihChingYu/CLIN_SKAT", force=T) and executed by loading the package into R using library(CLIN_SKAT). All outputs (tabular and graphical) can be downloaded in simple, publishable formats.

Conclusions

Statistical association analysis is often underpowered due to low sample sizes and high numbers of variants to be tested, limiting detection of causal ones. Therefore, retaining a subset of variants that are biologically meaningful seems to be a more effective strategy for identifying explainable associations while reducing the degrees of freedom. CLIN_SKAT offers users a one-stop R package that identifies disease risk variants with improved power via a series of tailor-made procedures that allows dimension reduction, by retaining functionally relevant variants, and incorporating ethnicity based priors. Furthermore, it also eliminates the requirement for high computational resources and bioinformatics expertise.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04987-2.

PathwayMultiomics: An R Package for Efficient Integrative Analysis of Multi-Omics Datasets With Matched or Un-matched Samples

Recent advances in technology have made multi-omics datasets increasingly available to researchers. To leverage the wealth of information in multi-omics data, a number of integrative analysis strategies have been proposed recently. However, effectively extracting biological insights from these large, complex datasets remains challenging. In particular, matched samples with multiple types of omics data measured on each sample are often required for multi-omics analysis tools, which can significantly reduce the sample size. Another challenge is that analysis techniques such as dimension reductions, which extract association signals in high dimensional datasets by estimating a few variables that explain most of the variations in the samples, are typically applied to whole-genome data, which can be computationally demanding. Here we present pathwayMultiomics, a pathway-based approach for integrative analysis of multi-omics data with categorical, continuous, or survival outcome variables. The input of pathwayMultiomics is pathway p-values for individual omics data types, which are then integrated using a novel statistic, the MiniMax statistic, to prioritize pathways dysregulated in multiple types of omics datasets. Importantly, pathwayMultiomics is computationally efficient and does not require matched samples in multi-omics data. We performed a comprehensive simulation study to show that pathwayMultiomics significantly outperformed currently available multi-omics tools with improved power and well-controlled false-positive rates. In addition, we also analyzed real multi-omics datasets to show that pathwayMultiomics was able to recover known biology by nominating biologically meaningful pathways in complex diseases such as Alzheimer's disease.

Multi-Omics Model Applied to Cancer Genetics

In this review, we focus on bioinformatic oncology as an integrative discipline that incorporates knowledge from the mathematical, physical, and computational fields to further the biomedical understanding of cancer. Before providing a deeper insight into the bioinformatics approach and utilities involved in oncology, we must understand what is a system biology framework and the genetic connection, because of the high heterogenicity of the backgrounds of people approaching precision medicine. In fact, it is essential to providing general theoretical information on genomics, epigenomics, and transcriptomics to understand the phases of multi-omics approach. We consider how to create a multi-omics model. In the last section, we describe the new frontiers and future perspectives of this field.

leapR: An R Package for Multiomic Pathway Analysis

A generalized goal of many high-throughput data studies is to identify functional mechanisms that underlie observed biological phenomena, whether they be disease outcomes or metabolic output. Increasingly, studies that rely on multiple sources of high-throughput data (genomic, transcriptomic, proteomic, metabolomic) are faced with a challenge of summarizing the data to generate testable hypotheses. However, this requires a time-consuming process to evaluate numerous statistical methods across numerous data sources. Here, we introduce the leapR package, a framework to rapidly assess biological pathway activity using diverse statistical tests and data sources, allowing facile integration of multisource data. The leapR package with a user manual and example workflow is available for download from GitHub (https://github.com/biodataganache/leapR).

COCOA: coordinate covariation analysis of epigenetic heterogeneity

A key challenge in epigenetics is to determine the biological significance of epigenetic variation among individuals. We present Coordinate Covariation Analysis (COCOA), a computational framework that uses covariation of epigenetic signals across individuals and a database of region sets to annotate epigenetic heterogeneity. COCOA is the first such tool for DNA methylation data and can also analyze any epigenetic signal with genomic coordinates. We demonstrate COCOA’s utility by analyzing DNA methylation, ATAC-seq, and multi-omic data in supervised and unsupervised analyses, showing that COCOA provides new understanding of inter-sample epigenetic variation. COCOA is available on Bioconductor (http://bioconductor.org/packages/COCOA).