Genome-wide stochastic adaptive DNA amplification at direct and inverted DNA repeats in the parasite Leishmania.

Gene amplification of specific loci has been described in all kingdoms of life. In the protozoan parasite Leishmania, the product of amplification is usually part of extrachromosomal circular or linear amplicons that are formed at the level of direct or inverted repeated sequences. A bioinformatics screen revealed that repeated sequences are widely distributed in the Leishmania genome and the repeats are chromosome-specific, conserved among species, and generally present in low copy number. Using sensitive PCR assays, we provide evidence that the Leishmania genome is continuously being rearranged at the level of these repeated sequences, which serve as a functional platform for constitutive and stochastic amplification (and deletion) of genomic segments in the population. This process is adaptive as the copy number of advantageous extrachromosomal circular or linear elements increases upon selective pressure and is reversible when selection is removed. We also provide mechanistic insights on the formation of circular and linear amplicons through RAD51 recombinase-dependent and -independent mechanisms, respectively. The whole genome of Leishmania is thus stochastically rearranged at the level of repeated sequences, and the selection of parasite subpopulations with changes in the copy number of specific loci is used as a strategy to respond to a changing environment.

Roles of Rad51 paralogs for promoting homologous recombination in Leishmania infantum.

To achieve drug resistance Leishmania parasite alters gene copy number by using its repeated sequences widely distributed through the genome. Even though homologous recombination (HR) is ascribed to maintain genome stability, this eukaryote exploits this potent mechanism driven by the Rad51 recombinase to form beneficial extrachromosomal circular amplicons. Here, we provide insights on the formation of these circular amplicons by analyzing the functions of the Rad51 paralogs. We purified three Leishmania infantum Rad51 paralogs homologs (LiRad51-3, LiRad51-4 and LiRad51-6) all of which directly interact with LiRad51. LiRad51-3, LiRad51-4 and LiRad51-6 show differences in DNA binding and annealing capacities. Moreover, it is also noteworthy that LiRad51-3 and LiRad51-4 are able to stimulate Rad51-mediated D-loop formation. In addition, we succeed to inactivate the LiRad51-4 gene and report a decrease of circular amplicons in this mutant. The LiRad51-3 gene was found to be essential for cell viability. Thus, we propose that the LiRad51 paralogs play crucial functions in extrachromosomal circular DNA amplification to circumvent drug actions and preserve survival.

Interactions between BRCA2 and RAD51 for promoting homologous recombination in Leishmania infantum.

In most organisms, the primary function of homologous recombination (HR) is to allow genome protection by the faithful repair of DNA double-strand breaks. The vital step of HR is the search for sequence homology, mediated by the RAD51 recombinase, which is stimulated further by proteins mediators such as the tumor suppressor BRCA2. The biochemical interplay between RAD51 and BRCA2 is unknown in Leishmania or Trypanosoma. Here we show that the Leishmania infantum BRCA2 protein possesses several critical features important for the regulation of DNA recombination at the genetic and biochemical level. A BRCA2 null mutant, generated by gene disruption, displayed genomic instability and gene-targeting defects. Furthermore, cytological studies show that LiRAD51 can no longer localize to the nucleus in this mutant. The Leishmania RAD51 and BRCA2 interact together and the purified proteins bind single-strand DNA. Remarkably, LiBRCA2 is a recombination mediator that stimulates the invasion of a resected DNA double-strand break in an undamaged template by LiRAD51 to form a D-loop structure. Collectively, our data show that LiBRCA2 and LiRAD51 promote HR at the genetic and biochemical level in L. infantum, the causative agent of visceral leishmaniasis.

Genome sequencing of the lizard parasite Leishmania tarentolae reveals loss of genes associated to the intracellular stage of human pathogenic species.

The Leishmania tarentolae Parrot-TarII strain genome sequence was resolved to an average 16-fold mean coverage by next-generation DNA sequencing technologies. This is the first non-pathogenic to humans kinetoplastid protozoan genome to be described thus providing an opportunity for comparison with the completed genomes of pathogenic Leishmania species. A high synteny was observed between all sequenced Leishmania species. A limited number of chromosomal regions diverged between L. tarentolae and L. infantum, while remaining syntenic to L. major. Globally, >90% of the L. tarentolae gene content was shared with the other Leishmania species. We identified 95 predicted coding sequences unique to L. tarentolae and 250 genes that were absent from L. tarentolae. Interestingly, many of the latter genes were expressed in the intracellular amastigote stage of pathogenic species. In addition, genes coding for products involved in antioxidant defence or participating in vesicular-mediated protein transport were underrepresented in L. tarentolae. In contrast to other Leishmania genomes, two gene families were expanded in L. tarentolae, namely the zinc metallo-peptidase surface glycoprotein GP63 and the promastigote surface antigen PSA31C. Overall, L. tarentolae's gene content appears better adapted to the promastigote insect stage rather than the amastigote mammalian stage.

