The use of next generation sequencing has the power to expedite the identification of drug resistance determinants and biomarkers and was applied successfully to drug resistance studies in Leishmania. This allowed the identification of modulation in gene expression, gene dosage alterations, changes in chromosome copy numbers and single nucleotide polymorphisms that correlated with resistance in Leishmania strains derived from the laboratory and from the field. An impressive heterogeneity at the population level was also observed, individual clones within populations often differing in both genotypes and phenotypes, hence complicating the elucidation of resistance mechanisms. This review summarizes the most recent highlights that whole genome sequencing brought to our understanding of Leishmania drug resistance and likely new directions.

Drug resistance is a major problem in leishmaniasis chemotherapy. RNA expression profiling using DNA microarrays is a suitable approach to study simultaneous events leading to a drug-resistance phenotype. Genomic analysis has been performed primarily with Old World Leishmania species and here we investigate molecular alterations in antimony resistance in the New World species L. amazonensis. We selected populations of L. amazonensis promastigotes for resistance to antimony by step-wise drug pressure. Gene expression of highly resistant mutants was studied using DNA microarrays. RNA expression profiling of antimony-resistant L. amazonensis revealed the overexpression of genes involved in drug resistance including the ABC transporter MRPA and several genes related to thiol metabolism. The MRPA overexpression was validated by quantitative real-time RT-PCR and further analysis revealed that this increased expression was correlated to gene amplification as part of extrachromosomal linear amplicons in some mutants and as part of supernumerary chromosomes in other mutants. The expression of several other genes encoding hypothetical proteins but also nucleobase and glucose transporter encoding genes were found to be modulated. Mechanisms classically found in Old World antimony resistant Leishmania were also highlighted in New World antimony-resistant L. amazonensis. These studies were useful to the identification of resistance molecular markers. Leishmania are unicellular microorganisms that can be transmitted to humans by the bite of sandflies. They cause a spectrum of diseases called leishmaniasis, which are classified as neglected tropical diseases by the World Health Organization. The treatment of leishmaniasis is based on the administration of antimony-containing drugs. These drugs have been used since 1947 and still constitute the mainstay for leishmaniasis treatment in several countries. One of the problems with these compounds is the emergence of resistance. Our work seeks to understand how these parasites become resistant to the drug. We studied antimony-resistant Leishmania amazonensis mutants. We analyzed gene expression at the whole genome level in antimony-resistant parasites and identified mechanisms used by Leishmania for resistance. This work could help us in developing new strategies for treatment in endemic countries where people are unresponsive to antimony-based chemotherapy. The identification of common mechanisms among different species of resistant parasites may also contribute to the development of diagnostic kits to identify and monitor the spread of resistance.

Antimony resistance complicates the treatment of infections caused by the parasite Leishmania. Using next generation sequencing, we sequenced the genome of four independent Leishmania guyanensis antimony-resistant (SbR) mutants and found different chromosomal alterations including aneuploidy, intrachromosomal gene amplification and gene deletion. A segment covering 30 genes on chromosome 19 was amplified intrachromosomally in three of the four mutants. The gene coding for the multidrug resistance associated protein A involved in antimony resistance was also amplified in the four mutants, most likely through chromosomal translocation. All mutants also displayed a reduced accumulation of antimony mainly due to genomic alterations at the level of the subtelomeric region of chromosome 31 harboring the gene coding for the aquaglyceroporin 1 (LgAQP1). Resistance involved the loss of LgAQP1 through subtelomeric deletions in three mutants. Interestingly, the fourth mutant harbored a single G133D point mutation in LgAQP1 whose role in resistance was functionality confirmed through drug sensitivity and antimony accumulation assays. In contrast to the Leishmania subspecies that resort to extrachromosomal amplification, the Viannia strains studied here used intrachromosomal amplification and locus deletion. This is the first report of a naturally occurred point mutation in AQP1 in antimony resistant parasites. Drug resistance remains a major concern in leishmaniasis chemotherapy, a neglected tropical disease that causes 60,000 deaths around the world annually. To better understand the molecular mechanisms behind drug resistance, we selected L. guyanensis parasites resistant to antimony, the first-line drug against this disease in many countries. Through whole-genome sequencing we found variations in the copy number of chromosomes in addition to gene amplification and gene deletion events in antimony-resistant parasites. A marker previously related to antimony resistance, the gene coding for multidrug resistant protein A, was found to be amplified. Transport studies revealed a reduced antimony accumulation in resistant parasites that we correlated with the deletion of the gene coding for the aquaglyceroporin 1 (AQP1) responsible for antimony uptake in Leishmania. Additionally, a point mutation in AQP1 was found to be associated with antimony resistance. These findings may contribute to the development of new chemotherapy strategies against leishmaniasis.