Multi-Omics Factor Analysis (MOFA2) for unsupervised integration of multiple data modalities. Identifies shared and view-specific sources of variation. Use when integrating RNA-seq, proteomics, methylation, or other omics to discover latent factors driving biological variation across modalities.
Reference examples tested with: scanpy 1.10+
Before using code patterns, verify installed versions match. If versions differ:
packageVersion('<pkg>') then ?function_name to verify parametersIf code throws ImportError, AttributeError, or TypeError, introspect the installed package and adapt the example to match the actual API rather than retrying.
"Find shared variation across my omics layers" → Discover latent factors that capture shared and modality-specific sources of biological variation in an unsupervised manner.
MOFA2::create_mofa() → prepare_mofa() → run_mofa()mofapy2 for training, muon for downstreamGoal: Load and align multiple omics matrices into a consistent format for MOFA2 input.
Approach: Read each omics layer, intersect to common samples, transpose to features-by-samples orientation.
library(MOFA2)
library(MultiAssayExperiment)
# Load individual omics matrices (samples x features)
rna <- as.matrix(read.csv('rnaseq_matrix.csv', row.names = 1))
protein <- as.matrix(read.csv('proteomics_matrix.csv', row.names = 1))
methylation <- as.matrix(read.csv('methylation_matrix.csv', row.names = 1))
# Ensure consistent sample names across views
common_samples <- Reduce(intersect, list(rownames(rna), rownames(protein), rownames(methylation)))
rna <- rna[common_samples, ]
protein <- protein[common_samples, ]
methylation <- methylation[common_samples, ]
# Transpose to features x samples (MOFA format)
data_list <- list(
RNA = t(rna),
Protein = t(protein),
Methylation = t(methylation)
)
Goal: Configure and train a MOFA2 model to discover shared and view-specific latent factors.
Approach: Set model and training options, then run variational inference to learn factor decomposition.
# Create MOFA object
mofa <- create_mofa(data_list)
# View data overview
plot_data_overview(mofa)
# Set model options
model_opts <- get_default_model_options(mofa)
model_opts$num_factors <- 15 # Number of factors to learn
# Set training options
train_opts <- get_default_training_options(mofa)
train_opts$convergence_mode <- 'slow'
train_opts$seed <- 42
# Prepare and train
mofa <- prepare_mofa(mofa, model_options = model_opts, training_options = train_opts)
mofa <- run_mofa(mofa, outfile = 'mofa_model.hdf5')
Goal: Quantify how much variance each factor explains per omics view and extract factor scores and loadings.
Approach: Plot variance decomposition and retrieve factor values (sample scores) and weights (feature loadings) as data frames.
# Variance explained per factor per view
plot_variance_explained(mofa, max_r2 = 15)
plot_variance_explained(mofa, plot_total = TRUE)[[2]]
# Factor values (sample scores)
factors <- get_factors(mofa, as.data.frame = TRUE)
# Factor weights (feature loadings)
weights <- get_weights(mofa, as.data.frame = TRUE)
Goal: Generate publication-quality plots of factor values, feature weights, and factor correlations.
Approach: Use MOFA2 built-in plotting functions for scatter plots, heatmaps, and correlation matrices.
# Scatter plot of factor values
plot_factor(mofa, factor = 1, color_by = 'Group')
# Heatmap of factor weights
plot_weights(mofa, view = 'RNA', factor = 1, nfeatures = 20)
# Top features per factor
plot_top_weights(mofa, view = 'RNA', factor = 1, nfeatures = 10)
# Correlation between factors
plot_factor_cor(mofa)
Goal: Identify biological pathways and gene sets associated with each MOFA factor.
Approach: Extract top-weighted features per factor and run gene set enrichment on factor weights.
# Extract top features for pathway analysis
top_rna_factor1 <- get_weights(mofa, views = 'RNA', factors = 1, as.data.frame = TRUE)
top_rna_factor1 <- top_rna_factor1[order(abs(top_rna_factor1$value), decreasing = TRUE), ]
gene_list <- head(top_rna_factor1$feature, 100)
# Gene set enrichment on factor weights
library(MOFA2)
enrichment <- run_enrichment(mofa, feature.sets = msigdb_genesets,
view = 'RNA', factors = 1:5)
plot_enrichment(enrichment, factor = 1, max.pathways = 15)
Goal: Annotate MOFA factors with clinical or experimental metadata for colored visualizations.
Approach: Load sample annotations and attach to the MOFA object, then plot factors colored by metadata variables.
# Load sample annotations
metadata <- read.csv('sample_metadata.csv', row.names = 1)
# Add to MOFA object
samples_metadata(mofa) <- metadata[samples_names(mofa)[[1]], ]
# Color by metadata
plot_factor(mofa, factor = 1, color_by = 'Condition')
plot_factors(mofa, factors = 1:3, color_by = 'Condition')
Goal: Apply MOFA to batch- or group-structured multi-omics data and compare factor activity across groups.
Approach: Organize data as nested list by group, train multi-group MOFA, and visualize group-specific factor patterns.
# For batch/group-structured data
data_list_grouped <- list(
group1 = list(RNA = rna_g1, Protein = prot_g1),
group2 = list(RNA = rna_g2, Protein = prot_g2)
)
mofa_grouped <- create_mofa(data_list_grouped)
mofa_grouped <- prepare_mofa(mofa_grouped)
mofa_grouped <- run_mofa(mofa_grouped)
# Compare factor activity across groups
plot_factor(mofa_grouped, factor = 1, group_by = 'group', color_by = 'group')
Goal: Apply MOFA to single-cell multi-modal data (CITE-seq, Multiome) with stochastic inference for scalability.
Approach: Extract modality matrices from a Seurat object, create MOFA with stochastic training for large cell counts.
# For single-cell multi-omics (CITE-seq, Multiome)
library(Seurat)
# Extract modalities from Seurat object
rna_mat <- GetAssayData(seurat_obj, assay = 'RNA', layer = 'data')
adt_mat <- GetAssayData(seurat_obj, assay = 'ADT', layer = 'data')
# Create MOFA with single-cell settings
mofa_sc <- create_mofa(list(RNA = rna_mat, ADT = adt_mat))
model_opts <- get_default_model_options(mofa_sc)
model_opts$num_factors <- 10
# Use stochastic inference for large datasets
train_opts <- get_default_training_options(mofa_sc)
train_opts$stochastic <- TRUE
Goal: Save MOFA factor scores, feature weights, and variance explained to CSV for downstream use.
Approach: Extract each result type as a data frame and write to disk.
# Save factor values
factors_df <- get_factors(mofa, as.data.frame = TRUE)
write.csv(factors_df, 'mofa_factors.csv', row.names = FALSE)
# Save weights
weights_df <- get_weights(mofa, as.data.frame = TRUE)
write.csv(weights_df, 'mofa_weights.csv', row.names = FALSE)
# Save variance explained
var_exp <- get_variance_explained(mofa)
write.csv(var_exp$r2_per_factor, 'mofa_variance_explained.csv')