Decreased glutamate transport in acivicin resistant Leishmania tarentolae.

Studies of drug resistance in the protozoan parasites of the genus Leishmania have been helpful in revealing biochemical pathways as potential drug targets. The chlorinated glutamine analogue acivicin has shown good activity against Leishmania cells and was shown to target several enzymes containing amidotransferase domains. We selected a Leishmania tarentolae clone for acivicin resistance. The genome of this resistant strain was sequenced and the gene coding for the amidotransferase domain-containing GMP synthase was found to be amplified. Episomal expression of this gene in wild-type L. tarentolae revealed a modest role in acivicin resistance. The most prominent defect observed in the resistant mutant was reduced uptake of glutamate, and through competition experiments we determined that glutamate and acivicin, but not glutamine, share the same transporter. Several amino acid transporters (AATs) were either deleted or mutated in the resistant cells. Some contributed to the acivicin resistance phenotype although none corresponded to the main glutamate transporter. Through sequence analysis one AAT on chromosome 22 corresponded to the main glutamate transporter. Episomal expression of the gene coding for this transporter in the resistant mutant restored glutamate transport and acivicin susceptibility. Its genetic knockout led to reduced glutamate transport and acivicin resistance. We propose that acivicin binds covalently to this transporter and as such leads to decreased transport of glutamate and acivicin thus leading to acivicin resistance.

The gamma-glutamylcysteine synthetase gene of Leishmania is essential and involved in response to oxidants.

Gamma-glutamylcysteine synthetase, encoded by the GSH1 gene, is the rate-limiting enzyme in the biosynthesis of glutathione and of trypanothione in Leishmania. The importance of GSH1 was assessed by generating GSH1 null mutants in Leishmania infantum. Removal of even a single wild-type allelic copy of GSH1 invariably led to the generation of an extra copy of GSH1, maintaining two intact wild-type alleles. However, by first supplementing the parasites with a rescue plasmid, we succeeded in obtaining both a single and null chromosomal GSH1 mutants. Parasites with one intact GSH1 chromosomal allele lost the rescuing plasmid but not the double knockout, when grown in the absence of antibiotic, indicating the essentiality of the GSH1 gene in Leishmania. Heterozygous mutants with one allele-inactivated transcribed less GSH1 mRNA and synthesized less glutathione and trypanothione. These mutants were more susceptible to oxidative stresses in vitro as promastigotes and showed decreased survival inside activated macrophages producing reactive oxygen or nitrogen species. These mutants showed a significant decreased survival in the presence of antimony (SbV) compared with control cells. All phenotypes were reverted in the add-back mutant, thus proving the importance of thiols in dealing with oxidants including the action of antimonials.

MRPA-independent mechanisms of antimony resistance in Leishmania infantum.

Control of both human and canine leishmaniasis is based on a very short list of chemotherapeutic agents, headed by antimonial derivatives (Sb). The utility of these molecules is severely threatened by high rates of drug resistance. The ABC transporter MRPA is one of the few key Sb resistance proteins described to date, whose role in detoxification has been thoroughly studied in Leishmania parasites. Nonetheless, its rapid amplification during drug selection complicates the discovery of other mechanisms potentially involved in Sb resistance. In this study, stepwise drug-resistance selection and next-generation sequencing were combined in the search for novel Sb-resistance mechanisms deployed by parasites when MRPA is abolished by targeted gene disruption. The gene mrpA is not essential in L. infantum, and its disruption leads to an Sb hypersensitive phenotype in both promastigotes and amastigotes. Five independent mrpA-/- mutants were selected for antimony resistance. These mutants displayed major changes in their ploidy, as well as extrachromosomal linear amplifications of the subtelomeric region of chromosome 23, which includes the genes coding for ABCC1 and ABCC2. Overexpression of ABCC2, but not of ABCC1, resulted in increased Sb tolerance in the mrpA-/- mutant. SNP analyses revealed three different heterozygous mutations in the gene coding for a serine acetyltransferase (SAT) involved in de novo cysteine synthesis in Leishmania. Overexpression of satQ390K, satG321R and satG325R variants led to a 2-3.2 -fold increase in Sb resistance in mrpA-/- parasites. Only satG321R and satG325R induced increased Sb resistance in wild-type parasites. These results reinforce and expand knowledge on the complex nature of Sb resistance in Leishmania parasites.

