Search PubMedSearch

SEARCH · Search PubMed

Results for “single cell-type gene expression”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

Perplexity as a Metric for Isoform Diversity in the Human Transcriptome.

Long-read sequencing (LRS) has revealed a far greater diversity of RNA isoforms than earlier technologies, increasing the critical need to determine which, and how many, isoforms per gene are biologically meaningful. To define the space of relevant isoforms from LRS, many existing analysis pipelines rely on arbitrary expression cutoffs, but a single threshold cannot accommodate the broad variability in isoform complexity across genes, cell-types, and disease states captured by LRS. To address this, we propose using perplexity-an interpretable measure derived from entropy-that quantifies the effective number of isoforms per gene based on the full, unfiltered isoform ratio distribution. Calculating perplexity for 124 ENCODE4 PacBio LRS datasets spanning 55 human cell types, we show that it provides intuitive assessments of isoform diversity and captures uncertainty across genes with varying complexity. Perplexity can be calculated at multiple gene regulatory levels-from transcript to protein-to compare how isoform diversity is reduced across stages of gene expression. On average, genes have an ORF-level perplexity of 2.1, indicating production of two distinct protein isoforms. We extended this analysis to evaluate expression variation across tissues and identified 4,593 ORFs across 3,102 genes with moderate to extreme tissue-specificity. We propose perplexity as a consistent, quantitative metric for interpreting isoform diversity across genes, cell types, and disease states. All results are compiled into a community resource to enable cross-study comparisons of novel isoforms.

Journal Article

Exploring Endoplasmic Reticulum Stress-Related Genes in Cartilage Defects: Implications for Diagnosis and Therapy.

INTRODUCTION: Cartilage defects (CDs) are orthopedic conditions with limited regenerative potential. This study aimed to identify endoplasmic reticulum (ER) stress-related biomarkers and construct a diagnostic model to enhance the early detection of CD. METHODS: This study analyzed the transcriptomic dataset GSE129147 to identify ER stressrelated differentially expressed genes (ERSRDEGs) between CD and control tissues using the limma package (version 3.58.1). Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analyses were employed for functional enrichment. Immune infiltration was assessed using cell-type identification, which involved estimating the relative subsets of RNA transcripts and single-sample gene set enrichment analysis. Diagnostic models were constructed using logistic regression, support vector machine, and least absolute shrinkage and selection operator regression. RESULTS: Twenty ERSRDEGs were identified, with CYBB, ATP6V1A, and TNFRSF12A significantly upregulated in CD samples. GO and KEGG analyses highlighted oxidative stress response and extracellular matrix remodeling as key mechanisms in CD pathogenesis. Immune analysis revealed an increase in regulatory T cells and a reduction in CD8. T cells. TNFRSF12A showed strong immune associations and, together with TWIST1 and ATP6V1A, formed the final preliminary diagnostic model. The preliminary LASSO model achieved satisfactory predictive accuracy (AUC: 0.7-0.9). DISCUSSION: These findings suggest that ER stress and immune imbalance jointly contribute to cartilage degeneration. The identified genes, particularly TNFRSF12A, TWIST1, and ATP6V1A, not only serve as potential biomarkers but also provide preliminary evidence for new mechanistic insights into stress-immune crosstalk in CD. CONCLUSION: This study reveals the key roles of ER stress and immune dysregulation in CDs. Moreover, the ERSRDEG-based diagnostic model provides preliminary bioinformatics evidence and potential molecular indicators for targeted diagnostics and therapies.

Humans

A contextual activity score (CAS) for inferring ADAR-associated transcriptional activity across RNA-seq, single-cell, and spatial transcriptomics.

