  ## Build the ClusterTree from Seurat object
  expression_seurat <- BuildClusterTree(expression_seurat,
                                        do.plot=TRUE)
  
  ## Plot custom phylogenetic tree that is more easily interpretable
  png(filename="./with_Batch_Cycle_Correction/ClusterTree.png", width=1800, height=1700, bg = "white", res = 150)
  plot(ladderize(expression_seurat@cluster.tree[[1]]), adj=0, 
       label.offset=50,
       edge.width=1) 
  nodelabels(frame="circle",bg="grey")
  dev.off()