Gene Set Analysis: Challenges, Opportunities, and Future Research

DeepGSEA: explainable deep gene set enrichment analysis for single-cell transcriptomic data

MOTIVATION

Gene set enrichment (GSE) analysis allows for an interpretation of gene expression through pre-defined gene set databases and is a critical step in understanding different phenotypes. With the rapid development of single-cell RNA sequencing (scRNA-seq) technology, GSE analysis can be performed on fine-grained gene expression data to gain a nuanced understanding of phenotypes of interest. However, with the cellular heterogeneity in single-cell gene profiles, current statistical GSE analysis methods sometimes fail to identify enriched gene sets. Meanwhile, deep learning has gained traction in applications like clustering and trajectory inference in single-cell studies due to its prowess in capturing complex data patterns. However, its use in GSE analysis remains limited, due to interpretability challenges.

RESULTS

In this paper, we present DeepGSEA, an explainable deep gene set enrichment analysis approach which leverages the expressiveness of interpretable, prototype-based neural networks to provide an in-depth analysis of GSE. DeepGSEA learns the ability to capture GSE information through our designed classification tasks, and significance tests can be performed on each gene set, enabling the identification of enriched sets. The underlying distribution of a gene set learned by DeepGSEA can be explicitly visualized using the encoded cell and cellular prototype embeddings. We demonstrate the performance of DeepGSEA over commonly used GSE analysis methods by examining their sensitivity and specificity with four simulation studies. In addition, we test our model on three real scRNA-seq datasets and illustrate the interpretability of DeepGSEA by showing how its results can be explained.

AVAILABILITY AND IMPLEMENTATION

https://github.com/Teddy-XiongGZ/DeepGSEA.

A genome scale transcriptional regulatory model of the human placenta

Gene regulation is essential to placental function and fetal development. We built a genome-scale transcriptional regulatory network (TRN) of the human placenta using digital genomic footprinting and transcriptomic data. We integrated 475 transcriptomes and 12 DNase hypersensitivity datasets from placental samples to globally and quantitatively map transcription factor (TF)-target gene interactions. In an independent dataset, the TRN model predicted target gene expression with an out-of-sample R2 greater than 0.25 for 73% of target genes. We performed siRNA knockdowns of four TFs and achieved concordance between the predicted gene targets in our TRN and differences in expression of knockdowns with an accuracy of >0.7 for three of the four TFs. Our final model contained 113,158 interactions across 391 TFs and 7712 target genes and is publicly available. We identified 29 TFs which were significantly enriched as regulators for genes previously associated with preterm birth, and eight of these TFs were decreased in preterm placentas.

Elucidating immune-related gene transcriptional programs via factorization of large-scale RNA-profiles

Recent developments in immunotherapy, including immune checkpoint blockade (ICB) and adoptive cell therapy (ACT), have encountered challenges such as immune-related adverse events and resistance, especially in solid tumors. To advance the field, a deeper understanding of the molecular mechanisms behind treatment responses and resistance is essential. However, the lack of functionally characterized immune-related gene sets has limited data-driven immunological research. To address this gap, we adopted non-negative matrix factorization on 83 human bulk RNA sequencing (RNA-seq) datasets and constructed 28 immune-specific gene sets. After rigorous immunologist-led manual annotations and orthogonal validations across immunological contexts and functional omics data, we demonstrated that these gene sets can be applied to refine pan-cancer immune subtypes, improve ICB response prediction and functionally annotate spatial transcriptomic data. These functional gene sets, informing diverse immune states, will advance our understanding of immunology and cancer research.

GOAT: efficient and robust identification of gene set enrichment

Gene set enrichment analysis is foundational to the interpretation of high throughput biology. Identifying enriched Gene Ontology (GO) terms or disease-associated gene sets within a list of gene effect sizes that represent experimental outcomes is an everyday task in life science that crucially depends on robust and sensitive statistical tools. We here present GOAT, a parameter-free algorithm for gene set enrichment analysis of preranked gene lists. The algorithm can precompute null distributions from standardized gene scores, enabling enrichment testing of the GO database in one second. Validations using synthetic data show that estimated gene set p-values are well calibrated under the null hypothesis and invariant to gene list length and gene set size. Application to various real-world proteomics and gene expression studies demonstrates that GOAT identifies more significant GO terms as compared to current methods. GOAT is freely available as an R package and user-friendly online tool for gene set enrichment analyses that includes interactive data visualizations: https://ftwkoopmans.github.io/goat .

Tumour microenvironment characterisation to stratify patients for hyperthermic intraperitoneal chemotherapy in high-grade serous ovarian cancer (OVHIPEC-1)

BACKGROUND

Hyperthermic intraperitoneal chemotherapy (HIPEC) improves survival in patients with Stage III ovarian cancer following interval cytoreductive surgery (CRS). Optimising patient selection is essential to maximise treatment efficacy and avoid overtreatment. This study aimed to identify biomarkers that predict HIPEC benefit by analysing gene signatures and cellular composition of tumours from participants in the OVHIPEC-1 trial.

METHODS

Whole-transcriptome RNA sequencing data were retrieved from high-grade serous ovarian cancer (HGSOC) samples from 147 patients obtained during interval CRS. We performed differential gene expression analysis and applied deconvolution methods to estimate cell-type proportions in bulk mRNA data, validated by histological assessment. We tested the interaction between treatment and potential predictors on progression-free survival using Cox proportional hazards models.

RESULTS

While differential gene expression analysis did not yield any predictive biomarkers, the cellular composition, as characterised by deconvolution, indicated that the absence of macrophages and the presence of B cells in the tumour microenvironment are potential predictors of HIPEC benefit. The histological assessment confirmed the predictive value of macrophage absence.

CONCLUSION

Immune cell composition, in particular macrophages absence, may predict response to HIPEC in HGSOC and these hypothesis-generating findings warrant further investigation.

CLINICAL TRIAL REGISTRATION

NCT00426257.

Molecular characterisation of lupus low disease activity state (LLDAS) and DORIS remission by whole-blood transcriptome-based pathways in a pan-European systemic lupus erythematosus cohort

OBJECTIVES

To unveil biological milieus underlying low disease activity (LDA) and remission versus active systemic lupus erythematosus (SLE).

METHODS

We determined differentially expressed pathways (DEPs) in SLE patients from the PRECISESADS project (NTC02890121) stratified into patients fulfilling and not fulfilling the criteria of (1) Lupus LDA State (LLDAS), (2) Definitions of Remission in SLE remission, and (3) LLDAS exclusive of remission.

RESULTS

We analysed data from 321 patients; 40.8% were in LLDAS, and 17.4% in DORIS remission. After exclusion of patients in remission, 28.3% were in LLDAS. Overall, 604 pathways differed significantly in LLDAS versus non-LLDAS patients with an false-discovery rate-corrected p (q)<0.05 and a robust effect size (dr)≥0.36. Accordingly, 288 pathways differed significantly between DORIS remitters and non-remitters (q<0.05 and dr≥0.36). DEPs yielded distinct molecular clusters characterised by differential serological, musculoskeletal, and renal activity. Analysis of partially overlapping samples showed no DEPs between LLDAS and DORIS remission. Drug repurposing potentiality for treating SLE was unveiled, as were important pathways underlying active SLE whose modulation could aid attainment of LLDAS/remission, including toll-like receptor (TLR) cascades, Bruton tyrosine kinase (BTK) activity, the cytotoxic T lymphocyte antigen 4 (CTLA-4)-related inhibitory signalling, and the nucleotide-binding oligomerization domain leucine-rich repeat-containing protein 3 (NLRP3) inflammasome pathway.

CONCLUSIONS

We demonstrated for the first time molecular signalling pathways distinguishing LLDAS/remission from active SLE. LLDAS/remission was associated with reversal of biological processes related to SLE pathogenesis and specific clinical manifestations. DEP clustering by remission better grouped patients compared with LLDAS, substantiating remission as the ultimate treatment goal in SLE; however, the lack of substantial pathway differentiation between the two states justifies LLDAS as an acceptable goal from a biological perspective.

Bridging systems biology and tissue engineering: Unleashing the full potential of complex 3D in vitro tissue models of disease

