Computational solutions for spatial transcriptomics

Adhesion-Related Pathways and Functional Polarization of Astrocytes in Traumatic Brain Injury: Insights from Single-cell RNA Sequencing

Traumatic brain injury (TBI) induces profound functional heterogeneity in astrocytes, yet the regulatory mechanisms underlying this diversity remain poorly understood. In this study, we analyzed single-cell RNA sequencing data from the cortex and hippocampus of TBI mouse models to characterize astrocyte subtypes and their functional dynamics. We identified two major reactive subtypes: A1 astrocytes, enriched in inflammatory response, synaptic regulation, and neurodegenerative disease-related pathways; and A2 astrocytes, enriched in lipid metabolism, extracellular matrix (ECM) remodeling, and phagosome formation pathways. These functional differences were consistently observed across datasets with varying injury severities. Notably, adhesion-related pathways-including gap junctions, adherens junctions, and calcium-dependent adhesion-showed significant subtype-specific expression patterns and temporal shifts. Pseudotime trajectory analysis further suggested a potential transition between A1 and A2 states, accompanied by dynamic regulation of adhesion-related genes. Our findings highlight the complex and context-dependent roles of astrocytes in TBI and propose cell adhesion as a key modulator of astrocyte functional polarization.

Decoding immune cell interactions during cardiac allograft vasculopathy: insights derived from bioinformatic strategies

Chronic allograft vasculopathy (CAV) is a major cause of late graft failure in heart transplant recipients, characterized by progressive intimal thickening and diffuse narrowing of the coronary arteries. Unlike atherosclerosis, CAV exhibits a distinct cellular composition and lesion distribution, yet its pathogenesis remains incompletely understood. A major challenge in CAV research has been the limited application of advanced "-omics" technologies, which have revolutionized the study of other vascular diseases. Recent advancements in single-cell and spatial transcriptomics, proteomics, and metabolomics have begun to uncover the complex immune-endothelial-stromal interactions driving CAV progression. Notably, single-cell RNA sequencing has identified previously unrecognized immune cell populations and signaling pathways implicated in endothelial injury and vascular remodeling after heart transplantation. Despite these breakthroughs, studies applying these technologies to CAV remain sparse, limiting the translation of these insights into clinical practice. This review aims to bridge this gap by summarizing recent findings from single-cell and multi-omic approaches, highlighting key discoveries, and discussing their implications for understanding CAV pathogenesis.

SpaGRN: Investigating spatially informed regulatory paths for spatially resolved transcriptomics data

Cells spatially organize into distinct cell types or functional domains through localized gene regulatory networks. However, current spatially resolved transcriptomics analyses fail to integrate spatial constraints and proximal cell influences, limiting the mechanistic understanding of tissue organization. Here, we introduce SpaGRN, a statistical framework that reconstructs cell-type- or functional-domain-specific, dynamic, and spatial regulons by coupling intracellular spatial regulatory causality with extracellular signaling path information. Benchmarking across synthetic and real datasets demonstrates SpaGRN's superior precision over state-of-the-art tools in identifying context-dependent regulons. Applied to diverse spatially resolved transcriptomics platforms (Stereo-seq, STARmap, MERFISH, CosMx, Slide-seq, and 10x Visium), complex cancerous samples, and 3D datasets of developing Drosophila embryos and larvae, SpaGRN not only provides a versatile toolkit for decoding receptor-mediated spatial regulons but also reveals spatiotemporal regulatory mechanisms underlying organogenesis and inflammation.

CellRomeR: an R package for clustering cell migration phenotypes from microscopy data

Motivation

The analysis of cell migration using time-lapse microscopy typically focuses on track characteristics for classification and statistical evaluation of migration behaviour. However, considerable heterogeneity can be seen in cell morphology and microscope signal intensity features within the migrating cell populations.

Results

To utilize this information in cell migration analysis, we introduce here an R package CellRomeR, designed for the phenotypic clustering of cells based on their morphological and motility features from microscopy images. Utilizing machine learning techniques and building on an iterative clustering projection method, CellRomeR offers a new approach to identify heterogeneity in cell populations. The clustering of cells along the migration tracks allows association of distinct cellular phenotypes with different cell migration types and detection of migration patterns associated with stable and unstable cell phenotypes. The user-friendly interface of CellRomeR and multiple visualization options facilitate an in-depth understanding of cellular behaviour, addressing previous challenges in clustering cell trajectories using microscope cell tracking data.

Availability and implementation

CellRomeR is available as an R package from https://github.com/elolab/CellRomeR.

Machine learning-based transcriptmics analysis reveals BMX, GRB10, and GADD45A as crucial biomarkers and therapeutic targets in sepsis

Sepsis is a life-threatening condition characterized by a dysregulated host response to infection, resulting in high mortality rates and complex clinical management. This study leverages transcriptomics and machine learning (ML) to identify critical biomarkers and therapeutic targets in sepsis. Analyzing microarray data from the Gene Expression Omnibus (GEO) datasets GSE28750, GSE26440, GSE13205, and GSE9960, we discovered three pivotal biomarkers that BMX (bone marrow tyrosine kinase gene on chromosome X), GRB10 (growth factor receptor bound protein 10), and GADD45A (growth arrest and DNA damage inducible alpha), exhibiting exceptional diagnostic accuracy (AUC >0.9). Functional enrichment analyses revealed that these genes play key roles in reactive oxygen species metabolism and immune response regulation. Specifically, GADD45A was positively correlated with eosinophils and inversely associated with activated NK cells, CD8 T cells, and activated memory CD4 T cells. BMX showed positive correlations with eosinophils, mast cells, and neutrophils, while GRB10 was linked to eosinophils and M2 macrophages. Additionally, we constructed a comprehensive mRNA-miRNA-lncRNA regulatory network, identifying key interactions that may drive sepsis pathogenesis. Molecular docking and dynamics simulations validated Bendroflumethiazide, Cianidanol, and Hexamidine as promising therapeutic agents targeting these biomarkers. In conclusion, this integrated approach provides profound insights into the molecular mechanisms underlying sepsis, pinpointing BMX, GRB10, and GADD45A as pivotal biomarkers and therapeutic targets. These findings significantly enhance our understanding of sepsis pathophysiology and lay the groundwork for developing personalized diagnostic and therapeutic strategies aimed at improving patient outcomes.

Benchmarking spatial transcriptomics technologies with the multi-sample SpatialBenchVisium dataset

BACKGROUND

Spatial transcriptomics allows gene expression to be measured within complex tissue contexts. Among the array of spatial capture technologies available is 10x Genomics' Visium platform, a popular method which enables transcriptome-wide profiling of tissue sections. Visium offers a range of sample handling and library construction methods which introduces a need for benchmarking to compare data quality and assess how well the technology can recover expected tissue features and biological signatures.

RESULTS

Here we present SpatialBenchVisium, a unique reference dataset generated from spleen tissue of mice responding to malaria infection spanning several tissue preparation protocols (both fresh frozen and FFPE, with either manual or CytAssist tissue placement). We note better quality control metrics in reference samples prepared using probe-based capture methods, particularly those processed with CytAssist, validating the improvement in data quality produced with the platform. Our analysis of replicate samples extends to explore spatially variable gene detection, the outcomes of clustering and cell deconvolution using matched single-cell RNA-sequencing data and publicly available reference data to identify cell types and tissue regions expected in the spleen. Multi-sample differential expression analysis recovered known gene signatures related to biological sex or gene knockout.

