merge_samples/tax_glom后phyloseq的PCoA样本无法按genus着色问题

核心问题 #

合并样本后PCoA图按genus着色失效（无颜色输出且无报错），根源在于合并样本过程中，样本元数据（sam_data）的genus列未正确关联到合并后的样本——即便已解决NA和列名前缀问题，仍可能存在genus值缺失、映射错误的情况。

排查步骤 #

先确认合并+tax_glom后phyloseq对象的样本元数据状态：

# 查看最终对象的样本元数据
head(sample_data(your_final_ps))
# 检查genus列的NA情况
table(is.na(sample_data(your_final_ps)$genus))

若输出显示genus列存在NA，或该列根本不存在，就是着色失效的直接原因。

解决方案 #

1. 合并样本时手动绑定正确的genus信息

merge_samples默认仅保留分组列（此处为species），会丢失其他分类阶元。需提前建立species-genus对应关系，合并后手动映射：

# 从原始phyloseq对象提取唯一的species-genus对应表（确保同一species对应唯一genus）
sp_genus_mapping <- unique(as.data.frame(sample_data(original_ps))[, c("species", "genus")])

# 按species合并样本
merged_ps <- merge_samples(original_ps, group = "species", fun = sum)

# 将genus映射到合并后的样本（合并后的样本名即为species名）
sample_data(merged_ps)$genus <- sp_genus_mapping$genus[match(sample_names(merged_ps), sp_genus_mapping$species)]

2. 执行tax_glom后确认元数据完整性

tax_glom仅合并OTU分类，不会修改样本元数据，但操作后仍需验证：

# 按目标分类阶元合并OTU（示例为真菌属水平，可按需调整taxrank参数）
glommed_ps <- tax_glom(merged_ps, taxrank = "genus")

# 确认genus列无NA
any(is.na(sample_data(glommed_ps)$genus))

3. 绘图时明确指定映射变量

用plot_ordination绘图时，确保color参数精准指向sam_data中的genus列：

# 计算Bray-Curtis距离
bray_dist <- distance(glommed_ps, method = "bray")
# 执行PCoA分析
pcoa_result <- ordinate(glommed_ps, method = "PCoA", distance = bray_dist)
# 绘图并按genus着色
plot_ordination(glommed_ps, pcoa_result, color = "genus") +
  geom_point(size = 3) +
  theme_bw() +
  labs(color = "Moth Genus")

额外提示 #

若原始数据中存在同一species对应多个genus的情况（分类信息错误），需先清理sam_data的分类数据，确保一个species仅对应一个genus，否则映射后仍会出现NA或错误值。

内容的提问来源于stack exchange，提问作者sondre