Rapid advances in tissue engineering have resulted in more complex and physiologically relevant 3D in vitro tissue models with applications in fundamental biology and therapeutic development. However, the complexity provided by these models is often not leveraged fully due to the reductionist methods used to analyze them. Computational and mathematical models developed in the field of systems biology can address this issue. Yet, traditional systems biology has been mostly applied to simpler in vitro models with little physiological relevance and limited cellular complexity. Therefore, integrating these two inherently interdisciplinary fields can result in new insights and move both disciplines forward. In this review, we provide a systematic overview of how systems biology has been integrated with 3D in vitro tissue models and discuss key application areas where the synergies between both fields have led to important advances with potential translational impact. We then outline key directions for future research and discuss a framework for further integration between fields.

Recreating metabolic interactions of the tumour microenvironment

Tumours are heterogeneous tissues containing diverse populations of cells and an abundant extracellular matrix (ECM). This tumour microenvironment prompts cancer cells to adapt their metabolism to survive and grow. Besides epigenetic factors, the metabolism of cancer cells is shaped by crosstalk with stromal cells and extracellular components. To date, most experimental models neglect the complexity of the tumour microenvironment and its relevance in regulating the dynamics of the metabolism in cancer. We discuss emerging strategies to model cellular and extracellular aspects of cancer metabolism. We highlight cancer models based on bioengineering, animal, and mathematical approaches to recreate cell-cell and cell-matrix interactions and patient-specific metabolism. Combining these approaches will improve our understanding of cancer metabolism and support the development of metabolism-targeting therapies.

Unraveling the molecular relevance of brain phenotypes: A comparative analysis of null models and test statistics

Correlating transcriptional profiles with imaging-derived phenotypes has the potential to reveal possible molecular architectures associated with cognitive functions, brain development and disorders. Competitive null models built by resampling genes and self-contained null models built by spinning brain regions, along with varying test statistics, have been used to determine the significance of transcriptional associations. However, there has been no systematic evaluation of their performance in imaging transcriptomics analyses. Here, we evaluated the performance of eight different test statistics (mean, mean absolute value, mean squared value, max mean, median, Kolmogorov-Smirnov (KS), Weighted KS and the number of significant correlations) in both competitive null models and self-contained null models. Simulated brain maps (n = 1,000) and gene sets (n = 500) were used to calculate the probability of significance (Psig) for each statistical test. Our results suggested that competitive null models may result in false positive results driven by co-expression within gene sets. Furthermore, we demonstrated that the self-contained null models may fail to account for distribution characteristics (e.g., bimodality) of correlations between all available genes and brain phenotypes, leading to false positives. These two confounding factors interacted differently with test statistics, resulting in varying outcomes. Specifically, the sign-sensitive test statistics (i.e., mean, median, KS, Weighted KS) were influenced by co-expression bias in the competitive null models, while median and sign-insensitive test statistics were sensitive to the bimodality bias in the self-contained null models. Additionally, KS-based statistics produced conservative results in the self-contained null models, which increased the risk of false negatives. Comprehensive supplementary analyses with various configurations, including realistic scenarios, supported the results. These findings suggest utilizing sign-insensitive test statistics such as mean absolute value, max mean in the competitive null models and the mean as the test statistic for the self-contained null models. Additionally, adopting the confounder-matched (e.g., coexpression-matched) null models as an alternative to standard null models can be a viable strategy. Overall, the present study offers insights into the selection of statistical tests for imaging transcriptomics studies, highlighting areas for further investigation and refinement in the evaluation of novel and commonly used tests.

Network enrichment significance testing in brain-phenotype association studies

Functional networks often guide our interpretation of spatial maps of brain-phenotype associations. However, methods for assessing enrichment of associations within networks of interest have varied in terms of both scientific rigor and underlying assumptions. While some approaches have relied on subjective interpretations, others have made unrealistic assumptions about spatial properties of imaging data, leading to inflated false positive rates. We seek to address this gap in existing methodology by borrowing insight from a method widely used in genetics research for testing enrichment of associations between a set of genes and a phenotype of interest. We propose network enrichment significance testing (NEST), a flexible framework for testing the specificity of brain-phenotype associations to functional networks or other sub-regions of the brain. We apply NEST to study enrichment of associations with structural and functional brain imaging data from a large-scale neurodevelopmental cohort study.

Assessment of Gene Set Enrichment Analysis using curated RNA-seq-based benchmarks

Pathway enrichment analysis is a ubiquitous computational biology method to interpret a list of genes (typically derived from the association of large-scale omics data with phenotypes of interest) in terms of higher-level, predefined gene sets that share biological function, chromosomal location, or other common features. Among many tools developed so far, Gene Set Enrichment Analysis (GSEA) stands out as one of the pioneering and most widely used methods. Although originally developed for microarray data, GSEA is nowadays extensively utilized for RNA-seq data analysis. Here, we quantitatively assessed the performance of a variety of GSEA modalities and provide guidance in the practical use of GSEA in RNA-seq experiments. We leveraged harmonized RNA-seq datasets available from The Cancer Genome Atlas (TCGA) in combination with large, curated pathway collections from the Molecular Signatures Database to obtain cancer-type-specific target pathway lists across multiple cancer types. We carried out a detailed analysis of GSEA performance using both gene-set and phenotype permutations combined with four different choices for the Kolmogorov-Smirnov enrichment statistic. Based on our benchmarks, we conclude that the classic/unweighted gene-set permutation approach offered comparable or better sensitivity-vs-specificity tradeoffs across cancer types compared with other, more complex and computationally intensive permutation methods. Finally, we analyzed other large cohorts for thyroid cancer and hepatocellular carcinoma. We utilized a new consensus metric, the Enrichment Evidence Score (EES), which showed a remarkable agreement between pathways identified in TCGA and those from other sources, despite differences in cancer etiology. This finding suggests an EES-based strategy to identify a core set of pathways that may be complemented by an expanded set of pathways for downstream exploratory analysis. This work fills the existing gap in current guidelines and benchmarks for the use of GSEA with RNA-seq data and provides a framework to enable detailed benchmarking of other RNA-seq-based pathway analysis tools.

Elucidating immune-related gene transcriptional programs via factorization of large-scale RNA-profiles

Recent developments in immunotherapy, including immune checkpoint blockade (ICB) and adoptive cell therapy, have encountered challenges such as immune-related adverse events and resistance, especially in solid tumors. To advance the field, a deeper understanding of the molecular mechanisms behind treatment responses and resistance is essential. However, the lack of functionally characterized immune-related gene sets has limited data-driven immunological research. To address this gap, we adopted non-negative matrix factorization on 83 human bulk RNA-seq datasets and constructed 28 immune-specific gene sets. After rigorous immunologist-led manual annotations and orthogonal validations across immunological contexts and functional omics data, we demonstrated that these gene sets can be applied to refine pan-cancer immune subtypes, improve ICB response prediction and functionally annotate spatial transcriptomic data. These functional gene sets, informing diverse immune states, will advance our understanding of immunology and cancer research.

Integrative bioinformatics and experimental validation of hub genetic markers in acne vulgaris: Toward personalized diagnostic and therapeutic strategies

BACKGROUND

Acne vulgaris is a widespread chronic inflammatory dermatological condition. The precise molecular and genetic mechanisms of its pathogenesis remain incompletely understood. This research synthesizes existing databases, targeting a comprehensive exploration of core genetic markers.

METHODS

Gene expression datasets (GSE6475, GSE108110, and GSE53795) were retrieved from the GEO. Differentially expressed genes (DEGs) were identified using the limma package. Enrichment analyses were conducted using GSVA for pathway assessment and clusterProfiler for GO and KEGG analyses. PPI networks and immune cell infiltration were analyzed using the STRING database and ssGSEA, respectively. We investigated the correlation between hub gene biomarkers and immune cell infiltration using Spearman's rank analysis. ROC curve analysis validated the hub genes' diagnostic accuracy. miRNet, TarBase v8.0, and ChEA3 identified miRNA/transcription factor-gene interactions, while DrugBank delineated drug-gene interactions. Experiments utilized HaCaT cells stimulated with Propionibacterium acnes, treated with retinoic acid and methotrexate, and evaluated using RT-qPCR, ELISA, western blot, lentiviral transduction, CCK-8, wound-healing, and transwell assays.

RESULTS

There were 104 genes with consistent differences across the three datasets of paired acne and normal skin. Functional analyses emphasized the significant enrichment of these DEGs in immune-related pathways. PPI network analysis pinpointed hub genes PTPRC, CXCL8, ITGB2, and MMP9 as central players in acne pathogenesis. Elevated levels of specific immune cell infiltration in acne lesions corroborated the inflammatory nature of the disease. ROC curve analysis identified the acne diagnostic potential of four hub genes. Key miRNAs, particularly hsa-mir-124-3p, and central transcription factors like TFEC were noted as significant regulators. In vitro validation using HaCaT cells confirmed the upregulation of hub genes following Propionibacterium acnes exposure, while CXCL8 knockdown reduced pro-inflammatory cytokines, cell proliferation, and migration. DrugBank insights led to the exploration of retinoic acid and methotrexate, both of which mitigated gene expression upsurge and inflammatory mediator secretion.