BACKGROUND AND OBJECTIVE: Adenosine-to-inosine RNA editing, catalyzed by Adenosine Deaminases Acting on RNA (ADARs), is a widespread modification involved in neural function, immune regulation, and cancer. The Alu Editing Index (AEI) is the standard metric to estimate ADAR activity but requires raw sequencing reads and is poorly suited for single-cell and spatial transcriptomic data. This study aimed to develop an alternative framework for inferring ADAR-associated transcriptional activity from gene expression data across diverse transcriptomic technologies. METHODS: We developed the Contextual Activity Score (CAS), a framework based on transcriptional signatures from ADAR perturbation experiments. Context-specific signatures were generated for human neurons, mouse neurons, and cancer models to infer ADAR1 and ADAR2 activity. CAS was computed from normalized gene expression matrices using regulon-based enrichment analysis. Performance was evaluated by comparing with the Alu Editing Index across bulk RNA sequencing datasets, simulated sequencing depths, and library preparation protocols. RESULTS: CAS showed strong concordance with the Alu Editing Index across multiple datasets, while remaining robust to reduced sequencing depth and different library protocols. Unlike the Alu Editing Index, CAS can be applied to single-cell and spatial transcriptomic data and enables the independent assessment of ADAR2 activity. In cancer and neuronal contexts, CAS captured biologically meaningful variations in ADAR-associated transcriptional activity at sample, cell-type, and spatial levels. CONCLUSION: CAS provides a scalable approach applicable across multiple RNA-seq protocols for estimating ADAR-associated transcriptional activity using gene expression data. This method, implemented in an open-source R package for broad adoption, expands the ability to study ADAR-associated transcriptional activity across transcriptomic modalities where direct editing quantification is challenging, such as single-cell and spatial transcriptomics.

Adenosine Deaminase

Early oligodendrocyte dysfunction signature in Alzheimer's disease: Insights from DNA methylomics and transcriptomics.

Much research into the aetiology of Alzheimer's disease (AD) has focused on neuronal cell types, while studies on the contribution of glial cells, particularly oligodendrocytes (OLGs), are only starting to emerge. Altered brain DNA methylation, an epigenetic modification that provides the interplay between genetics and environmental cues to tightly regulate gene expression, is well documented in AD. Yet, cell-type-specific investigations remain limited. Here, we examine the role of DNA methylation and OLGs in AD, and how such changes may impact gene expression. We performed weighted-gene correlation network analysis (WGCNA) on multiple brain omics AD datasets across species: human DNA methylation data from 4 brain regions, human brain single-nuclei RNA sequencing data and mouse brain RNA sequencing data. We compared AD-associated network modules enriched for OLG genes across AD brain regions, as well as with other neurodegenerative disease DNA methylation datasets. We identified a DNA methylation signature associated with AD, enriched for OLGs, and preserved across brain regions representing early and late AD pathology stages. Genes within this signature showed altered expression in AD OLGs, confirming cell-type specificity and relevance to AD. This OLG signature was also preserved in transgenic mice with early Aβ pathology and in other neurodegenerative diseases without Aβ pathology. We reveal a consistent pattern of OLG dysfunction spanning early to late stages of AD, across DNA methylation and gene expression. Our findings highlight OLG-associated DNA methylation changes as important in AD pathogenesis, and possibly in other neurodegenerative diseases, opening new avenues for therapeutic development.

Alzheimer Disease

Sparse deconvolution of cell type medleys in spatial transcriptomics.

Mapping cell distributions across spatial locations with whole-genome coverage is essential for understanding cellular responses and signaling However, current deconvolution models aim to estimate the proportions of distinct cell types in each spatial transcriptomics spot by integrating reference single-cell data. These models often assume strong overlap between the reference and spatial datasets, neglecting biology-grounded constraints such as sparsity and cell-type variations, as well as technical sparsity. As a result, these methods rely on over-permissive algorithms that ignore given constraints leading to inaccurate predictions, particularly in heterogeneous or unmatched datasets. We introduce Weight-Induced Sparse Regression (WISpR), a machine learning algorithm that integrates spot-specific hyperparameters and sparsity-driven modeling. Unlike conventional approaches that neglect biology-grounded constraints, WISpR accurately predicts cell-type distributions while preserving biological coherence, i.e., spatially and functionally consistent cell-type localization, even in unmatched datasets. Benchmarking against five alternative methods across ten datasets, WISpR consistently outperformed competitors and predicted cellular landscapes in both normal and cancerous tissues. By leveraging sparse cell-type arrangements, WISpR provides biologically informed, high-resolution cellular maps. Its ability to decode tissue organization in both healthy and diseased states highlights WISpR's practical utility for spatial transcriptomics, particularly in challenging settings involving noise, sparsity, or reference mismatches.