Galangin ameliorates PTU-induced vitiligo in zebrafish and B16F10 cells by increasing melanogenesis through activation of the p38/JNK MAPK pathway

Objective

Alpinia officinarum Hance is a traditional herb in Xinjiang for the treatment of vitiligo, and galangin (GA) is a flavonoid isolated from its roots. However, its therapeutic mechanism remains unclear.

Methods

In this study, 1-phenyl-2-thiourea (PTU) was used to establish a vitiligo model in zebrafish. After successful modeling, different concentrations of GA (1 and 2 μM) were administered, and the distribution of melanin granules was observed by assaying the melanin content, masson-fontana staining and tyrosinase activity. Transcriptomic analysis and molecular docking were used to identify potential GA-related pathways and targets for improving vitiligo. In addition, we evaluated the proliferation of B16F10 cells by PTU induction and also observed cellular melanin distribution using masson-fontana staining. Finally, Western blot was performed to detect the proteins of the relevant pathways.

Results

The results showed that GA significantly increased melanin production and tyrosinase activity in depigmented zebrafish. In addition, we found that GA decreased ROS and MDA levels and increased the expression of GSH, CAT and T-SOD. In addition, transcriptome analysis indicated that GA likely acts through the mitogen-activated protein kinase (MAPK) signaling pathway. GA has a strong binding affinity for important targets.GA significantly increased the expression of genes such as mapk8b, mapk14a, mapk3, mitf, tyr, tyrp1b, tyrp1a, dct, and oca2, and decreased the expression of genes such as expression of genes such as raf1 and egfr. In addition, GA enhanced the viability of B16F10 cells, increased intracellular melanin content, and increased the expression of proteins such as p38, JNK1/2/3, TYR, MITF, TRP1, TRP2, and so on.

Conclusion

GA increases melanin production and distribution, improves tyrosinase activity, upregulates the expression of related genes and proteins through activation of MAPK and tyrosine metabolic pathways, downregulates oxidative stress, and then regulates changes in melanin synthesis to improve vitiligo.

How Can Spatial Transcriptomic Profiling Advance Our Understanding of Skin Diseases?

Spatial transcriptomic (ST) profiling is the mapping of gene expression within cell populations with preservation of positional context and represents an exciting new approach to develop our understanding of local and regional influences upon skin biology in health and disease. With the ability to probe from a few hundred transcripts to the entire transcriptome, multiple ST approaches are now widely available. In this paper, we review the ST field and discuss its application to dermatology. Its potential to advance our understanding of skin biology in health and disease is highlighted through the illustrative examples of 3 research areas: cutaneous aging, tumorigenesis, and psoriasis.

A best practices framework for spatial biology studies in drug discovery and development: enabling successful cohort studies using digital spatial profiling

The discovery of biomarkers, essential for successful drug development, is often hindered by the limited availability of tissue samples, typically obtained through core needle biopsies. Standard 'omics platforms can consume significant amounts of tissue, forcing scientist to trade off spatial context for high-plex assays, such as genome-wide assays. While bulk gene expression approaches and standard single-cell transcriptomics have been valuable in defining various molecular and cellular mechanisms, they do not retain spatial context. As such, they have limited power in resolving tissue heterogeneity and cell-cell interactions. Current spatial transcriptomics platforms offer limited transcriptome coverage and have low throughput, restricting the number of samples that can be analyzed daily or even weekly. While the Digital Spatial Profiling (DSP) method does not provide single-cell resolution, it presents a significant advancement by enabling scalable whole transcriptome and ultrahigh-plex protein analysis from distinct tissue compartments and structures using a single tissue slide. These capabilities overcome significant constraints in biomarker analysis in solid tissue specimens. These advancements in tissue profiling play a crucial role in deepening our understanding of disease biology and in identifying potential therapeutic targets and biomarkers. To enhance the use of spatial biology tools in drug discovery and development, the DSP Scientific Consortium has created best practices guidelines. These guidelines, built on digital spatial profiling data and expertise, offer a practical framework for designing spatial studies and using current and future spatial biology platforms. The aim is to improve tissue analysis in all research areas supporting drug discovery and development.

Desflurane attenuates renal ischemia-reperfusion injury by modulating ITGB1/CD9 and reducing oxidative stress in tubular epithelial cells

OBJECTIVE

Renal ischemia-reperfusion (I/R) injury triggers significant oxidative stress and inflammation, leading to tubular epithelial cell (TEC) damage. This study investigates the protective role of Desflurane (DFE) in renal I/R by modulating the ITGB1/CD9 signaling pathway and mitigating oxidative damage.

METHODS

Single-cell RNA sequencing (scRNA-seq) and transcriptome analysis identified ITGB1 as a key regulatory gene in TECs during renal I/R. The effects of DFE on ITGB1/CD9 expression were evaluated through in vitro experiments using RT-qPCR, Western blot, and TUNEL assays. A mouse model of renal I/R was employed to assess renal function and oxidative stress markers under DFE treatment.

RESULTS

DFE reduced ITGB1 and CD9 expression, resulting in decreased TEC apoptosis and enhanced proliferation. In vivo, DFE-treated mice exhibited improved renal function, with significantly lower serum creatinine and blood urea nitrogen (BUN) levels. Additionally, DFE reduced oxidative stress, as indicated by decreased malondialdehyde (MDA) and myeloperoxidase (MPO) activity, alongside increased superoxide dismutase (SOD) activity.

CONCLUSION

DFE confers renal protection against I/R by modulating the ITGB1/CD9 axis and reducing oxidative stress, offering a promising therapeutic strategy for mitigating kidney damage.

Spatial biology - unravelling complexity within the glioblastoma microenvironment

The advent and refinement of state-of-the-art spatial biology technologies have facilitated analysis that combines the advantages of high-throughput single cell analysis with techniques that preserve tissue architecture. This combination of cellular phenotyping with retained spatial context provides a much greater understanding of cellular interactions within the tumour microenvironment (TME). For glioblastoma, with its significant intra-tumoural heterogeneity, cellular plasticity, and complex TME, appreciating and understanding these spatial patterns may prove key to improving patient outcomes. This review examines the advances in spatial biology techniques, discusses how these methodologies are being applied to study glioblastoma, and explores how spatial information improves understanding of the TME. Ultimately, it is this spatial context that will accelerate the identification of more effective treatments for glioblastoma.

Single-cell sequencing in diabetic retinopathy: progress and prospects

Diabetic retinopathy is a major ocular complication of diabetes, characterized by progressive retinal microvascular damage and significant visual impairment in working-age adults. Traditional bulk RNA sequencing offers overall gene expression profiles but does not account for cellular heterogeneity. Single-cell RNA sequencing overcomes this limitation by providing transcriptomic data at the individual cell level and distinguishing novel cell subtypes, developmental trajectories, and intercellular communications. Researchers can use single-cell sequencing to draw retinal cell atlases and identify the transcriptomic features of retinal cells, enhancing our understanding of the pathogenesis and pathological changes in diabetic retinopathy. Additionally, single-cell sequencing is widely employed to analyze retinal organoids and single extracellular vesicles. Single-cell multi-omics sequencing integrates omics information, whereas stereo-sequencing analyzes gene expression and spatiotemporal data simultaneously. This review discusses the protocols of single-cell sequencing for obtaining single cells from retina and accurate sequencing data. It highlights the applications and advancements of single-cell sequencing in the study of normal retinas and the pathological changes associated with diabetic retinopathy. This underscores the potential of these technologies to deepen our understanding of the pathogenesis of diabetic retinopathy that may lead to the introduction of new therapeutic strategies.