CONCLUSION

This comprehensive study elucidated pivotal genes associated with acne pathogenesis, notably PTPRC, CXCL8, ITGB2, and MMP9. The findings underscore potential biomarkers, therapeutic targets, and the therapeutic potential of agents like retinoic acid and methotrexate. The congruence between bioinformatics and experimental validations suggests promising avenues for personalized acne treatments.

Meta-analysis of differential gene expression in lower motor neurons isolated by laser capture microdissection from post-mortem ALS spinal cords

Introduction

ALS is a fatal neurodegenerative disease for which underlying mechanisms are incompletely understood. The motor neuron is a central player in ALS pathogenesis but different transcriptome signatures have been derived from bulk analysis of post-mortem tissue and iPSC-derived motor neurons (iPSC-MNs).

Methods

This study performed a meta-analysis of six gene expression studies (microarray and RNA-seq) in which laser capture microdissection (LCM) was used to isolate lower motor neurons from post-mortem spinal cords of ALS and control (CTL) subjects. Differentially expressed genes (DEGs) with consistent ALS versus CTL expression differences across studies were identified.

Results

The analysis identified 222 ALS-increased DEGs (FDR <0.10, SMD >0.80) and 278 ALS-decreased DEGs (FDR <0.10, SMD < -0.80). ALS-increased DEGs were linked to PI3K-AKT signaling, innate immunity, inflammation, motor neuron differentiation and extracellular matrix. ALS-decreased DEGs were associated with the ubiquitin-proteosome system, microtubules, axon growth, RNA-binding proteins and synaptic membrane. ALS-decreased DEG mRNAs frequently interacted with RNA-binding proteins (e.g., FUS, HuR). The complete set of DEGs (increased and decreased) overlapped significantly with genes near ALS-associated SNP loci (p < 0.01). Transcription factor target motifs with increased proximity to ALS-increased DEGs were identified, most notably DNA elements predicted to interact with forkhead transcription factors (e.g., FOXP1) and motor neuron and pancreas homeobox 1 (MNX1). Some of these DNA elements overlie ALS-associated SNPs within known enhancers and are predicted to have genotype-dependent MNX1 interactions. DEGs were compared to those identified from SOD1-G93A mice and bulk spinal cord segments or iPSC-MNs from ALS patients. There was good correspondence with transcriptome changes from SOD1-G93A mice (r ≤ 0.408) but most DEGs were not differentially expressed in bulk spinal cords or iPSC-MNs and transcriptome-wide effect size correlations were weak (bulk tissue: r ≤ 0.207, iPSC-MN: r ≤ 0.037).

Conclusion

This study defines a robust transcriptome signature from LCM-based motor neuron studies of post-mortem tissue from ALS and CTL subjects. This signature differs from those obtained from analysis of bulk spinal cord segments and iPSC-MNs. Results provide insight into mechanisms underlying gene dysregulation in ALS and highlight connections between these mechanisms, ALS genetics, and motor neuron biology.

Finding new analgesics: Computational pharmacology faces drug discovery challenges

Despite the worldwide prevalence and huge burden of pain, pain is an undertreated phenomenon. Currently used analgesics have several limitations regarding their efficacy and safety. The discovery of analgesics possessing a novel mechanism of action has faced multiple challenges, including a limited understanding of biological processes underpinning pain and analgesia and poor animal-to-human translation. Computational pharmacology is currently employed to face these challenges. In this review, we discuss the theory, methods, and applications of computational pharmacology in pain research. Computational pharmacology encompasses a wide variety of theoretical concepts and practical methodological approaches, with the overall aim of gaining biological insight through data acquisition and analysis. Data are acquired from patients or animal models with pain or analgesic treatment, at different levels of biological organization (molecular, cellular, physiological, and behavioral). Distinct methodological algorithms can then be used to analyze and integrate data. This helps to facilitate the identification of biological molecules and processes associated with pain phenotype, build quantitative models of pain signaling, and extract translatable features between humans and animals. However, computational pharmacology has several limitations, and its predictions can provide false positive and negative findings. Therefore, computational predictions are required to be validated experimentally before drawing solid conclusions. In this review, we discuss several case study examples of combining and integrating computational tools with experimental pain research tools to meet drug discovery challenges.

RNA-sequencing of Human Kidney Allografts and Delineation of T-Cell Genes, Gene Sets, and Pathways Associated With Acute T Cell-mediated Rejection

BACKGROUND

Delineation of T-cell genes, gene sets, pathways, and T-cell subtypes associated with acute T cell-mediated rejection (TCMR) may improve its management.

METHODS

We performed bulk RNA-sequencing of 34 kidney allograft biopsies (16 Banff TCMR and 18 no rejection [NR] biopsies) from 34 adult recipients of human kidneys. Computational analysis was performed to determine the differential intragraft expression of T-cell genes at the level of single-gene, gene set, and pathways.

RESULTS

T-cell signaling pathway gene sets for plenary T-cell activation were overrepresented in TCMR biopsies compared with NR biopsies. Heightened expression of T-cell signaling genes was validated using external TCMR biopsies. Pro- and anti-inflammatory immune gene sets were enriched, and metabolism gene sets were depleted in TCMR biopsies compared with NR biopsies. Gene signatures of regulatory T cells, Th1 cells, Th2 cells, Th17 cells, T follicular helper cells, CD4 tissue-resident memory T cells, and CD8 tissue-resident memory T cells were enriched in TCMR biopsies compared with NR biopsies. T-cell exhaustion and anergy were also molecular attributes of TCMR. Gene sets associated with antigen processing and presentation, and leukocyte transendothelial migration were overexpressed in TCMR biopsies compared with NR biopsies. Cellular deconvolution of graft infiltrating cells by gene expression patterns identified CD8 T cell to be the most abundant T-cell subtype infiltrating the allograft during TCMR.

CONCLUSIONS

Our delineation of intragraft T-cell gene expression patterns, in addition to yielding new biological insights, may help prioritize T-cell genes and T-cell subtypes for therapeutic targeting.

Comorbidity Genes of Alzheimer's Disease and Type 2 Diabetes Associated with Memory and Cognitive Function

Alzheimer's disease (AD) and type 2 diabetes mellitus (T2DM) are comorbidities that result from the sharing of common genes. The molecular background of comorbidities can provide clues for the development of treatment and management strategies. Here, the common genes involved in the development of the two diseases and in memory and cognitive function are reviewed. Network clustering based on protein-protein interaction network identified tightly connected gene clusters that have an impact on memory and cognition among the comorbidity genes of AD and T2DM. Genes with functional implications were intensively reviewed and relevant evidence summarized. Gene information will be useful in the discovery of biomarkers and the identification of tentative therapeutic targets for AD and T2DM.

Gene set correlation enrichment analysis for interpreting and annotating gene expression profiles

Pathway analysis, including nontopology-based (non-TB) and topology-based (TB) methods, is widely used to interpret the biological phenomena underlying differences in expression data between two phenotypes. By considering dependencies and interactions between genes, TB methods usually perform better than non-TB methods in identifying pathways that include closely relevant or directly causative genes for a given phenotype. However, most TB methods may be limited by incomplete pathway data used as the reference network or by difficulties in selecting appropriate reference networks for different research topics. Here, we propose a gene set correlation enrichment analysis method, Gscore, based on an expression dataset-derived coexpression network to examine whether a differentially expressed gene (DEG) list (or each of its DEGs) is associated with a known gene set. Gscore is better able to identify target pathways in 89 human disease expression datasets than eight other state-of-the-art methods and offers insight into how disease-wide and pathway-wide associations reflect clinical outcomes. When applied to RNA-seq data from COVID-19-related cells and patient samples, Gscore provided a means for studying how DEGs are implicated in COVID-19-related pathways. In summary, Gscore offers a powerful analytical approach for annotating individual DEGs, DEG lists, and genome-wide expression profiles based on existing biological knowledge.

Robustness evaluations of pathway activity inference methods on gene expression data

BACKGROUND

