Abstract
Increasing drug resistance towards first line antimony-derived compounds has forced the introduction of novel therapies in leishmaniasis endemic areas including amphotericin B and miltefosine. However, their use is threatened by the emergence and spread of drug-resistant strains. In order to discover stage-dependent resistance genes, we have adapted the Cos-Seq approach through the introduction of macrophage infections in the pipeline. A L. infantum intracellular amastigote population complemented with a L. infantum cosmid library was submitted to increasing concentrations of miltefosine, amphotericin B and pentavalent antimonials in experimental infections of THP-1 cells. For each step of selection, amastigotes were extracted and cosmids were isolated and submitted to next-generation sequencing, followed by subsequent gene-enrichment analyses. Cos-Seq screen in amastigotes revealed four highly enriched loci for antimony, five for miltefosine and one for amphotericin B. Of these, a total of seven cosmids were recovered and tested for resistance in both promastigotes and amastigotes. Candidate genes within the pinpointed genomic regions were validated using single gene overexpression in wild-type parasites and/or gene disruption by means of a CRISPR-Cas9-based approach. This led to the identification and validation of a stage-independent antimony-resistance gene (LinJ.06.1010) coding for a putative leucine rich repeat protein and a novel amastigote-specific miltefosine-resistance gene (LinJ.32.0050) coding for a member of the SEC13 family of WD-repeat proteins. This study further reinforces the power of Cos-Seq approach to discover novel drug-resistance genes, some of which are life-stages specific.