A single-cell and spatial RNA-seq database for Alzheimer's disease (ssREAD).

Alzheimer's Disease (AD) pathology has been increasingly explored through single-cell and single-nucleus RNA-sequencing (scRNA-seq & snRNA-seq) and spatial transcriptomics (ST). However, the surge in data demands a comprehensive, user-friendly repository. Addressing this, we introduce a single-cell and spatial RNA-seq database for Alzheimer's disease (ssREAD). It offers a broader spectrum of AD-related datasets, an optimized analytical pipeline, and improved usability. The database encompasses 1,053 samples (277 integrated datasets) from 67 AD-related scRNA-seq & snRNA-seq studies, totaling 7,332,202 cells. Additionally, it archives 381 ST datasets from 18 human and mouse brain studies. Each dataset is annotated with details such as species, gender, brain region, disease/control status, age, and AD Braak stages. ssREAD also provides an analysis suite for cell clustering, identification of differentially expressed and spatially variable genes, cell-type-specific marker genes and regulons, and spot deconvolution for integrative analysis. ssREAD is freely available at https://bmblx.bmi.osumc.edu/ssread/ .

A weighted two-stage sequence alignment framework to identify motifs from ChIP-exo data.

In this study, we introduce TESA (weighted two-stage alignment), an innovative motif prediction tool that refines the identification of DNA-binding protein motifs, essential for deciphering transcriptional regulatory mechanisms. Unlike traditional algorithms that rely solely on sequence data, TESA integrates the high-resolution chromatin immunoprecipitation (ChIP) signal, specifically from ChIP-exonuclease (ChIP-exo), by assigning weights to sequence positions, thereby enhancing motif discovery. TESA employs a nuanced approach combining a binomial distribution model with a graph model, further supported by a "bookend" model, to improve the accuracy of predicting motifs of varying lengths. Our evaluation, utilizing an extensive compilation of 90 prokaryotic ChIP-exo datasets from proChIPdb and 167 H. sapiens datasets, compared TESA's performance against seven established tools. The results indicate TESA's improved precision in motif identification, suggesting its valuable contribution to the field of genomic research.

Targeting metabolic sensing switch GPR84 on macrophages for cancer immunotherapy.

INTRODUCTION: As one of the major components of the tumor microenvironment, tumor-associated macrophages (TAMs) possess profound inhibitory activity against T cells and facilitate tumor escape from immune checkpoint blockade therapy. Converting this pro-tumorigenic toward the anti-tumorigenic phenotype thus is an important strategy for enhancing adaptive immunity against cancer. However, a plethora of mechanisms have been described for pro-tumorigenic differentiation in cancer, metabolic switches to program the anti-tumorigenic property of TAMs are elusive. MATERIALS AND METHODS: From an unbiased analysis of single-cell transcriptome data from multiple tumor models, we discovered that anti-tumorigenic TAMs uniquely express elevated levels of a specific fatty acid receptor, G-protein-coupled receptor 84 (GPR84). Genetic ablation of GPR84 in mice leads to impaired pro-inflammatory polarization of macrophages, while enhancing their anti-inflammatory phenotype. By contrast, GPR84 activation by its agonist, 6-n-octylaminouracil (6-OAU), potentiates pro-inflammatory phenotype via the enhanced STAT1 pathway. Moreover, 6-OAU treatment significantly retards tumor growth and increases the anti-tumor efficacy of anti-PD-1 therapy. CONCLUSION: Overall, we report a previously unappreciated fatty acid receptor, GPR84, that serves as an important metabolic sensing switch for orchestrating anti-tumorigenic macrophage polarization. Pharmacological agonists of GPR84 hold promise to reshape and reverse the immunosuppressive TME, and thereby restore responsiveness of cancer to overcome resistance to immune checkpoint blockade.

Graph Fourier transform for spatial omics representation and analyses of complex organs.

