Cited by 785Single-cell and spatial transcriptomicsBioinformatics and Genomic NetworksGene expression and cancer classification

The DOI matches its Crossref record; 92 of the 92 references it cites were found in a registry or an index; 89 of the 95 sources its claims cite were checked for support.

hdWGCNA identifies co-expression networks in high-dimensional transcriptomics data

Open access

Published version of a preprint

Abstract

Biological systems are immensely complex, organized into a multi-scale hierarchy of functional units based on tightly regulated interactions between distinct molecules, cells, organs, and organisms. While experimental methods enable transcriptome-wide measurements across millions of cells, popular bioinformatic tools do not support systems-level analysis. Here we present hdWGCNA, a comprehensive framework for analyzing co-expression networks in high-dimensional transcriptomics data such as single-cell and spatial RNA sequencing (RNA-seq). hdWGCNA provides functions for network inference, gene module identification, gene enrichment analysis, statistical tests, and data visualization. Beyond conventional single-cell RNA-seq, hdWGCNA is capable of performing isoform-level network analysis using long-read single-cell data. We showcase hdWGCNA using data from autism spectrum disorder and Alzheimer's disease brain samples, identifying disease-relevant co-expression network modules. hdWGCNA is directly compatible with Seurat, a widely used R package for single-cell and spatial transcriptomics analysis, and we demonstrate the scalability of hdWGCNA by analyzing a dataset containing nearly 1 million cells.

Overview

AI snapshot

Written from the paper’s full text.

TL;DR

hdWGCNA is an R package/framework for co-expression network analysis in single-cell and spatial transcriptomics data, extending Seurat and adding metacell/metaspot aggregation, network inference, module detection, enrichment, and visualization. The authors demonstrate it across human brain, mouse brain spatial, long-read isoform, ASD, AD, and million-cell PBMC datasets, identifying disease- and cell-type-relevant modules.

Key findings

  • hdWGCNA identified 96 co-expression modules across six major cell types in 36,671 human PFC snRNA-seq cells from 11 cognitively normal donors; ODC modules were reproducible in an independent dataset, with all ODC-specific modules significantly preserved (Z summary preservation ≥5).
  • In mouse brain Visium ST, hdWGCNA found 12 spatial modules from 2,696 anterior and 3,353 posterior spots, including white-matter myelination module SM1 and ventricle/BBB-associated SM4 enriched for hemoglobin subunits.
  • Isoform-level analysis of 2,190 radial glia lineage cells from mouse hippocampus ScISOrSeq yielded 11 modules; 61.5% of retained genes had a single isoform, 18.2% had multiple isoforms all in the same module, and 20.4% had isoforms spread across several modules.
  • In ASD inhibitory neurons (20,249 INH nuclei), hdWGCNA found 14 modules; modules INH-M11, INH-M13, and INH-M3 were enriched for ASD-associated genes, and all modules except INH-M4 showed significant ME differences in at least one INH subpopulation. In AD microglia consensus analysis, four modules were identified, with MG-M1 containing DAM genes and MG-M2 containing homeostatic markers.

Methods

  • Design: hdWGCNA builds on WGCNA-style co-expression network analysis, using metacells for scRNA-seq or metaspots for spatial transcriptomics to reduce sparsity, then signed adjacency, soft-power thresholding, signed topological overlap matrix (TOM), and Dynamic Tree Cut for module detection.
  • Data: applications included human PFC snRNA-seq (36,671 cells), mouse brain Visium ST (6,049 spots total), mouse hippocampus ScISOrSeq (6,832 nuclei; 2,190 radial glia lineage cells), ASD snRNA-seq (20,249 inhibitory neurons), AD microglia from three snRNA-seq datasets (9,904 microglia nuclei from 189,127 total nuclei), and a PBMC scRNA-seq dataset (965,363 cells retained).
  • Analysis: module eigengenes (MEs) computed by SVD with optional Harmony correction; downstream tests included differential module eigengene analysis, module preservation with permutations, EGAD functional coherence, XGBoost hub-gene prediction, Monocle3 pseudotime, scDRS polygenic risk scoring, and projection of bulk RNA-seq modules into snRNA-seq.