Humans

Combined somatic mutation and transcriptome analysis reveals region-specific differences in clonal architecture in human cortex.

The human cerebral cortex is specialized into regions, but little is known about how human cellular lineages shape cortical regional variation and neuronal cell-type distribution during development. Here, we map single-cell lineages of human cortical regions and neuronal subtypes using >1,000 somatic single-nucleotide variants (sSNVs) identified from deep bulk whole-genome sequencing and analyzed over 25 regions and >72,000 single cells. In the fronto-parietal cortex, sSNVs are rarely restricted, marking neuron-generating clones that disperse into neighboring regions. In contrast, the primary visual cortex harbors 30%-70% more sSNVs than the neighboring secondary visual cortex. Clones at this border exhibit more restricted dispersion, suggesting late developmental lineage segregation. Single-nucleus sSNV and whole-transcriptome analysis reveal glutamatergic neuron clones with modest regional restrictions that share low-mosaic sSNVs with some GABAergic neurons, suggesting a recent dorsal cortical progenitor. Our analysis reveals human-specific cortical lineage patterns, regional differences in clonal patterns, and late divergence of some glutamatergic/GABAergic lineages.

Humans

Shared genetic architecture and neurobiological pathways of problematic alcohol use and anxiety disorders.

Problematic alcohol use (PAU) and anxiety disorders (ANX) frequently co-occur, implying shared genetic and neurobiological foundations. However, the directionality of potential causal relationships and the specific mechanisms underlying the overlap remain unclear. Thus, we investigated the shared genetic architecture and neurobiological pathways between PAU and ANX using a multimethod genomic approach. We analyzed summary statistics from genome-wide association studies (GWAS) of PAU and ANX using Mendelian Randomization to assess causal associations between ANX and PAU. We used MiXeR to assess the overall shared genomic architecture, Local Analysis of (co)Variant Association to estimate regional genetic correlations, and conjunctional false discovery rate (conjFDR) to identify individual overlapping loci. We used FUMA to map single-nucleotide polymorphisms (SNPs) to independent loci, conduct differential gene expression analyses across 30 general and 54 specific tissue types, and perform cell-type specificity analyses using a human brain cell atlas. Druggability of identified targets was also evaluated. Mendelian Randomization analyses indicated bidirectional causal associations between ANX and PAU. MiXeR identified moderate polygenic overlap (52.5%) and genetic correlation (rg = 0.44) between the traits, with high effect direction concordance among shared estimated causal variants (86.4%). ConjFDR identified 97 shared lead SNPs, of which 89 had concordant and 8 discordant effects on PAU and ANX. These loci mapped to 97 genes, including DRD2 and PDE4B, genes linked to dopaminergic and cAMP signaling pathways, respectively. Concordant gene expression was enriched in brain, nerve, adrenal gland, esophagus, stomach, and colon, with enriched expression specifically in the prefrontal cortex, anterior cingulate cortex, hippocampus, hypothalamus, substantia nigra and amygdala. FUMA cell-type enrichment analysis identified associations predominantly in neurons from the cerebral cortex, hippocampus, and thalamus. We found substantial genetic and neurobiological overlap between PAU and ANX, highlighting reciprocal, causal relationships between the traits, with differentially expressed genes enriched in addiction- and anxiety-relevant brain regions. These findings support shared genetic and neurobiological mechanisms linking PAU and ANX, while acknowledging that some signals may reflect broader internalizing or psychiatric liability.

Journal Article

Comparative cellular analysis of motor cortex in human, marmoset and mouse.