Trends and Challenges of the Modern Pathology Laboratory for Biopharmaceutical Research Excellence

Pathology, a fundamental discipline that bridges basic scientific discovery to the clinic, is integral to successful drug development. Intrinsically multimodal and multidimensional, anatomic pathology continues to be empowered by advancements in molecular and digital technologies enabling the spatial tissue detection of biomolecules such as genes, transcripts, and proteins. Over the past two decades, breakthroughs in spatial molecular biology technologies and advancements in automation and digitization of laboratory processes have enabled the implementation of higher throughput assays and the generation of extensive molecular data sets from tissue sections in biopharmaceutical research and development research units. It is our goal to provide readers with some rationale, advice, and ideas to help establish a modern molecular pathology laboratory to meet the emerging needs of biopharmaceutical research. This manuscript provides (1) a high-level overview of the current state and future vision for excellence in research pathology practice and (2) shared perspectives on how to optimally leverage the expertise of discovery, toxicologic, and translational pathologists to provide effective spatial, molecular, and digital pathology data to support modern drug discovery. It captures insights from the experiences, challenges, and solutions from pathology laboratories of various biopharmaceutical organizations, including their approaches to troubleshooting and adopting new technologies.

Molecular profiling unveils pyroptosis markers in preterm birth

Through a comprehensive examination of pyroptosis-related differential expressed genes (PRDEGs), this work investigates the molecular complexities of spontaneous preterm birth (SPTB), also known as premature delivery, before the due date. Through the process of merging and correcting batch effects in the GSE120480 and GSE73714 datasets, we were able to identify 36 PRDEGs that exhibited significant expression differentiation in SPTB. Through functional enrichment and pathway analysis, their importance in amino acid transport and cytokine receptor interaction has been highlighted. Among the genes that have emerged as crucial, CEBPA, APOA1, and CEP55 have been identified. The relevance of these molecules was demonstrated using experimental knockdowns, which also suggested that they could be used as molecular biomarkers and therapeutic targets for SPTB.

Multi-omics in exploring the pathophysiology of diabetic retinopathy

Diabetic retinopathy (DR) is a leading global cause of vision impairment, with its prevalence increasing alongside the rising rates of diabetes mellitus (DM). Despite the retina's complex structure, the underlying pathology of DR remains incompletely understood. Single-cell RNA sequencing (scRNA-seq) and recent advancements in multi-omics analyses have revolutionized molecular profiling, enabling high-throughput analysis and comprehensive characterization of complex biological systems. This review highlights the significant contributions of scRNA-seq, in conjunction with other multi-omics technologies, to DR research. Integrated scRNA-seq and transcriptomic analyses have revealed novel insights into DR pathogenesis, including alternative transcription start site events, fluctuations in cell populations, altered gene expression profiles, and critical signaling pathways within retinal cells. Furthermore, by integrating scRNA-seq with genetic association studies and multi-omics analyses, researchers have identified novel biomarkers, susceptibility genes, and potential therapeutic targets for DR, emphasizing the importance of specific retinal cell types in disease progression. The integration of scRNA-seq with metabolomics has also been instrumental in identifying specific metabolites and dysregulated pathways associated with DR. It is highly conceivable that the continued synergy between scRNA-seq and other multi-omics approaches will accelerate the discovery of underlying mechanisms and the development of novel therapeutic interventions for DR.

VGAE-CCI: variational graph autoencoder-based construction of 3D spatial cell-cell communication network

Cell-cell communication plays a critical role in maintaining normal biological functions, regulating development and differentiation, and controlling immune responses. The rapid development of single-cell RNA sequencing and spatial transcriptomics sequencing (ST-seq) technologies provides essential data support for in-depth and comprehensive analysis of cell-cell communication. However, ST-seq data often contain incomplete data and systematic biases, which may reduce the accuracy and reliability of predicting cell-cell communication. Furthermore, other methods for analyzing cell-cell communication mainly focus on individual tissue sections, neglecting cell-cell communication across multiple tissue layers, and fail to comprehensively elucidate cell-cell communication networks within three-dimensional tissues. To address the aforementioned issues, we propose VGAE-CCI, a deep learning framework based on the Variational Graph Autoencoder, capable of identifying cell-cell communication across multiple tissue layers. Additionally, this model can be applied to spatial transcriptomics data with missing or partially incomplete data and can clustered cells at single-cell resolution based on spatial encoding information within complex tissues, thereby enabling more accurate inference of cell-cell communication. Finally, we tested our method on six datasets and compared it with other state of art methods for predicting cell-cell communication. Our method outperformed other methods across multiple metrics, demonstrating its efficiency and reliability in predicting cell-cell communication.

Application of spatial-omics to the classification of kidney biopsy samples in transplantation

Improvement of long-term outcomes through targeted treatment is a primary concern in kidney transplant medicine. Currently, the validation of a rejection diagnosis and subsequent treatment depends on the histological assessment of allograft biopsy samples, according to the Banff classification system. However, the lack of (early) disease-specific tissue markers hinders accurate diagnosis and thus timely intervention. This challenge mainly results from an incomplete understanding of the pathophysiological processes underlying late allograft failure. Integration of large-scale multimodal approaches for investigating allograft biopsy samples might offer new insights into this pathophysiology, which are necessary for the identification of novel therapeutic targets and the development of tailored immunotherapeutic interventions. Several omics technologies - including transcriptomic, proteomic, lipidomic and metabolomic tools (and multimodal data analysis strategies) - can be applied to allograft biopsy investigation. However, despite their successful application in research settings and their potential clinical value, several barriers limit the broad implementation of many of these tools into clinical practice. Among spatial-omics technologies, mass spectrometry imaging, which is under-represented in the transplant field, has the potential to enable multi-omics investigations that might expand the insights gained with current clinical analysis technologies.

Rediscovering the human thymus through cutting-edge technologies

Recent technological advances have transformed our understanding of the human thymus. Innovations such as high-resolution imaging, single-cell omics, and organoid cultures, including thymic epithelial cell (TEC) differentiation and culture, and improvements in biomaterials, have further elucidated the thymus architecture, cellular dynamics, and molecular mechanisms underlying T cell development, and have unraveled previously unrecognized levels of stromal cell heterogeneity. These advancements offer unprecedented insights into thymic biology and hold promise for the development of novel therapeutic strategies for immune-related disorders.

Hsp90aa1/JUN/Ccl2 regulatory axis mediates migration and differentiation of NSPCs, promoting the onset and progression of early post-ischemic stroke epilepsy

Early-onset epilepsy following ischemic stroke is a severe neurological condition, the pathogenesis of which remains incompletely understood. Recent studies suggest that Neural stem/progenitor cells (NSPCs) play a crucial role in the disease process, yet the precise molecular mechanisms regulating NSPCs have not been thoroughly investigated. This study utilized single-cell transcriptome sequencing and bioinformatics analysis to identify disease-related genes, which were subsequently validated in both in vitro and in vivo experiments. The findings revealed that Hsp90aa1 (heat shock protein 90 kDa alpha, class A member 1), Jun proto-oncogene (JUN), and CC Motif Ligation 2 (Ccl2) constitute an important regulatory axis influencing the migration and differentiation of NSPCs, potentially impacting the onset and progression of early-onset epilepsy post-ischemic stroke. Additionally, the expression of Hsp90aa1 was found to influence the likelihood of seizure occurrence and the severity of brain ischemia.

