A Telomeric Cluster of Antimony Resistance Genes on Chromosome 34 of Leishmania infantum

CRISPR-Cas9 high-throughput screening to study drug resistance in Leishmania infantum

Parasites of the genus Leishmania pose a global health threat with limited treatment options. New drugs are urgently needed, and genomic screens have the potential to accelerate target discovery, mode of action, and resistance mechanisms against these new drugs. We describe here our effort in developing a genome-wide CRISPR-Cas9 screen in Leishmania, an organism lacking a functional nonhomologous end joining system that must rely on microhomology-mediated end joining, single-strand annealing, or homologous recombination for repairing Cas9-induced double-stranded DNA breaks. A new vector for cloning and expressing single guide RNAs (sgRNAs) was designed and proven to be effective in a small pilot project while enriching specific sgRNAs during drug selection. We then developed a whole-genome library of 49,754 sgRNAs, targeting all the genes of Leishmania infantum. This library was transfected in L. infantum expressing Cas9, and these cells were selected for resistance to two antileishmanials, miltefosine and amphotericin B. The sgRNAs the most enriched in the miltefosine screen targeted the miltefosine transporter gene, but sgRNAs targeting genes coding for a RING-variant protein and a transmembrane protein were also enriched. The sgRNAs the most enriched by amphotericin B targeted the sterol 24 C methyltransferase genes and a hypothetical gene. Through gene disruption experiments, we proved that loss of function of these genes was associated with resistance. This study describes the feasibility of carrying out whole-genome CRISPR-Cas9 screens in Leishmania provided that a strong selective pressure is applied. Such a screen can be used for accelerating the development of urgently needed antileishmanial drugs.IMPORTANCELeishmaniasis, a global health threat, lacks adequate treatment options and drug resistance exacerbates the challenge. This study introduces a CRISPR-Cas9 screening approach in Leishmania infantum, unraveling mechanisms of drug resistance at a genome-wide scale. Our screen was applied against two main antileishmanial drugs, and guides were enriched upon drug selection. These guides targeted known and new targets, hence validating the use of this screen against Leishmania. This strategy provides a powerful tool to expedite drug discovery as well as potential therapeutic targets against this neglected tropical disease.

Unlocking the power of AI models: exploring protein folding prediction through comparative analysis

Protein structure determination has made progress with the aid of deep learning models, enabling the prediction of protein folding from protein sequences. However, obtaining accurate predictions becomes essential in certain cases where the protein structure remains undescribed. This is particularly challenging when dealing with rare, diverse structures and complex sample preparation. Different metrics assess prediction reliability and offer insights into result strength, providing a comprehensive understanding of protein structure by combining different models. In a previous study, two proteins named ARM58 and ARM56 were investigated. These proteins contain four domains of unknown function and are present in Leishmania spp. ARM refers to an antimony resistance marker. The study's main objective is to assess the accuracy of the model's predictions, thereby providing insights into the complexities and supporting metrics underlying these findings. The analysis also extends to the comparison of predictions obtained from other species and organisms. Notably, one of these proteins shares an ortholog with Trypanosoma cruzi and Trypanosoma brucei, leading further significance to our analysis. This attempt underscored the importance of evaluating the diverse outputs from deep learning models, facilitating comparisons across different organisms and proteins. This becomes particularly pertinent in cases where no previous structural information is available.

Extracellular vesicles in parasitic diseases - from pathogenesis to future diagnostic tools

Parasitic diseases are still a major public health problem especially among individuals of low socioeconomic status in underdeveloped countries. In recent years it has been demonstrated that parasites can release extracellular vesicles that participate in the host-parasite communication, immune evasion, and in governing processes associated with host infection. Extracellular vesicles are membrane-bound structures released into the extracellular space that can carry several types of biomolecules, including proteins, lipids, nucleic acids, and metabolites, which directly impact the target cells. Extracellular vesicles have attracted wide attention due to their relevance in host-parasite communication and for their potential value in applications such as in the diagnostic biomarker discovery. This review of the literature aimed to join the current knowledge on the role of extracellular vesicles in host-parasite interaction and summarize its molecular content, providing information for the acquisition of new tools that can be used in the diagnosis of parasitic diseases. These findings shed light to the potential of extracellular vesicle cargo derived from protozoan parasites as novel diagnostic tools.

Antimony resistance and gene expression in Leishmania: spotlight on molecular and proteomic aspects