The primary motor cortex (M1) is essential for voluntary fine-motor control and is functionally conserved across mammals1. Here, using high-throughput transcriptomic and epigenomic profiling of more than 450,000 single nuclei in humans, marmoset monkeys and mice, we demonstrate a broadly conserved cellular makeup of this region, with similarities that mirror evolutionary distance and are consistent between the transcriptome and epigenome. The core conserved molecular identities of neuronal and non-neuronal cell types allow us to generate a cross-species consensus classification of cell types, and to infer conserved properties of cell types across species. Despite the overall conservation, however, many species-dependent specializations are apparent, including differences in cell-type proportions, gene expression, DNA methylation and chromatin state. Few cell-type marker genes are conserved across species, revealing a short list of candidate genes and regulatory mechanisms that are responsible for conserved features of homologous cell types, such as the GABAergic chandelier cells. This consensus transcriptomic classification allows us to use patch-seq (a combination of whole-cell patch-clamp recordings, RNA sequencing and morphological characterization) to identify corticospinal Betz cells from layer 5 in non-human primates and humans, and to characterize their highly specialized physiology and anatomy. These findings highlight the robust molecular underpinnings of cell-type diversity in M1 across mammals, and point to the genes and regulatory pathways responsible for the functional identity of cell types and their species-specific adaptations.

Animals

Cell-type signatures of Alzheimer's disease shared across population groups.

Genomic studies at single-cell resolution have identified several cell types associated with clinical and pathological traits in Alzheimer's disease1-9, but have not examined associations that are shared across populations. To bridge this gap, here we use single-nucleus RNA sequencing and assay for transposase-accessible chromatin with sequencing to profile cortical and subcortical regions in post-mortem brain-tissue samples from Latin, white (excluding Latin) and African American (excluding Latin) individuals. Using discrete and continuous dissections of molecular programs, we identify cell-type-specific clusters associated with Alzheimer's disease in a region-specific manner across all three population groups, including microglial (GPNMB+ and CD74+ subgroups), astrocytic (SERPINH1+, CD44+ and WIF1+ subgroups) and neuronal (SST+ GABAergic and superficial-layer glutamatergic) signatures. We also report continuous gene-expression factors in astrocytes and oligodendrocytes that are not captured by discrete cluster assignments, but which show strong associations with disease phenotypes; these factors are enriched for genes associated with annotated functions such as lipid processing and neurotransmitter reuptake. Finally, we find that molecular programs reveal six distinct subgroups of individuals with cognitive impairment that span all three populations, are not captured by neuropathology, and are instead distinguished by molecular signatures that are not universally present but are nonetheless associated with ante-mortem impairment. Overall, our study identifies key cell types and gene programs implicated in Alzheimer's disease that are shared across population groups, and underscores how representative sampling can capture both shared signatures and disease heterogeneity, thereby enabling better prioritization of key cell types for further investigation.

Female

A single-cell atlas of multiple myeloma defines malignant archetypes and proliferative states.

Multiple myeloma (MM) is a plasma-cell malignancy with extensive genomic and transcriptional heterogeneity, limiting disease classification and precision therapy. Here we generated a clinically annotated, population-scale, single-cell atlas of MM from 341 individuals spanning the disease and treatment continuum. We identified five recurrent malignant transcriptional archetypes and an orthogonal proliferative program associated with genomic features, therapeutic resistance and clinical outcomes. Validation in the independent CoMMpass cohort demonstrated robustness, prognostic relevance and portability across platforms. We developed a single-cell, target-discovery pipeline prioritizing malignant enrichment, cell-type specificity and tissue restriction, identifying FCRL2 as a plasma-restricted or B cell-lineage-restricted surface target expressed by malignant plasma cells. FCRL2-targeted chimeric antigen receptor T cells demonstrated antigen-specific activity in vitro and survival benefit in vivo. Together, these data provide a clinically actionable blueprint for patient stratification and precision target nomination in plasma-cell malignancies.

Multiple Myeloma

Spatial isoform sequencing at single-cell resolution reveals cell-type-specific spatial isoform variability in multiple brain cell types.