Next-generation spatial transcriptomics: unleashing the power to gear up translational oncology

The growing advances in spatial transcriptomics (ST) stand as the new frontier bringing unprecedented influences in the realm of translational oncology. This has triggered systemic experimental design, analytical scope, and depth alongside with thorough bioinformatics approaches being constantly developed in the last few years. However, harnessing the power of spatial biology and streamlining an array of ST tools to achieve designated research goals are fundamental and require real-world experiences. We present a systemic review by updating the technical scope of ST across different principal basis in a timeline manner hinting on the generally adopted ST techniques used within the community. We also review the current progress of bioinformatic tools and propose in a pipelined workflow with a toolbox available for ST data exploration. With particular interests in tumor microenvironment where ST is being broadly utilized, we summarize the up-to-date progress made via ST-based technologies by narrating studies categorized into either mechanistic elucidation or biomarker profiling (translational oncology) across multiple cancer types and their ways of deploying the research through ST. This updated review offers as a guidance with forward-looking viewpoints endorsed by many high-resolution ST tools being utilized to disentangle biological questions that may lead to clinical significance in the future.

lncRNAs as prognostic markers and therapeutic targets in cuproptosis-mediated cancer

Long non-coding RNAs (lncRNAs) have emerged as crucial regulators in various cellular processes, including cancer progression and stress response. Recent studies have demonstrated that copper accumulation induces a unique form of cell death known as cuproptosis, with lncRNAs playing a key role in regulating cuproptosis-associated pathways. These lncRNAs may trigger cell-specific responses to copper stress, presenting new opportunities as prognostic markers and therapeutic targets. This paper delves into the role of lncRNAs in cuproptosis-mediated cancer, underscoring their potential as biomarkers and targets for innovative therapeutic strategies. A thorough review of scientific literature was conducted, utilizing databases such as PubMed, Google Scholar, and ScienceDirect, with search terms like 'lncRNAs,' 'cuproptosis,' and 'cancer.' Studies were selected based on their relevance to lncRNA regulation of cuproptosis pathways and their implications for cancer prognosis and treatment. The review highlights the significant contribution of lncRNAs in regulating cuproptosis-related genes and pathways, impacting copper metabolism, mitochondrial stress responses, and apoptotic signaling. Specific lncRNAs are potential prognostic markers in breast, lung, liver, ovarian, pancreatic, and gastric cancers. The objective of this article is to explore the role of lncRNAs as potential prognostic markers and therapeutic targets in cancers mediated by cuproptosis.

Computational Strategies and Algorithms for Inferring Cellular Composition of Spatial Transcriptomics Data

Spatial transcriptomics technology has been an essential and powerful method for delineating tissue architecture at the molecular level. However, due to the limitations of the current spatial techniques, the cellular information cannot be directly measured but instead spatial spots typically varying from a diameter of 0.2 to 100 µm are characterized. Therefore, it is vital to apply computational strategies for inferring the cellular composition within each spatial spot. The main objective of this review is to summarize the most recent progresses in estimating the exact cellular proportions for each spatial spot, and to prospect the future directions of this field.

Gene regulatory networks in disease and ageing

The precise control of gene expression is required for the maintenance of cellular homeostasis and proper cellular function, and the declining control of gene expression with age is considered a major contributor to age-associated changes in cellular physiology and disease. The coordination of gene expression can be represented through models of the molecular interactions that govern gene expression levels, so-called gene regulatory networks. Gene regulatory networks can represent interactions that occur through signal transduction, those that involve regulatory transcription factors, or statistical models of gene-gene relationships based on the premise that certain sets of genes tend to be coexpressed across a range of conditions and cell types. Advances in experimental and computational technologies have enabled the inference of these networks on an unprecedented scale and at unprecedented precision. Here, we delineate different types of gene regulatory networks and their cell-biological interpretation. We describe methods for inferring such networks from large-scale, multi-omics datasets and present applications that have aided our understanding of cellular ageing and disease mechanisms.

MICAL1 promotes the proliferation in acute myeloid leukemia and is associated with clinical prognosis and immune infiltration

Acute myeloid leukemia (AML) is one of the most common hematopoietic malignancies that has a poor prognosis and a high rate of relapse. Dysregulated metabolism plays an important role in AML progression. This study aimed to conduct a comprehensive analysis of MRGs using TCGA and GEO datasets and further explore the potential function of critical MRGs in AML progression. In this study, we identified 17 survival-related differentially expressed MRGs in AML using TCGA and GEO datasets. The 150 AML samples were divided into three molecular subtypes using 17 MRGs, and we found that three molecular subtypes exhibited a different association with ferroptosis, cuproptosis and m6A related genes. Moreover, a prognostic signature that comprised nine MRGs and had good predictive capacity was established by LASSO-Cox stepwise regression analysis. Among the 17 MRGs, our attention focused on MICAL1 which was highly expressed in many types of tumors, including AML and its overexpression was also confirmed in several AML cell lines. We also found that the expression of MICAL1 was associated with several immune cells. Moreover, functional experiments revealed that knockdown of MICAL1 distinctly suppressed the proliferation of AML cells. Overall, this study not only contributes to a deeper understanding of the molecular mechanisms underlying AML but also provides potential targets and prognostic markers for AML treatment. These findings offer robust support for further research into therapeutic strategies and mechanisms related to AML, with the potential to improve the prognosis and quality of life for AML patients. Nevertheless, further research is needed to validate these findings and explore more in-depth molecular mechanisms.

Single-cell sequencing, spatial transcriptome ad periodontitis: Rethink pathogenesis and classification

OBJECTIVE

This narrative review illuminates on the application of single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics (ST) in periodontitis and highlights the probability of relating cell population and gene signatures to the pathogenesis of the disease for a better diagnosis.

METHODS

An electronic search of the literature in the PubMed database for the keywords, "single cell sequencing" OR "spatial transcriptomics" and "periodontitis" OR "gingiva" OR "oral mucosa" yielded 486 research articles and reviews. After filtering duplicates and careful curation, 22 papers conducted in humans were retained.

RESULTS

The molecular mechanisms underlying periodontitis are complex and involve the interaction of multiple cells and various gene expressions. Most residing cells in periodontal tissues participate in maintaining homeostasis and health, while in addition to infiltrating immune cells contribute to the fight against the bacterial insult.

CONCLUSION

scRNA-seq and ST have provided new insights into the cellular and molecular changes associated with periodontitis for a better diagnosis and clinical outcome. New functions of cells and genes are revealed with these techniques; however, no cells or gene signatures are attributed to periodontitis so far.

Deciphering brain cellular and behavioral mechanisms: Insights from single-cell and spatial RNA sequencing

