Research Appraisals
Evidence-based critical appraisals of the latest medical research, systematically evaluated using Oxford CEBM methodology.
Showing 16 appraisals
Biological psychiatry
Systematic Review and Mega-Analysis of the Peripheral Blood Transcriptome in Depression Implicates Dysregulation of Lymphoid Cells and Histones
BACKGROUND: Depression has been associated with transcriptomic changes in peripheral blood. However, the contribution of specific immune cell subsets or pathways remains unclear, and findings have been variable across previous studies, which have not tended to account for sample cellular composition. METHODS: We performed a systematic review of peripheral blood transcriptome studies in depression. For the 5 datasets meeting criteria (total N = 6011), we performed harmonized reprocessing and cell composition-adjusted differential gene and transcript analyses, followed by a bias- and inflation-adjusted weighted z score mega-analysis. We investigated the biological pathways and cell subsets implicated by the results. We also performed a sex-stratified gene network mega-analysis using consensus weighted gene coexpression network analysis (WGCNA). RESULTS: Few genes showed robust differential gene expression (DGE) in depression. Depression was reproducibly associated with decreases in replication-dependent histones and with a decrease in oxidative phosphorylation pathways in females only. Cell source analyses implicated lymphoid cells (T cells and natural killer cells) as likely contributors to the depression differential expression signature. WGCNA mega-analysis revealed multiple consensus modules associated with depression, with a PUF60-related module upregulated in both female and male depression in sex-stratified analyses. Two genes predicted to be causally relevant to depression by transcriptome-wide association studies (GPX4 and GYPE) showed significant DGE. CONCLUSIONS: These results are convergent with immunogenetic evidence implicating lymphoid cell dysregulation in depression, while also highlighting histone alterations as a key molecular signature in depression. They also indicate the importance of large-scale datasets for biomarker discovery in the context of heterogeneous disorders such as depression.
16 June 2026
Read appraisal →Cell systems
GFETM: Genome foundation-based embedded topic model for scATAC-seq modeling
Single-cell assay for transposase-accessible chromatin with sequencing (scATAC-seq) enables investigation of open chromatin landscapes at single-cell resolution. However, analyzing scATAC-seq data remains challenging due to inherent sparsity and noise. Genome foundation models (GFMs), pre-trained on extensive DNA sequence datasets, have demonstrated effectiveness in genome analysis. Because open chromatin regions (OCRs) harbor salient sequence features, we hypothesized that leveraging GFMs' sequence embeddings could enhance scATAC-seq modeling accuracy and generalizability. We introduce the genome foundation embedded topic model (GFETM), an interpretable deep learning framework combining GFMs with the embedded topic model (ETM) for scATAC-seq analysis. By integrating DNA sequence embeddings extracted by a GFM from OCRs, GFETM demonstrates superior accuracy and generalizability and captures cell-state-specific transcription factor (TF) activity with both zero-shot inference and attention-mechanism analysis. Finally, the topic mixtures inferred by GFETM reveal biologically meaningful epigenomic signatures of kidney diabetes. A record of this paper's transparent peer review process is included in the supplemental information.
22 May 2026
Read appraisal →Cell systems
Learning multi-cellular representations of single-cell transcriptomics data enables characterization of patient-level disease states
Single-cell RNA sequencing (scRNA-seq) is a prominent tool for studying human disease biology. The availability of massive scRNA-seq datasets and advanced machine learning has driven the development of single-cell foundation models that provide versatile cell representations based on expression. To understand disease states, we need to account for entire tissue ecosystems, simultaneously considering many different interacting cells. Here, we tackle this challenge by generating patient-level representations derived from multi-cellular expression contexts measured with scRNA-seq of tissues. We develop PaSCient, a machine learning model that employs a multi-level representation learning paradigm and provides importance scores at the cell and gene levels for fine-grained analysis of disease characteristics. We apply PaSCient to a large-scale scRNA-seq atlas of 12.5 million cells from over 2,700 patients. Comprehensive benchmarking demonstrates the superiority of PaSCient in disease classification and multiple downstream applications, including dimensionality reduction, feature prioritization, and patient subgroup discovery.
22 May 2026
Read appraisal →Molecular genetics and genomics : MGG
Comparison between a conventional tool and deep learning models for RNA velocity analysis of scRNA-Seq data
Single-cell RNA sequencing (scRNA-Seq) enables analysis of gene expression at single-cell resolution. RNA velocity analysis infers the temporal dynamics of transcriptional states from the relative abundances of spliced/unspliced mRNA quantified via scRNA-Seq. Classical RNA velocity approaches, such as scVelo, implement gene-specific kinetic modeling. Deep learning methods including DeepVelo, VeloVI, LatentVelo, SymVelo, and scTour are based on variational autoencoders (VAEs), which allow to enhance the robustness and accuracy by leveraging nonlinear latent representations. Here, we systematically evaluated the performance of deep learning RNA velocity tools by comparing with the scVelo dynamical model to access the possible advantages of VAE-base methods. For this purpose, public datasets (GSE149689 and GSE203233) were initially processed using a standard scRNA-Seq pipeline. Comparisons among results of selected velocity tools were conducted using cosine similarity of velocity vectors to assess directional concordance, and by mean squared error analysis of trajectory continuity for the deep learning models. Overall, VAE methods produced significant, richer, and more directionally coherent and consistent velocity fields than the classical model. Our findings indicate that deep learning models provide more consistent and biologically plausible cell-state trajectories, although at the expense of higher computational demands and reliance on accurate splicing quantification. Altogether, our results underscore the relevance of VAE-based frameworks to advance RNA velocity analysis while highlighting the need for careful preprocessing.
14 May 2026
Read appraisal →Journal of chemical information and modeling
scFNSA: A Factorized Node-Set Attentive Framework for Single-Cell Multi-Omics Integration
Single-cell multi-omics technologies enable simultaneous interrogation of transcriptional and epigenomic states at single-cell resolution, providing powerful means to dissect cellular heterogeneity and regulatory mechanisms. However, effective integration of single-cell multi-omics data remains challenging due to extreme sparsity, high dimensionality, and cross-modal heterogeneity, which often lead to distorted similarity structures, information loss, and limited biological interpretability in existing methods. Here, we propose scFNSA, a factorized node-set attentive framework for integrative analysis of single-cell multi-omics data that explicitly decouples modality-specific feature learning from graph-based relational modeling. scFNSA first learns aligned low-dimensional representations for individual modalities through multi-view variational inference and subsequently models complex cellular dependencies via attention-guided graph learning with structured node masking. This design enables robust cross-modal structural alignment while effectively mitigating noise propagation and oversmoothing commonly encountered in deep graph neural networks. Across multiple single-cell multi-omics datasets, scFNSA consistently improves integrative representation quality and cell type resolution compared with state-of-the-art approaches. By decoupling feature extraction from relational inference, scFNSA provides a robust and interpretable framework for single-cell multi-omics integration, facilitating more accurate characterization of cellular states and underlying regulatory landscapes.
13 May 2026
Read appraisal →Genome research
Incorporating valuable prior knowledge to improve deep learning prediction of genetic perturbation responses
Genetic perturbation response prediction plays a critical role in virtual cell research, yet the performance of current deep learning models still leaves room for improvement. In this study, we present a prior-guided response inference model (PRIM) that leverages a valuable priori knowledge of gene expression in control cells to model the effects of perturbations at the cellular level. This allows PRIM to predict the amount of change in gene expression after perturbation to better approximate the real situation. Compared with existing deep learning approaches, PRIM achieves superior performance across multiple data sets and notable advantages in predicting combinatorial perturbation responses. Moreover, it is more lightweight than current deep learning models and enables faster forward inference. Importantly, PRIM effectively captures nonadditive genetic interactions and shows the potential to uncover associations between combinatorial perturbations and new biologically meaningful phenotypes. These findings provide new insights into the application of deep learning for predicting cellular responses to genetic perturbations.
11 May 2026
Read appraisal →Nature reviews. Genetics
Interpretation, extrapolation and perturbation of single cells.
Single-cell analyses have transitioned from descriptive atlasing towards inferring causal effects and mechanistic relationships that capture cellular logic. Technological advances and the growing scale of observational and interventional datasets have fuelled the development of machine learning methods aimed at identifying such dependencies and extrapolating perturbation effects. Here, we review and connect these approaches according to their modelling concepts (including representation learning, causal inference, mechanistic discovery, disentanglement and population tracing), underlying assumptions and downstream tasks. We propose a unifying ontology to guide practitioners in selecting the most suitable methods for a given biological question, with detailed technical descriptions provided in an online resource . Finally, we identify promising computational directions and underexplored data properties that could pave the way for future developments.
4 May 2026
Read appraisal →Scientific reports
scDecorr: feature decorrelation based representation learning enables self-supervised alignment of multiple single-cell experiments
Single-cell RNA sequencing (scRNA-seq) has revolutionized our understanding of cellular heterogeneity in complex biological systems. However, analyzing and integrating scRNA-seq data poses unique computational challenges due to sparsity, high variability, and technical batch effects. Here, we propose a novel framework called scDecorr for robust representation learning and data integration for scRNA-seq analysis. Our approach leverages the idea of feature decorrelation-based self-supervised learning (SSL) to obtain efficient low-dimensional representations of individual cells without relying on cell-type annotations. By maximizing similarity among distorted embeddings while decorrelating their components, scDecorr captures the biological signature while eliminating technical noise. Furthermore, scDecorr incorporates unsupervised domain adaptation to bridge the gap between batches with different distributions, enabling effective integration of scRNA-seq data from diverse sources. Our framework achieves domain-invariant representations by learning cell embeddings independently across domains and employing domain-specific batch normalization. We evaluate scDecorr on a variety of single-cell datasets and demonstrate its ability to integrate batches without losing the inherent biological variance, thereby facilitating optimal clustering. The representations generated by scDecorr also exhibit robustness in label transfer tasks, allowing for effective transfer of cell-type labels from reference to query datasets. Overall, scDecorr offers a powerful tool for efficient analysis and integration of large and complex scRNA-seq datasets, advancing our understanding of cellular processes and disease mechanisms. The code is available here https://github.com/hayatlab/scdecorr .
3 May 2026
Read appraisal →Analytical chemistry
scAClc: A Multi-Objective Adaptive Clustering Framework for Single-Cell Transcriptomics via Contrastive and Resolution-Aware Representation Learning
Single-cell RNA sequencing (scRNA-seq) enables whole-transcriptomic profiling at single-cell resolution, facilitating the construction of virtual cell representations that capture the full spectrum of cellular identities. Realizing this goal hinges on accurate clustering, which remains challenging due to data sparsity, high dimensionality, noise, and the need to specify cluster numbers a priori. We propose scAClc, a novel clustering framework featuring multiobjective optimization and adaptive resolution discovery, designed to address these limitations through three key innovations. First, a Hierarchical Gene Relevance Module integrates global gene variability with local neighborhood-specific signals to eliminate redundancy while retaining biologically informative features. Second, an Anchor-Centered Contrastive Learning Module adaptively selects representative anchors to guide embedding learning, promoting compact intracluster structure and clear intercluster separation. Third, based on the robust low-dimensional embedding, we propose a Self-Adaptive Resolution Discovery Module to automatically infer the number of clusters by jointly modeling intra- and intercluster distances. Extensive experiments on 15 real scRNA-seq data sets demonstrate that scAClc consistently outperforms six state-of-the-art methods across multiple evaluation metrics. Ablation studies further confirmed the complementary contributions of each module. In addition, interpretability analysis effectively mitigates the "black box" nature of clustering models and sheds light on the biological mechanisms underlying cell clustering. The source code is publicly available at https://github.com/scAClc/scAClc.
23 Apr 2026
Read appraisal →Cell systems
Context-informed subgraph foundation models enable interpretable protein-function prediction
Protein-function prediction is crucial for elucidating molecular mechanisms driving biological processes and therapeutics development. Despite numerous computational tools demonstrating promising performance, they fall short when predicting rare, uncharacterized functions or indirect activities. Here, we present COSMOS, a context-aware Gene Ontology (GO) subgraph mining system for protein-function prediction. By leveraging inductive subgraph foundation models and an enriched knowledge graph of protein-GO relationships, COSMOS performs zero-shot, few-shot, and low-homology protein-function prediction. Built on 7,923,952 functional semantic relationships, COSMOS demonstrates robust capabilities to (1) generate state-of-the-art predictions for GO classes with sparse or no experimental annotations, (2) provide interpretable functional subgraphs for transparent rationale analysis, and (3) deliver complementary benefits when integrated with existing embedding-based prediction methods. We anticipate that COSMOS will serve as a complementary approach to conventional protein annotation methods and an interpretable tool for predicting protein functions within underexplored GO classes, thereby advancing genomics and therapeutic research.
18 Apr 2026
Read appraisal →Proceedings of the National Academy of Sciences of the United States of America
Coalescence and translation: A language model for population genetics.
Probabilistic models such as the sequentially Markovian coalescent have long provided a powerful framework for population genetic inference, enabling reconstruction of demographic history and ancestral relationships from genomic data. However, these methods are inherently specialized, relying on predefined assumptions and/or limited scalability. Recent advances in simulation and deep learning provide an alternative approach: learning directly to generalize from synthetic genetic data to infer hidden evolutionary processes. Here we reframe the inference of coalescence times as a problem of translation between two biological languages: the sparse, observable patterns of mutation along the genome and the unobservable ancestral recombination graph that gave rise to them. Inspired by large language models, we develop cxt, a decoder-only transformer that autoregressively predicts coalescent events conditioned on local mutational context. We show that cxt performs competitively with state-of-the-art Markov Chain Monte Carlo-based likelihood models across a broad range of demographic scenarios, matching their accuracy in-distribution and approaching it in out-of-distribution, with the potential for improvement via fine-tuning. Trained on simulations spanning the stdpopsim catalog [J. R. Adrion et al., eLife 9, e54967 (2020); M. E. Lauterbur et al., eLife 12, e84874 (2023); G. Gower et al., Accessible, realistic genome simulation with selection using stdpopsim. bioRxiv [Preprint] (2025). https://doi.org/10.1101/2025.03.23.644823.], the model generalizes robustly and enables efficient inference at scale, producing over a million coalescence predictions in minutes. cxt also produces well-calibrated approximate posteriors, enabling principled uncertainty quantification. We apply cxt to population genomic data from both humans and mosquitoes, highlighting the model's ability to deal with the complexities of empirical data.
17 Apr 2026
Read appraisal →Physical chemistry chemical physics : PCCP
TransG4: an interpretable deep-learning approach for sequence-based G-quadruplex prediction
G-quadruplexes (G4) are non-canonical nucleic acid secondary structures formed in guanine-rich regions and have been shown to regulate diverse cellular processes such as gene expression, DNA replication, and telomere maintenance, with increasing evidence linking G4 to cancer and other human diseases. G4 predominantly emerge in guanine-rich regions and are implicated in a spectrum of molecular interactions and disease phenotypes, thus researchers are interested in the formation of G4. However, predicting the formation of G4 from nucleotide sequences is a persistent problem. Existing computational tools for G4 prediction are either rule-based on domain knowledge or rely on a single neural network model like a convolutional neural network (CNN), which lacks interpretability and struggles to capture long-range dependencies among bases. Here, we introduce TransG4, a novel neural network architecture that integrates a CNN, a transformer, and bidirectional gated recurrent units (BiGRUs) to identify potential G4 structures. TransG4 demonstrates strong predictive performance on both G4-seq and rG4-seq datasets, accurately predicting DNA mismatch scores and consistently outperforming existing methods in RNA RSR-ratio prediction. Attention-based interpretations further show that TransG4 captures biologically meaningful motifs consistent with canonical G4 structures, providing an interpretable and generalizable framework and representing a novel and impactful contribution to sequence-based G4 propensity prediction.
17 Apr 2026
Read appraisal →PloS one
Exploration of schizophrenia-associated gene modules using graph theory, co-expression networks, and dimensionality reduction.
The PsychENCODE consortium has generated a comprehensive multi-omic dataset from human brain samples, spanning healthy individuals and those with neuropsychiatric disorders. In this study, we focus exclusively on schizophrenia and analyze PsychENCODE transcriptomic data to construct gene regulatory and co-expression networks, aiming to uncover biologically and clinically relevant gene modules. We apply three independent analytical approaches to the same dataset. First, we preprocess the data using Surrogate Variable Analysis (SVA) to correct for latent variation and normalize gene expression. For graph-based analysis, we use Pearson correlation and the igraph package to construct gene co-expression networks and apply Prim's algorithm to generate minimum spanning trees (MST) and compute centrality measures. Next, we implement Weighted Gene Co-expression Network Analysis (WGCNA), assuming a scale-free topology, to identify modules associated with schizophrenia traits. Dimensionality reduction is performed using Principal Component Analysis (PCA), with visualization aided by t-distributed Stochastic Neighbor Embedding (t-SNE) and Multidimensional Scaling (MDS). Functional enrichment is carried out using Gene Ontology (GO) and KEGG pathway databases. Each method reveals distinct yet complementary biological signatures associated with schizophrenia. The igraph-based approach highlights differentially expressed genes (DEGs) with high centrality, uncovering hub genes involved in mRNA export, synaptic signaling, and ion transport. WGCNA identifies co-expression modules strongly correlated with schizophrenia diagnosis, enriched for immune response, histone modification, and mRNA surveillance pathways. PCA isolates key genes contributing to diagnostic variance, with enrichment in neurotransmitter release cycles and cytokine signaling. Collectively, these results underscore the involvement of immune, synaptic, and epigenetic processes in schizophrenia and demonstrate the power of using multiple, orthogonal analytical lenses.
17 Apr 2026
Read appraisal →Journal of chemical information and modeling
ncProFormer: A CNN-enhanced Transformer for ncRNA Coding-Potential Prediction
Recent studies have revealed that some noncoding RNAs (ncRNAs) bear translational potential, and their encoded micropeptides have essential functions in multiple biological processes. However, accurate identification of coding-capable ncRNAs remains challenging due to weak translation signals, low conservation, and heterogeneous data distributions. Herein, we propose ncProFormer, a deep learning framework tailored for ncRNA coding-potential prediction. ncProFormer integrates the nucleic-acid language model GENA-LM to obtain contextual sequence embeddings, adopts an all-token representation strategy, and employs a convolutional neural network (CNN)-enhanced transformer encoder to jointly capture local nucleotide patterns and long-range dependencies. ncProFormer consistently outperformed the existing methods across the in-house human data set, the external validation data set, the public CPPred benchmark data set. More importantly, this study presents the first cross-species evaluation in ncRNA coding-potential prediction. Without retraining, ncProFormer maintained its strong predictive performance on mouse and rat data sets, showing that the learned biological representations are transferable and it is robust under the distributional shift and cross-species conditions. Collectively, these findings establish ncProFormer as an effective and generalizable framework for uncovering the coding potential of ncRNAs, thus offering a promising computational tool for characterizing ncRNA functions across diverse transcriptomic contexts.
14 Apr 2026
Read appraisal →Genome research
spRefine denoises and imputes spatial transcriptomic data with a reference-free framework powered by genomic language model
The analysis of spatial transcriptomic data is hindered by high noise levels and missing gene measurements, challenges that are further compounded by the higher cost of spatial data compared to traditional single-cell data. To overcome this challenge, we introduce spRefine, a deep learning framework that leverages genomic language models to jointly denoise and impute spatial transcriptomic data. Our results demonstrate that spRefine yields more robust cell- and spot-level representations after denoising and imputation, substantially improving data integration. In addition, spRefine serves as a strong framework for model pretraining and the discovery of novel biological signals, as highlighted by multiple downstream applications across data sets of varying scales. Notably, spRefine enhances the accuracy of spatial aging clock estimations and uncovers new aging-related relationships associated with key biological processes, such as neuronal function loss, which offers new insights for analyzing aging effect with spatial transcriptomics.
8 Apr 2026
Read appraisal →Genes (Basel)
The Evolving Role of Artificial Intelligence in Medical Genetics: Advancing Healthcare, Research, and Biosafety Management
The integration of artificial intelligence (AI) with medical genetics is transforming healthcare by addressing the analytical challenges posed by the vast complexity of multi-omics data. This review explores the synergistic convergence of these fields, highlighting AI's transformative role in enhancing diagnostic precision, enabling non-invasive molecular profiling through imaging-genetics, and advancing predictive and personalized medicine via polygenic risk scores and pharmacogenomics. AI is also emerging as a powerful generative tool in therapeutic design, accelerating drug discovery, protein engineering, and precision gene editing. However, this powerful synergy introduces significant ethical, regulatory, and biosecurity challenges, including data privacy, algorithmic bias, and the dual-use risks of AI-enabled genetic engineering. The future envisions a responsible co-evolution, with multimodal AI and the concept of the Digital Twin driving precision medicine, underpinned by interdisciplinary collaboration to ensure fairness, transparency, and societal trust. This article charts the current landscape and proposes actionable directions, emphasizing the need for robust governance to harness AI's potential while mitigating its risks for the benefit of human health.
1 Dec 2025
Read appraisal →