Autoencoder to Identify Sex-Specific Sub-phenotypes in Alzheimer's Disease Progression Using Longitudinal Electronic Health Records.

Alzheimer's Disease (AD) is a complex neurodegenerative disorder significantly influenced by sex differences, with approximately two-thirds of AD patients being women. Characterizing the sex-specific AD progression and identifying its progression trajectory is a crucial step to developing effective risk stratification and prevention strategies. In this study, we developed an autoencoder to uncover sex-specific sub-phenotypes in AD progression leveraging longitudinal electronic health record (EHR) data from OneFlorida+ Clinical Research Consortium. Specifically, we first constructed temporal patient representation using longitudinal EHRs from a sex-stratified AD cohort. We used a long short-term memory (LSTM)-based autoencoder to extract and generate latent representation embeddings from sequential clinical records of patients. We then applied hierarchical agglomerative clustering to the learned representations, grouping patients based on their progression sub-phenotypes. The experimental results show we successfully identified five primary sex-based AD sub-phenotypes with corresponding progression pathways with high confidence. These sex-specific sub-phenotypes not only illustrated distinct AD progression patterns but also revealed differences in clinical characteristics and comorbidities between females and males in AD development. These findings could provide valuable insights for advancing personalized AD intervention and treatment strategies.

Enhancer-driven gene regulatory networks inference from single-cell RNA-seq and ATAC-seq data.

Deciphering the intricate relationships between transcription factors (TFs), enhancers, and genes through the inference of enhancer-driven gene regulatory networks (eGRNs) is crucial in understanding gene regulatory programs in a complex biological system. This study introduces STREAM, a novel method that leverages a Steiner forest problem model, a hybrid biclustering pipeline, and submodular optimization to infer eGRNs from jointly profiled single-cell transcriptome and chromatin accessibility data. Compared to existing methods, STREAM demonstrates enhanced performance in terms of TF recovery, TF-enhancer linkage prediction, and enhancer-gene relation discovery. Application of STREAM to an Alzheimer's disease dataset and a diffuse small lymphocytic lymphoma dataset reveals its ability to identify TF-enhancer-gene relations associated with pseudotime, as well as key TF-enhancer-gene relations and TF cooperation underlying tumor cells.

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.

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.

Research on property and burning behavior of flammable casing for underground coal gasification.

In this work, the comprehensive properties of flammable casing for underground coal gasification is systematically investigated, including flammable casing material physical, chemical and mechanical properties and full-size flammable casing mechanical properties and burning behavior. The flammable casing material consists of magnesium alloy matrix and rare earth particles, thermal conductivity and expansion property of which are weak. Results of high-temperature tensile test reveal that flammable casing material has good high temperature strength which declines by 30 % at 300 °C. Corrosion rate of flammable casing material is relatively high without extra protection. The full-size flammable casing possesses considerable mechanical property, thread property and high temperature collapse resistance. Burning of flammable casing is safe and stable. Burning rate of flammable casing material can be effectively controlled by water flow. Combustion product of flammable casing presents powder condition, which has no risk of blocking the gasification channel. To sum up, flammable casing is necessary to the realization of underground coal gasifying process, which plays the significant role of the development and application of underground coal gasification technology.

Distinctive patterns of sulfide- and butyrate-metabolizing bacteria after bariatric surgery: potential implications for colorectal cancer risk.

Despite improved cardiometabolic outcomes following bariatric surgery, its long-term impact on colorectal cancer (CRC) risk remains uncertain. In parallel, the influence of bariatric surgery on the host microbiome and relationships with disease outcomes is beginning to be appreciated. Therefore, we investigated the impact of Roux-en-Y gastric bypass (RYGB) and vertical sleeve gastrectomy (VSG) on the patterns of sulfide-reducing and butyrate-producing bacteria, which are hypothesized to modulate CRC risk after bariatric surgery. In this single-center, cross-sectional study, we included 15 pre-surgery subjects with severe obesity and patients who are at a median (range) of 25.6 (9.9-46.5) months after RYGB (n = 16) or VSG (n = 10). The DNA abundance of fecal bacteria and enzymes involved in butyrate and sulfide metabolism were identified using metagenomic sequencing. Differences between pre-surgery and post-RYGB or post-VSG cohorts were quantified using the linear discriminant analysis (LDA) effect size (LEfSe) method. Our sample was predominantly female (87%) with a median (range) age of 46 (23-71) years. Post-RYGB and post-VSG patients had a higher DNA abundance of fecal sulfide-reducing bacteria than pre-surgery controls (LDA = 1.3-4.4, p < .05). The most significant enrichments were for fecal E. coli, Acidaminococcus and A. finegoldii after RYGB, and for A. finegoldii, S. vestibularis, V. parvula after VSG. As for butyrate-producing bacteria, R. faecis was more abundant, whereas B. dentium and A. hardus were lower post-RYGB vs. pre-surgery. B. dentium was also lower in post-VSG vs. pre-surgery. Consistent with these findings, our analysis showed a greater enrichment of sulfide-reducing enzymes after bariatric surgery, especially RYGB, vs. pre-surgery. The DNA abundance of butyrate-producing enzymes was lower post-RYGB. In conclusion, the two most used bariatric surgeries, RYGB and VSG, are associated with microbiome patterns that are potentially implicated in CRC risk. Future studies are needed to validate and understand the impact of these microbiome changes on CRC risk after bariatric surgery.

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.

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.