The brain is a complex computing system composed of a multitude of interacting neurons. The computational outputs of this system determine the behavior and perception of every individual. Each brain cell expresses thousands of genes that dictate the cell's function and physiological properties. Therefore, deciphering the molecular expression of each cell is of great significance for understanding its characteristics and role in brain function. Additionally, the positional information of each cell can provide crucial insights into their involvement in local brain circuits. In this review, we briefly overview the principles of single-cell RNA sequencing and spatial transcriptomics, the potential issues and challenges in their data processing, and their applications in brain research. We further outline several promising directions in neuroscience that could be integrated with single-cell RNA sequencing, including neurodevelopment, the identification of novel brain microstructures, cognition and behavior, neuronal cell positioning, molecules and cells related to advanced brain functions, sleep-wake cycles/circadian rhythms, and computational modeling of brain function. We believe that the deep integration of these directions with single-cell and spatial RNA sequencing can contribute significantly to understanding the roles of individual cells or cell types in these specific functions, thereby making important contributions to addressing critical questions in those fields. This article is categorized under: RNA Evolution and Genomics > Computational Analyses of RNA RNA in Disease and Development > RNA in Development RNA in Disease and Development > RNA in Disease.

Gene count normalization in single-cell imaging-based spatially resolved transcriptomics

BACKGROUND

Recent advances in imaging-based spatially resolved transcriptomics (im-SRT) technologies now enable high-throughput profiling of targeted genes and their locations in fixed tissues. Normalization of gene expression data is often needed to account for technical factors that may confound underlying biological signals.

RESULTS

Here, we investigate the potential impact of different gene count normalization methods with different targeted gene panels in the analysis and interpretation of im-SRT data. Using different simulated gene panels that overrepresent genes expressed in specific tissue regions or cell types, we demonstrate how normalization methods based on detected gene counts per cell differentially impact normalized gene expression magnitudes in a region- or cell type-specific manner. We show that these normalization-induced effects may reduce the reliability of downstream analyses including differential gene expression, gene fold change, and spatially variable gene analysis, introducing false positive and false negative results when compared to results obtained from gene panels that are more representative of the gene expression of the tissue's component cell types. These effects are not observed with normalization approaches that do not use detected gene counts for gene expression magnitude adjustment, such as with cell volume or cell area normalization.

CONCLUSIONS

We recommend using non-gene count-based normalization approaches when feasible and evaluating gene panel representativeness before using gene count-based normalization methods if necessary. Overall, we caution that the choice of normalization method and gene panel may impact the biological interpretation of the im-SRT data.

Elucidating the molecular landscape of tendinitis: the role of inflammasome-related genes and immune interactions

Tendinitis, characterized by the inflammation of tendons, poses significant challenges in both diagnosis and treatment due to its multifaceted etiology and complex pathophysiology. This study aimed to dissect the molecular mechanisms underlying tendinitis, with a particular focus on inflammasome-related genes and their interactions with the immune system. Through comprehensive gene expression analysis and bioinformatics approaches, we identified distinct expression profiles of inflammasome genes, such as NLRP6, NLRP1, and MEFV, which showed significant correlations with immune checkpoint molecules, indicating a pivotal role in the inflammatory cascade of tendinitis. Additionally, MYD88 and CD36 were found to be closely associated with HLA family molecules, underscoring their involvement in immune response modulation. Contrary to expectations, chemokines exhibited minimal correlation with inflammasome genes, suggesting an unconventional inflammatory pathway in tendinitis. Transcription factors like SP110 and CREB5 emerged as key regulators of inflammasome genes, providing insight into the transcriptional control mechanisms in tendinitis. Furthermore, potential therapeutic targets were identified through the DGidb database, highlighting drugs that could modulate the activity of inflammasome genes, offering new avenues for targeted tendinitis therapy. Our findings elucidate the complex molecular landscape of tendinitis, emphasizing the significant role of inflammasomes and immune interactions, and pave the way for the development of novel diagnostic and therapeutic strategies.

SpotSweeper: spatially-aware quality control for spatial transcriptomics

Quality control (QC) is a crucial step to ensure the reliability and accuracy of the data obtained from RNA sequencing experiments, including spatially-resolved transcriptomics (SRT). Existing QC approaches for SRT that have been adopted from single-nucleus RNA sequencing (snRNA-seq) methods are confounded by spatial biology and are inappropriate for SRT data. In addition, no methods currently exist for identifying histological tissue artifacts unique to SRT. Here, we introduce SpotSweeper, spatially-aware QC methods for identifying local outliers and regional artifacts in SRT. SpotSweeper evaluates the quality of individual spots relative to their local neighborhood, thus minimizing bias due to biological heterogeneity, and uses multiscale methods to detect regional artifacts. Using SpotSweeper on publicly available data, we identified a consistent set of Visium barcodes/spots as systematically low quality and demonstrate that SpotSweeper accurately identifies two distinct types of regional artifacts, resulting in improved downstream clustering and marker gene detection for spatial domains.

Spatial transcriptomics in cancer research and potential clinical impact: a narrative review

Spatial transcriptomics (ST) provides novel insights into the tumor microenvironment (TME). ST allows the quantification and illustration of gene expression profiles in the spatial context of tissues, including both the cancer cells and the microenvironment in which they are found. In cancer research, ST has already provided novel insights into cancer metastasis, prognosis, and immunotherapy responsiveness. The clinical precision oncology application of next-generation sequencing (NGS) and RNA profiling of tumors relies on bulk methods that lack spatial context. The ability to preserve spatial information is now possible, as it allows us to capture tumor heterogeneity and multifocality. In this narrative review, we summarize precision oncology, discuss tumor sequencing in the clinic, and review the available ST research methods, including seqFISH, MERFISH (Vizgen), CosMx SMI (NanoString), Xenium (10x), Visium (10x), Stereo-seq (STOmics), and GeoMx DSP (NanoString). We then review the current ST literature with a focus on solid tumors organized by tumor type. Finally, we conclude by addressing an important question: how will spatial transcriptomics ultimately help patients with cancer?

HE2Gene: image-to-RNA translation via multi-task learning for spatial transcriptomics data

MOTIVATION

Tissue context and molecular profiling are commonly used measures in understanding normal development and disease pathology. In recent years, the development of spatial molecular profiling technologies (e.g. spatial resolved transcriptomics) has enabled the exploration of quantitative links between tissue morphology and gene expression. However, these technologies remain expensive and time-consuming, with subsequent analyses necessitating high-throughput pathological annotations. On the other hand, existing computational tools are limited to predicting only a few dozen to several hundred genes, and the majority of the methods are designed for bulk RNA-seq.

RESULTS

In this context, we propose HE2Gene, the first multi-task learning-based method capable of predicting tens of thousands of spot-level gene expressions along with pathological annotations from H&E-stained images. Experimental results demonstrate that HE2Gene is comparable to state-of-the-art methods and generalizes well on an external dataset without the need for re-training. Moreover, HE2Gene preserves the annotated spatial domains and has the potential to identify biomarkers. This capability facilitates cancer diagnosis and broadens its applicability to investigate gene-disease associations.

AVAILABILITY AND IMPLEMENTATION

The source code and data information has been deposited at https://github.com/Microbiods/HE2Gene.

Enhanced clustering-based differential expression analysis method for RNA-seq data