Leishmaniasis is a vector-borne parasitic disease caused by Leishmania parasites with a spectrum of clinical manifestations, ranging from skin lesions to severe visceral complications. Treatment of this infection has been extremely challenging with the concurrent emergence of drug resistance. The differential gene expression and the discrepancies in protein functions contribute to the appearance of 2 distinct phenotypes: resistant and sensitive, but the current diagnostic tools fail to differentiate between them. The identification of gene expression patterns and molecular mechanisms coupled with antimony (Sb) resistance can be leveraged to prompt diagnosis and select the most effective treatment methods. The present study attempts to use comparative expression of Sb resistance-associated genes in resistant and sensitive Leishmania, to disclose their relative abundance in clinical or in vitro selected isolates to gain an understanding of the molecular mechanisms of Sb response/resistance. Data suggest that the analysis of resistance gene expression would verify the Sb resistance or susceptibility only to a certain extent; however, none of the individual expression patterns of the studied genes was diagnostic as a biomarker of Sb response of Leishmania. The findings highlighted will be useful in bridging the knowledge gap and discovering innovative diagnostic tools and novel therapeutic targets.

The paradigm of intracellular parasite survival and drug resistance in leishmanial parasite through genome plasticity and epigenetics: Perception and future perspective

Leishmania is an intracellular, zoonotic, kinetoplastid eukaryote with more than 1.2 million cases all over the world. The leishmanial chromosomes are divided into polymorphic chromosomal ends, conserved central domains, and antigen-encoding genes found in telomere-proximal regions. The genome flexibility of chromosomal ends of the leishmanial parasite is known to cause drug resistance and intracellular survival through the evasion of host defense mechanisms. Therefore, in this review, we discuss the plasticity of Leishmania genome organization which is the primary cause of drug resistance and parasite survival. Moreover, we have not only elucidated the causes of such genome plasticity which includes aneuploidy, epigenetic factors, copy number variation (CNV), and post-translation modification (PTM) but also highlighted their impact on drug resistance and parasite survival.

Global distribution of treatment resistance gene markers for leishmaniasis

Background

Pentavalent antimonials (Sb(V)) such as meglumine antimoniate (Glucantime®) and sodium stibogluconate (Pentostam®) are used as first‐line treatments for leishmaniasis, either alone or in combination with second‐line drugs such as amphotericin B (Amp B), miltefosine (MIL), methotrexate (MTX), or cryotherapy. Therapeutic aspects of these drugs are now challenged because of clinical resistance worldwide.

Methods

We reviewedthe recent original studies were assessed by searching in electronic databases such as Scopus, Pubmed, Embase, and Web of Science.

Results

Studies on molecular biomarkers involved in drug resistance are essential for monitoring the disease. We reviewed genes and mechanisms of resistance to leishmaniasis, and the geographical distribution of these biomarkers in each country has also been thoroughly investigated.

Conclusion

Due to the emergence of resistant genes mainly in anthroponotic Leishmania species such as L. donovani and L. tropica, as the causative agents of ACL and AVL, respectively, selection of an appropriate treatment modality is essential. Physicians should be aware of the presence of such resistance for the selection of proper treatment modalities in endemic countries.

Experimental Strategies to Explore Drug Action and Resistance in Kinetoplastid Parasites

Kinetoplastids are the causative agents of leishmaniasis, human African trypanosomiasis, and American trypanosomiasis. They are responsible for high mortality and morbidity in (sub)tropical regions. Adequate treatment options are limited and have several drawbacks, such as toxicity, need for parenteral administration, and occurrence of treatment failure and drug resistance. Therefore, there is an urgency for the development of new drugs. Phenotypic screening already allowed the identification of promising new chemical entities with anti-kinetoplastid activity potential, but knowledge on their mode-of-action (MoA) is lacking due to the generally applied whole-cell based approach. However, identification of the drug target is essential to steer further drug discovery and development. Multiple complementary techniques have indeed been used for MoA elucidation. In this review, the different 'omics' approaches employed to define the MoA or mode-of-resistance of current reference drugs and some new anti-kinetoplastid compounds are discussed.

Coupling chemical mutagenesis to next generation sequencing for the identification of drug resistance mutations in Leishmania

Current genome-wide screens allow system-wide study of drug resistance but detecting small nucleotide variants (SNVs) is challenging. Here, we use chemical mutagenesis, drug selection and next generation sequencing to characterize miltefosine and paromomycin resistant clones of the parasite Leishmania. We highlight several genes involved in drug resistance by sequencing the genomes of 41 resistant clones and by concentrating on recurrent SNVs. We associate genes linked to lipid metabolism or to ribosome/translation functions with miltefosine or paromomycin resistance, respectively. We prove by allelic replacement and CRISPR-Cas9 gene-editing that the essential protein kinase CDPK1 is crucial for paromomycin resistance. We have linked CDPK1 in translation by functional interactome analysis, and provide evidence that CDPK1 contributes to antimonial resistance in the parasite. This screen is powerful in exploring networks of drug resistance in an organism with diploid to mosaic aneuploid genome, hence widening the scope of its applicability.

