Comparison of Whole Genome Sequencing versus Standard Molecular Diagnostics for Species Identification in the Leishmania Viannia Subgenus

The prognosis and treatment of New World tegumentary leishmaniasis is dependent on the infecting species, yet such species identification in the Leishmania Viannia subgenus poses a diagnostic challenge. Currently, speciation relies on standard molecular techniques such as restriction fragment length polymorphism (RFLP) analysis, and Sanger sequencing (SS). Whole-genome sequencing (WGS) is a robust and increasingly cost-efficient tool that may improve Leishmania species identification. We evaluated WGS versus standard RFLP-SS for species identification in three reference and five clinical strains of Leishmania Viannia spp. Internal transcribed spacer1 (its1), cysteine proteinase b (cpb), and heat shock protein 70 (hsp70) polymerase chain reaction-restriction fragment length polymorphism (RFLP) was performed, followed by SS of the its2, cpb, hsp70, and mannose phosphate isomerase (mpi) loci. After de novo assembly, sequences were mapped, and homology compared with both reference strains and reference genomes on National Center for Biotechnology Information. All American Type Culture Collection strains were confirmed to be single-species of L. V. braziliensis, L. V. guyanensis, or L. V. panamensis by WGS. Conversely, RFLP-SS was able to definitively identify one of three isolates to the species level. Clinical samples were identified as either single-species (N = 3), mixed (N = 1), or hybrid (N = 1) infections by WGS, while standard molecular diagnosis required multi-target composite analysis for identification due to loci-dependent results by RFLP-SS. We have corroborated the utility of WGS as a diagnostic tool to speciate members of the L. Viannia subgenus and to discriminate between mixed and hybrid infections. WGS is a potentially useful complement to multistaged RFLP-SS for species identification in Leishmania infections.