RNA-seq is a tool for measuring gene expression and is commonly used to identify differentially expressed genes (DEGs). Gene clustering has been widely used to classify DEGs with similar expression patterns, but rarely used to identify DEGs themselves. We recently reported that the clustering-based method (called MBCdeg1 and 2) for identifying DEGs has great potential. However, these methods left room for improvement. This study reports on the improvement (named MBCdeg3). We compared a total of six competing methods: three conventional R packages (edgeR, DESeq2, and TCC) and three versions of MBCdeg (i.e., MBCdeg1, 2, and 3) corresponding to three different normalization algorithms. As MBCdeg3 performs well in many simulation scenarios of RNA-seq count data, MBCdeg3 replaces MBCdeg1 and 2 in our previous report. •MBCdeg3 is a method for both identification and classification of DEGs from RNA-seq count data.•MBCdeg3 is available as a function of R, which is common in the field of expression analysis.•MBCdeg3 performs well in a variety of simulation scenarios for RNA-seq count data.

SpatialPrompt: spatially aware scalable and accurate tool for spot deconvolution and domain identification in spatial transcriptomics

Efficiently mapping of cell types in situ remains a major challenge in spatial transcriptomics. Most spot deconvolution tools ignore spatial coordinate information and perform extremely slow on large datasets. Here, we introduce SpatialPrompt, a spatially aware and scalable tool for spot deconvolution and domain identification. SpatialPrompt integrates gene expression, spatial location, and single-cell RNA sequencing (scRNA-seq) dataset as reference to accurately infer cell-type proportions of spatial spots. SpatialPrompt uses non-negative ridge regression and graph neural network to efficiently capture local microenvironment information. Our extensive benchmarking analysis on Visium, Slide-seq, and MERFISH datasets demonstrated superior performance of SpatialPrompt over 15 existing tools. On mouse hippocampus dataset, SpatialPrompt achieves spot deconvolution and domain identification within 2 minutes for 50,000 spots. Overall, domain identification using SpatialPrompt was 44 to 150 times faster than existing methods. We build a database housing 40 plus curated scRNA-seq datasets for seamless integration with SpatialPrompt for spot deconvolution.

Cellular diversity through space and time: adding new dimensions to GBM therapeutic development

The current median survival for glioblastoma (GBM) patients is only about 16 months, with many patients succumbing to the disease in just a matter of months, making it the most common and aggressive primary brain cancer in adults. This poor outcome is, in part, due to the lack of new treatment options with only one FDA-approved treatment in the last decade. Advances in sequencing techniques and transcriptomic analyses have revealed a vast degree of heterogeneity in GBM, from inter-patient diversity to intra-tumoral cellular variability. These cutting-edge approaches are providing new molecular insights highlighting a critical role for the tumor microenvironment (TME) as a driver of cellular plasticity and phenotypic heterogeneity. With this expanded molecular toolbox, the influence of TME factors, including endogenous (e.g., oxygen and nutrient availability and interactions with non-malignant cells) and iatrogenically induced (e.g., post-therapeutic intervention) stimuli, on tumor cell states can be explored to a greater depth. There exists a critical need for interrogating the temporal and spatial aspects of patient tumors at a high, cell-level resolution to identify therapeutically targetable states, interactions and mechanisms. In this review, we discuss advancements in our understanding of spatiotemporal diversity in GBM with an emphasis on the influence of hypoxia and immune cell interactions on tumor cell heterogeneity. Additionally, we describe specific high-resolution spatially resolved methodologies and their potential to expand the impact of pre-clinical GBM studies. Finally, we highlight clinical attempts at targeting hypoxia- and immune-related mechanisms of malignancy and the potential therapeutic opportunities afforded by single-cell and spatial exploration of GBM patient specimens.

Heterogeneity of immune cells and their communications unveiled by transcriptome profiling in acute inflammatory lung injury

Background

Acute Respiratory Distress Syndrome (ARDS) or its earlier stage Acute lung injury (ALI), is a worldwide health concern that jeopardizes human well-being. Currently, the treatment strategies to mitigate the incidence and mortality of ARDS are severely restricted. This limitation can be attributed, at least in part, to the substantial variations in immunity observed in individuals with this syndrome.

Methods

Bulk and single cell RNA sequencing from ALI mice and single cell RNA sequencing from ARDS patients were analyzed. We utilized the Seurat program package in R and cellmarker 2.0 to cluster and annotate the data. The differential, enrichment, protein interaction, and cell-cell communication analysis were conducted.

Results

The mice with ALI caused by pulmonary and extrapulmonary factors demonstrated differential expression including Clec4e, Retnlg, S100a9, Coro1a, and Lars2. We have determined that inflammatory factors have a greater significance in extrapulmonary ALI, while multiple pathways collaborate in the development of pulmonary ALI. Clustering analysis revealed significant heterogeneity in the relative abundance of immune cells in different ALI models. The autocrine action of neutrophils plays a crucial role in pulmonary ALI. Additionally, there was a significant increase in signaling intensity between B cells and M1 macrophages, NKT cells and M1 macrophages in extrapulmonary ALI. The CXCL, CSF3 and MIF, TGFβ signaling pathways play a vital role in pulmonary and extrapulmonary ALI, respectively. Moreover, the analysis of human single-cell revealed DCs signaling to monocytes and neutrophils in COVID-19-associated ARDS is stronger compared to sepsis-related ARDS. In sepsis-related ARDS, CD8+ T and Th cells exhibit more prominent signaling to B-cell nucleated DCs. Meanwhile, both MIF and CXCL signaling pathways are specific to sepsis-related ARDS.

Conclusion

This study has identified specific gene signatures and signaling pathways in animal models and human samples that facilitate the interaction between immune cells, which could be targeted therapeutically in ARDS patients of various etiologies.

spVC for the detection and interpretation of spatial gene expression variation

Spatially resolved transcriptomics technologies have opened new avenues for understanding gene expression heterogeneity in spatial contexts. However, existing methods for identifying spatially variable genes often focus solely on statistical significance, limiting their ability to capture continuous expression patterns and integrate spot-level covariates. To address these challenges, we introduce spVC, a statistical method based on a generalized Poisson model. spVC seamlessly integrates constant and spatially varying effects of covariates, facilitating comprehensive exploration of gene expression variability and enhancing interpretability. Simulation and real data applications confirm spVC's accuracy in these tasks, highlighting its versatility in spatial transcriptomics analysis.

SPADE: spatial deconvolution for domain specific cell-type estimation

Understanding gene expression in different cell types within their spatial context is a key goal in genomics research. SPADE (SPAtial DEconvolution), our proposed method, addresses this by integrating spatial patterns into the analysis of cell type composition. This approach uses a combination of single-cell RNA sequencing, spatial transcriptomics, and histological data to accurately estimate the proportions of cell types in various locations. Our analyses of synthetic data have demonstrated SPADE's capability to discern cell type-specific spatial patterns effectively. When applied to real-life datasets, SPADE provides insights into cellular dynamics and the composition of tumor tissues. This enhances our comprehension of complex biological systems and aids in exploring cellular diversity. SPADE represents a significant advancement in deciphering spatial gene expression patterns, offering a powerful tool for the detailed investigation of cell types in spatial transcriptomics.

GFAP-isoforms in the nervous system: Understanding the need for diversity

Glial fibrillary acidic protein (GFAP) is an intermediate filament (IF) protein expressed in specific types of glial cells in the nervous system. The expression of GFAP is highly regulated during brain development and in neurological diseases. The presence of distinct GFAP-isoforms in various cell types, developmental stages, and diseases indicates that GFAP (post-)transcriptional regulation has a role in glial cell physiology and pathology. GFAP-isoforms differ in sub-cellular localisation, IF-network assembly properties, and IF-dynamics which results in distinct molecular interactions and mechanical properties of the IF-network. Therefore, GFAP (post-)transcriptional regulation is likely a mechanism by which radial glia, astrocytes, and glioma cells can modulate cellular function.