Here, Bhattacharya et al. chemically mutagenize Leishmania and identify genes associated with resistance to miltefosine and paromomycin by next generation sequencing. The study shows that a protein kinase (CDPK1) can mediate resistance to paromomycin by affecting translation.

Genomewide Analysis of Mode of Action of the S-Adenosylmethionine Analogue Sinefungin in Leishmania infantum

To further our understanding of one-carbon metabolism in the protozoan parasite Leishmania, we conducted genomic screens to study how the parasite responded to sinefungin (SNF) selection. SNF is a structural analogue of S-adenosylmethionine (AdoMet), a key methyl group donor to a number of biomolecules. One screen consisted of sequencing SNF-resistant mutants generated by stepwise selection with gradually increasing drug concentrations. These studies demonstrated deletion of the AdoMet transporter (AdoMetT1) by intergenic recombination as a crucial loss-of-function marker for SNF resistance. The second screen consisted of Cos-seq, a gain-of-function cosmid-based genomewide functional screen with increasing SNF concentration coupled to next-generation sequencing. Cosmids enriched in that screen and sequenced led to the identification of (i) the AdoMet synthetase (METK) as the major SNF target, (ii) an mRNA [(guanine-N7)-methyltransferase (CMT1)], (iii) a leucine carboxyl methyltransferase (LCMT), (iv) two tryparedoxin genes, and (v) two protein phosphatase regulatory genes. Further functional exploration indicated that LCMT interacts with one phosphatase catalytic subunit (PP2AC) and that mutation of the C-terminal leucine residue of PP2AC affects sinefungin susceptibility. These holistic screens led to the identification of transporters, biosynthetic genes, RNA and protein methyltransferases, as well as phosphatases linked to AdoMet-mediated functions in Leishmania IMPORTANCE The two main cellular metabolic one-carbon donors are reduced folates and S-adenosylmethionine, whose biosynthetic pathways have proven highly effective in chemotherapeutic interventions in various cell types. Sinefungin, a nucleoside analogue of S-adenosylmethionine, was shown to have potent activity against the protozoan parasite Leishmania Here, we studied resistance to sinefungin using whole-genome approaches as a way to further our understanding of the role of S-adenosylmethionine in this parasite and to reveal novel potential drug targets. These approaches allowed the characterization of novel features related to S-adenosylmethionine function in Leishmania which could further help in the development of sinefungin-like compounds against this pathogenic parasite.

Antimony resistance in Leishmania (Viannia) braziliensis clinical isolates from atypical lesions associates with increased ARM56/ARM58 transcripts and reduced drug uptake

BACKGROUND

In addition to the limited therapeutic arsenal and the side effects of antileishmanial agents, drug resistance hinders disease control. In Brazil, Leishmania braziliensis causes atypical (AT) tegumentary leishmaniasis lesions, frequently refractory to treatment.

OBJECTIVES

The main goal of this study was to characterise antimony (Sb)-resistant (SbR) L. braziliensis strains obtained from patients living in Xakriabá indigenous community, Minas Gerais, Brazil.

METHODS

The aquaglyceroporin 1-encoding gene (AQP1) from L. braziliensis clinical isolates was sequenced, and its function was evaluated by hypo-osmotic shock. mRNA levels of genes associated with Sb resistance were measured by quantitative reverse transcription polymerase chain reaction (qRT-PCR). Atomic absorption was used to measure Sb uptake.

FINDINGS

Although clinical isolates presented delayed recovery time in hypo-osmotic shock, AQP1 function was maintained. Isolate 340 accumulated less Sb than all other isolates, supporting the 65-fold downregulation of AQP1 mRNA levels. Both 330 and 340 isolates upregulated antimony resistance marker (ARM) 56/ARM58 and multidrug resistant protein A (MRPA); however, only ARM58 upregulation was an exclusive feature of SbR field isolates. CA7AE seemed to increase drug uptake in L. braziliensis and represented a tool to study the role of glycoconjugates in Sb transport.

MAIN CONCLUSIONS

There is a clear correlation between ARM56/58 upregulation and Sb resistance in AT-harbouring patients, suggesting the use of these markers as potential indicators to help the treatment choice and outcome, preventing therapeutic failure.

A single amino acid substitution (H451Y) in Leishmania calcium-dependent kinase SCAMK confers high tolerance and resistance to antimony

Background

