Coefficient of Variation (cV) was calculated using FinCal v0.6.3 (55 ). Differential gene expression analysis was performed with limma package (49 (link)). Permutational multivariate analysis of variance (PERMANOVA) was used to compare distances between gene expression profiles of CL lesions and healthy skin. Receiver operating characteristic (ROC) curves were made with plotROC package to calculate area under the curve (AUC) for potential markers of treatment outcome (56 (link)). Mann-Whitney U tests and Spearman’s rho (ρ) correlations were performed in the R statistical environment or GraphPrism v7, and P<0.05 was considered statistically significant.