With the exponential growth of high-throughput technologies, multiple pathway analysis methods have been proposed to estimate pathway activities from gene expression profiles. These pathway activity inference methods can be divided into two main categories: non-Topology-Based (non-TB) and Pathway Topology-Based (PTB) methods. Although some review and survey articles discussed the topic from different aspects, there is a lack of systematic assessment and comparisons on the robustness of these approaches.

RESULTS

Thus, this study presents comprehensive robustness evaluations of seven widely used pathway activity inference methods using six cancer datasets based on two assessments. The first assessment seeks to investigate the robustness of pathway activity in pathway activity inference methods, while the second assessment aims to assess the robustness of risk-active pathways and genes predicted by these methods. The mean reproducibility power and total number of identified informative pathways and genes were evaluated. Based on the first assessment, the mean reproducibility power of pathway activity inference methods generally decreased as the number of pathway selections increased. Entropy-based Directed Random Walk (e-DRW) distinctly outperformed other methods in exhibiting the greatest reproducibility power across all cancer datasets. On the other hand, the second assessment shows that no methods provide satisfactory results across datasets.

CONCLUSION

However, PTB methods generally appear to perform better in producing greater reproducibility power and identifying potential cancer markers compared to non-TB methods.

Insights from metabolomics and transcriptomics studies on Perfluorohexanesulfonic acid (PFHxS) exposed zebrafish embryos

Perfluorohexanesulfonic acid (PFHxS) is a short-chain perfluoroalkyl substance widely used to replace the banned perfluorooctanesulfonic acid (PFOS) in various industrial and household products. It can be found in the environment and human bodies; however, its potential toxicities are not well studied. Zebrafish have been extensively used as a model for studying toxicants, and currently, two studies have reported on the toxicity of PFHxS in zebrafish from different approaches. Ulhaq and Tse (J Hazard Mater. 2023; 457: 131722) conducted general biological experiments and applied transcriptomics to demonstrate that PFHxS at a concentration of 5 μM could affect glucose and fatty acid metabolism, leading to oxidative stress, developmental defects, and cell cycle arrest. Xu et al. (Sci Total Environ. 2023; 887: 163770) employed metabolomics and showed that concentrations of various metabolites changed after exposure to 3 and 10 μM PFHxS. As we observed a match between the metabolomics data and our biochemistry experimental findings, we integrated the two studies, which enabled us to unfold the possible mechanism of the deregulated metabolites. We identified 22 differential expressed genes (DEGs) in the tricarboxylic acid (TCA) cycle, 17 DEGs in glcyolytic process, including the critical glucokinase under the carbon metabolism. Besides, genes likes aldehyde dehydrogenases, and histone-lysine N-methyltransferases that participate in lipid peroxidation and amino metabolism respectively were spotted. Lastly, we further strengthen our discoveries by undergoing the gene set enrichment analysis. This article could provide insights into the toxicity of PFHxS, as well as prospects for environmental studies.

Deciphering the genetic landscape of obesity: a data-driven approach to identifying plausible causal genes and therapeutic targets

OBJECTIVES

Genome-wide association studies (GWAS) have successfully revealed numerous susceptibility loci for obesity. However, identifying the causal genes, pathways, and tissues/cell types responsible for these associations remains a challenge, and standardized analysis workflows are lacking. Additionally, due to limited treatment options for obesity, there is a need for the development of new pharmacological therapies. This study aimed to address these issues by performing step-wise utilization of knowledgebase for gene prioritization and assessing the potential relevance of key obesity genes as therapeutic targets.

METHODS AND RESULTS

First, we generated a list of 28,787 obesity-associated SNPs from the publicly available GWAS dataset (approximately 800,000 individuals in the GIANT meta-analysis). Then, we prioritized 1372 genes with significant in silico evidence against genomic and transcriptomic data, including transcriptionally regulated genes in the brain from transcriptome-wide association studies. In further narrowing down the gene list, we selected key genes, which we found to be useful for the discovery of potential drug seeds as demonstrated in lipid GWAS separately. We thus identified 74 key genes for obesity, which are highly interconnected and enriched in several biological processes that contribute to obesity, including energy expenditure and homeostasis. Of 74 key genes, 37 had not been reported for the pathophysiology of obesity. Finally, by drug-gene interaction analysis, we detected 23 (of 74) key genes that are potential targets for 78 approved and marketed drugs.

CONCLUSIONS

Our results provide valuable insights into new treatment options for obesity through a data-driven approach that integrates multiple up-to-date knowledgebases.

Reactome pathway analysis from whole-blood transcriptome reveals unique characteristics of systemic sclerosis patients at the preclinical stage

Objective

This study aims to characterize differential expressed pathways (DEP) in subjects with preclinical systemic sclerosis (PreSSc) characterized uniquely by Raynaud phenomenon, specific autoantibodies, and/or capillaroscopy positive for scleroderma pattern.

Methods

Whole-blood samples from 33 PreSSc with clinical prospective data (baseline and after 4 years of follow-up) and 16 matched healthy controls (HC) were analyzed for global gene expression transcriptome analysis via RNA sequencing. Functional Analysis of Individual Microarray Expression method annotated Reactome individualized pathways. ANOVA analysis identified DEP whose predictive capability were tested in logistic regression models after extensive internal validation.

Results

At 4 years, 42.4% subjects progressed (evolving PreSSc), while the others kept stable PreSSc clinical features (stable PreSSc). At baseline, out of 831 pathways, 541 DEP were significant at a false discovery rate <0.05, differentiating PreSSc versus HC with an AUROC = 0.792 ± 0.242 in regression models. Four clinical groups were identified via unsupervised clustering (HC, HC and PreSSc with HC-like features, PreSSc and HC with PreSSc-like features, and PreSSc). Biological signatures changed with disease progression while remaining unchanged in stable subjects. The magnitude of change was related to the baseline cluster, yet no DEP at baseline was predictive of progression. Disease progression was mostly related to changes in signal transduction pathways especially linked to calcium-related events and inositol 1,4,5-triphosphate metabolism.

Conclusion

PreSSc had distinguished Reactome pathway signatures compared to HC. Progression to definite SSc was characterized by a shift in biological fingertips. Calcium-related events promoting endothelial damage and vasculopathy may be relevant to disease progression.

Addressing erroneous scale assumptions in microbe and gene set enrichment analysis

By applying Differential Set Analysis (DSA) to sequence count data, researchers can determine whether groups of microbes or genes are differentially enriched. Yet sequence count data suffer from a scale limitation: these data lack information about the scale (i.e., size) of the biological system under study, leading some authors to call these data compositional (i.e., proportional). In this article, we show that commonly used DSA methods that rely on normalization make strong, implicit assumptions about the unmeasured system scale. We show that even small errors in these scale assumptions can lead to positive predictive values as low as 9%. To address this problem, we take three novel approaches. First, we introduce a sensitivity analysis framework to identify when modeling results are robust to such errors and when they are suspect. Unlike standard benchmarking studies, this framework does not require ground-truth knowledge and can therefore be applied to both simulated and real data. Second, we introduce a statistical test that provably controls Type-I error at a nominal rate despite errors in scale assumptions. Finally, we discuss how the impact of scale limitations depends on a researcher's scientific goals and provide tools that researchers can use to evaluate whether their goals are at risk from erroneous scale assumptions. Overall, the goal of this article is to catalyze future research into the impact of scale limitations in analyses of sequence count data; to illustrate that scale limitations can lead to inferential errors in practice; yet to also show that rigorous and reproducible scale reliant inference is possible if done carefully.

Challenges and opportunities in sharing microbiome data and analyses

Microbiome data, metadata and analytical workflows have become 'big' in terms of volume and complexity. Although the infrastructure and technologies to share data have been established, the interdisciplinary and multi-omic nature of the field can make resources difficult to identify and use. Following best practices for data deposition requires substantial effort, with sometimes little obvious reward. Gaps remain where microbiome-specific resources for data sharing or reproducibility do not yet exist. We outline available best practices, challenges to their adoption and opportunities in data sharing in microbiome research. We showcase examples of best practices and advocate for their enforcement and incentivization for data sharing. This includes recognition of data curation and sharing endeavours by individuals, institutions, journals and funders. Opportunities for progress include enabling microbiome-specific databases to incorporate future methods for data analysis, integration and reuse.

An interpretable machine learning pipeline based on transcriptomics predicts phenotypes of lupus patients

