## Keep the significant markers and sort them by significance

sign_markers <- names(which(sort(de_out4$adjp[, "adjp_BCRXLvsRef"]) < FDR_cutoff))

sign_adjp <- de_out4$adjp[sign_markers , "adjp_BCRXLvsRef"]

## Normalize expression to mean = 0 and sd = 1

plot_differential_heatmap_wrapper(expr_norm = expr_median_sample_norm, 
  sign_adjp = sign_adjp, condition = md$condition[mm], 
  color_conditions = color_conditions)