Limitations

  • Stated: single-cell transcriptomic measurements are noisy, with technical noise from dropout events and experimental protocols; hdWGCNA attempts to address this by constructing networks in metacell or metaspot profiles and by using module preservation testing for reproducibility.
  • Evident from the text: the metacell algorithm requires a dimensional reduction of the input matrix, and technical artifacts in that reduction can influence metacell construction; the optimal number of cells per metacell may differ across cell types and tissues.
  • Evident from the text: pseudo-bulk aggregation is suggested as likely suitable for co-expression network analysis only given sufficient sample size, and several analyses rely on external datasets and parameter choices for validation.

Model-written summary. Check it against the paper before relying on it.

Integrity signals

Registry

Crossref · 10.1016/j.crmeth.2023.100498

Cites

92 of 92 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Figures

Graphical abstract
Graphical abstract CC BY · PubMed Central
Figure 1 Overview of the hdWGCNA workflow and application in the human prefrontal cortex (A) Schematic overview of the standard hdWGCNA workflow on a scRNA-seq dataset. UMAP plot shows 36,671 cells from 11 cognitively normal donors in the Zhou et al. human prefrontal cortex (PFC) dataset. ASC, astrocytes; EX, excitatory neurons; INH, inhibitory neurons; MG, microglia; ODC, oligodendrocytes; OPC, oligodendrocyte progenitor cells. (B) Density plot showing the distribution of pairwise Pearson correlations between genes from the single-cell (sc) expression matrix and metacell expression matrices with varying values of the K-nearest neighbors parameter K. (C) Expression matrix density (1, sparsity) for the sc, pseudo-bulk (pb), and metacell matrices with varying values of K in each cell type. (D) Heatmap of scaled gene expression for the top five hub genes by kME in INH-M6, EX-M2, ODC-M3, OPC-M2, ASC-M18, and MG-M14. (E) snRNA-seq UMAP colored by module eigengene (ME) for selected modules as in (D). (F) UMAP plot of the ODC co-expression network. Each node represents a single gene, and edges represent co-expression links between genes and module hub genes. Point size is scaled by kME. Nodes are colored by co-expression module assignment. The top two hub genes per module are labeled. Network edges were downsampled for visual clarity. (G) snRNA-seq UMAP as in (A) colored by MEs for the 10 ODC co-expression modules as in (F). (H) Module preservation analysis of the ODC modules in the Morabito et al.12 human PFC dataset. The module’s size versus the preservation statistic (Z preservation) is shown for each module. Z<5, not preserved; 10>Z≥5, moderately preserved; Z≥10, highly preserved.
Figure 1 Overview of the hdWGCNA workflow and application in the human prefrontal cortex (A) Schematic overview of the standard hdWGCNA workflow on a scRNA-seq dataset. UMAP plot shows 36,671 cells from 11 cognitively normal donors in the Zhou et al. human prefrontal cortex (PFC) dataset. ASC, astrocytes; EX, excitatory neurons; INH, inhibitory neurons; MG, microglia; ODC, oligodendrocytes; OPC, oligodendrocyte progenitor cells. (B) Density plot showing the distribution of pairwise Pearson correlations between genes from the single-cell (sc) expression matrix and metacell expression matrices with varying values of the K-nearest neighbors parameter K. (C) Expression matrix density (1, sparsity) for the sc, pseudo-bulk (pb), and metacell matrices with varying values of K in each cell type. (D) Heatmap of scaled gene expression for the top five hub genes by kME in INH-M6, EX-M2, ODC-M3, OPC-M2, ASC-M18, and MG-M14. (E) snRNA-seq UMAP colored by module eigengene (ME) for selected modules as in (D). (F) UMAP plot of the ODC co-expression network. Each node represents a single gene, and edges represent co-expression links between genes and module hub genes. Point size is scaled by kME. Nodes are colored by co-expression module assignment. The top two hub genes per module are labeled. Network edges were downsampled for visual clarity. (G) snRNA-seq UMAP as in (A) colored by MEs for the 10 ODC co-expression modules as in (F). (H) Module preservation analysis of the ODC modules in the Morabito et al.12 human PFC dataset. The module’s size versus the preservation statistic (Z preservation) is shown for each module. Z<5, not preserved; 10>Z≥5, moderately preserved; Z≥10, highly preserved.CC BY · PubMed Central
Figure 2 Runtime, memory usage, and performance of hdWGCNA (A and B) We ran the main co-expression network analysis functions of the hdWGCNA R package on 65,415 neuronal cells in a human brain dataset9 from 54 samples, and tracked the runtime (A) and memory usage upper bound (B) for different-sized subsets of the data ranging from 1,000 through 50,000 cells. (C) Violin plots showing distributions of EGAD23 neighbor-voting area under the receiver operating characteristic curve (AUC) scores in each of the cell-type-specific co-expression networks from the human PFC dataset.11 (D) Violin plots showing distributions of multifunctionality AUC scores in each of the cell-type-specific co-expression networks from the human PFC dataset. ASC, astrocytes; EX, excitatory neurons; INH, inhibitory neurons; MG, microglia; ODC, oligodendrocytes; OPC, oligodendrocyte progenitor cells. (E) Performance of the XGBoost regularized regression models used to predict gene expression based on the expression of the top 10 module hub genes for all 96 co-expression modules from the Zhou et al.11 human PFC dataset. Violin plots showing the test set root-mean-square error (RMSE) comparing the predicted expression with observed for each gene, split by each co-expression module. Modules are ordered within each cell type from lowest mean RMSE to highest.
Figure 2 Runtime, memory usage, and performance of hdWGCNA (A and B) We ran the main co-expression network analysis functions of the hdWGCNA R package on 65,415 neuronal cells in a human brain dataset9 from 54 samples, and tracked the runtime (A) and memory usage upper bound (B) for different-sized subsets of the data ranging from 1,000 through 50,000 cells. (C) Violin plots showing distributions of EGAD23 neighbor-voting area under the receiver operating characteristic curve (AUC) scores in each of the cell-type-specific co-expression networks from the human PFC dataset.11 (D) Violin plots showing distributions of multifunctionality AUC scores in each of the cell-type-specific co-expression networks from the human PFC dataset. ASC, astrocytes; EX, excitatory neurons; INH, inhibitory neurons; MG, microglia; ODC, oligodendrocytes; OPC, oligodendrocyte progenitor cells. (E) Performance of the XGBoost regularized regression models used to predict gene expression based on the expression of the top 10 module hub genes for all 96 co-expression modules from the Zhou et al.11 human PFC dataset. Violin plots showing the test set root-mean-square error (RMSE) comparing the predicted expression with observed for each gene, split by each co-expression module. Modules are ordered within each cell type from lowest mean RMSE to highest.CC BY · PubMed Central
Figure 3 Spatial co-expression networks represent regional expression patterns in the mouse brain (A) Visium spatial transcriptomics (ST) in anterior (left, 2,696 spots) and posterior (right, 3,353 spots) mouse brain sections, colored by Louvain clusters annotated by anatomical regions. (B) UMAP plot of the mouse brain ST co-expression network. Each node represents a single gene, and edges represent co-expression links between genes and module hub genes. Point size is scaled by kME. Nodes are colored by co-expression module assignment. The top five hub genes per module are labeled. Network edges were downsampled for visual clarity. (C) ST samples colored by MEs for the 12 spatial co-expression modules. Gray color indicates an ME value less than zero.
Figure 3 Spatial co-expression networks represent regional expression patterns in the mouse brain (A) Visium spatial transcriptomics (ST) in anterior (left, 2,696 spots) and posterior (right, 3,353 spots) mouse brain sections, colored by Louvain clusters annotated by anatomical regions. (B) UMAP plot of the mouse brain ST co-expression network. Each node represents a single gene, and edges represent co-expression links between genes and module hub genes. Point size is scaled by kME. Nodes are colored by co-expression module assignment. The top five hub genes per module are labeled. Network edges were downsampled for visual clarity. (C) ST samples colored by MEs for the 12 spatial co-expression modules. Gray color indicates an ME value less than zero.CC BY · PubMed Central
Figure 4 Isoform co-expression network analysis reveals fate-specific expression programs in the hippocampal radial glia lineage (A) UMAP plot of cells from the mouse hippocampus ScISOrSeq dataset.8 Major cell types are labeled and the cells used for co-expression network analysis are colored. This dataset contains expression information for 96,093 isoforms and 31,053 genes in 6,832 cells from one mouse brain sample. ASC, astrocytes; CPX, choroid plexus epithelial cells; EPD, ependymal cells; EX, excitatory neurons; GRN, granule neurons; INH, inhibitory neurons; MAC, macrophages; NPC, neuronal intermediate progenitor cells; MG, microglia; OPC, oligodendrocyte progenitor cells; RGL, radial glia; VASC, vasculature cells. (B) UMAP plot of the radial glia lineage, colored by Monocle 337 pseudotime assignment. Top left, ependymal (EPD) trajectory; top right, astrocyte (ASC) trajectory; bottom left, neuronal intermediate progenitor cell (NPC) trajectory. (C) UMAP plot of the radial glia lineage isoform co-expression network. Each node represents a single isoform, and edges represent co-expression links between isoforms and module hub isoforms. Point size is scaled by kMEiso. Nodes are colored by co-expression module assignment. Network edges were downsampled for visual clarity. (D) Donut chart showing the percentage of genes with one isoform, with multiple isoforms that are all assigned to the same module, and with multiple isoforms that are spread across more than one module. (E) Module eigenisoforms (MEiso) as a function of pseudotime for each co-expression module. For each module, a separate locally estimated scatterplot smoothing (LOESS) regression line is shown for each developmental trajectory. (F) Dot plot showing selected GO term enrichment results for each co-expression module. (G) Gene models for selected isoforms of Gfap, colored by co-expression module assignment. (H) Gene models for selected isoforms of H3f3b, colored by co-expression module assignment. (I) Top: gene models for selected isoforms of Cd9, colored by co-expression module assignment. Bottom: Swan38 graphical representation of Cd9 alternative splicing isoforms. Splice sites and transcript start/end sites are represented as nodes; introns and exons are represented as connections between nodes. These two isoforms are distinguished by alternative TSS usage. Gene models from the GENCODE VM23 comprehensive transcript set are shown below transcripts in panels (G)–(I).
Figure 4 Isoform co-expression network analysis reveals fate-specific expression programs in the hippocampal radial glia lineage (A) UMAP plot of cells from the mouse hippocampus ScISOrSeq dataset.8 Major cell types are labeled and the cells used for co-expression network analysis are colored. This dataset contains expression information for 96,093 isoforms and 31,053 genes in 6,832 cells from one mouse brain sample. ASC, astrocytes; CPX, choroid plexus epithelial cells; EPD, ependymal cells; EX, excitatory neurons; GRN, granule neurons; INH, inhibitory neurons; MAC, macrophages; NPC, neuronal intermediate progenitor cells; MG, microglia; OPC, oligodendrocyte progenitor cells; RGL, radial glia; VASC, vasculature cells. (B) UMAP plot of the radial glia lineage, colored by Monocle 337 pseudotime assignment. Top left, ependymal (EPD) trajectory; top right, astrocyte (ASC) trajectory; bottom left, neuronal intermediate progenitor cell (NPC) trajectory. (C) UMAP plot of the radial glia lineage isoform co-expression network. Each node represents a single isoform, and edges represent co-expression links between isoforms and module hub isoforms. Point size is scaled by kMEiso. Nodes are colored by co-expression module assignment. Network edges were downsampled for visual clarity. (D) Donut chart showing the percentage of genes with one isoform, with multiple isoforms that are all assigned to the same module, and with multiple isoforms that are spread across more than one module. (E) Module eigenisoforms (MEiso) as a function of pseudotime for each co-expression module. For each module, a separate locally estimated scatterplot smoothing (LOESS) regression line is shown for each developmental trajectory. (F) Dot plot showing selected GO term enrichment results for each co-expression module. (G) Gene models for selected isoforms of Gfap, colored by co-expression module assignment. (H) Gene models for selected isoforms of H3f3b, colored by co-expression module assignment. (I) Top: gene models for selected isoforms of Cd9, colored by co-expression module assignment. Bottom: Swan38 graphical representation of Cd9 alternative splicing isoforms. Splice sites and transcript start/end sites are represented as nodes; introns and exons are represented as connections between nodes. These two isoforms are distinguished by alternative TSS usage. Gene models from the GENCODE VM23 comprehensive transcript set are shown below transcripts in panels (G)–(I).CC BY · PubMed Central
Figure 5 Co-expression network analysis of inhibitory neurons in Autism spectrum disorder (A) UMAP plot of 121,451 nuclei from the cortex of 22 ASD donors, 24 controls, and eight epilepsy donors profiled with snRNA-seq. Inhibitory neuron subtypes are highlighted. ASC, astrocytes; EX, excitatory neurons; INH, inhibitory neurons; MG, microglia; ODC, oligodendrocytes; OPC, oligodendrocyte progenitor cells. (B) Gene co-expression network derived from inhibitory neurons, represented as a two-dimensional UMAP embedding of the TOM. Nodes represent genes, colored by module assignment. Module hub genes with prior evidence of ASD association from SFARI are labeled. Edges represent co-expression relationships between genes and module hub genes. Network edges were downsampled for visual clarity. (C) Gene overlap analysis comparing ASD-associated genes from SFARI and INH co-expression modules, using Fisher’s exact test. × indicates that the overlap was not significant (FDR >0.05). (D) snRNA-seq UMAP plots as in (A) colored by MEs for INH co-expression modules. (E) Violin plots showing MEs in each INH cluster. Two-sided Wilcoxon test was used to compare ASD versus control samples. Nuclei from epilepsy donors were excluded in this comparison. Not significant (ns), p > 0.05; ^(∗)p≤0.05; ∗∗p≤0.01; ∗∗∗p≤0.001; ∗∗∗∗p≤0.0001. (F) Selected GO enrichment results for each co-expression module.
Figure 5 Co-expression network analysis of inhibitory neurons in Autism spectrum disorder (A) UMAP plot of 121,451 nuclei from the cortex of 22 ASD donors, 24 controls, and eight epilepsy donors profiled with snRNA-seq. Inhibitory neuron subtypes are highlighted. ASC, astrocytes; EX, excitatory neurons; INH, inhibitory neurons; MG, microglia; ODC, oligodendrocytes; OPC, oligodendrocyte progenitor cells. (B) Gene co-expression network derived from inhibitory neurons, represented as a two-dimensional UMAP embedding of the TOM. Nodes represent genes, colored by module assignment. Module hub genes with prior evidence of ASD association from SFARI are labeled. Edges represent co-expression relationships between genes and module hub genes. Network edges were downsampled for visual clarity. (C) Gene overlap analysis comparing ASD-associated genes from SFARI and INH co-expression modules, using Fisher’s exact test. × indicates that the overlap was not significant (FDR >0.05). (D) snRNA-seq UMAP plots as in (A) colored by MEs for INH co-expression modules. (E) Violin plots showing MEs in each INH cluster. Two-sided Wilcoxon test was used to compare ASD versus control samples. Nuclei from epilepsy donors were excluded in this comparison. Not significant (ns), p > 0.05; ^(∗)p≤0.05; ∗∗p≤0.01; ∗∗∗p≤0.001; ∗∗∗∗p≤0.0001. (F) Selected GO enrichment results for each co-expression module.CC BY · PubMed Central
Figure 6 Consensus network analysis of microglia in AD (A) Left: table showing the number of samples and the number of microglia nuclei from published AD snRNA-seq datasets used for co-expression network analysis. Right: integrated UMAP plot of nuclei from three snRNA-seq datasets. (B) UMAP plot of microglia, colored by Monocle 337 pseudotime assignment. (C) UMAP plot of the microglia co-expression network. Each node represents a single gene, and edges represent co-expression links between genes and module hub genes. Point size is scaled by kME. Nodes are colored by co-expression module assignment. The top 10 hub genes per module are labeled, as well as additional genes of interest. Network edges were downsampled for visual clarity. (D) Selected Gene Ontology (GO) terms enriched in co-expression modules. Bar plots show the log-scaled enrichment of each term. (E) MEs as a function of pseudotime; points are averaged MEs in 50 pseudotime bins of equal size. Line represents LOESS regression with a 95% confidence interval. (F) Microglia UMAP colored by ME. (G) Differential module eigengene (DME) results in 10 pseudotime bins of equal size. For each pseudotime bin, we performed DME analysis between cells from AD (positive fold change) and control samples. × symbol indicates that the test did not reach significance (Wilcoxon rank-sum test Bonferroni-adjusted p value > 0.05). (H) Top: microglia UMAP colored by AD single-cell disease relevance score (scDRS)52Z score. Bottom: scDRS Z score as a function of pseudotime, points are averaged scDRS Z scores in 50 pseudotime bins of equal size. Line represents linear regression with a 95% confidence interval. (I) Heatmap of Pearson correlations of MEs and scDRS Z scores, split by cells from AD and control samples. (J) Abbreviations denote the following brain regions: SFG, superior frontal gyrus; EC, entorhinal cortex; OC, occipital cortex; OTC, occipitotemporal cortex. ∗∗Highly preserved (Z≥10); ^(∗)moderately preserved (10>Z≥5); x, not preserved (Z<5).
Figure 6 Consensus network analysis of microglia in AD (A) Left: table showing the number of samples and the number of microglia nuclei from published AD snRNA-seq datasets used for co-expression network analysis. Right: integrated UMAP plot of nuclei from three snRNA-seq datasets. (B) UMAP plot of microglia, colored by Monocle 337 pseudotime assignment. (C) UMAP plot of the microglia co-expression network. Each node represents a single gene, and edges represent co-expression links between genes and module hub genes. Point size is scaled by kME. Nodes are colored by co-expression module assignment. The top 10 hub genes per module are labeled, as well as additional genes of interest. Network edges were downsampled for visual clarity. (D) Selected Gene Ontology (GO) terms enriched in co-expression modules. Bar plots show the log-scaled enrichment of each term. (E) MEs as a function of pseudotime; points are averaged MEs in 50 pseudotime bins of equal size. Line represents LOESS regression with a 95% confidence interval. (F) Microglia UMAP colored by ME. (G) Differential module eigengene (DME) results in 10 pseudotime bins of equal size. For each pseudotime bin, we performed DME analysis between cells from AD (positive fold change) and control samples. × symbol indicates that the test did not reach significance (Wilcoxon rank-sum test Bonferroni-adjusted p value > 0.05). (H) Top: microglia UMAP colored by AD single-cell disease relevance score (scDRS)52Z score. Bottom: scDRS Z score as a function of pseudotime, points are averaged scDRS Z scores in 50 pseudotime bins of equal size. Line represents linear regression with a 95% confidence interval. (I) Heatmap of Pearson correlations of MEs and scDRS Z scores, split by cells from AD and control samples. (J) Abbreviations denote the following brain regions: SFG, superior frontal gyrus; EC, entorhinal cortex; OC, occipital cortex; OTC, occipitotemporal cortex. ∗∗Highly preserved (Z≥10); ^(∗)moderately preserved (10>Z≥5); x, not preserved (Z<5).CC BY · PubMed Central
Figure 7 Projecting bulk RNA-seq co-expression modules into a single-cell dataset (A) UMAP plot of 57,950 nuclei from an snRNA-seq dataset of the human PFC from AD (N = 11) and control (N = 7) PFC samples.12 Cells are colored by major cell type assignment. (B) Multi-region consensus co-expression modules from Morabito et al.56 bulk RNA-seq analysis projected into the snRNA-seq dataset as in (A). (C) Co-expression modules from the AMP-AD bulk RNA-seq dataset63 projected into the snRNA-seq dataset as in (A). CBE, cerebellum; DLPFC, dorsolateral PFC; FP, frontal pole; IFG, inferior frontal gyrus; PHG, parahippocampal gyrus; STG, superior temporal gyrus; TCX, temporal cortex; ASC, astrocytes; EX, excitatory neurons; INH, inhibitory neurons; MG, microglia; ODC, oligodendrocytes; OPC, oligodendrocyte progenitor cells; VASC, vascular cells.
Figure 7 Projecting bulk RNA-seq co-expression modules into a single-cell dataset (A) UMAP plot of 57,950 nuclei from an snRNA-seq dataset of the human PFC from AD (N = 11) and control (N = 7) PFC samples.12 Cells are colored by major cell type assignment. (B) Multi-region consensus co-expression modules from Morabito et al.56 bulk RNA-seq analysis projected into the snRNA-seq dataset as in (A). (C) Co-expression modules from the AMP-AD bulk RNA-seq dataset63 projected into the snRNA-seq dataset as in (A). CBE, cerebellum; DLPFC, dorsolateral PFC; FP, frontal pole; IFG, inferior frontal gyrus; PHG, parahippocampal gyrus; STG, superior temporal gyrus; TCX, temporal cortex; ASC, astrocytes; EX, excitatory neurons; INH, inhibitory neurons; MG, microglia; ODC, oligodendrocytes; OPC, oligodendrocyte progenitor cells; VASC, vascular cells.CC BY · PubMed Central