Spatial omics technologies are capable of deciphering detailed components of complex organs or tissue in cellular and subcellular resolution. A robust, interpretable, and unbiased representation method for spatial omics is necessary to illuminate novel investigations into biological functions, whereas a mathematical theory deficiency still exists. We present SpaGFT (Spatial Graph Fourier Transform), which provides a unique analytical feature representation of spatial omics data and elucidates molecular signatures linked to critical biological processes within tissues and cells. It outperformed existing tools in spatially variable gene prediction and gene expression imputation across human/mouse Visium data. Integrating SpaGFT representation into existing machine learning frameworks can enhance up to 40% accuracy of spatial domain identification, cell type annotation, cell-to-spot alignment, and subcellular hallmark inference. SpaGFT identified immunological regions for B cell maturation in human lymph node Visium data, characterized secondary follicle variations from in-house human tonsil CODEX data, and detected extremely rare subcellular organelles such as Cajal body and Set1/COMPASS. This new method lays the groundwork for a new theoretical model in explainable AI, advancing our understanding of tissue organization and function.

MarsGT: Multi-omics analysis for rare population inference using single-cell graph transformer.

Rare cell populations are key in neoplastic progression and therapeutic response, offering potential intervention targets. However, their computational identification and analysis often lag behind major cell types. To fill this gap, we introduce MarsGT: Multi-omics Analysis for Rare population inference using a Single-cell Graph Transformer. It identifies rare cell populations using a probability-based heterogeneous graph transformer on single-cell multi-omics data. MarsGT outperforms existing tools in identifying rare cells across 550 simulated and four real human datasets. In mouse retina data, it reveals unique subpopulations of rare bipolar cells and a Müller glia cell subpopulation. In human lymph node data, MarsGT detects an intermediate B cell population potentially acting as lymphoma precursors. In human melanoma data, it identifies a rare MAIT-like population impacted by a high IFN-I response and reveals the mechanism of immunotherapy. Hence, MarsGT offers biological insights and suggests potential strategies for early detection and therapeutic intervention of disease.

A Bioinformatics Tool for Identifying Intratumoral Microbes from the ORIEN Dataset.

Evidence supports significant interactions among microbes, immune cells, and tumor cells in at least 10%-20% of human cancers, emphasizing the importance of further investigating these complex relationships. However, the implications and significance of tumor-related microbes remain largely unknown. Studies have demonstrated the critical roles of host microbes in cancer prevention and treatment responses. Understanding interactions between host microbes and cancer can drive cancer diagnosis and microbial therapeutics (bugs as drugs). Computational identification of cancer-specific microbes and their associations is still challenging due to the high dimensionality and high sparsity of intratumoral microbiome data, which requires large datasets containing sufficient event observations to identify relationships, and the interactions within microbial communities, the heterogeneity in microbial composition, and other confounding effects that can lead to spurious associations. To solve these issues, we present a bioinformatics tool, microbial graph attention (MEGA), to identify the microbes most strongly associated with 12 cancer types. We demonstrate its utility on a dataset from a consortium of nine cancer centers in the Oncology Research Information Exchange Network. This package has three unique features: species-sample relations are represented in a heterogeneous graph and learned by a graph attention network; it incorporates metabolic and phylogenetic information to reflect intricate relationships within microbial communities; and it provides multiple functionalities for association interpretations and visualizations. We analyzed 2,704 tumor RNA sequencing samples and MEGA interpreted the tissue-resident microbial signatures of each of 12 cancer types. MEGA can effectively identify cancer-associated microbial signatures and refine their interactions with tumors. SIGNIFICANCE: Studying the tumor microbiome in high-throughput sequencing data is challenging because of the extremely sparse data matrices, heterogeneity, and high likelihood of contamination. We present a new deep learning tool, MEGA, to refine the organisms that interact with tumors.

Inference of disease-associated microbial gene modules based on metagenomic and metatranscriptomic data.

The identification of microbial characteristics associated with diseases is crucial for disease diagnosis and therapy. However, the presence of heterogeneity, high dimensionality, and large amounts of microbial data presents tremendous challenges in discovering key microbial features. In this paper, we present IDAM, a novel computational method for inferring disease-associated gene modules from metagenomic and metatranscriptomic data. This method integrates gene context conservation (uber-operons) and regulatory mechanisms (gene co-expression patterns) within a mathematical graph model to explore gene modules associated with specific diseases. It alleviates reliance on prior meta-data. We applied IDAM to publicly available datasets from inflammatory bowel disease, melanoma, type 1 diabetes mellitus, and irritable bowel syndrome. The results demonstrated the superior performance of IDAM in inferring disease-associated characteristics compared to existing popular tools. Furthermore, we showcased the high reproducibility of the gene modules inferred by IDAM using independent cohorts with inflammatory bowel disease. We believe that IDAM can be a highly advantageous method for exploring disease-associated microbial characteristics. The source code of IDAM is freely available at https://github.com/OSU-BMBL/IDAM, and the web server can be accessed at https://bmblx.bmi.osumc.edu/idam/.