Machine learning (ML) has the potential to identify subsets of patients with distinct phenotypes from gene expression data. However, phenotype prediction using ML has often relied on identifying important genes without a systems biology context. To address this, we created an interpretable ML approach based on blood transcriptomics to predict phenotype in systemic lupus erythematosus (SLE), a heterogeneous autoimmune disease. We employed a sequential grouped feature importance algorithm to assess the performance of gene sets, including immune and metabolic pathways and cell types, known to be abnormal in SLE in predicting disease activity and organ involvement. Gene sets related to interferon, tumor necrosis factor, the mitoribosome, and T cell activation were the best predictors of phenotype with excellent performance. These results suggest potential relationships between the molecular pathways identified in each model and manifestations of SLE. This ML approach to phenotype prediction can be applied to other diseases and tissues.

Single-cell gene set scoring with nearest neighbor graph smoothed data (gssnng)

Summary

Gene set scoring (or enrichment) is a common dimension reduction task in bioinformatics that can be focused on the differences between groups or at the single sample level. Gene sets can represent biological functions, molecular pathways, cell identities, and more. Gene set scores are context dependent values that are useful for interpreting biological changes following experiments or perturbations. Single sample scoring produces a set of scores, one for each member of a group, which can be analyzed with statistical models that can include additional clinically important factors such as gender or age. However, the sparsity and technical noise of single-cell expression measures create difficulties for these methods, which were originally designed for bulk expression profiling (microarrays, RNAseq). This can be greatly remedied by first applying a smoothing transformation that shares gene measure information within transcriptomic neighborhoods. In this work, we use the nearest neighbor graph of cells for matrix smoothing to produce high quality gene set scores on a per-cell, per-group, level which is useful for visualization and statistical analysis.

Availability and implementation

The gssnng software is available using the python package index (PyPI) and works with Scanpy AnnData objects. It can be installed using "pip install gssnng." More information and demo notebooks: see https://github.com/IlyaLab/gssnng.

Delving into gene-set multiplex networks facilitated by a k-nearest neighbor-based measure of similarity

Gene sets are functional units for living cells. Previously, limited studies investigated the complex relations among gene sets, but documents about their altering patterns across biological conditions still need to be prepared. In this study, we adopted and modified a classical k-nearest neighbor-based association function to detect inter-gene-set similarities. Based on this method, we built multiplex networks of gene sets for the first time; these networks contain layers of gene sets corresponding to different populations of cells. The context-based multiplex networks can capture meaningful biological variation and have considerable differences from knowledge-based networks of gene sets built on Jaccard similarity, as demonstrated in this study. Furthermore, at the scale of individual gene sets, the structural coefficients of gene sets (multiplex PageRank centrality, clustering coefficient, and participation coefficient) disclose the diversity of gene sets from the perspective of structural properties and make it easier to identify unique gene sets. In gene set enrichment analysis (GSEA), each gene set is treated independently, and its contextual and relational attributes are ignored. The structural coefficients of gene sets can supplement GSEA with information about the overall picture of gene sets, promoting the constructive reorganization of the enriched terms and helping researchers better prioritize and select gene sets.

Essential procedures of single-cell RNA sequencing in multiple myeloma and its translational value

Multiple myeloma (MM) is a malignant neoplasm characterized by clonal proliferation of abnormal plasma cells. In many countries, it ranks as the second most prevalent malignant neoplasm of the hematopoietic system. Although treatment methods for MM have been continuously improved and the survival of patients has been dramatically prolonged, MM remains an incurable disease with a high probability of recurrence. As such, there are still many challenges to be addressed. One promising approach is single-cell RNA sequencing (scRNA-seq), which can elucidate the transcriptome heterogeneity of individual cells and reveal previously unknown cell types or states in complex tissues. In this review, we outlined the experimental workflow of scRNA-seq in MM, listed some commonly used scRNA-seq platforms and analytical tools. In addition, with the advent of scRNA-seq, many studies have made new progress in the key molecular mechanisms during MM clonal evolution, cell interactions and molecular regulation in the microenvironment, and drug resistance mechanisms in target therapy. We summarized the main findings and sequencing platforms for applying scRNA-seq to MM research and proposed broad directions for targeted therapies based on these findings.

Ferroptosis regulator NOS2 is closely associated with the prognosis and cell malignant behaviors of hepatoblastoma: a bioinformatic and in vitro study

Background

Hepatoblastoma (HB) is the most common liver tumor in children with easy metastasis. The emergence of ferroptosis as a novel form of cell death has gained increased attention in various human cancers. However, the roles of ferroptosis-related (FR) genes in HB remain elusive.

Methods

The GSE133039, GSE131329, and GSE81928 datasets were utilized for screening core FR genes in HB. Through Lasso regression analysis and using the support vector machine recursive feature elimination (SVM-RFE) algorithm, three candidate FR genes were obtained for characterizing HB. Their expression patterns and their clinical associations were explored through the 'Limma' R package, and their diagnostic potential was evaluated using ROC curves. Nitric oxide synthase 2 (NOS2) emerged as a candidate for further analyses. The CIBERSORT algorithm and GSEA dataset were used to respectively investigate the immune and metabolism effects of NOS2; the former was validated through immunofluorescence. The GSDC database was employed to analyze the correlation between NOS2 expression and the therapeutic efficacy of multiple drugs. PCR, Western blotting, colony formation assays, and Transwell experiments, were used to determine biological functions of NOS2 in HB cells. Potential upstream transcription factors of NOS2 were predicted through the TRRUST, hTFtarget, GeneCards, and JASPAR databases.

Results

NQO1, SLC1A4, and NOS2 were identified as potential genes in HB and found to be significantly upregulated in tumor samples. Nevertheless, only NOS2 was closely associated with HB clinicopathological characteristics; high NOS2 expression indicated poor prognosis, metastatic tendency, and late clinical stage. Immune analyses indicated that high NOS2 expression was concomitant with decreased infiltration levels of CD8+ T cells but increased infiltration levels of macrophages. GSEA revealed that NOS2 failed to affect the enrichments of glycolysis, fatty acid metabolism, and cholesterol biosynthesis in HB. Moreover, NOS2 was positively correlated with the IC50 values of trametinib, lapatinib, and cisplatin. NOS2 overexpression promoted the proliferation, migration and invasion of HepG2 and HuH-6 cells. JUND was identified as a potential transcriptional regulator of NOS2 by binding to its promoter (5'-TTCTGACTCTTTT-3').

Conclusion

NOS2 plays a significant role in HB clinical assessments and holds promise as a novel therapeutic target.

PSF toolkit: an R package for pathway curation and topology-aware analysis

Most high throughput genomic data analysis pipelines currently rely on over-representation or gene set enrichment analysis (ORA/GSEA) approaches for functional analysis. In contrast, topology-based pathway analysis methods, which offer a more biologically informed perspective by incorporating interaction and topology information, have remained underutilized and inaccessible due to various limiting factors. These methods heavily rely on the quality of pathway topologies and often utilize predefined topologies from databases without assessing their correctness. To address these issues and make topology-aware pathway analysis more accessible and flexible, we introduce the PSF (Pathway Signal Flow) toolkit R package. Our toolkit integrates pathway curation and topology-based analysis, providing interactive and command-line tools that facilitate pathway importation, correction, and modification from diverse sources. This enables users to perform topology-based pathway signal flow analysis in both interactive and command-line modes. To showcase the toolkit's usability, we curated 36 KEGG signaling pathways and conducted several use-case studies, comparing our method with ORA and the topology-based signaling pathway impact analysis (SPIA) method. The results demonstrate that the algorithm can effectively identify ORA enriched pathways while providing more detailed branch-level information. Moreover, in contrast to the SPIA method, it offers the advantage of being cut-off free and less susceptible to the variability caused by selection thresholds. By combining pathway curation and topology-based analysis, the PSF toolkit enhances the quality, flexibility, and accessibility of topology-aware pathway analysis. Researchers can now easily import pathways from various sources, correct and modify them as needed, and perform detailed topology-based pathway signal flow analysis. In summary, our PSF toolkit offers an integrated solution that addresses the limitations of current topology-based pathway analysis methods. By providing interactive and command-line tools for pathway curation and topology-based analysis, we empower researchers to conduct comprehensive pathway analyses across a wide range of applications.

Pathway-level multi-omics analysis of the molecular mechanisms underlying the toxicity of long-term tacrolimus exposure

