WGCNA with eigengene expression plots of ortholog expression for the E5 timeseries molecular project

Today I updated a previously generated WGCNA to include plots that show module eigengene expression against factors of interest.

Overview

Please see my last post for all the details on the WGCNA. Find my previous heatmap with correlations for each module in my past post here. The code can be found here.

Module eigengene expression

As a reminder, here is the heatmap showing module correlations with factors of interest and physiological metrics.

I then generated plots that show module eigengene expression values for each module.

The module eigengene is the first principal component of the expression values of all genes within a module. It represents the dominant shared expression pattern of that module across samples. In other words, the eigengene value for each sample summarizes how “high” or “low” the genes in that module are expressed collectively, allowing us to compare module-level expression across species, timepoints, sites, or physiological states rather than examining individual genes.

Here are some notes on these plots.

  • They are faceted by WGCNA module (15 total modules)
  • The metric of interest is on the x-axis
  • Each panel with a red star shows modules that are significantly correlated with that factor or metric

Modules by site

Modules by time

Modules by species

Modules by maximal photosynthesis rate

Modules by photosynthetic efficiency

Modules by calcification

Modules by symbiont cell density

Modules by antioxidant capacity

Modules by host biomass

Modules by compensation irradiance

Modules by saturating irradiance

Modules by host protein content

Modules by host:symbiont biomass ratio

Modules by respiration

Modules by symbiont biomass

Modules by total chlorophyll per symbiont cell

Modules by total chlorophyll per unit biomass

Source Code

The complete analysis is documented in the R Markdown notebook:

Written on November 5, 2025