ssREAD: A Single-cell and Spatial RNA-seq Database for Alzheimer's Disease.

Alzheimer's Disease (AD) is a neurodegenerative malady predominantly affecting the elderly and exhibits its debilitating effects on a dementia-prone population. Recently, the advent of innovative technologies, such as single-cell and single-nucleus RNA-sequencing (scRNA-seq & snRNA-seq) and spatial transcriptomics (ST), has reformed our investigative approaches toward comprehending AD's neuropathological intricacies and underpinning regulatory mechanisms, encompassing sub-cellular, cellular, and spatial dimensions. In light of the overwhelming proliferation of single-cell and ST data associated with AD, the imperative for a comprehensive, user-friendly database that addresses the scientific community's analytical demands has never been more paramount. Introduced initially in 2020, scREAD presented itself as a pioneering repository that systematized publicly available scRNA-seq and snRNA-seq datasets derived from post-mortem human brain tissues and mouse models mirroring AD pathology. Here, we introduce ssREAD, a substantial upgrade over scREAD, enriching the platform with a broader spectrum of datasets, an optimized analytical pipeline, and enhanced usability and visibility. Specifically, ssREAD amalgamates an impressive portfolio of over 189 datasets extracted from 35 distinct AD-related scRNA-seq and snRNA-seq studies, encompassing a staggering 2,572,355 cells. In addition, we have diligently curated and archived 300 ST datasets, originating from 12 human and mouse brain studies, which include two focused on AD and ten control studies. Every dataset within our repository is meticulously annotated, bearing critical identifiers including species, gender, brain region, disease/control status, age, and AD stages. Besides the collection of above datasets in ssREAD, it delivers an exhaustive analysis suite offering cell clustering and annotation, inference of differentially expressed and spatially variable genes, identification of cell-type-specific marker genes and regulons, and spot deconvolution for integrative analysis of ST and scRNA-seq & snRNA-seq data from public domains. All these resources are freely accessible through a user-friendly, consolidated web portal available at https://bmblx.bmi.osumc.edu/ssread/.

Identification of HPV oncogene and host cell differentiation associated cellular heterogeneity in cervical cancer via single-cell transcriptomic analysis.

Human Papillomaviruses (HPVs) are associated with around 5%-10% of human cancer, notably nearly 99% of cervical cancer. The mechanisms HPV interacts with stratified epithelium (differentiated layers) during the viral life cycle, and oncogenesis remain unclear. In this study, we used single-cell transcriptome analysis to study viral gene and host cell differentiation-associated heterogeneity of HPV-positive cervical cancer tissue. We examined the HPV16 genes-E1, E6, and E7, and found they expressed differently across nine epithelial clusters. We found that three epithelial clusters had the highest proportion of HPV-positive cells (33.6%, 37.5%, and 32.4%, respectively), while two exhibited the lowest proportions (7.21% and 5.63%, respectively). Notably, the cluster with the most HPV-positive cells deviated significantly from normal epithelial layer markers, exhibiting functional heterogeneity and altered epithelial structuring, indicating that significant molecular heterogeneity existed in cancer tissues and that these cells exhibited unique/different gene signatures compared with normal epithelial cells. These HPV-positive cells, compared to HPV-negative, showed different gene expressions related to the extracellular matrix, cell adhesion, proliferation, and apoptosis. Further, the viral oncogenes E6 and E7 appeared to modify epithelial function via distinct pathways, thus contributing to cervical cancer progression. We investigated the HPV and host transcripts from a novel viewpoint focusing on layer heterogeneity. Our results indicated varied HPV expression across epithelial clusters and epithelial heterogeneity associated with viral oncogenes, contributing biological insights to this critical field of study.

Identification of HPV oncogene and host cell differentiation associated cellular heterogeneity in cervical cancer via single-cell transcriptomic analysis.