Tacrolimus (TAC)-based treatment is associated with nephrotoxicity and hepatotoxicity; however, the underlying molecular mechanisms responsible for this toxicity have not been fully explored. This study elucidated the molecular processes underlying the toxic effects of TAC using an integrative omics approach. Rats were sacrificed after 4 weeks of daily oral TAC administration at a dose of 5 mg/kg. The liver and kidney underwent genome-wide gene expression profiling and untargeted metabolomics assays. Molecular alterations were identified using individual data profiling modalities and further characterized by pathway-level transcriptomics-metabolomics integration analysis. Metabolic disturbances were mainly related to an imbalance in oxidant-antioxidant status, as well as in lipid and amino acid metabolism in the liver and kidney. Gene expression profiles also indicated profound molecular alterations, including in genes associated with a dysregulated immune response, proinflammatory signals, and programmed cell death in the liver and kidney. Joint-pathway analysis indicated that the toxicity of TAC was associated with DNA synthesis disruption, oxidative stress, and cell membrane permeabilization, as well as lipid and glucose metabolism. In conclusion, our pathway-level integration of transcriptome and metabolome and conventional analyses of individual omics profiles, provided a more comprehensive picture of the molecular changes resulting from TAC toxicity. This study also serves as a valuable resource for subsequent investigations aiming to understand the mechanism underlying the molecular toxicology of TAC.

Molecular Mechanisms of Endometriosis Revealed Using Omics Data

Endometriosis is a gynecological disorder prevalent in women of reproductive age. The primary symptoms include dysmenorrhea, irregular menstruation, and infertility. However, the pathogenesis of endometriosis remains unclear. With the advent of high-throughput technologies, various omics experiments have been conducted to identify genes related to the pathophysiology of endometriosis. This review highlights the molecular mechanisms underlying endometriosis using omics. When genes identified in omics experiments were compared with endometriosis disease genes identified in independent studies, the number of overlapping genes was moderate. However, the characteristics of these genes were found to be equivalent when functional gene set enrichment analysis was performed using gene ontology and biological pathway information. These findings indicate that omics technology provides invaluable information regarding the pathophysiology of endometriosis. Moreover, the functional characteristics revealed using enrichment analysis provide important clues for discovering endometriosis disease genes in future research.

Flame (v2.0): advanced integration and interpretation of functional enrichment results from multiple sources

Functional enrichment is the process of identifying implicated functional terms from a given input list of genes or proteins. In this article, we present Flame (v2.0), a web tool which offers a combinatorial approach through merging and visualizing results from widely used functional enrichment applications while also allowing various flexible input options. In this version, Flame utilizes the aGOtool, g: Profiler, WebGestalt, and Enrichr pipelines and presents their outputs separately or in combination following a visual analytics approach. For intuitive representations and easier interpretation, it uses interactive plots such as parameterizable networks, heatmaps, barcharts, and scatter plots. Users can also: (i) handle multiple protein/gene lists and analyse union and intersection sets simultaneously through interactive UpSet plots, (ii) automatically extract genes and proteins from free text through text-mining and Named Entity Recognition (NER) techniques, (iii) upload single nucleotide polymorphisms (SNPs) and extract their relative genes, or (iv) analyse multiple lists of differentially expressed proteins/genes after selecting them interactively from a parameterizable volcano plot. Compared to the previous version of 197 supported organisms, Flame (v2.0) currently allows enrichment for 14 436 organisms.

AVAILABILITY AND IMPLEMENTATION

Web Application: http://flame.pavlopouloslab.info. Code: https://github.com/PavlopoulosLab/Flame. Docker: https://hub.docker.com/r/pavlopouloslab/flame.

Coherent pathway enrichment estimation by modeling inter-pathway dependencies using regularized regression

MOTIVATION

Gene set enrichment methods are a common tool to improve the interpretability of gene lists as obtained, for example, from differential gene expression analyses. They are based on computing whether dysregulated genes are located in certain biological pathways more often than expected by chance. Gene set enrichment tools rely on pre-existing pathway databases such as KEGG, Reactome, or the Gene Ontology. These databases are increasing in size and in the number of redundancies between pathways, which complicates the statistical enrichment computation.

RESULTS

We address this problem and develop a novel gene set enrichment method, called pareg, which is based on a regularized generalized linear model and directly incorporates dependencies between gene sets related to certain biological functions, for example, due to shared genes, in the enrichment computation. We show that pareg is more robust to noise than competing methods. Additionally, we demonstrate the ability of our method to recover known pathways as well as to suggest novel treatment targets in an exploratory analysis using breast cancer samples from TCGA.

AVAILABILITY AND IMPLEMENTATION