Spatial long-read technologies are increasingly common but usually lack single-cell resolution. This leaves unanswered whether spatially variable isoforms reflect variability within one cell type or differences in region-specific cell-type composition. Here, we developed Spl-ISO-Seq2 (500-nm resolution) and accompanying software, Spl-IsoQuant-2 and Spl-IsoFind, enabling long-read sequencing of >450 million barcodes versus 80,000 previously. Applying this to the adult mouse brain, we compared differential isoform abundance between known regions and spatial isoform patterns independent of predefined regions. Both identified overlapping hits, for example, Rps24 in oligodendrocytes. For known Snap25 spatial isoform variation, we show that it occurs in excitatory neurons. The region-agnostic approach also uncovered patterns missed by region-based comparisons, for example, for Ighm. Notably, many spatial isoform signals are not driven by cell-type composition alone. Finally, our software is applicable to many spatial and single-cell protocols, demonstrating reproducibility between platforms (for example, Visium HD/Stereo-seq). Overall, our experimental/analytical methods enable a submicron-resolution-isoform view and open avenues for spatial isoform disease research.

Animals

Decoding regional keratinization in human oral mucosa through high-resolution spatial transcriptomics.

Oral mucosa exhibits region-specific keratinization, essential for periodontal health, yet the spatial and molecular mechanisms driving these differences remain poorly understood. This study aimed to generate a high-resolution spatial transcriptomic atlas of the human oral mucosa around the mucogingival junction, to reveal stromal-epithelial interactions, that distinguish keratinized from non-keratinized programs. Formalin-fixed paraffin-embedded specimens from the mucogingival junction area of two healthy donors were analyzed with the 10 × Genomics Visium HD platform, yielding two keratinized and two non-keratinized regions. Spatial clustering, pseudotime trajectory inference, cell-type integration with a single-cell reference, and ligand-receptor network analysis were applied to delineate epithelial and stromal compartments. Sixteen reproducible clusters, recapitulating tissue architecture, were identified and revealed distinct transcriptional signatures, distinguishing gingiva from lining mucosa. Pseudotime analysis revealed bifurcating epithelial lineages, originating from a shared basal progenitor layer into keratinized and non-keratinized programs. Gingival keratinization was driven by stromal collagen ligands (COL1A1, COL1A2, COL6A1, COL6A2) engaging epithelial receptors (CD44, SDC1), further reinforced within the epithelium by desmosomal adhesion via DSG1-DSC2/3. Gingival keratinization emerges from integrated stromal collagen signaling and epithelial adhesion. This spatially resolved framework advances understanding of oral mucosal specialization and provides a foundation for biologically guided regenerative therapies.

Humans

CeLLTra: aligning cell names with gene expression via a pathway-informed transformer.