For almost a century, antimonials have remained the first-line drugs for the treatment of leishmaniasis. However, little is known about their mode of action and clinical resistance mechanisms.

Objectives

We have previously shown that Leishmania nicotinamidase (PNC1) is an essential enzyme for parasite NAD+ homeostasis and virulence in vivo. Here, we found that parasites lacking the pnc1 gene (Δpnc1) are hypersusceptible to the active form of antimony (SbIII) and used these mutant parasites to better understand antimony’s mode of action and the mechanisms leading to resistance.

Methods

SbIII-resistant WT and Δpnc1 parasites were selected in vitro by a stepwise selection method. NAD(H)/NADP(H) dosages and quantitative RT–PCR experiments were performed to explain the susceptibility differences observed between strains. WGS and a marker-free CRISPR/Cas9 base-editing approach were used to identify and validate the role of a new resistance mutation.

Results

NAD+-depleted Δpnc1 parasites were highly susceptible to SbIII and this phenotype could be rescued by NAD+ precursor or trypanothione precursor supplementation. Δpnc1 parasites could become resistant to SbIII by an unknown mechanism. WGS revealed a unique amino acid substitution (H451Y) in an EF-hand domain of an orphan calcium-dependent kinase, recently named SCAMK. When introduced into a WT reference strain by base editing, the H451Y mutation allowed Leishmania parasites to survive at extreme concentrations of SbIII, potentiating the rapid emergence of resistant parasites.

Conclusions

These results establish that Leishmania SCAMK is a new central hub of antimony’s mode of action and resistance development, and uncover the importance of drug tolerance mutations in the evolution of parasite drug resistance.

Gene Replacement by Homologous Recombination

While homologous recombination-based gene replacement is about to be supplanted by more modern approaches, it is still retaining usefulness for genes that prove to be poor targets for CRISPR/cas-based approaches. Homologous recombination has proven to be relatively robust to minor sequence mismatches between GOI-flanking sequences and the gene replacement constructs, and the faithfulness of recombination events is easily verified by whole-genome sequencing. Moreover, the availability of custom synthetic gene production by numerous service providers should allow for a relatively quick generation of null mutants without the need to introduce additional protein-coding genes beyond the selection markers.

Extracellular Vesicle-Mediated Communication Within Host-Parasite Interactions

Extracellular vesicles (EVs) are small membrane-surrounded structures released by different kinds of cells (normal, diseased, and transformed cells) in vivo and in vitro that contain large amounts of important substances (such as lipids, proteins, metabolites, DNA, RNA, and non-coding RNA (ncRNA), including miRNA, lncRNA, tRNA, rRNA, snoRNA, and scaRNA) in an evolutionarily conserved manner. EVs, including exosomes, play a role in the transmission of information, and substances between cells that is increasingly being recognized as important. In some infectious diseases such as parasitic diseases, EVs have emerged as a ubiquitous mechanism for mediating communication during host-parasite interactions. EVs can enable multiple modes to transfer virulence factors and effector molecules from parasites to hosts, thereby regulating host gene expression, and immune responses and, consequently, mediating the pathogenic process, which has made us rethink our understanding of the host-parasite interface. Thus, here, we review the present findings regarding EVs (especially exosomes) and recognize the role of EVs in host-parasite interactions. We hope that a better understanding of the mechanisms of parasite-derived EVs may provide new insights for further diagnostic biomarker, vaccine, and therapeutic development.

Reverse Epidemiology: An Experimental Framework to Drive Leishmania Biomarker Discovery in situ by Functional Genetic Screening Using Relevant Animal Models

Leishmania biomarker discovery remains an important challenge that needs to be revisited in light of our increasing knowledge on parasite-specific biology, notably its genome instability. In the absence of classical transcriptional regulation in these early-branching eukaryotes, fluctuations in transcript abundance can be generated by gene and chromosome amplifications, which have been linked to parasite phenotypic variability with respect to virulence, tissue tropism, and drug resistance. Conducting in vitro evolutionary experiments to study mechanisms of Leishmania environmental adaptation, we recently validated the link between parasite genetic amplification and fitness gain, thus defining gene and chromosome copy number variations (CNVs) as important Leishmania biomarkers. These experiments also demonstrated that long-term Leishmania culture adaptation can strongly interfere with epidemiologically relevant, genetic signals, which challenges current protocols for biomarker discovery, all of which rely on in vitro expansion of clinical isolates. Here we propose an experimental framework independent of long-term culture termed “reverse” epidemiology, which applies established protocols for functional genetic screening of cosmid-transfected parasites in animal models for the identification of clinically relevant genetic loci that then inform targeted field studies for their validation as Leishmania biomarkers.