Down regulation of KMP-11 in Leishmania infantum axenic antimony resistant amastigotes as revealed by a proteomic screen.

The therapeutic mainstay against the protozoan parasite Leishmania is still based on the antiquated pentavalent antimonials, but resistance is increasing in several parts of the world. Resistance is now partly understood in laboratory promastigote isolates, but the mechanism leading to drug resistance in amastigote isolates is lagging behind. Here we describe a comparative proteomic analysis of a genetically related pair of antimonial-sensitive and -resistant Leishmania infantum axenic amastigote strains. The proteomics screen has highlighted a number of proteins differentially expressed in the resistant parasite. The expression of the protein argininosuccinate synthetase (ARGG) was increased in the drug resistant mutant while a decrease in the expression of the kinetoplastid membrane protein (KMP-11) correlated with the drug resistance phenotype. This proteomic screen highlighted several novel proteins that are putatively involved in resistance to antimonials.

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.

A protein of the leucine-rich repeats (LRRs) superfamily is implicated in antimony resistance in Leishmania infantum amastigotes.

Pentavalent antimonial containing drugs (SbV) are the mainstay for the control of the protozoan parasite Leishmania but resistance to this class of drug is now prevalent in several endemic areas. We describe here the use of functional cloning where an expression cosmid bank derived from Leishmania infantum was transfected in L. infantum axenic amastigotes and selected for potassium antimonyl tartrate (SbIII) resistance. This strategy allowed the isolation of a cosmid encoding for a novel resistance protein, LinJ34.0570, which belongs to the superfamily of leucine-rich repeat (LRR) proteins. Parasites overexpressing this LRR protein, which is part of the LRR_CC subfamily, were resistant to SbIII as axenic amastigotes and to SbV as intracellular parasites. This work pinpoints a novel protein that can contribute to antimonial resistance in Leishmania.

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.

Growth arrested live-attenuated Leishmania infantum KHARON1 null mutants display cytokinesis defect and protective immunity in mice.

There is no safe and efficacious vaccine against human leishmaniasis available and live attenuated vaccines have been used as a prophylactic alternative against the disease. In order to obtain an attenuated Leishmania parasite for vaccine purposes, we generated L. infantum KHARON1 (KH1) null mutants (ΔLikh1). This gene was previously associated with growth defects in L. mexicana. ΔLikh1 was obtained and confirmed by PCR, qPCR and Southern blot. We also generate a KH1 complemented line with the introduction of episomal copies of KH1. Although ΔLikh1 promastigote forms exhibited a growth pattern similar to the wild-type line, they differ in morphology without affecting parasite viability. L. infantum KH1-deficient amastigotes were unable to sustain experimental infection in macrophages, forming multinucleate cells which was confirmed by in vivo attenuation phenotype. The cell cycle analysis of ΔLikh1 amastigotes showed arrested cells at G2/M phase. ΔLikh1-immunized mice presented reduced parasite burden upon challenging with virulent L. infantum, when compared to naïve mice. An effect associated with increased Li SLA-specific IgG serum levels and IL-17 production. Thus, ΔLikh1 parasites present an infective-attenuated phenotype due to a cytokinesis defect, whereas it induces immunity against visceral leishmaniasis in mouse model, being a candidate for antileishmanial vaccine purposes.

Different Mutations in a P-type ATPase Transporter in Leishmania Parasites are Associated with Cross-resistance to Two Leading Drugs by Distinct Mechanisms.

Leishmania infantum is an etiological agent of the life-threatening visceral form of leishmaniasis. Liposomal amphotericin B (AmB) followed by a short administration of miltefosine (MF) is a drug combination effective for treating visceral leishmaniasis in endemic regions of India. Resistance to MF can be due to point mutations in the miltefosine transporter (MT). Here we show that mutations in MT are also observed in Leishmania AmB-resistant mutants. The MF-induced MT mutations, but not the AmB induced mutations in MT, alter the translocation/uptake of MF. Moreover, mutations in the MT selected by AmB or MF have a major impact on lipid species that is linked to cross-resistance between both drugs. These alterations include changes of specific phospholipids, some of which are enriched with cyclopropanated fatty acids, as well as an increase in inositolphosphoceramide species. Collectively these results provide evidence of the risk of cross-resistance emergence derived from current AmB-MF sequential or co-treatments for visceral leishmaniasis.

Inactivation of the cytosolic and mitochondrial serine hydroxymethyl transferase genes in Leishmania major.