pareg is freely available as an R package on Bioconductor (https://bioconductor.org/packages/release/bioc/html/pareg.html) as well as on https://github.com/cbg-ethz/pareg. The GitHub repository also contains the Snakemake workflows needed to reproduce all results presented here.

Exploratory Assessment of Proteomic Network Changes in Cerebrospinal Fluid of Mild Cognitive Impairment Patients: A Pilot Study

(1) Background: Despite the existence of well-established, CSF-based biomarkers such as amyloid-β and phosphorylated-tau, the pathways involved in the pathophysiology of Alzheimer's disease (AD) remain an active area of research. (2) Methods: We measured 3072 proteins in CSF samples of AD-biomarker positive mild cognitive impairment (MCI) participants (n = 38) and controls (n = 48), using the Explore panel of the Olink proximity extension assay (PEA). We performed group comparisons, association studies with diagnosis, age, and APOE ε4 status, overrepresentation analysis (ORA), and gene set enrichment analysis (GSEA) to determine differentially expressed proteins and dysregulated pathways. (3) Results: GSEA results demonstrated an enrichment of granulocyte-related and chemotactic pathways (core enrichment proteins: ITGB2, ITGAM, ICAM1, SELL, SELP, C5, IL1A). Moreover, some of the well-replicated, differentially expressed proteins in CSF included: ITGAM, ITGB2, C1QA, TREM2, GFAP, NEFL, MMP-10, and a novel tau-related marker, SCRN1. (4) Conclusion: Our results highlight the upregulation of neuroinflammatory pathways, especially chemotactic and granulocyte recruitment in CSF of early AD patients.

Towards artificial intelligence-enabled extracellular vesicle precision drug delivery

Extracellular Vesicles (EVs), particularly exosomes, recently exploded into nanomedicine as an emerging drug delivery approach due to their superior biocompatibility, circulating stability, and bioavailability in vivo. However, EV heterogeneity makes molecular targeting precision a critical challenge. Deciphering key molecular drivers for controlling EV tissue targeting specificity is in great need. Artificial intelligence (AI) brings powerful prediction ability for guiding the rational design of engineered EVs in precision control for drug delivery. This review focuses on cutting-edge nano-delivery via integrating large-scale EV data with AI to develop AI-directed EV therapies and illuminate the clinical translation potential. We briefly review the current status of EVs in drug delivery, including the current frontier, limitations, and considerations to advance the field. Subsequently, we detail the future of AI in drug delivery and its impact on precision EV delivery. Our review discusses the current universal challenge of standardization and critical considerations when using AI combined with EVs for precision drug delivery. Finally, we will conclude this review with a perspective on future clinical translation led by a combined effort of AI and EV research.

An innovative index to incorporate transcriptomic data into weight of evidence approaches for environmental risk assessment

The sharp decrease in the cost of RNA-sequencing and the rapid improvement in computational analysis of eco-toxicogenomic data have brought new insights into the adverse effects of chemicals on aquatic organisms. Yet, transcriptomics is generally applied qualitatively in environmental risk assessments, hampering more effective exploitation of this evidence through multidisciplinary studies. In view of this limitation, a methodology is here presented to quantitatively elaborate transcriptional data in support to environmental risk assessment. The proposed methodology makes use of results from the application of Gene Set Enrichment Analysis to recent studies investigating the response of Mytilus galloprovincialis and Ruditapes philippinarum exposed to contaminants of emerging concern. The degree of changes in gene sets and the relevance of physiological reactions are integrated in the calculation of a hazard index. The outcome is then classified according to five hazard classes (from absent to severe), providing an evaluation of whole-transcriptome effects of chemical exposure. The application to experimental and simulated datasets proved that the method can effectively discriminate different levels of altered transcriptomic responses when compared to expert judgement (Spearman correlation coefficient of 0.96). A further application to data collected in two independent studies of Salmo trutta and Xenopus tropicalis exposed to contaminants confirmed the potential extension of the methodology to other aquatic species. This methodology can serve as a proof of concept for the integration of "genomic tools" in environmental risk assessment based on multidisciplinary investigations. To this end, the proposed transcriptomic hazard index can now be incorporated into quantitative Weight of Evidence approaches and weighed, with results from other types of analysis, to elucidate the role of chemicals in adverse ecological effects.

Transcriptomic and Proteomic of Gastrocnemius Muscle in Peripheral Artery Disease

BACKGROUND

Few effective therapies exist to improve lower extremity muscle pathology and mobility loss due to peripheral artery disease (PAD), in part because mechanisms associated with functional impairment remain unclear.

METHODS

To better understand mechanisms of muscle impairment in PAD, we performed in-depth transcriptomic and proteomic analyses on gastrocnemius muscle biopsies from 31 PAD participants (mean age, 69.9 years) and 29 age- and sex-matched non-PAD controls (mean age, 70.0 years) free of diabetes or limb-threatening ischemia.

RESULTS

Transcriptomic and proteomic analyses suggested activation of hypoxia-compensatory mechanisms in PAD muscle, including inflammation, fibrosis, apoptosis, angiogenesis, unfolded protein response, and nerve and muscle repair. Stoichiometric proportions of mitochondrial respiratory proteins were aberrant in PAD compared to non-PAD, suggesting that respiratory proteins not in complete functional units are not removed by mitophagy, likely contributing to abnormal mitochondrial activity. Supporting this hypothesis, greater mitochondrial respiratory protein abundance was significantly associated with greater complex II and complex IV respiratory activity in non-PAD but not in PAD. Rate-limiting glycolytic enzymes, such as hexokinase and pyruvate kinase, were less abundant in muscle of people with PAD compared with non-PAD participants, suggesting diminished glucose metabolism.

CONCLUSIONS

In PAD muscle, hypoxia induces accumulation of mitochondria respiratory proteins, reduced activity of rate-limiting glycolytic enzymes, and an enhanced integrated stress response that modulates protein translation. These mechanisms may serve as targets for disease modification.

Integrating genetics and transcriptomics to study major depressive disorder: a conceptual framework, bioinformatic approaches, and recent findings

Major depressive disorder (MDD) is a complex and heterogeneous psychiatric syndrome with genetic and environmental influences. In addition to neuroanatomical and circuit-level disturbances, dysregulation of the brain transcriptome is a key phenotypic signature of MDD. Postmortem brain gene expression data are uniquely valuable resources for identifying this signature and key genomic drivers in human depression; however, the scarcity of brain tissue limits our capacity to observe the dynamic transcriptional landscape of MDD. It is therefore crucial to explore and integrate depression and stress transcriptomic data from numerous, complementary perspectives to construct a richer understanding of the pathophysiology of depression. In this review, we discuss multiple approaches for exploring the brain transcriptome reflecting dynamic stages of MDD: predisposition, onset, and illness. We next highlight bioinformatic approaches for hypothesis-free, genome-wide analyses of genomic and transcriptomic data and their integration. Last, we summarize the findings of recent genetic and transcriptomic studies within this conceptual framework.

MYB insufficiency disrupts proteostasis in hematopoietic stem cells, leading to age-related neoplasia

MYB plays a key role in gene regulation throughout the hematopoietic hierarchy and is critical for the maintenance of normal hematopoietic stem cells (HSC). Acquired genetic dysregulation of MYB is involved in the etiology of a number of leukemias, although inherited noncoding variants of the MYB gene are a susceptibility factor for many hematological conditions, including myeloproliferative neoplasms (MPN). The mechanisms that connect variations in MYB levels to disease predisposition, especially concerning age dependency in disease initiation, are completely unknown. Here, we describe a model of Myb insufficiency in mice that leads to MPN, myelodysplasia, and leukemia in later life, mirroring the age profile of equivalent human diseases. We show that this age dependency is intrinsic to HSC, involving a combination of an initial defective cellular state resulting from small effects on the expression of multiple genes and a progressive accumulation of further subtle changes. Similar to previous studies showing the importance of proteostasis in HSC maintenance, we observed altered proteasomal activity and elevated proliferation indicators, followed by elevated ribosome activity in young Myb-insufficient mice. We propose that these alterations combine to cause an imbalance in proteostasis, potentially creating a cellular milieu favoring disease initiation.

Interpreting omics data with pathway enrichment analysis

Pathway enrichment analysis is indispensable for interpreting omics datasets and generating hypotheses. However, the foundations of enrichment analysis remain elusive to many biologists. Here, we discuss best practices in interpreting different types of omics data using pathway enrichment analysis and highlight the importance of considering intrinsic features of various types of omics data. We further explain major components that influence the outcomes of a pathway enrichment analysis, including defining background sets and choosing reference annotation databases. To improve reproducibility, we describe how to standardize reporting methodological details in publications. This article aims to serve as a primer for biologists to leverage the wealth of omics resources and motivate bioinformatics tool developers to enhance the power of pathway enrichment analysis.

NaRnEA: An Information Theoretic Framework for Gene Set Analysis

Gene sets are being increasingly leveraged to make high-level biological inferences from transcriptomic data; however, existing gene set analysis methods rely on overly conservative, heuristic approaches for quantifying the statistical significance of gene set enrichment. We created Nonparametric analytical-Rank-based Enrichment Analysis (NaRnEA) to facilitate accurate and robust gene set analysis with an optimal null model derived using the information theoretic Principle of Maximum Entropy. By measuring the differential activity of ~2500 transcriptional regulatory proteins based on the differential expression of each protein's transcriptional targets between primary tumors and normal tissue samples in three cohorts from The Cancer Genome Atlas (TCGA), we demonstrate that NaRnEA critically improves in two widely used gene set analysis methods: Gene Set Enrichment Analysis (GSEA) and analytical-Rank-based Enrichment Analysis (aREA). We show that the NaRnEA-inferred differential protein activity is significantly correlated with differential protein abundance inferred from independent, phenotype-matched mass spectrometry data in the Clinical Proteomic Tumor Analysis Consortium (CPTAC), confirming the statistical and biological accuracy of our approach. Additionally, our analysis crucially demonstrates that the sample-shuffling empirical null models leveraged by GSEA and aREA for gene set analysis are overly conservative, a shortcoming that is avoided by the newly developed Maximum Entropy analytical null model employed by NaRnEA.

Abolishing storage lipids induces protein misfolding and stress responses in Yarrowia lipolytica

Yarrowia lipolytica naturally saves excess carbon as storage lipids. Engineering efforts allow redirecting the high precursor flux required for lipid synthesis toward added-value chemicals such as polyketides, flavonoids, and terpenoids. To redirect precursor flux from storage lipids to other products, four genes involved in triacylglycerol and sterol ester synthesis (DGA1, DGA2, LRO1, and ARE1) can be deleted. To elucidate the effect of the deletions on cell physiology and regulation, we performed chemostat cultivations under carbon and nitrogen limitations, followed by transcriptome analysis. We found that storage lipid-free cells show an enrichment of the unfolded protein response, and several biological processes related to protein refolding and degradation are enriched. Additionally, storage lipid-free cells show an altered lipid class distribution with an abundance of potentially cytotoxic free fatty acids under nitrogen limitation. Our findings not only highlight the importance of lipid metabolism on cell physiology and proteostasis, but can also aid the development of improved chassy strains of Y. lipolytica for commodity chemical production.

Genome-Wide Gene-Set Analysis Identifies Molecular Mechanisms Associated with ALS

Amyotrophic lateral sclerosis (ALS) is a fatal late-onset motor neuron disease characterized by the loss of the upper and lower motor neurons. Our understanding of the molecular basis of ALS pathology remains elusive, complicating the development of efficient treatment. Gene-set analyses of genome-wide data have offered insight into the biological processes and pathways of complex diseases and can suggest new hypotheses regarding causal mechanisms. Our aim in this study was to identify and explore biological pathways and other gene sets having genomic association to ALS. Two cohorts of genomic data from the dbGaP repository were combined: (a) the largest available ALS individual-level genotype dataset (N = 12,319), and (b) a similarly sized control cohort (N = 13,210). Following comprehensive quality control pipelines, imputation and meta-analysis, we assembled a large European descent ALS-control cohort of 9244 ALS cases and 12,795 healthy controls represented by genetic variants of 19,242 genes. Multi-marker analysis of genomic annotation (MAGMA) gene-set analysis was applied to an extensive collection of 31,454 gene sets from the molecular signatures database (MSigDB). Statistically significant associations were observed for gene sets related to immune response, apoptosis, lipid metabolism, neuron differentiation, muscle cell function, synaptic plasticity and development. We also report novel interactions between gene sets, suggestive of mechanistic overlaps. A manual meta-categorization and enrichment mapping approach is used to explore the overlap of gene membership between significant gene sets, revealing a number of shared mechanisms.

Cortical profiles of numerous psychiatric disorders and normal development share a common pattern

The neurobiological bases of the association between development and psychopathology remain poorly understood. Here, we identify a shared spatial pattern of cortical thickness (CT) in normative development and several psychiatric and neurological disorders. Principal component analysis (PCA) was applied to CT of 68 regions in the Desikan-Killiany atlas derived from three large-scale datasets comprising a total of 41,075 neurotypical participants. PCA produced a spatially broad first principal component (PC1) that was reproducible across datasets. Then PC1 derived from healthy adult participants was compared to the pattern of CT differences associated with psychiatric and neurological disorders comprising a total of 14,886 cases and 20,962 controls from seven ENIGMA disease-related working groups, normative maturation and aging comprising a total of 17,697 scans from the ABCD Study® and the IMAGEN developmental study, and 17,075 participants from the ENIGMA Lifespan working group, as well as gene expression maps from the Allen Human Brain Atlas. Results revealed substantial spatial correspondences between PC1 and widespread lower CT observed in numerous psychiatric disorders. Moreover, the PC1 pattern was also correlated with the spatial pattern of normative maturation and aging. The transcriptional analysis identified a set of genes including KCNA2, KCNS1 and KCNS2 with expression patterns closely related to the spatial pattern of PC1. The gene category enrichment analysis indicated that the transcriptional correlations of PC1 were enriched to multiple gene ontology categories and were specifically over-represented starting at late childhood, coinciding with the onset of significant cortical maturation and emergence of psychopathology during the prepubertal-to-pubertal transition. Collectively, the present study reports a reproducible latent pattern of CT that captures interregional profiles of cortical changes in both normative brain maturation and a spectrum of psychiatric disorders. The pubertal timing of the expression of PC1-related genes implicates disrupted neurodevelopment in the pathogenesis of the spectrum of psychiatric diseases emerging during adolescence.

Overview of gene expression techniques with an emphasis on vitamin D related studies

Each cell controls when and how its genes must be expressed for proper function. Every function in a cell is driven by signaling molecules through various regulatory cascades. Different cells in a multicellular organism may express very different sets of genes, even though they contain the same DNA. The set of genes expressed in a cell determines the set of proteins and functional RNAs it contains, giving it its unique properties. Malfunction in gene expression harms the cell and can lead to the development of various disease conditions. The use of rapid high-throughput gene expression profiling unravels the complexity of human disease at various levels. Peripheral blood mononuclear cells (PBMC) have been used frequently to understand gene expression homeostasis in various disease conditions. However, more studies are required to validate whether PBMC gene expression patterns accurately reflect the expression of other cells or tissues. Vitamin D, which is responsible for a multitude of health consequences, is also an immune modulatory hormone with major biological activities in the innate and adaptive immune systems. Vitamin D exerts its diverse biological effects in target tissues by regulating gene expression and its deficiency, is recognized as a public health problem worldwide. Understanding the genetic factors that affect vitamin D has the potential benefit that it will make it easier to identify individuals who require supplementation. Different technological advances in gene expression can be used to identify and assess the severity of disease and aid in the development of novel therapeutic interventions. This review focuses on different gene expression approaches and various clinical studies of vitamin D to investigate the role of gene expression in identifying the molecular signature of the disease.

GSCA: an integrated platform for gene set cancer analysis at genomic, pharmacogenomic and immunogenomic levels

Cancer initiation and progression are likely caused by the dysregulation of biological pathways. Gene set analysis (GSA) could improve the signal-to-noise ratio and identify potential biological insights on the gene set level. However, platforms exploring cancer multi-omics data using GSA methods are lacking. In this study, we upgraded our GSCALite to GSCA (gene set cancer analysis, http://bioinfo.life.hust.edu.cn/GSCA) for cancer GSA at genomic, pharmacogenomic and immunogenomic levels. In this improved GSCA, we integrated expression, mutation, drug sensitivity and clinical data from four public data sources for 33 cancer types. We introduced useful features to GSCA, including associations between immune infiltration with gene expression and genomic variations, and associations between gene set expression/mutation and clinical outcomes. GSCA has four main functional modules for cancer GSA to explore, analyze and visualize expression, genomic variations, tumor immune infiltration, drug sensitivity and their associations with clinical outcomes. We used case studies of three gene sets: (i) seven cell cycle genes, (ii) tumor suppressor genes of PI3K pathway and (iii) oncogenes of PI3K pathway to prove the advantage of GSCA over single gene analysis. We found novel associations of gene set expression and mutation with clinical outcomes in different cancer types on gene set level, while on single gene analysis level, they are not significant associations. In conclusion, GSCA is a user-friendly web server and a useful resource for conducting hypothesis tests by using GSA methods at genomic, pharmacogenomic and immunogenomic levels.

From SNP to pathway-based GWAS meta-analysis: do current meta-analysis approaches resolve power and replication in genetic association studies?

Genome-wide association studies (GWAS) have benefited greatly from enhanced high-throughput technology in recent decades. GWAS meta-analysis has become increasingly popular to highlight the genetic architecture of complex traits, informing about the replicability and variability of effect estimations across human ancestries. A wealth of GWAS meta-analysis methodologies have been developed depending on the input data and the outcome information of interest. We present a survey of current approaches from SNP to pathway-based meta-analysis by acknowledging the range of resources and methodologies in the field, and we provide a comprehensive review of different categories of Genome-Wide Meta-analysis methods employed. These methods highlight different levels at which GWAS meta-analysis may be done, including Single Nucleotide Polymorphisms, Genes and Pathways, for which we describe their framework outline. We also discuss the strengths and pitfalls of each approach and make suggestions regarding each of them.

The hitchhikers' guide to RNA sequencing and functional analysis

DNA and RNA sequencing technologies have revolutionized biology and biomedical sciences, sequencing full genomes and transcriptomes at very high speeds and reasonably low costs. RNA sequencing (RNA-Seq) enables transcript identification and quantification, but once sequencing has concluded researchers can be easily overwhelmed with questions such as how to go from raw data to differential expression (DE), pathway analysis and interpretation. Several pipelines and procedures have been developed to this effect. Even though there is no unique way to perform RNA-Seq analysis, it usually follows these steps: 1) raw reads quality check, 2) alignment of reads to a reference genome, 3) aligned reads' summarization according to an annotation file, 4) DE analysis and 5) gene set analysis and/or functional enrichment analysis. Each step requires researchers to make decisions, and the wide variety of options and resulting large volumes of data often lead to interpretation challenges. There also seems to be insufficient guidance on how best to obtain relevant information and derive actionable knowledge from transcription experiments. In this paper, we explain RNA-Seq steps in detail and outline differences and similarities of different popular options, as well as advantages and disadvantages. We also discuss non-coding RNA analysis, multi-omics, meta-transcriptomics and the use of artificial intelligence methods complementing the arsenal of tools available to researchers. Lastly, we perform a complete analysis from raw reads to DE and functional enrichment analysis, visually illustrating how results are not absolute truths and how algorithmic decisions can greatly impact results and interpretation.