Human Papillomaviruses (HPVs) are associated with around 5-10% of human cancer, notably nearly 99% of cervical cancer. The mechanisms HPV interacts with stratified epithelium (differentiated layers) during the viral life cycle, and oncogenesis remain unclear. In this study, we used single-cell transcriptome analysis to study viral gene and host cell differentiation-associated heterogeneity of HPV-positive cervical cancer tissue. We examined the HPV16 genes - E1, E6, and E7, and found they expressed differently across nine epithelial clusters. We found that three epithelial clusters had the highest proportion of HPV-positive cells (33.6%, 37.5%, and 32.4%, respectively), while two exhibited the lowest proportions (7.21% and 5.63%, respectively). Notably, the cluster with the most HPV-positive cells deviated significantly from normal epithelial layer markers, exhibiting functional heterogeneity and altered epithelial structuring, indicating that significant molecular heterogeneity existed in cancer tissues and that these cells exhibited unique/different gene signatures compared with normal epithelial cells. These HPV-positive cells, compared to HPV-negative, showed different gene expressions related to the extracellular matrix, cell adhesion, proliferation, and apoptosis. Further, the viral oncogenes E6 and E7 appeared to modify epithelial function via distinct pathways, thus contributing to cervical cancer progression. We investigated the HPV and host transcripts from a novel viewpoint focusing on layer heterogeneity. Our results indicated varied HPV expression across epithelial clusters and epithelial heterogeneity associated with viral oncogenes, contributing biological insights to this critical field of study.

MarsGT: Multi-omics analysis for rare population inference using single-cell graph transformer.

Rare cell populations are key in neoplastic progression and therapeutic response, offering potential intervention targets. However, their computational identification and analysis often lag behind major cell types. To fill this gap, we introduced MarsGT: Multi-omics Analysis for Rare population inference using Single-cell Graph Transformer. It identifies rare cell populations using a probability-based heterogeneous graph transformer on single-cell multi-omics data. MarsGT outperformed existing tools in identifying rare cells across 400 simulated and four real human datasets. In mouse retina data, it revealed unique subpopulations of rare bipolar cells and a Müller glia cell subpopulation. In human lymph node data, MarsGT detected an intermediate B cell population potentially acting as lymphoma precursors. In human melanoma data, it identified a rare MAIT-like population impacted by a high IFN-I response and revealed the mechanism of immunotherapy. Hence, MarsGT offers biological insights and suggests potential strategies for early detection and therapeutic intervention of disease.

Repression of ferroptotic cell death by mitochondrial calcium signaling.

The uptake of Ca2+ into and extrusion of calcium from the mitochondrial matrix, regulated by the mitochondrial Ca2+ uniporter (MCU), is a fundamental biological process that has crucial impacts on cellular metabolism, signaling, growth and survival. Herein, we report that the embryonic lethality of Mcu-deficient mice is fully rescued by orally supplementing ferroptosis inhibitor lipophilic antioxidant vitamin E and ubiquinol. Mechanistically, we found MCU promotes acetyl-CoA-mediated GPX4 acetylation at K90 residue, and K90R mutation impaired the GPX4 enzymatic activity, a step that is crucial for ferroptosis. Structural analysis supports the possibility that GPX4 K90R mutation alters the conformational state of the molecule, resulting in disruption of a salt bridge formation with D23, which was confirmed by mutagenesis studies. Finally, we report that deletion of MCU in cancer cells caused a marked reduction in tumor growth in multiple cancer models. In summary, our study provides a first direct link between mitochondrial calcium level and sustained GPX4 enzymatic activity to regulate ferroptosis, which consequently protects cancer cells from ferroptosis.

A bioinformatics tool for identifying intratumoral microbes from the ORIEN dataset.