Impaired Human Cardiac Cell Development due to NOTCH1 Deficiency.

BACKGROUND: NOTCH1 pathogenic variants are implicated in multiple types of congenital heart defects including hypoplastic left heart syndrome, where the left ventricle is underdeveloped. It is unknown how NOTCH1 regulates human cardiac cell lineage determination and cardiomyocyte proliferation. In addition, mechanisms by which NOTCH1 pathogenic variants lead to ventricular hypoplasia in hypoplastic left heart syndrome remain elusive. METHODS: CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats)/Cas9 genome editing was utilized to delete NOTCH1 in human induced pluripotent stem cells. Cardiac differentiation was carried out by sequential modulation of WNT signaling, and NOTCH1 knockout and wild-type differentiating cells were collected at day 0, 2, 5, 10, 14, and 30 for single-cell RNA-seq. RESULTS: Human NOTCH1 knockout induced pluripotent stem cells are able to generate functional cardiomyocytes and endothelial cells, suggesting that NOTCH1 is not required for mesoderm differentiation and cardiovascular development in vitro. However, disruption of NOTCH1 blocks human ventricular-like cardiomyocyte differentiation but promotes atrial-like cardiomyocyte generation through shortening the action potential duration. NOTCH1 deficiency leads to defective proliferation of early human cardiomyocytes, and transcriptomic analysis indicates that pathways involved in cell cycle progression and mitosis are downregulated in NOTCH1 knockout cardiomyocytes. Single-cell transcriptomic analysis reveals abnormal cell lineage determination of cardiac mesoderm, which is manifested by the biased differentiation toward epicardial and second heart field progenitors at the expense of first heart field progenitors in NOTCH1 knockout cell populations. CONCLUSIONS: NOTCH1 is essential for human ventricular-like cardiomyocyte differentiation and proliferation through balancing cell fate determination of cardiac mesoderm and modulating cell cycle progression. Because first heart field progenitors primarily contribute to the left ventricle, we speculate that pathogenic NOTCH1 variants lead to biased differentiation of first heart field progenitors, blocked ventricular-like cardiomyocyte differentiation, and defective cardiomyocyte proliferation, which collaboratively contribute to left ventricular hypoplasia in hypoplastic left heart syndrome.

CEMIG: prediction of the cis-regulatory motif using the de Bruijn graph from ATAC-seq.

Sequence motif discovery algorithms enhance the identification of novel deoxyribonucleic acid sequences with pivotal biological significance, especially transcription factor (TF)-binding motifs. The advent of assay for transposase-accessible chromatin using sequencing (ATAC-seq) has broadened the toolkit for motif characterization. Nonetheless, prevailing computational approaches have focused on delineating TF-binding footprints, with motif discovery receiving less attention. Herein, we present Cis rEgulatory Motif Influence using de Bruijn Graph (CEMIG), an algorithm leveraging de Bruijn and Hamming distance graph paradigms to predict and map motif sites. Assessment on 129 ATAC-seq datasets from the Cistrome Data Browser demonstrates CEMIG's exceptional performance, surpassing three established methodologies on four evaluative metrics. CEMIG accurately identifies both cell-type-specific and common TF motifs within GM12878 and K562 cell lines, demonstrating its comparative genomic capabilities in the identification of evolutionary conservation and cell-type specificity. In-depth transcriptional and functional genomic studies have validated the functional relevance of CEMIG-identified motifs across various cell types. CEMIG is available at https://github.com/OSU-BMBL/CEMIG, developed in C++ to ensure cross-platform compatibility with Linux, macOS and Windows operating systems.

Spatially resolved transcriptomics reveals genes associated with the vulnerability of middle temporal gyrus in Alzheimer's disease.