Verification

Claims the paper makes

68 Own result · 38 Supported · 31 Own interpretation · 22 No basis given · 15 Partly supported · 10 Own method · 5 Couldn’t confirm · 2 Sources not found · 2 No full text of the sources · 1 Supported so far · 1 Not checked yet · accuracy not yet measured

Show all 105 claims

  • Probably not a substantial claim. 71
  • It may state two separate claims as one. 35
90 of 90
  • The paper relies on this without support.

    Popular bioinformatic tools do not support systems-level analysis.

    “While experimental methods enable transcriptome-wide measurements across millions of cells, popular bioinformatic tools do not support systems-level analysis.”
  • The paper relies on this without support.

    Sparsity and noise inherent in single-cell data can lead to spurious gene-gene correlations, thereby complicating co-expression network analysis.

    “The sparsity and noise inherent in single-cell data can lead to spurious gene-gene correlations, thereby complicating co-expression network analysis.”
  • The paper relies on this without support.

    The correlation structure of single-cell or spatial transcriptomic data varies greatly for different subsets (cell types, cell states, anatomical regions).

    “Additionally, the correlation structure of single-cell or spatial transcriptomic data varies greatly for different subsets (cell types, cell states, anatomical regions).”
  • The paper relies on this without support.

    Data sparsity is inherent to the current generation of sequencing-based spatial transcriptomics technologies.

    “Data sparsity is also inherent to the current generation of these technologies, therefore we propose a metaspot aggregation approach prior to network analysis (Figure S6).”
  • The authors state this as an assumption.

    Co-expressed genes are involved in similar biological processes, so co-expression network structures should be predictive of biological pathway membership.

    “In principle, we expect that co-expressed genes are involved in similar biological processes, and therefore co-expression network structures should be predictive of biological pathway membership.”
  • The authors state this as an assumption.

    Genes within the same co-expression modules derived from specific cell types should be functionally related or co-regulated.

    “In principle, genes within the same co-expression modules derived from specific cell types should be functionally related or co-regulated.”
  • The paper relies on this without support.

    Co-expression networks behave as functional biological units.

    “Co-expression networks behave as functional biological units; therefore, we reason that the hub genes and other members of MG-M1 represent candidates for an expanded set of human DAM genes including ACTB, TPT1, and EEF1A1.”
  • The paper relies on this without support.

    It remains cost-prohibitive for most researchers to perform high-dimensional -omics studies of large patient cohorts.

    “To date, it remains cost-prohibitive for most researchers to perform high-dimensional -omics studies of large patient cohorts, but there are numerous large-scale disease-relevant bulk RNA-seq datasets containing thousands of samples from consortia such as the Encyclopedia of DNA Elements (ENCODE),60 the Genotype-Tissue Expression (GTEx) project,61 and The Cancer Genome Atlas (TCGA).62”
  • The paper relies on this without support.

    Classical bioinformatic approaches such as differential gene expression do not provide information about the broader context of genes in specific pathways or regulatory regimes.

    “Classical bioinformatic approaches for transcriptomics analysis such as differential gene expression are useful for finding individual genes that are altered in a particular disease or condition of interest, but they do not provide information about the broader context of these genes in specific pathways or regulatory regimes.”
  • The paper relies on this without support.

    Biological processes such as development or regeneration require coordination of distinct sets of genes in certain cell types with spatial specificity.

    “For example, biological processes such as development or regeneration require coordination of distinct sets of genes in certain cell types with spatial specificity.”
  • The paper relies on this without support.

    Single-cell gene expression datasets typically contain many more zero-valued entries than non-zero-valued entries, meaning that these datasets are sparse.

    “Single-cell gene expression datasets typically contain many more zero valued entries than non-zero valued entries, meaning that these datasets are sparse.”
  • The paper relies on this without support.

    Conventional single-cell gene expression assays yield sparse gene expression matrices.

    “Conventional single-cell gene expression assays yield sparse gene expression matrices.”
  • The paper relies on this without support.

    Correlations of sparse vectors may lead to downstream conclusions that are not robust or reproducible.

    “In general, correlations of sparse vectors may lead to downstream conclusions that are not robust or reproducible.”
  • The paper relies on this without support.

    Transcriptomic measurements of single cells are generally noisy, imposing challenges and limitations in the analysis of these datasets.

    “Transcriptomic measurements of single cells are generally noisy, imposing challenges and limitations in the analysis of these datasets.”
  • The paper relies on this without support.

    Technical noise may arise from dropout events or from various steps in the experimental protocols, potentially making downstream data analysis and interpretation more difficult.

    “Technical noise may arise from dropout events or from various steps in the experimental protocols, potentially making downstream data analysis and interpretation more difficult.”
  • The paper relies on this without support.

    Sequencing-based spatial transcriptomics approaches such as 10x Genomics Visium yield sparse transcriptomic profiles, introducing the same potential pitfalls as single-cell data for co-expression network analysis.

    “Sequencing-based ST approaches such as the 10× Genomics Visium platform also yield sparse transcriptomic profiles, thus introducing the same potential pitfalls as single-cell data for co-expression network analysis.”
  • The paper relies on this without support.

    The widespread adoption of single-cell genomics has led to many biologists running their own computational analysis.

    “The widespread adoption of single-cell genomics has led to many biologists running their own computational analysis, and we designed the hdWGCNA R package with these individuals in mind through our various step-by-step tutorials and detailed documentation.”