Leishmania has two serine hydroxylmethyl transferase (SHMT) genes, one coding for a cytosolic and the other for a mitochondrial enzyme. Trypanosoma cruzi has only the gene coding for the cytosolic enzyme and Trypanosoma brucei has no SHMT. We tested whether these genes were dispensable for growth in Leishmania major. By gene inactivation we succeeded in generating three cells lines one without the cytosolic cSHMT, one without the mitochondrial mSHMT, and finally one L. major line without any SHMT. SHMT is thus dispensable for growth of Leishmania in rich medium. The ability of the various shmt null mutants to grow in defined medium was tested and the growth of the shmt null mutant was dependent on the presence of serine. Overall this work has shown that SHMT is dispensable for Leishmania growth but it may be necessary when growing in environments poor in serine.

Pterin transport and metabolism in Leishmania and related trypanosomatid parasites.

The folate metabolic pathway has been exploited successfully for the development of antimicrobial and antineoplasic agents. Inhibitors of this pathway, however, are not useful against Leishmania and other trypanosomatids. Work on the mechanism of methotrexate resistance in Leishmania has dramatically increased our understanding of folate and pterin metabolism in this organism. The metabolic and cellular functions of the reduced form of folates and pterins are beginning to be established and this work has led to several unexpected findings. Moreover, the currently ongoing sequencing efforts on trypanosomatid genomes are suggesting the presence of several gene products that are likely to require folates and pterins. A number of the properties of folate and pterin metabolism are unique suggesting that these pathways are valid and worthwhile targets for drug development.

Genetic polymorphisms and drug susceptibility in four isolates of Leishmania tropica obtained from Canadian soldiers returning from Afghanistan.

BACKGROUND: Cutaneous leishmaniasis (CL) is a vector-borne parasitic disease characterized by the presence of one or more lesions on the skin that usually heal spontaneously after a few months. Most cases of CL worldwide occur in Southwest Asia, Africa and South America, and a number of cases have been reported among troops deployed to Afghanistan. No vaccines are available against this disease, and its treatment relies on chemotherapy. The aim of this study was to characterize parasites isolated from Canadian soldiers at the molecular level and to determine their susceptibility profile against a panel of antileishmanials to identify appropriate therapies. METHODOLOGY/PRINCIPAL FINDINGS: Parasites were isolated from skin lesions and characterized as Leishmania tropica based on their pulsed field gel electrophoresis profiles and pteridine reductase 1 (PTR1) sequences. Unusually high allelic polymorphisms were observed at several genetic loci for the L. tropica isolates that were characterized. The drug susceptibility profile of intracellular amastigote parasites was determined using an established macrophage assay. All isolates were sensitive to miltefosine, amphotericin B, sodium stibogluconate (Pentostam) and paromomycin, but were not susceptible to fluconazole. Variable levels of susceptibility were observed for the antimalarial agent atovaquone/proguanil (Malarone). Three Canadian soldiers from this study were successfully treated with miltefosine. CONCLUSIONS/SIGNIFICANCE: This study shows high heterogeneity between the two L. tropica allelic versions of a gene but despite this, L. tropica isolated from Afghanistan are susceptible to several of the antileishmanial drugs available.

Discovery of factors linked to antimony resistance in Leishmania panamensis through differential proteome analysis.

The rate of treatment failure to antileishmanial chemotherapy in Latin America is up to 64%. Parasite drug resistance contributes to an unknown proportion of treatment failures. Identification of clinically relevant molecular mechanisms responsible for parasite drug resistance is critical to the conservation of available drugs and to the discovery of novel targets to reverse the resistant phenotype. We conducted comparative proteomic-based analysis of Leishmania (Viannia) panamensis lines selected in vitro for resistance to trivalent antimony (Sb(III)) to identify factors associated with antimony resistance. Using 2-dimensional gel electrophoresis, two distinct sub-proteomes (soluble in NP-40/urea and Triton X-114, respectively) of promastigotes of WT and Sb(III)-resistant lines were generated. Overall, 9 differentially expressed putative Sb-resistance factors were detected and identified by mass spectrometry. These constituted two major groups: (a) proteins involved in general stress responses and (b) proteins with highly specific metabolic and transport functions, potentially directly contributing to the Sb-resistance mechanism. Notably, the sulfur amino acid-metabolizing enzymes S-adenosylmethionine synthetase (SAMS) and S-adenosylhomocysteine hydrolase (SAHH) were over-expressed in Sb(III)-resistant lines and Sb(III)-resistant clinical isolates. These enzymes play a central role in the upstream synthesis of precursors of trypanothione, a key molecule involved in Sb-resistance in Leishmania parasites, and suggest involvement of epigenetic regulation in response to drug exposure. These data re-enforce the importance of thiol metabolism in Leishmania Sb resistance, reveal previously unrecognized steps in the mechanism(s) of Sb tolerance, and suggest a cross-talk between drug resistance, metabolism and virulence.