Single-Cell Transcriptome Identifies the Renal Cell Type Tropism of Human BK Polyomavirus

BK polyomavirus (BKPyV) infection is the main factor affecting the prognosis of kidney transplant recipients, as no antiviral agent is yet available. A better understanding of the renal-cell-type tropism of BKPyV can serve to develop new treatment strategies. In this study, the single-cell transcriptomic analysis demonstrated that the ranking of BKPyV tropism for the kidney was proximal tubule cells (PT), collecting duct cells (CD), and glomerular endothelial cells (GEC) according to the signature of renal cell type and immune microenvironment. In normal kidneys, we found that BKPyV infection-related transcription factors P65 and CEBPB were PT-specific transcription factors, and PT showed higher glycolysis/gluconeogenesis activities than CD and GEC. Furthermore, in the BKPyV-infected kidneys, the percentage of late viral transcripts in PT was significantly higher than in CD and GEC. In addition, PT had the smallest cell-cell interactions with immune cells compared to CD and GEC in both normal and BKPyV-infected kidneys. Subsequently, we indirectly demonstrated the ranking of BKPyV tropism via the clinical observation of sequential biopsies. Together, our results provided in-depth insights into the renal cell-type tropism of BKPyV in vivo at single-cell resolution and proposed a novel antiviral target.