High-throughput Cos-Seq screen with intracellular Leishmania infantum for the discovery of novel drug-resistance mechanisms

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.

•

The Cos-Seq led to the discovery of several new genomic regions selected with drugs.

•

This work led to the validation of novel drug-resistance genes in Leishmania.

•

Gene LinJ.06.1010 is involved in antimony resistance in both life stages of the parasite.

•

Gene LinJ.32.0050 is involved in miltefosine resistance in amastigotes.

•

The amastigote screen is labour intensive but complements screens in promastigotes.

Hsp90 inhibitors radicicol and geldanamycin have opposing effects on Leishmania Aha1-dependent proliferation

Hsp90 and its co-chaperones are essential for the medically important parasite Leishmania donovani, facilitating life cycle control and intracellular survival. Activity of Hsp90 is regulated by co-chaperones of the Aha1 and P23 families. In this paper, we studied the expression of L. donovani Aha1 in two life cycle stages, its interaction with Hsp90 and the phenotype of Aha1 null mutants during the insect stage and inside infected macrophages. This study provides a detailed in vitro analysis of the function of Aha1 in Leishmania parasites and the first instance of a reverse genetic analysis of Aha1 in a protozoan parasite. While Aha1 is non-essential under standard growth conditions and at elevated temperature, Aha1 protects against ethanol stress. However, both overexpression and lack of Aha1 affected parasite growth in the presence of the Hsp90 inhibitors radicicol (RAD) and geldanamycin (GA). Under RAD pressure, P23 and Aha1 act in an antagonistic way. By contrast, expression levels of both co-chaperones have similar effects under GA treatment, indicating different inhibition mechanisms by the two compounds. Aha1 is also secreted in virulence-enhancing exosomes. This may explain why the loss of Aha1 reduces the infectivity of L. donovani in ex vivo mouse macrophages, indicating a role during the intracellular mammalian stage.

Biophysical and Pharmacological Characterization of Energy-Dependent Efflux of Sb in Laboratory-Selected Resistant Strains of Leishmania (Viannia) Subgenus

The growing resistance of leishmaniasis to first-line drugs like antimonials in some regions limits the control of this parasitic disease. The precise mechanisms involved in Leishmania antimony resistance are still subject to debate. The reduction of intracellular Sb^III accumulation is a common change observed in both laboratory-selected and field isolated resistant Leishmania strains, but the exact transport pathways involved in antimony resistance have not yet been elucidated. In order to functionally characterize the antimony transport routes responsible for resistance, we performed systematic transport studies of Sb^III in wild-type and resistant strains of L. (Viannia) guyanensis and L. (V.) braziliensis. Those include influx and efflux assays and the influence of ABC transporters and metabolism inhibitors: prochlorperazine, probenecid, verapamil, BSO, and sodium azide. The mRNA levels of genes associated with antimony resistance (MRPA, GSH1, ODC, AQP1, ABCI4, and ARM58) were also investigated in addition to intracellular thiol levels. A strong reduction of Sb influx was observed in L. guyanensis resistant mutant (LgSbR), but not in L. braziliensis (LbSbR). Both mutants showed increased energy-dependent efflux of Sb^III, when compared to their respective parental strains. In LgSbR, BSO and prochlorperazine inhibited antimony efflux and resistance was associated with increased MRPA and GSH1 mRNA levels, while in LbSbR antimony efflux was inhibited by probenicid and prochlorperazine in absence of resistance-associated gene modulation. Intracellular thiol levels were increased in both Sb-resistant mutants. An energy-dependent Sb^III efflux pathway sensitive to prochlorperazine was clearly evidenced in both Sb-resistant mutants. In conclusion, the present study allowed the biophysical and pharmacological characterization of energy-dependent Sb efflux pathway apparently independent of MRPA, ABCI4, and ARM58 upregulation, in Leishmania (Vianna) mutant selected in vitro for resistance to Sb^III. Prochlorperazine has also been identified as an effective chemosensitizer in both Sb resistant mutants, which acts through inhibition of the active efflux of Sb.

Plasticity of the Leishmania genome leading to gene copy number variations and drug resistance

Leishmania has a plastic genome, and drug pressure can select for gene copy number variation (CNV). CNVs can apply either to whole chromosomes, leading to aneuploidy, or to specific genomic regions. For the latter, the amplification of chromosomal regions occurs at the level of homologous direct or inverted repeated sequences leading to extrachromosomal circular or linear amplified DNAs. This ability of Leishmania to respond to drug pressure by CNVs has led to the development of genomic screens such as Cos-Seq, which has the potential of expediting the discovery of drug targets for novel promising drug candidates.