Modulation of gene expression in human macrophages treated with the anti-leishmania pentavalent antimonial drug sodium stibogluconate.

Within the mammalian host, Leishmania donovani is an obligatory intracellular protozoan parasite that resides and multiplies exclusively in the phagolysosomes of macrophages. Leishmania control relies primarily on chemotherapy, with the mainstay being pentavalent antimony (SbV) complexed to carbohydrates in the form of sodium stibogluconate (Pentostam) or meglumine antimoniate (Glucantime). The mode of action of SbV is still not known precisely. To explore the effect of SbV on macrophage gene expression, a microarray analysis was performed using Affymetrix focus arrays to compare gene expression profiles in noninfected and L. donovani-infected THP-1 monocytic cells treated or not treated with sodium stibogluconate. Under our experimental conditions, SbV changed the expression of a few host genes, and this was independent of whether cells were infected or not infected with Leishmania. Leishmania infection had a greater effect on the modulation of host gene expression. Statistical analyses have indicated that the expression of eight genes was modified by at least twofold upon SbV treatment, with six genes upregulated and two genes downregulated. One gene whose expression was affected by SbV was the heme oxygenase gene HMOX-1, and this change was observed both in the monocytic cell line THP-1 and in primary human monocyte-derived macrophages. Another pathway that was affected was the glutathione biosynthesis pathway, where the expression of the glutamate-cysteine ligase modifier subunit was increased upon SbV treatment. Our analysis has suggested that, under our experimental conditions, the expression of a few genes is altered upon SbV treatment, and some of these encoded proteins may be implicated in the yet-to-be-defined mode of action of SbV.

Role of ABC transporter MRPA, gamma-glutamylcysteine synthetase and ornithine decarboxylase in natural antimony-resistant isolates of Leishmania donovani.

OBJECTIVES: The resistance of clinical isolates of Leishmania donovani to sodium antimony gluconate (SAG), the mainstay of treatment in Indian visceral leishmaniasis, has become a critical issue in India. The present work investigates the mechanism of resistance to SAG in parasites isolated from patients who are unresponsive to SAG. METHODS AND RESULTS: Susceptibility to SAG as determined in vitro with intracellular amastigotes correlated well with the clinical response. The ABC transporter gene MRPA was amplified in resistant field isolates as part of an extrachromosomal circle. Co-amplification of the pterin reductase gene (PTR1) and MRPA suggests amplification of the H locus in SAG-resistant isolates. Amplification of MRPA was correlated to increased RNA as determined by real-time PCR. MRPA is an ABC-thiol transporter, and cysteine and glutathione were increased in the resistant isolates. Ornithine decarboxylase (a rate limiting enzyme in polyamine biosynthesis), and gamma-glutamylcysteine synthetase (a rate limiting enzyme in glutathione biosynthesis), the two building blocks of the main cellular thiol trypanothione, were overexpressed in some of the resistant isolates. CONCLUSIONS: A variety of resistance mechanisms to SAG, most of them consistent with a model based on the study of resistance in vitro, were present in clinical isolates from the same geographical region.

Role of the ABC transporter MRPA (PGPA) in antimony resistance in Leishmania infantum axenic and intracellular amastigotes.

Antimonial compounds are the mainstay for the treatment of infections with the protozoan parasite Leishmania. We present our studies on Leishmania infantum amastigote parasites selected for resistance to potassium antimonyl tartrate [Sb(III)]. Inside macrophages, the Sb(III)-selected cells are cross-resistant to sodium stibogluconate (Pentostam), the main drug used against Leishmania. Putative alterations in the level of expression of more than 40 genes were compared between susceptible and resistant axenic amastigotes using customized DNA microarrays. The expression of three genes coding for the ABC transporter MRPA (PGPA), S-adenosylhomocysteine hydrolase, and folylpolyglutamate synthase was found to be consistently increased. The levels of cysteine were found to be increased in the mutant. Transfection of the MRPA gene was shown to confer sodium stibogluconate resistance in intracellular parasites. This MRPA-mediated resistance could be reverted by using the glutathione biosynthesis-specific inhibitor buthionine sulfoximine. These results highlight for the first time the role of MRPA in antimony resistance in the amastigote stage of the parasite and suggest a strategy for reversing resistance.