Human middle temporal gyrus (MTG) is a vulnerable brain region in early Alzheimer's disease (AD), but little is known about the molecular mechanisms underlying this regional vulnerability. Here we utilize the 10 × Visium platform to define the spatial transcriptomic profile in both AD and control (CT) MTG. We identify unique marker genes for cortical layers and the white matter, and layer-specific differentially expressed genes (DEGs) in human AD compared to CT. Deconvolution of the Visium spots showcases the significant difference in particular cell types among cortical layers and the white matter. Gene co-expression analyses reveal eight gene modules, four of which have significantly altered co-expression patterns in the presence of AD pathology. The co-expression patterns of hub genes and enriched pathways in the presence of AD pathology indicate an important role of cell-cell-communications among microglia, oligodendrocytes, astrocytes, and neurons, which may contribute to the cellular and regional vulnerability in early AD. Using single-molecule fluorescent in situ hybridization, we validated the cell-type-specific expression of three novel DEGs (e.g., KIF5A, PAQR6, and SLC1A3) and eleven previously reported DEGs associated with AD pathology (i.e., amyloid beta plaques and intraneuronal neurofibrillary tangles or neuropil threads) at the single cell level. Our results may contribute to the understanding of the complex architecture and neuronal and glial response to AD pathology of this vulnerable brain region.

Molecular Regulation of Heme Oxygenase-1 Expression by E2F Transcription Factor 2 in Lung Fibroblast Cells: Relevance to Idiopathic Pulmonary Fibrosis.

Idiopathic pulmonary fibrosis (IPF) is a fatal chronic lung disease. Heme oxygenase-1 (HMOX1/HO-1) is an enzyme that catalyzes the degradation of heme. The role of HO-1 in the pathogenesis of IPF has been studied; however, the molecular regulation of HO-1 and its role in IPF are still unclear. In this study, we found that HO-1 protein levels significantly increased in lung myofibroblasts in IPF patients and in lungs in a murine model of bleomycin-induced lung fibrosis. In addition, we observed that administration of a E2F transcription factor inhibitor elevated HO-1 mRNA and protein levels in lung fibroblasts. Downregulation of E2F2 by siRNA transfection increased HO-1 mRNA and protein levels, while overexpression of E2F2 reduced HO-1 levels. However, overexpression of E2F2 did not alter hemin-induced HO-1 protein levels. Furthermore, modulation of HO-1 levels regulated TGF-ß1-induced myofibroblast differentiation without altering the phosphorylation of Smad2/3 in lung fibroblast cells. Moreover, the phosphorylation of protein kinase B (Akt) was significantly upregulated in HO-1-depleted lung fibroblast cells. In summary, this study demonstrated that E2F2 regulates the baseline expression of HO-1, but has no effect on modulating HO-1 expression by hemin. Finally, elevated HO-1 expression contributes to the TGF-ß1-induced lung myofibroblast differentiation through the activation of the serine/threonine kinase AKT pathway. Overall, our findings suggest that targeting E2F2/HO-1 might be a new therapeutic strategy to treat fibrotic diseases such as IPF.

DESSO-DB: A web database for sequence and shape motif analyses and identification.

Cis-regulatory motif (motif for short) identification and analyses are essential steps in detecting gene regulatory mechanisms. Deep learning (DL) models have shown substantial advances in motif prediction. In parallel, intuitive and integrative web databases are needed to make effective use of DL models and ensure easy access to the identified motifs. Here, we present DESSO-DB, a web database developed to allow efficient access to the identified motifs and diverse motif analyses. DESSO-DB provides motif prediction results and visualizations of 690 ENCODE human Chromatin Immunoprecipitation sequencing (ChIP-seq) data (including 161 transcription factors (TFs) in 91 cell lines) and 1,677 human ChIP-seq data (including 547 TFs in 359 cell lines) from Cistrome DB using DESSO, which is an in-house developed DL tool for motif prediction. It also provides online motif finding and scanning functions for new ChIP-seq/ATAC-seq datasets and downloadable motif results of the above 690 DECODE datasets, 126 cancer ChIP-seq, 55 RNA Crosslinking-Immunoprecipitation and high-throughput sequencing (CLIP-seq) data. DESSO-DB is deployed on the Google Cloud Platform, providing stabilized and efficient resources freely to the public. DESSO-DB is free and available at http://cloud.osubmi.com/DESSO/.

Microglia coordinate cellular interactions during spinal cord repair in mice.

Traumatic spinal cord injury (SCI) triggers a neuro-inflammatory response dominated by tissue-resident microglia and monocyte derived macrophages (MDMs). Since activated microglia and MDMs are morphologically identical and express similar phenotypic markers in vivo, identifying injury responses specifically coordinated by microglia has historically been challenging. Here, we pharmacologically depleted microglia and use anatomical, histopathological, tract tracing, bulk and single cell RNA sequencing to reveal the cellular and molecular responses to SCI controlled by microglia. We show that microglia are vital for SCI recovery and coordinate injury responses in CNS-resident glia and infiltrating leukocytes. Depleting microglia exacerbates tissue damage and worsens functional recovery. Conversely, restoring select microglia-dependent signaling axes, identified through sequencing data, in microglia depleted mice prevents secondary damage and promotes recovery. Additional bioinformatics analyses reveal that optimal repair after SCI might be achieved by co-opting key ligand-receptor interactions between microglia, astrocytes and MDMs.