Navigating the landscapes of spatial transcriptomics: How computational methods guide the way

Spatially resolved transcriptomics has been dramatically transforming biological and medical research in various fields. It enables transcriptome profiling at single-cell, multi-cellular, or sub-cellular resolution, while retaining the information of geometric localizations of cells in complex tissues. The coupling of cell spatial information and its molecular characteristics generates a novel multi-modal high-throughput data source, which poses new challenges for the development of analytical methods for data-mining. Spatial transcriptomic data are often highly complex, noisy, and biased, presenting a series of difficulties, many unresolved, for data analysis and generation of biological insights. In addition, to keep pace with the ever-evolving spatial transcriptomic experimental technologies, the existing analytical theories and tools need to be updated and reformed accordingly. In this review, we provide an overview and discussion of the current computational approaches for mining of spatial transcriptomics data. Future directions and perspectives of methodology design are proposed to stimulate further discussions and advances in new analytical models and algorithms. This article is categorized under: RNA Methods > RNA Analyses in Cells RNA Evolution and Genomics > Computational Analyses of RNA RNA Export and Localization > RNA Localization.

Spatial transcriptomics reveals prognosis-associated cellular heterogeneity in the papillary thyroid carcinoma microenvironment

BACKGROUND

Papillary thyroid carcinoma (PTC) is the most common malignant endocrine tumour, and its incidence and prevalence are increasing considerably. Cellular heterogeneity in the tumour microenvironment is important for PTC prognosis. Spatial transcriptomics is a powerful technique for cellular heterogeneity study.

METHODS

In conjunction with a clinical pathologist identification method, spatial transcriptomics was employed to characterise the spatial location and RNA profiles of PTC-associated cells within the tissue sections. The spatial RNA-clinical signature genes for each cell type were extracted and applied to outlining the distribution regions of specific cells on the entire section. The cellular heterogeneity of each cell type was further revealed by ContourPlot analysis, monocle analysis, trajectory analysis, ligand-receptor analysis and Gene Ontology enrichment analysis.

RESULTS

The spatial distribution region of tumour cells, typical and atypical follicular cells (FCs and AFCs) and immune cells were accurately and comprehensively identified in all five PTC tissue sections. AFCs were identified as a transitional state between FCs and tumour cells, exhibiting a higher resemblance to the latter. Three tumour foci were shared among all patients out of the 13 observed. Notably, tumour foci No. 2 displayed elevated expression levels of genes associated with lower relapse-free survival in PTC patients. We discovered key ligand-receptor interactions, including LAMB3-ITGA2, FN1-ITGA3 and FN1-SDC4, involved in the transition of PTC cells from FCs to AFCs and eventually to tumour cells. High expression of these patterns correlated with reduced relapse-free survival. In the tumour immune microenvironment, reduced interaction between myeloid-derived TGFB1 and TGFBR1 in tumour focus No. 2 contributed to tumourigenesis and increased heterogeneity. The spatial RNA-clinical analysis method developed here revealed prognosis-associated cellular heterogeneity in the PTC microenvironment.

CONCLUSIONS

The occurrence of tumour foci No. 2 and three enhanced ligand-receptor interactions in the AFC area/tumour foci reduced the relapse-free survival of PTC patients, potentially leading to improved prognostic strategies and targeted therapies for PTC patients.

Comparative analysis of multiplexed in situ gene expression profiling technologies

The burgeoning interest in in situ multiplexed gene expression profiling technologies has opened new avenues for understanding cellular behavior and interactions. In this study, we present a comparative benchmark analysis of six in situ gene expression profiling methods, including both commercially available and academically developed methods, using publicly accessible mouse brain datasets. We find that standard sensitivity metrics, such as the number of unique molecules detected per cell, are not directly comparable across datasets due to substantial differences in the incidence of off-target molecular artifacts impacting specificity. To address these challenges, we explored various potential sources of molecular artifacts, developed novel metrics to control for them, and utilized these metrics to evaluate and compare different in situ technologies. Finally, we demonstrate how molecular false positives can seriously confound spatially-aware differential expression analysis, requiring caution in the interpretation of downstream results. Our analysis provides guidance for the selection, processing, and interpretation of in situ spatial technologies.

Spatial multi-omics: novel tools to study the complexity of cardiovascular diseases

Spatial multi-omic studies have emerged as a promising approach to comprehensively analyze cells in tissues, enabling the joint analysis of multiple data modalities like transcriptome, epigenome, proteome, and metabolome in parallel or even the same tissue section. This review focuses on the recent advancements in spatial multi-omics technologies, including novel data modalities and computational approaches. We discuss the advancements in low-resolution and high-resolution spatial multi-omics methods which can resolve up to 10,000 of individual molecules at subcellular level. By applying and integrating these techniques, researchers have recently gained valuable insights into the molecular circuits and mechanisms which govern cell biology along the cardiovascular disease spectrum. We provide an overview of current data analysis approaches, with a focus on data integration of multi-omic datasets, highlighting strengths and weaknesses of various computational pipelines. These tools play a crucial role in analyzing and interpreting spatial multi-omics datasets, facilitating the discovery of new findings, and enhancing translational cardiovascular research. Despite nontrivial challenges, such as the need for standardization of experimental setups, data analysis, and improved computational tools, the application of spatial multi-omics holds tremendous potential in revolutionizing our understanding of human disease processes and the identification of novel biomarkers and therapeutic targets. Exciting opportunities lie ahead for the spatial multi-omics field and will likely contribute to the advancement of personalized medicine for cardiovascular diseases.

Spatial omics: An innovative frontier in aging research

Disentangling the impact of aging on health and disease has become critical as population aging progresses rapidly. Studying aging at the molecular level is complicated by the diverse aging profiles and dynamics. However, the examination of cellular states within aging tissues in situ is hampered by the lack of high-resolution spatial data. Emerging spatial omics technologies facilitate molecular and spatial analysis of tissues, providing direct access to precise information on various functional regions and serving as a favorable tool for unraveling the heterogeneity of aging. In this review, we summarize the recent advances in spatial omics application in multi-organ aging research, which has enhanced the understanding of aging mechanisms from multiple standpoints. We also discuss the main challenges in spatial omics research to date, the opportunities for further developing the technology, and the potential applications of spatial omics in aging and aging-related diseases.

Ursa: A Comprehensive Multiomics Toolbox for High-Throughput Single-Cell Analysis

The burgeoning amount of single-cell data has been accompanied by revolutionary changes to computational methods to map, quantify, and analyze the outputs of these cutting-edge technologies. Many are still unable to reap the benefits of these advancements due to the lack of bioinformatics expertise. To address this issue, we present Ursa, an automated single-cell multiomics R package containing 6 automated single-cell omics and spatial transcriptomics workflows. Ursa allows scientists to carry out post-quantification single or multiomics analyses in genomics, transcriptomics, epigenetics, proteomics, and immunomics at the single-cell level. It serves as a 1-stop analytic solution by providing users with outcomes to quality control assessments, multidimensional analyses such as dimension reduction and clustering, and extended analyses such as pseudotime trajectory and gene-set enrichment analyses. Ursa aims bridge the gap between those with bioinformatics expertise and those without by providing an easy-to-use bioinformatics package for scientists in hoping to accelerate their research potential. Ursa is freely available at https://github.com/singlecellomics/ursa.

