Selective whole-genome amplification reveals population genetics ofLeishmania braziliensisdirectly from patient skin biopsies

Abstract

ABSTRACTIn Brazil,Leishmania braziliensisis the main causative agent of the neglected tropical disease, cutaneous leishmaniasis (CL). CL presents on a spectrum of disease severity with a high rate of treatment failure. Yet the parasite factors that contribute to disease presentation and treatment outcome are not well understood, in part because successfully isolating and culturing parasites from patient lesions remains a major technical challenge. Here we describe the development of selective whole genome amplification (SWGA) forLeishmaniaand show that this method enables culture-independent analysis of parasite genomes obtained directly from primary patient skin samples, allowing us to circumvent artifacts associated with adaptation to culture. We show that SWGA can be applied to multipleLeishmaniaspecies residing in different host species, suggesting that this method is broadly useful in both experimental infection models and clinical studies. SWGA carried out directly on skin biopsies collected from patients in Corte de Pedra, Bahia, Brazil, showed extensive genomic diversity. Finally, as a proof-of-concept, we demonstrated that SWGA data can be integrated with published whole genome data from cultured parasite isolates to identify variants unique to specific geographic regions in Brazil where treatment failure rates are known to be high. SWGA provides a relatively simple method to generateLeishmaniagenomes directly from patient samples, unlocking the potential to link parasite genetics with host clinical phenotypes.AUTHOR SUMMARYLeishmania braziliensisis the main cause of cutaneous leishmaniasis in Brazil. Due to limitations in culturing, it is important to study the parasite in a culture-independent manner. We use selective whole genome amplification (SWGA) to explore parasite genomic diversity directly from patient biopsies. This method is inexpensive and can be broadly used to generate parasite genome sequence data sampled from differentLeishmaniaspecies infecting different mammalian hosts. We found high diversity among theL. braziliensisgenomes from Bahia, Brazil, which correlated with geographic location. By integrating these data with publicly available genome sequences from other studies spanning four countries in South America, we identified variants unique to Northeast Brazil that may be linked to high regional rates of treatment failure.