Evidence supports significant interactions among microbes, immune cells, and tumor cells in at least 10-20% of human cancers, emphasizing the importance of further investigating these complex relationships. However, the implications and significance of tumor-related microbes remain largely unknown. Studies have demonstrated the critical roles of host microbes in cancer prevention and treatment responses. Understanding interactions between host microbes and cancer can drive cancer diagnosis and microbial therapeutics (bugs as drugs). Computational identification of cancer-specific microbes and their associations is still challenging due to the high dimensionality and high sparsity of intratumoral microbiome data, which requires large datasets containing sufficient event observations to identify relationships, and the interactions within microbial communities, the heterogeneity in microbial composition, and other confounding effects that can lead to spurious associations. To solve these issues, we present a bioinformatics tool, MEGA, to identify the microbes most strongly associated with 12 cancer types. We demonstrate its utility on a dataset from a consortium of 9 cancer centers in the Oncology Research Information Exchange Network (ORIEN). This package has 3 unique features: species-sample relations are represented in a heterogeneous graph and learned by a graph attention network; it incorporates metabolic and phylogenetic information to reflect intricate relationships within microbial communities; and it provides multiple functionalities for association interpretations and visualizations. We analyzed 2704 tumor RNA-seq samples and MEGA interpreted the tissue-resident microbial signatures of each of 12 cancer types. MEGA can effectively identify cancer-associated microbial signatures and refine their interactions with tumors.

Single-cell biological network inference using a heterogeneous graph transformer.

Single-cell multi-omics (scMulti-omics) allows the quantification of multiple modalities simultaneously to capture the intricacy of complex molecular mechanisms and cellular heterogeneity. Existing tools cannot effectively infer the active biological networks in diverse cell types and the response of these networks to external stimuli. Here we present DeepMAPS for biological network inference from scMulti-omics. It models scMulti-omics in a heterogeneous graph and learns relations among cells and genes within both local and global contexts in a robust manner using a multi-head graph transformer. Benchmarking results indicate DeepMAPS performs better than existing tools in cell clustering and biological network construction. It also showcases competitive capability in deriving cell-type-specific biological networks in lung tumor leukocyte CITE-seq data and matched diffuse small lymphocytic lymphoma scRNA-seq and scATAC-seq data. In addition, we deploy a DeepMAPS webserver equipped with multiple functionalities and visualizations to improve the usability and reproducibility of scMulti-omics data analysis.

Explainable Deep Hypergraph Learning Modeling the Peptide Secondary Structure Prediction.

Accurately predicting peptide secondary structures remains a challenging task due to the lack of discriminative information in short peptides. In this study, PHAT is proposed, a deep hypergraph learning framework for the prediction of peptide secondary structures and the exploration of downstream tasks. The framework includes a novel interpretable deep hypergraph multi-head attention network that uses residue-based reasoning for structure prediction. The algorithm can incorporate sequential semantic information from large-scale biological corpus and structural semantic information from multi-scale structural segmentation, leading to better accuracy and interpretability even with extremely short peptides. The interpretable models are able to highlight the reasoning of structural feature representations and the classification of secondary substructures. The importance of secondary structures in peptide tertiary structure reconstruction and downstream functional analysis is further demonstrated, highlighting the versatility of our models. To facilitate the use of the model, an online server is established which is accessible via http://inner.wei-group.net/PHAT/. The work is expected to assist in the design of functional peptides and contribute to the advancement of structural biology research.

Computational methods and challenges in analyzing intratumoral microbiome data.

The human microbiome is intimately related to cancer biology and plays a vital role in the efficacy of cancer treatments, including immunotherapy. Extraordinary evidence has revealed that several microbes influence tumor development through interaction with the host immune system, that is, immuno-oncology-microbiome (IOM). This review focuses on the intratumoral microbiome in IOM and describes the available data and computational methods for discovering biological insights of microbial profiling from host bulk, single-cell, and spatial sequencing data. Critical challenges in data analysis and integration are discussed. Specifically, the microorganisms associated with cancer and cancer treatment in the context of IOM are collected and integrated from the literature. Lastly, we provide our perspectives for future directions in IOM research.

scGNN 2.0: a graph neural network tool for imputation and clustering of single-cell RNA-Seq data.

