Bulk WGCNA analysis with omicverse

Perform bulk WGCNA on expression data to identify gene co-expression modules.

1.2k|145|Updated Mar 22, 2021
One-click install
npx skills add https://github.com/Starlitnightly/omicverse --skill bulk-wgcna-analysis-with-omicverse
Or copy as Structured Prompt for Agent
Please help me install this Agent Skill.
Skill: Bulk WGCNA analysis with omicverse
Source: https://github.com/Starlitnightly/omicverse/tree/main/.claude/skills/bulk-wgcna-analysis
Command: npx skills add https://github.com/Starlitnightly/omicverse --skill bulk-wgcna-analysis-with-omicverse

SYSTEM DOCUMENTATION & REQUIREMENTS

💡 This Skill requires omicverse, pandas, statsmodels, scanpy, matplotlib, and includes references (resource) components.

What problem does it solve?

Identifying biologically meaningful gene co-expression modules and hub genes from bulk RNA-seq data can be computationally intensive and require specialized algorithms. This Skill streamlines the Weighted Gene Co-expression Network Analysis (WGCNA) workflow.

Core Features & Use Cases

  • Expression Data Preprocessing: Prepare and filter bulk expression data for WGCNA.
  • Co-expression Network Construction: Build adjacency and topological overlap matrices to detect gene modules.
  • Hub Gene Identification: Extract and visualize key hub genes within specific modules of interest.
  • Use Case: Analyze bulk RNA-seq data from a disease cohort to find gene modules associated with disease progression, identify key hub genes within these modules, and correlate them with clinical traits.

Quick Start

Perform WGCNA on my 5xFAD expression data, visualize the resulting modules, and extract the top 10 hub genes from the 'lightgreen' module.

Frequently Asked Questions about Bulk WGCNA analysis with omicverse

High-intent search queries and answers about installing and using this skill.

FAQPage Schema
How do I identify gene co-expression modules from bulk RNA-seq data?

WGCNA (Weighted Gene Co-expression Network Analysis) identifies gene co-expression modules by constructing adjacency and topological overlap matrices from expression data, then detecting modules through dendrogram clustering. This reveals groups of genes with coordinated expression patterns across your samples.

What preprocessing steps are needed before running WGCNA on my expression data?

WGCNA requires filtering expression data using MAD (median absolute deviation), selecting variable genes, and ensuring data quality. These preprocessing steps remove noise and focus analysis on biologically informative genes before network construction.

Can I use WGCNA to correlate gene modules with clinical traits?

Yes, WGCNA enables module–trait analysis by calculating eigengenes (representative expression profiles for each module) and correlating them with sample metadata. This reveals which gene modules associate with clinical phenotypes or disease progression.

How do I extract and visualize hub genes from a specific co-expression module?

Hub gene extraction ranks genes within a module by connectivity and correlation strength, then sub-network plotting visualizes the top hub genes and their relationships. This identifies the most central genes driving module function.

Is WGCNA suitable for large-scale transcriptomic datasets?

Yes, WGCNA is designed for high-dimensional bulk RNA-seq data with many genes and samples. It efficiently handles large datasets and is implemented here via omicverse for streamlined analysis on complex transcriptomic studies.

What distinguishes soft-thresholding and adjacency matrix construction in WGCNA?

Soft-thresholding applies weighted power transformation to correlation values, preserving continuous information. This constructs the adjacency matrix, which feeds into topological overlap matrix (TOM) calculation to improve module detection reliability over hard thresholding.