17 of 17

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

Each contradiction between two passages listed here was confirmed by an independent check.

5 of 5

Citations

Cites

Show all 92 references

Research works printed but not linked: 1

  1. Fleming, 2019, CellBender remove-background: a deep generative model for unsupervised removal of background noise from scRNA-seq datasets, bioRxiv

Cited by

Sorted by citations

Metadata

Authors
Morabito, S., Reese, F., Rahimzadeh, N., Miyoshi, E., & Swarup, V.
Affiliations
University of California, Irvine (US)
Corresponding author
Swarup, Vivek
Venue
Cell Reports Methods
Publisher
Elsevier BV
Published
2023-06
Volume
3
Issue
6
Pages
100498
ISSN
2667-2375 (ISSN-L)
OpenAlex
W4380370818
Cited by
773
References
92
Citations by year
2026: 267 · 2025: 285 · 2024: 180 · 2023: 40
FWCI
58.73
Type
journal article
Language
en
Open access
Locations
  • Cell Reports Methods (publishedVersion, cc-by)
  • bioRxiv (Cold Spring Harbor Laboratory) (acceptedVersion, cc-by-nc-nd)
  • PubMed (publishedVersion)
  • PubMed Central (submittedVersion, cc-by)
  • eScholarship (California Digital Library) (submittedVersion, cc-by)
Funders
  • Uganda Cancer Institute
  • NIDA (1U01DA053826)
  • National Institute on Aging (1RF1AG071683, 3U19AG068054-02S, U54 AG054349-06)
  • NINDS (P01NS084974- 06A1)
  • Dr. Miriam and Sheldon G. Adelson Medical Research Foundation
  • University of California, Irvine
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Single-cell and spatial transcriptomics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Bioinformatics and Genomic Networks
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Gene expression and cancer classification
Keywords
co-expression network, single-cell RNA sequencing, spatial transcriptomics, autism spectrum disorder, Alzheimer's disease, brain transcriptomics
MeSH terms
Autism Spectrum Disorder, Alzheimer Disease, Humans, Gene Expression Profiling, Gene Regulatory Networks, Transcriptome

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

A language model drew the claims from the paper on 2026-10-13, under version 3 of the claim definition (model deepseek-v4.1-flash, prompt extract-v3).

A second model labelled each claim on 2026-10-13 (model jev-1.13-20260917, definition revision gold-guideline-v0.2).

Cited sources are checked with a language model’s help (gpt-6-sol).

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review