MOTIVATION: Single-cell RNA sequencing (scRNA-Seq) technology enables detailed exploration of gene expression at the individual cell level, crucial for annotating cell types and understanding cellular diversity. Traditional methods for cell type annotation often rely on marker genes and manual labeling, posing challenges due to low data quality and incomplete reference datasets. RESULTS: We developed CeLLTra, a novel contrastive learning framework that leverages a Transformer-based model integrating biological pathway information to group genes into super tokens, effectively capturing comprehensive gene expression from scRNA-Seq data. By combining this pathway-informed Transformer with a pretrained domain-specific language model, CeLLTra accurately aligns cell-type annotations with gene expression profiles. Evaluations on a large-scale human scRNA-Seq dataset showed that CeLLTra significantly outperformed state-of-the-art methods in supervised and zero-shot cell-type prediction. Additionally, CeLLTra generalized well to external datasets, improving clustering performance and enabling better characterization of cancerous cell states in tumor-infiltrating myeloid cells from non-small cell lung cancer patients. AVAILABILITY AND IMPLEMENTATION: CeLLTra is freely available on GitHub (https://github.com/WJZheng-group/CeLLTra) and Zenodo (https://doi.org/10.5281/zenodo.17666735). The datasets underlying this article are the following: GSE201333 and GSE127465. All these datasets are publicly available and can be freely accessed on the Gene Expression Omnibus repository.

Humans

scSNViz: visualization and analysis of cell-specific expressed SNVs.

MOTIVATION: Accurately characterizing expressed genetic variation at the single-cell level is essential for understanding transcriptional heterogeneity, allelic regulation, and mutational dynamics within complex tissues. However, few tools enable comprehensive visualization and quantitative analysis of expressed variants across individual cells. RESULTS: scSNViz is an R package for the exploration, quantification, and visualization of expressed single-nucleotide variants (SNVs) from cell-barcoded single-cell RNA sequencing (scRNA-seq) data. The software supports estimation of variant allele fractions, clustering of SNV expression profiles, and 2D and 3D visualization of individual SNVs or user-defined SNV groups. Beyond visualization, scSNViz facilitates investigation of cell-, cluster-, or lineage-specific variant expression patterns, as well as allelic dynamics including imprinting, random allele inactivation, and transcriptional bursting. It interoperates seamlessly with established single-cell frameworks-Seurat for clustering, Slingshot for trajectory inference, scType for cell-type annotation, and CopyKat for copy-number profiling-enabling integrative multi-omic analyses of expressed variation. AVAILABILITY AND IMPLEMENTATION: scSNViz is implemented in R and freely available at https://github.com/HorvathLab/scSNViz (DOI: 10.5281/zenodo.17307516). The package includes comprehensive documentation and example workflows designed for users with limited bioinformatics experience.

Software

ARCADIA reveals spatially dependent transcriptional programs through integration of scRNA-seq and spatial proteomics.

MOTIVATION: Cellular states are strongly influenced by spatial context, but single-cell RNA sequencing (scRNA-seq) loses information about local tissue organization, while spatial proteomic assays capture limited marker panels that constrain transcriptomic inference. Integrating these modalities can elucidate how spatial niches shape transcriptional programs, yet existing approaches depend on either feature-level correspondence such as gene-protein linkage or cell-level barcode pairing, which is often unavailable. RESULTS: We present ARCADIA (ARchetype-based Clustering and Alignment with Dual Integrative Autoencoders), a generative framework for cross-modal integration that operates without cell barcode pairing and does not assume direct feature-to-feature correspondence. ARCADIA identifies modality-specific archetypes, that is, convex combinations of cells representing extreme phenotypic states, and aligns these anchors across modalities by minimizing the discrepancy between their cell-type composition profiles. The aligned archetypes define a shared coordinate system that anchors dual variational autoencoders (VAEs) trained with cross-modal geometric regularization, preserving archetype structure and spatial neighborhood information while enabling bidirectional translation between modalities. On semi-synthetic CITE-seq data, ARCADIA outperforms existing weak-linkage methods. Applied to independent human tonsil scRNA-seq and CODEX data, ARCADIA reconstructs known tissue architecture and reveals spatially dependent transcriptional programs linking B-cell maturation and T-cell activation or exhaustion to microenvironmental niches. AVAILABILITY AND IMPLEMENTATION: Source code is accessible at https://github.com/azizilab/ARCADIA_public. Reproducibility scripts and data are available at https://github.com/azizilab/arcadia_reproducibility.

Proteomics

GBMdeconvoluteR accurately infers proportions of neoplastic and immune cell populations from bulk glioblastoma transcriptomics data.

BACKGROUND: Characterizing and quantifying cell types within glioblastoma (GBM) tumors at scale will facilitate a better understanding of the association between the cellular landscape and tumor phenotypes or clinical correlates. We aimed to develop a tool that deconvolutes immune and neoplastic cells within the GBM tumor microenvironment from bulk RNA sequencing data. METHODS: We developed an IDH wild-type (IDHwt) GBM-specific single immune cell reference consisting of B cells, T-cells, NK-cells, microglia, tumor associated macrophages, monocytes, mast and DC cells. We used this alongside an existing neoplastic single cell-type reference for astrocyte-like, oligodendrocyte- and neuronal progenitor-like and mesenchymal GBM cancer cells to create both marker and gene signature matrix-based deconvolution tools. We applied single-cell resolution imaging mass cytometry (IMC) to ten IDHwt GBM samples, five paired primary and recurrent tumors, to determine which deconvolution approach performed best. RESULTS: Marker-based deconvolution using GBM-tissue specific markers was most accurate for both immune cells and cancer cells, so we packaged this approach as GBMdeconvoluteR. We applied GBMdeconvoluteR to bulk GBM RNAseq data from The Cancer Genome Atlas and recapitulated recent findings from multi-omics single cell studies with regards associations between mesenchymal GBM cancer cells and both lymphoid and myeloid cells. Furthermore, we expanded upon this to show that these associations are stronger in patients with worse prognosis. CONCLUSIONS: GBMdeconvoluteR accurately quantifies immune and neoplastic cell proportions in IDHwt GBM bulk RNA sequencing data and is accessible here: https://gbmdeconvoluter.leeds.ac.uk.

Humans

Multi-omics analysis reveals distinct spatial compartmentalization of lung repair niches in pediatric ARDS.

BACKGROUND: Pediatric acute respiratory distress syndrome (PARDS), often triggered by viral infections, is a life-threatening condition. Despite its severity, children demonstrate significantly better survival rates and superior lung repair compared to adults. However, the mechanisms underlying this age-specific advantage remain incompletely understood. PATIENTS AND METHODS: We conducted a pilot multi-omics study of influenza-associated PARDS integrating single-cell RNA sequencing (scRNA-seq) of pediatric lung tissue and bronchoalveolar lavage fluid (BALF), spatial transcriptomics, and plasma proteomics. Analyses were harmonized with the Human Lung Cell Atlas (HLCA) reference, reanalysis of public pediatric PARDS airway scRNA-seq, and contextual comparisons to adult lethal COVID-19 lung. RESULTS: Tissue scRNA-seq and spatial data indicated outcome-linked divergence in PARDS. Survivor showed spatially restricted repair with preserved alveolar type II (AT2) cells, AT2-to-alveolar type I (AT1) differentiation signatures, and higher KRT17, whereas fatal case and adults exhibited diffuse immune activation with pro-fibrotic and pro-apoptotic signaling. In BALF, KRT17-positive airway stress–repair epithelial cells (hillock-like) increased from the acute to recovery phase, and plasma proteomics showed higher circulating KRT17 in survivors. HLCA-based label transfer strengthened cell-type definitions and enabled pediatric–adult comparisons suggesting biological and developmental differences; the adult lethal COVID-19 atlas provided a benchmark with attenuated epithelial repair and prominent collagen CTHRC1-pathologic fibroblasts. Fibroblast programs were regionally compartmentalized, with injury-enriched CTHRC1+ states versus alveolar fibroblasts in preserved areas, and showed stronger injury–homeostasis anti-correlation in fatalities. Myeloid remodeling included BALF transitions from FCN1-high inflammatory states toward FABP4-positive resident-like states, consistent with public pediatric datasets showing reduced inflammatory and interferon-stimulated gene (ISG) modules and severity-linked increases in aged neutrophils. CONCLUSIONS: This pilot multi-omics case series outlines putative pediatric lung repair niches in influenza-associated PARDS. KRT17-positive transitional epithelium, preserved AT2 differentiation, and restoration of resident-like macrophages may align with recovery, whereas diffuse immune activation and CTHRC1-enriched fibroblast programs may accompany worse outcomes. HLCA-guided annotations and adult benchmarks indicate possible age-related differences, warranting validation in larger multi-center cohorts.

Humans

OmnibusX: A unified platform for accessible multi-omics analysis.

OmnibusX is an integrated, privacy-centric platform that enables code-free multi-omics data analysis by bridging computational methodologies with user-friendly interfaces. Designed to overcome challenges posed by fragmented analytical tools and high computational barriers, OmnibusX consolidates workflows for diverse technologies - including bulk RNA-seq, single-cell RNA-seq, single-cell ATAC-seq, and spatial transcriptomics - into a single, cohesive application. The application integrates established open-source tools such as Scanpy, DESeq2, SciPy, and scikit-learn into transparent, reproducible pipelines, offering users control over analytical parameters. Additionally, OmnibusX features proprietary modules, including a highly accurate cell-type prediction engine and an interactive plotting editor for generating publication-quality visualizations. Available as a standalone desktop application and an enterprise edition for centralized server deployment, OmnibusX ensures all data processing is conducted locally, eliminating external data transfer and usage tracking. By lowering technical barriers and enhancing reproducibility, OmnibusX aims to accelerate biological discovery and foster robust, data-driven collaborations. A fully documented trial version is accessible at: https://omnibusx.com/apps.

Computational Biology