Single cell RNA sequencing – a valuable tool for cancer immunotherapy: a mini review

Single-cell RNA sequencing (scRNA-seq) technology has made great strides in research over the last decade. Data analysis has been aided by developments in bioinformatics tools and artificial intelligence, allowing biological and clinical researchers to get a deeper understanding of the different cell clusters and their dynamics within tumours. Combining conventional treatment modalities like chemotherapy and radiation with immunotherapy is a growing trend in cancer treatment. Hence, knowledge of the tumour microenvironment and the effect of each treatment modality on the TME, at a single cell level can provide treating clinicians with better clues for patient stratification and prognostication. With this knowledge, immunotherapy could become successful in treating a wide range of cancers, opening the path for the creation of even more effective treatment strategies. Despite the widespread availability of scRNA-seq technology, computational analysis and data interpretation are still challenges. Worldwide, such challenges are being addressed by various researchers, strengthening the contribution of this technology towards cancer elimination. In this mini-review, we primarily focus on the technique, its workflow, and the computational aspects of scRNA technology, along with an overview of the current challenges in the analysis and interpretation of the data generated.

Spatial Transcriptomics: Emerging Technologies in Tissue Gene Expression Profiling

In this Perspective, we discuss the current status and advances in spatial transcriptomics technologies, which allow high-resolution mapping of gene expression in intact cell and tissue samples. Spatial transcriptomics enables the creation of high-resolution maps of gene expression patterns within their native spatial context, adding an extra layer of information to the bulk sequencing data. Spatial transcriptomics has expanded significantly in recent years and is making a notable impact on a range of fields, including tissue architecture, developmental biology, cancer, and neurodegenerative and infectious diseases. The latest advancements in spatial transcriptomics have resulted in the development of highly multiplexed methods, transcriptomic-wide analysis, and single-cell resolution utilizing diverse technological approaches. In this Perspective, we provide a detailed analysis of the molecular foundations behind the main spatial transcriptomics technologies, including methods based on microdissection, in situ sequencing, single-molecule FISH, spatial capturing, selection of regions of interest, and single-cell or nuclei dissociation. We contextualize the detection and capturing efficiency, strengths, limitations, tissue compatibility, and applications of these techniques as well as provide information on data analysis. In addition, this Perspective discusses future directions and potential applications of spatial transcriptomics, highlighting the importance of the continued development to promote widespread adoption of these techniques within the research community.

Bayesian network models identify co-operative GPCR:G protein interactions that contribute to G protein coupling

Cooperative interactions in protein-protein interfaces demonstrate the interdependency or the linked network-like behavior of interface interactions and their effect on the coupling of proteins. Cooperative interactions also could cause ripple or allosteric effects at a distance in protein-protein interfaces. Although they are critically important in protein-protein interfaces it is challenging to determine which amino acid pair interactions are cooperative. In this work we have used Bayesian network modeling, an interpretable machine learning method, combined with molecular dynamics trajectories to identify the residue pairs that show high cooperativity and their allosteric effect in the interface of G protein-coupled receptor (GPCR) complexes with G proteins. Our results reveal a strong co-dependency in the formation of interface GPCR:G protein contacts. This observation indicates that cooperativity of GPCR:G protein interactions is necessary for the coupling and selectivity of G proteins and is thus critical for receptor function. We have identified subnetworks containing polar and hydrophobic interactions that are common among multiple GPCRs coupling to different G protein subtypes (Gs, Gi and Gq). These common subnetworks along with G protein-specific subnetworks together confer selectivity to the G protein coupling. This work underscores the potential of data-driven Bayesian network modeling in elucidating the intricate dependencies and selectivity determinants in GPCR:G protein complexes, offering valuable insights into the dynamic nature of these essential cellular signaling components.

Promise of spatially resolved omics for tumor research

Tumors are spatially heterogeneous tissues that comprise numerous cell types with intricate structures. By interacting with the microenvironment, tumor cells undergo dynamic changes in gene expression and metabolism, resulting in spatiotemporal variations in their capacity for proliferation and metastasis. In recent years, the rapid development of histological techniques has enabled efficient and high-throughput biomolecule analysis. By preserving location information while obtaining a large number of gene and molecular data, spatially resolved metabolomics (SRM) and spatially resolved transcriptomics (SRT) approaches can offer new ideas and reliable tools for the in-depth study of tumors. This review provides a comprehensive introduction and summary of the fundamental principles and research methods used for SRM and SRT techniques, as well as a review of their applications in cancer-related fields.

Spatial mapping of the DNA adducts in cancer

DNA adducts and strand breaks are induced by various exogenous and endogenous agents. Accumulation of DNA damage is implicated in many disease processes, including cancer, aging, and neurodegeneration. The continuous acquisition of DNA damage from exogenous and endogenous stressors coupled with defects in DNA repair pathways contribute to the accumulation of DNA damage within the genome and genomic instability. While mutational burden offers some insight into the level of DNA damage a cell may have experienced and subsequently repaired, it does not quantify DNA adducts and strand breaks. Mutational burden also infers the identity of the DNA damage. With advances in DNA adduct detection and quantification methods, there is an opportunity to identify DNA adducts driving mutagenesis and correlate with a known exposome. However, most DNA adduct detection methods require isolation or separation of the DNA and its adducts from the context of the nuclei. Mass spectrometry, comet assays, and other techniques precisely quantify lesion types but lose the nuclear context and even tissue context of the DNA damage. The growth in spatial analysis technologies offers a novel opportunity to leverage DNA damage detection with nuclear and tissue context. However, we lack a wealth of techniques capable of detecting DNA damage in situ. Here, we review the limited existing in situ DNA damage detection methods and examine their potential to offer spatial analysis of DNA adducts in tumors or other tissues. We also offer a perspective on the need for spatial analysis of DNA damage in situ and highlight Repair Assisted Damage Detection (RADD) as an in situ DNA adduct technique with the potential to integrate with spatial analysis and the challenges to be addressed.

JUMP: replicability analysis of high-throughput experiments with applications to spatial transcriptomic studies

MOTIVATION

Replicability is the cornerstone of scientific research. The current statistical method for high-dimensional replicability analysis either cannot control the false discovery rate (FDR) or is too conservative.

RESULTS

We propose a statistical method, JUMP, for the high-dimensional replicability analysis of two studies. The input is a high-dimensional paired sequence of p-values from two studies and the test statistic is the maximum of p-values of the pair. JUMP uses four states of the p-value pairs to indicate whether they are null or non-null. Conditional on the hidden states, JUMP computes the cumulative distribution function of the maximum of p-values for each state to conservatively approximate the probability of rejection under the composite null of replicability. JUMP estimates unknown parameters and uses a step-up procedure to control FDR. By incorporating different states of composite null, JUMP achieves a substantial power gain over existing methods while controlling the FDR. Analyzing two pairs of spatially resolved transcriptomic datasets, JUMP makes biological discoveries that otherwise cannot be obtained by using existing methods.

AVAILABILITY AND IMPLEMENTATION

An R package JUMP implementing the JUMP method is available on CRAN (https://CRAN.R-project.org/package=JUMP).