MOTIVATION: Gene expression imputation has been an essential step of the single-cell RNA-Seq data analysis workflow. Among several deep-learning methods, the debut of scGNN gained substantial recognition in 2021 for its superior performance and the ability to produce a cell-cell graph. However, the implementation of scGNN was relatively time-consuming and its performance could still be optimized. RESULTS: The implementation of scGNN 2.0 is significantly faster than scGNN thanks to a simplified close-loop architecture. For all eight datasets, cell clustering performance was increased by 85.02% on average in terms of adjusted rand index, and the imputation Median L1 Error was reduced by 67.94% on average. With the built-in visualizations, users can quickly assess the imputation and cell clustering results, compare against benchmarks and interpret the cell-cell interaction. The expanded input and output formats also pave the way for custom workflows that integrate scGNN 2.0 with other scRNA-Seq toolkits on both Python and R platforms. AVAILABILITY AND IMPLEMENTATION: scGNN 2.0 is implemented in Python (as of version 3.8) with the source code available at https://github.com/OSU-BMBL/scGNN2.0. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Deep transfer learning of cancer drug responses by integrating bulk and single-cell RNA-seq data.

Drug screening data from massive bulk gene expression databases can be analyzed to determine the optimal clinical application of cancer drugs. The growing amount of single-cell RNA sequencing (scRNA-seq) data also provides insights into improving therapeutic effectiveness by helping to study the heterogeneity of drug responses for cancer cell subpopulations. Developing computational approaches to predict and interpret cancer drug response in single-cell data collected from clinical samples can be very useful. We propose scDEAL, a deep transfer learning framework for cancer drug response prediction at the single-cell level by integrating large-scale bulk cell-line data. The highlight in scDEAL involves harmonizing drug-related bulk RNA-seq data with scRNA-seq data and transferring the model trained on bulk RNA-seq data to predict drug responses in scRNA-seq. Another feature of scDEAL is the integrated gradient feature interpretation to infer the signature genes of drug resistance mechanisms. We benchmark scDEAL on six scRNA-seq datasets and demonstrate its model interpretability via three case studies focusing on drug response label prediction, gene signature identification, and pseudotime analysis. We believe that scDEAL could help study cell reprogramming, drug selection, and repurposing for improving therapeutic efficacy.

SUSD2 suppresses CD8+ T cell antitumor immunity by targeting IL-2 receptor signaling.

Dysfunctional CD8+ T cells, which have defective production of antitumor effectors, represent a major mediator of immunosuppression in the tumor microenvironment. Here, we show that SUSD2 is a negative regulator of CD8+ T cell antitumor function. Susd2-/- effector CD8+ T cells showed enhanced production of antitumor molecules, which consequently blunted tumor growth in multiple syngeneic mouse tumor models. Through a quantitative mass spectrometry assay, we found that SUSD2 interacted with interleukin (IL)-2 receptor α through sushi domain-dependent protein interactions and that this interaction suppressed the binding of IL-2, an essential cytokine for the effector functions of CD8+ T cells, to IL-2 receptor α. SUSD2 was not expressed on regulatory CD4+ T cells and did not affect the inhibitory function of these cells. Adoptive transfer of Susd2-/- chimeric antigen receptor T cells induced a robust antitumor response in mice, highlighting the potential of SUSD2 as an immunotherapy target for cancer.

Define and visualize pathological architectures of human tissues from spatially resolved transcriptomics using deep learning.

Spatially resolved transcriptomics provides a new way to define spatial contexts and understand the pathogenesis of complex human diseases. Although some computational frameworks can characterize spatial context via various clustering methods, the detailed spatial architectures and functional zonation often cannot be revealed and localized due to the limited capacities of associating spatial information. We present RESEPT, a deep-learning framework for characterizing and visualizing tissue architecture from spatially resolved transcriptomics. Given inputs such as gene expression or RNA velocity, RESEPT learns a three-dimensional embedding with a spatial retained graph neural network from spatial transcriptomics. The embedding is then visualized by mapping into color channels in an RGB image and segmented with a supervised convolutional neural network model. Based on a benchmark of 10x Genomics Visium spatial transcriptomics datasets on the human and mouse cortex, RESEPT infers and visualizes the tissue architecture accurately. It is noteworthy that, for the in-house AD samples, RESEPT can localize cortex layers and cell types based on pre-defined region- or cell-type-enriched genes and furthermore provide critical insights into the identification of amyloid-beta plaques in Alzheimer's disease. Interestingly, in a glioblastoma sample analysis, RESEPT distinguishes tumor-enriched, non-tumor, and regions of neuropil with infiltrating tumor cells in support of clinical and prognostic cancer applications.