Decreased glutamate transport in acivicin resistant Leishmania tarentolae

Gaétan Roy; Arijit Bhattacharya; Philippe Leprohon; Marc Ouellette

doi:10.1371/journal.pntd.0010046

Abstract

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.

Author summary

Studies of drug resistance in the protozoan parasites of the genus Leishmania have been helpful in revealing biochemical pathways as potential drug targets. Here we report on the characterization at the genomics and metabolomics levels of a L. tarentolae strain made resistant to acivicin, an analogue of glutamine with activity against this parasite. We found that resistance to acivicin is accompanied by a reduced uptake and intracellular levels of glutamate and that both are expected to share the same transporter. Through gene overexpression and disruption studies we identified the main amino acid transporter responsible for glutamate uptake.

Citation: Roy G, Bhattacharya A, Leprohon P, Ouellette M (2021) Decreased glutamate transport in acivicin resistant Leishmania tarentolae. PLoS Negl Trop Dis 15(12): e0010046. https://doi.org/10.1371/journal.pntd.0010046

Editor: Andrew Paul Jackson, University of Liverpool, UNITED KINGDOM

Received: July 9, 2021; Accepted: December 2, 2021; Published: December 16, 2021

Copyright: © 2021 Roy et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: The data underlying the results presented in the study are available from the Sequence Reads Archive under accession SRX10910691.

Funding: This work was supported by a Canadian Institutes of Health Research (www.cihr-irsc.ca) Foundation grant to M.O (grant number FND167283). M.O holds a Canada Research Chair in Antimicrobial Resistance. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

The protozoan parasites of the genus Leishmania are distributed worldwide and are responsible for a number of clinical manifestations. A handful number of drugs are available which are unfortunately associated with toxicity, difficulty in administration, cost or drug resistance. There is a need for novel drugs, and there has been many studies investigating different metabolic pathways in Leishmania that could eventually be exploited for drug development (reviewed in [1–4]).

While Leishmania species are purine auxotroph (reviewed in [5]), they can synthesize de novo and salvage pyrimidines [6]. Recent studies of Leishmania cells selected for resistance to the pyrimidine analogue 5-fluorouracil revealed a number of salvage enzymes and the capacity of Leishmania to transport uracil [7]. In T. brucei, pyrimidine starvation is lethal but can be rescued by salvage, suggesting that it may not be an essential function for T. brucei [8]. Recent work, however, has shown that one of the salvage enzymes, thymidine kinase, is essential for parasite viability [9,10] thus suggesting that pyrimidine metabolism may finally indeed be a valid pathway for drug development. The use of drugs has been helpful for understanding biochemical pathways and for example the studies of resistance to the model drugs methotrexate and sinefungin were helpful in our understanding of one-carbon metabolism in Leishmania [11–13].

In this study we selected L. tarentolae for resistance to acivicin. Acivicin is a chlorinated glutamine analogue, first developed as an anti-cancer drug but shown to be highly active against promastigotes and intracellular L. donovani [14]. Acivicin was shown to target carbamyl phosphate synthase (CPS), the first enzyme of the de novo pyrimidine biosynthetic pathway [14]. Protozoan parasites have a unique bifunctional enzyme with a N-terminal glutamine amidotransferase domain fused to the C-terminal CPS domain and the activity of acivicin is due to the inhibition of the glutamine amidotransferase domain [15]. Acivicin was shown to be active against other enzymes containing an amidotransferase domain such as the pyrimidine metabolic enzyme CTP synthase [16] or the purine metabolic enzyme GMP synthase (GMPS) [17]. Acivicin has been shown to be active against T. brucei by targeting its CTP synthase [16].

We describe here the characterization of L. tarentolae selected for acivicin resistance by a combination of genomics and metabolomics approaches which revealed a link between acivicin and glutamate accumulation.

Methods

Cells lines, culture conditions

L. tarentolae TarII [18] was grown as promastigotes at 25°C in SDM-79 medium supplemented with 10% heat inactivated fetal bovine serum and 10 μg/mL hemin pH 7.0. The TarII strain was selected for acivicin (Sigma-Aldrich) resistance using a stepwise selection up to 150 μM resulting in the mutant L. tarentolae A150.1. L. tarentolae A150.1 was grown without acivicin for 5, 10, and 20 passages to obtain partial revertants.

DNA manipulation

The genes of interest were amplified from genomic DNA with Phusion Taq polymerase (NEB). After sequencing, to ascertain their correct sequence, the fragments were cloned in the XbaI- HindIII sites of the Leishmania expression vectors pSPαHYGα or pSPαNEOα [19]. Gene transfection was carried out by electroporation and transfectants were selected with 1 mM of hygromycin or 0.1 mM of G418. Gene knockout work was performed essentially as described [20] by the integration of NEO and HYG inactivation cassettes facilitated by the CRISPR-Cas9 system. Southern blot of PFGE and total DNA gels were done by hybridization with [α-³²P] dCTP-labeled DNA probes according to standard protocols.

Pulse-field gel electrophoresis (PFGE)

Intact chromosomes were prepared from exponential phase culture of L. tarentolae and lysed in 1% low melting agarose plugs. The parasites were lysed in MLM buffer (0.5M EDTA pH 9.5, 1% SDS and 350 mg/mL proteinase K) overnight at 50°C. The chromosome were separated in TBE 0.5× (from TBE 10×: 1M Tris, 1M Boric Acid, 0.02 M EDTA) by Pulsed-Field Gel Electrophoresis (PFGE) using a Bio-Rad CHEF-DR III apparatus at 5 V/cm and 120° separation angle as described previously [21].

Whole genome sequencing and analysis

Paired-ends sequencing libraries were prepared from L. tarentolae genomic DNA with the Nextera DNA sample prep kit and sequenced on an Illumina HiSeq platform with 101-nucleotides paired-ends reads. Sequence reads were aligned to the L. tarentolae Parrot TarII PacBio genome assembly [22] using bwa-mem [23]. Read duplicates were marked using Picard and GATK was applied for single nucleotide variants and small insertions or deletions discovery [24]. Copy numbers variations were derived from read depth coverage as described earlier [25].

Cellular transport assays

Transport assay were performed as previously described [26] with minor modifications. Briefly, for each time point 5 × 10⁷ parasites in exponential growth-phase were washed 3 times with phosphate-buffered saline (PBS) and pelleted. The cell pellet was resuspended in 100 uL of PBS, mixed with an equal volume of different dilutions of cold glutamate or glutamine traced with Tritium-labeled L-[3,4³H]-glutamic acid (PerkinElmer) or L-[3,4³H(N)]-glutamine (Perkin Elmer) in PBS at room temperature and layered over 100 uL of dibutylphthalate (Sigma) in a 1.5 mL tube. Transport assays were stopped at various time points by spinning the cells through the inert dibutylphthalate layer for 20 seconds at maximum speed in a microcentrifuge. Unincorporated glutamate or glutamine was removed by aspiration and cells were washed twice with 1 mL of cold PBS. Competition studies were carried out using 0.1 mM of cold glutamate traced with 1.0 μCi/assay of L-[3,4-³H]- glutamate and competed with 1mM glutamate or glutamine or acivicin. Glutamate and glutamine uptake was counted in a liquid scintillation counter (Beckman) and normalized to cell number. The background transport value was removed by subtracting the accumulation value obtained on ice. Incorporated glutamate or glutamine was calculated as follows: Glu(or Gln)_i = cpm_i × [Glu (or Gln)] × ν × cpm_st^-1 × t^-1 where Glu_i (or Gln_i) is the incorporated glutamate or glutamine, [Glu] or [Gln] is the nanomolar concentration of radiolabeled glutamate or glutamine, ν is the volume of radiolabeled glutamate or glutamine (100 uL in all cases), cpm_st is the total cpm measured for the added radiolabeled glutamate or glutamine and t is the time of incubation measured in minutes.

Glutamate and glutamine measurements in Leishmania by LC-MS

The amino acids glutamate and glutamine were quantitated essentially as described [27] by liquid chromatography (LC) on an Acquity UPLC I-Class binary pump system and mass spectrometry (MS) on a Q-TOF Synapt G2-Si system with an electrospray ionization source in positive ionization mode. The relative concentrations of intracellular glutamine and glutamate in metabolite extracts were determined against calibration curves.

Statistical analysis

Growth curves were performed in biological triplicates (unless specified otherwise), each performed with technical triplicates. For statistical analysis a two-tailed unpaired t-test with GraphPad Prism software was performed unless mentioned otherwise.

Results

Leishmania tarentolae acivicin resistant mutant and gene copy number variation

L. tarentolae promastigotes were susceptible to acivicin in the sub μM range (Table 1). One L. tarentolae clone was selected step by step for acivicin resistance up to 150 μM (L. tarentolae A150.1) and was close to 100-fold more resistant than the wild-type (WT) cell (Table 1). We sequenced the genome of L. tarentolae A150.1 that was compared to the L. tarentolae WT genome sequence [28] re-sequenced by long-reads sequencing [22]. A total of 20.7M reads were obtained for the mutant, producing a genome assembly of 31.2 Mb with a coverage depth of ≥30 fold for more than 90% of the nucleotides. Copy number variations (CNVs) were identified by comparing the coverage of uniquely mapped reads between L. tarentolae A150.1 and the WT line. Several cases of supernumerary chromosomes were observed in L. tarentolae A150.1 (Fig 1). Most of these had a gain of 1.5 chromosome copies compared to WT parasites, with the exception of chromosome 20 which is tetraploid in L. tarentolae A150.1 and thus gained 2 copies compared to L. tarentolae WT. The L. tarentolae A150.1 mutant also lost one allele for a few chromosomes (Fig 1).

Download:

Fig 1. Polyploidy of chromosomes in L. tarentolae wild-type and A150.1 cells.

The number of alleles per chromosome were inferred from the sequencing data.

https://doi.org/10.1371/journal.pntd.0010046.g001

Download:

Table 1. Susceptibility to acivicin for the L. tarentolae strains used in this study.

https://doi.org/10.1371/journal.pntd.0010046.t001

Normalized read depth coverage allowed the identification of genomic loci that were amplified or deleted in L. tarentolae A150.1. These are characterized by stretches of genes whose normalized read coverage abruptly varies compared to the WT (S1 Fig and S1 Table). Among the most striking observations were a subtelomeric region of ~55 kb on chromosome 22 that was amplified ~16-fold in L. tarentolae A150.1 based on the normalized read counts (Fig 2A). Amplification of smaller regions on chromosomes 35 (~32 fold) (Fig 2B) and 15 (~4-fold) (S1 Fig) and deletions on chromosome 31 (~2-fold) (Fig 2C) and at the end of chromosome 29 (~2 fold) (S1 Fig) were also observed. The profile of the different chromosomes suggests several other genomics events but these are either of a lesser fold difference or covering a single gene (S1 Fig), the latter of which should be considered with caution.

Download:

Fig 2. Gene amplification and deletion in L. tarentolae A150.1.

Log2-transformed L. tarentolae A150.1/WT reads ratios revealed amplification of genomic segments on chromosome 22 (A) and 35 (B), as well as gene deletion on chromosome 31 (C). Southern blot hybridizations of chromosomes separated by pulse-field gel electrophoresis using 32P-labelled DNA probes derived from chromosome 22 (D) and 35 (E). Southern blot of genomic DNAs digested with NheI and hybridized to a LtaPh_3106161 probe (F). The bands at 7.3 and 7.5 Kb correspond to LtaPh_3106161 while the band at 15 Kb corresponds to at least one AAT homologue that is not deleted. The ethidium bromide gel on the left shows equal lane loading. For panels D to F, the probes used for hybridization correspond to the complete open reading frame for the genes indicated under the blots. The gene sequence of LtaPh_2200100 in panel D shares a high homology with the gene LtaPh_0613000 which likely explains the two bands observed for both A150.1 and WT. Lane 1, L. tarentolae WT; Lane 2, L. tarentolae A150.1.

https://doi.org/10.1371/journal.pntd.0010046.g002

The change in copy number of specific genomic regions was confirmed by Southern blots. When we migrated the chromosomes of L. tarentolae cells in a pulse field gel (PFGE) we could readily see a low molecular weight band, smaller than the smallest chromosome, in L. tarentolae A150.1 which was absent in the WT cells (Fig 2D). The migration profile of this extra band is a signature for linear amplification [29]. It corresponds to the CNV observed for chromosome 22 (Fig 2A), since this amplicon hybridized with a number of probes spanning chromosome 22 (Fig 2D). The hybridization data was consistent with the breakpoint observed with sequence reads along chromosome 22, the linear amplicon extending from the telomeric sequences up to LtaPh_2201300 (Fig 2D and S1 Table).

Our CNV analysis also highlighted a small region amplified on chromosome 35 and hybridization with probes derived from the region indeed confirmed the amplification of a locus containing 3 genes as part of extrachromosomal circles, as suggested by the signature hybridization smear in L. tarentolae A150.1 (Fig 2E). The deletion of a genomic segment covering 10 genes on chromosome 31 revealed by the sequencing data (Fig 2C) was also supported by a Southern blot of NheI digested DNAs hybridized to a LtaPh_3106161 probe (Fig 2F). The product of this gene is part of a family of amino acid transporters (AATs) and would fall into the AAT1 subfamily according to the phylogenetic analysis and nomenclature of Jackson [30].

Single nucleotide polymorphisms in L. tarentolae acivicin resistant A150.1 mutant

Outside CNVs, single nucleotide polymorphisms (SNPs) are also important contributors of drug resistance in Leishmania [31]. The genome of L. tarentolae A150.1 had a total of 2336 SNPs compared to L. tarentolae WT, of which 416 were in genes and non-synonymous. Of these, 51 were homozygous and 365 were heterozygous. A total of 2519 small insertions/deletions (indels) were also detected, of which 5 were in genes and homozygous and 85 were in genes and heterozygous. The list of homozygous and heterozygous SNPs and indels in coding sequences is provided in S2 Table. The list of SNPs and indels within intergenic regions can be found in S3 Table. Since acivicin is a glutamine analogue we noticed SNPs in three AATs, including the homozygous change in LtaPh_3646100 and LtaPh_3106161, and the heterozygous substitutions in LtaPh_0712500 (S2 Table), belonging respectively to the AAT14, AAT1 and AAT19 subfamilies [30]. Another mutation potentially linked to acivicin resistance (see discussion) consists in the homozygous change in the glucosamine-fructose-6-phosphate aminotransferase LtaPh_0609200 (S2 Table).

Genes involved in acivicin resistance

We first concentrated on genes with altered CNVs and that could produce resistance to acivicin. The linear episome derived from chromosome 22 contains 14 genes (LtaPh_2200051 to LtaPh_2201300) and one of them, LtaPh_2201100, encodes for the Leishmania GMPS orthologue. This is one known target of acivicin and we directly tested its role by cloning it into an expression vector that was transfected in L. tarentolae WT cells. The L. tarentolae GMPS transfectant was 1.8 fold more resistant to acivicin than the control cells (Table 1). This was a minor albeit statistically significant contribution to resistance. The identity of the other genes co-amplified with GMPS did not suggest a role for any of them in acivicin resistance. There were only three genes (LtaPh_3528000-LtaPh_3528200) on the circular amplicon derived from chromosome 35 but transfection into L. tarentolae WT cells did not confer increased resistance to acivicin (Table 1).

The sequence data revealed two deletions on chromosome 31 that covered putative AATs (Fig 2C and S1 Table). Chromosome 31 is polyploid (Fig 1) and the deletion signals are coherent with the genes being loss from at least half of chromosome 31 alleles in L. tarentolae A150.1 (Fig 2C and 2F and S1 Table). A decrease in copy number of putative amino acid transporter genes prompted us to measure the transport of glutamate and glutamine since these are close analogues to acivicin. We first confirmed that L. tarentolae accumulates both glutamate (Fig 3A) and glutamine (Fig 3B). Different gene products possibly transport these two amino acids since acivicin competed with glutamate (Fig 3A) but not glutamine (Fig 3B) accumulation. Acivicin seems to compete equally well as cold glutamate (Fig 3A). Interestingly, the accumulation of glutamate, but not glutamine, was greatly diminished in L. tarentolae A150.1 (Fig 3C and 3D). The intracellular levels of glutamate and glutamine were also calculated in L. tarentolae WT and A150.1 by quantitative metabolomics using liquid-chromatography mass spectrometry (LC-MS) [27]. The level of glutamate was decreased by 2-fold in L. tarentolae A150.1 (Fig 3E) while the level of glutamine was 5-fold higher in the mutant compared to L. tarentolae WT (Fig 3F). We also monitored the stability of the resistance phenotype and we found that part of the resistance was lost with the number of passages in the absence of acivicin (Fig 4A). Interestingly this partial loss in resistance was paralleled with an increased capacity of transporting glutamate (Fig 4B).

Download:

Fig 3. Accumulation of glutamine and glutamate in L. tarentolae WT and A150.1.

(A) Time-dependent uptake of 0.1mM glutamate traced with L-[3,43H]-glutamate in L. tarentolae WT pSPαHYGα alone (open squares) and in the presence of 1mM acivicin (filled circles) or 1mM unlabelled glutamate (filled squares). (B) Time-dependent uptake of 0.1mM glutamine traced with L-[3,43H]-glutamine in L. tarentolae WT pSPαHYGα alone (open squares) and in the presence of 1mM acivicin (filled circles) or 1mM unlabelled glutamine (filled squares). (C) Time-dependent uptake of 0.1mM glutamate traced with L-[3,43H]-glutamate in L. tarentolae WT (open squares) and A150.1 (open triangles). (D) Time-dependent uptake of 0.1mM glutamine traced with L-[3,43H(N)]-glutamine in L. tarentolae WT (open squares) and A150.1 (open triangles). Levels of intracellular glutamate (E) and glutamine (F) in L. tarentolae WT and A150.1 as measured by liquid-chromatography coupled to mass spectrometry. ** p<0.01; ***p<0.001.

https://doi.org/10.1371/journal.pntd.0010046.g003

Download:

Fig 4. Gradual increase in sensitivity to acivicin and glutamate transport in L. tarentolae A150.1 grown in absence of acivicin.

The L. tarentolae A150.1 mutant was grown for 5, 10 and 20 passages in absence of acivicin. At those passages, the level of acivicin resistance relative to A150.1 (A) and glutamate accumulation (B) were measured. For the latter, we measured the time-dependent uptake of 0.1mM glutamate traced with L-[3,43H]-glutamate.

https://doi.org/10.1371/journal.pntd.0010046.g004

Overexpression of the AAT1 LtaPh_3106161, whose gene copies are partly deleted in A150.1 (Fig 2C and 2F and S1 Table), in L. tarentolae A150.1 did not alter glutamine accumulation but exhibited a modest but significant increase in glutamate accumulation (Fig 5A). Overexpression of LtaPh_3106161 in A150.1 did revert acivicin resistance by 3-fold (Table 1), supporting the phenotypic consequences of its deletion in A150.1. This re-sensitization is specific since overexpression of the close homologue LtaPh_3103700 did not revert resistance (Table 1). Mutations (SNPs) were found in two additional putative amino acids transporters (LtaPh_0712500 and LtaPh_3646100) (S2 Table). Transfection of LtaPh_0712500 in L. tarentolae A150.1 did not impact glutamine transport but similarly to the AAT1 LtaPh_3106161 it did increase slightly glutamate accumulation (Fig 5B). It did not however impacted acivicin resistance (Table 1). Transfection of LtaPh_3646100 in A150.1 was intriguingly associated with a small decrease in glutamine accumulation but in this transfectant the accumulation of glutamate was significantly increased (Fig 5C). While the transfection of the WT versions of these three AATs (deleted or mutated in A150.1) in L. tarentolae A150.1 reverted partly the defect in glutamate transport in A150.1, none of these AATs appear to qualify as the main glutamate transporter.

Download:

Fig 5. Accumulation of glutamate and glutamine in L. tarentolae A150.1 and in its transfectants with amino acid transporter (AAT) genes.

Time-dependent uptake of 0.1mM glutamate (left panels) or glutamine (right panels) traced with respectively L-[3,43H]-glutamate or L-[3,43H(N)]-glutamine in the L. tarentolae A150.1 mutant transfected with either a control vector (empty circles) or with LtaPh_3106161 (A), LtaPh_07.12500 (B), LtaPh_36.46100 (C), LmjF22.0230 (D) (filled circles). For the latter (D) we also show the accumulation of glutamate in wild-type cells transfected with LmjF22.0230 (filled squares) or with the control vector (empty squares). * p<0.05; ** p<0.01; ***p<0.001. For glutamate in panel D the level of significance indicated is for the comparison between mock-transfected and LmjF.22.0230-transfected cells (empty vs filled symbols) of the same strain (circles or squares). We noticed that the uptake of glutamate and glutamine in panel D were slightly lower in the wild-type control in comparison to what shown in Fig 3A and 3B, but overall glutamate uptake is active in the wild-type and not in the mutant A150.1.

https://doi.org/10.1371/journal.pntd.0010046.g005

Members of the solute carrier (SLC) proteins are structurally and functionally diverse and some subclasses were shown to transport amino acids using H(+) as the major driving force. Active transport of glutamate in Leishmania appears to be sensitive to extracellular concentration of H(+) but not to Na(+) [32]. When we aligned sequences of Leishmania putative amino acid permeases of unknown function we found that LtaPh_2202300 (LmjF22.0230) (AAT22 according to [30]) had 27% identity to SLC38A6 (SNAT6) (S2 Fig), a transporter with properties that could qualify it as a transporter of glutamate and which is conserved among Leishmania species (S3 Fig). We studied the L. major orthologue, as there were gaps with the L. tarentolae sequence genome draft when this work was started. We thus cloned the L. major orthologue LmjF22.0230 AAT22 in an expression vector that was transfected into L. tarentolae WT and A150.1. While episomal expression of LmjF22.0230 did not change the transport properties to glutamine in L. tarentolae WT (Fig 5D), the transport of glutamate was higher (Fig 5D). Transfection of LmjF22.0230 also restored glutamate transport in L. tarentolae A150.1 close to wild-type levels (Fig 5D). This regained capacity of transporting glutamate was associated with partial reversion of acivicin resistance (Table 1). In light of this result we reanalysed the sequence of AAT22 in A150.1 but could not detect either SNPs or CNVs, confirming our initial analysis that this gene is unaltered. To further the role of AAT22 in glutamate transport we inactivated LtaPh_2202300 using the CRISPR-Cas9 genome editing technology. Two inactivation cassettes with the NEO and PUR selectable markers were designed (Fig 6A) and transfected along with an episome driving the expression of a guide RNA targeting LtaPh_2202300 in a L. tarentolae WT cell line expressing the Cas9 nuclease [20]. A null mutant of LtaPh_2202300 was obtained as illustrated by a Southern blot analysis and the absence of an AAT22 WT allele (Fig 6B). The accumulation of glutamate but not glutamine was greatly diminished in that recombinant strain (Fig 6C). This strain was 1.6-fold more resistant to acivicin (Table 1).

Download:

Fig 6. Inactivation of the glutamate transporter LtaPh_2202300.

(A) Genomic organization of the LtaPh_2202300 locus. LtaPh_2202300 upstream and downstream fragments used for the integration of the NEO and PUR cassettes are shown. The integration of these cassettes was mediated by a CRIPR-Cas9 strategy [20]. (B) Southern blots of L. tarentolae genomic DNA digested with ApaI and XbaI and hybridized to a LtaPh_2202300 5’UTR probe. Lane 1, L. tarentolae wild type; lane 2, L. tarentolae LtaPh_2202300 null mutant with integrated NEO and PUR cassettes. (C) Time dependence uptake of 0.1mM glutamate traced with L-[3,43H]-glutamate (left) and of 0.1mM glutamine traced with L-[3,43H(N)]-glutamine (right) in the LtaPh_2202300 null mutant (filled circles) or the control strain (empty circles). ***p<0.001.

https://doi.org/10.1371/journal.pntd.0010046.g006

Discussion

Our characterization of a L. tarentolae mutant selected for acivicin resistance indicates that resistance is complex with many factors contributing modestly to the resistance phenotype. The most prominent defect was a decreased accumulation of glutamate in this mutant. Indeed, acivicin appears to be transported by the glutamate transporter and not by the glutamine transporter. Very little is known about glutamate transport in Leishmania and kinetoplastid parasites. It was characterized biochemically in L. amazonensis [32] and in T. cruzi [33] but the identity of the transporters remains unknown. One well studied AAT is the L. donovani AAP24 protein (AAT20 according to [30]) responsible for proline uptake [34] and while it regulates homeostasis of glutamate, the latter is not a substrate. Next generation sequencing is now a useful technique for analysing Leishmania drug resistant mutants but in the case of L. tar A150.1 genomics could not explain the decreased glutamate, and hence most likely acivicin, accumulation. While the genome sequence of the A150.1 mutant revealed many candidates that could explain this uptake defect (e.g. three AAT genes were either deleted or mutated) none appears to play a major role in glutamate accumulation. To gain insight into the main transporter of glutamate (and acivicin) we then resorted to phylogenetic analyses of AAT orthologs. This led to the gene LtaPh_2202300 (AAT22), which is neither mutated nor changed in copy number in A150.1, as the main glutamate transporter. Indeed, transfecting AAT22 in A150.1 restored glutamate transport while glutamate transport was diminished significantly in a L. tarentolae LtaPh_2202300 knockout cell line. Yet, we still observed glutamate transport in the AAT22 null mutant and thus there are other routes by which glutamate can enter into the cell. The 3 AAT genes found mutated or deleted in A150.1 were associated with low glutamate transport (Fig 4A–4C) and this likely explains why the transport of glutamate is lower in the mutant A150.1 in comparison to the AAT22 knockout.

How can we explain that AAT22 is not mutated in the acivicin resistant mutant A150.1? Acivicin is known to bind covalently to many of its enzymatic targets [35–37] and it was even advocated as a model molecule for the design of covalent enzymatic inhibitors [38]. Thus acivicin may covalently binds the transporter hence reducing its activity. Our observation that glutamate uptake and intracellular levels are diminished in cells resistant to acivicin (Fig 3) is consistent with this hypothesis. Cells can live without the main glutamate transporter, the covalent binding of acivicin to the transporter is likely to be sufficient to reduce its activity, and therefore there is no need for mutation or deletion of this gene. It is possible that acivicin does not bind as strongly to other AAT transporters and thus they need to be mutated or deleted to reduce the uptake of acivicin and by ricochet of glutamate. The resistance reversion work (Fig 4) is providing further support for the covalent binding hypothesis of a drug to a transport protein as a novel mechanism of drug resistance in Leishmania. Indeed, with increased passages in the absence of acivicin we observed a decreased resistance to acivicin which was correlated with increased accumulation of glutamate (and most likely acivicin) (Fig 4). Since the main glutamate transporter is neither deleted nor mutated in this mutant it suggests that with gradual passages more of the AAT transporter is acivicin free and its function is gradually restored. Transport is not to wild-type levels because other transporters are altered. Modeling and docking experiments with acivicin and glutamate could be done to support our hypothesis but the lack of a protein structure sufficiently similar to AAT22 is preventing for now such analysis.

In addition to the defect in glutamate transport we found that one potential target, GMPS, is amplified in A150.1. Its transfection only show low level resistance, however. While another locus encoded on chromosome 35 was also amplified as part of an extrachromosomal circle, we could not find a role in acivicin resistance for any of the three genes tested. Usually gene amplification is linked to a phenotype in drug resistant Leishmania cells [39] and one gene part of the chromosome 35 locus may need another partner for producing resistance. One other mutation potentially linked to acivicin resistance include the homozygous change in the glucosamine-fructose-6-phosphate aminotransferase LtaPh_0609200 (S2 Table) which if non-functional would divert glutamine away from the synthesis of glutamate and D-glucosamine-6-P, leading to the accumulation of glutamine and to a shortage of glutamate. Our measure of glutamate and glutamine levels by LC-MS (Fig 3E and 3F) are consistent with such hypothesis. We were not able to reconstruct the level of resistance found in the mutant A150.1. Possibly a number of mutations are acting synergistically. This would need to be tested, and the sequencing of the genome of cells at intermediate steps of acivicin selection could help in finding the right order of combination of genes to be transfected to reconstruct resistance.

Our study of acivicin resistance in L. tarentolae has highlighted the complexity of mutations involved in resistance to a single drug and has allowed the discovery that acivicin uptake is linked to glutamate uptake which led to the discovery of the Leishmania glutamate transporter. We propose that acivicin binds covalently to the glutamate transporter and this binding contribute to the reduction in glutamate (and acivicin) uptake and to decreased susceptibility to acivicin.

Supporting information

S1 Fig. Chromosomal amplification and deletion in L. tarentolae A150.1.

Log₂-transformed L. tarentolae A150.1/WT gene reads ratios normalized for gene length and total read counts for the 36 chromosomes. Genes along the chromosome are shown as dots on the X axis while the log₂-transformed L. tarentolae A150.1/WT reads ratios are indicated by the y axis.

https://doi.org/10.1371/journal.pntd.0010046.s001

(PDF)

S2 Fig. Alignment of SNAT6 homologs.

Alignment of the Homo sapiens solute carrier family 38 member 6 (SLC38A6), Trypanosoma cruzi strain CL Brener gene TcCLB.507811.100, Rattus norvegicus solute carrier family 38 member 6 (Slc38a6) and L. major strain Friedlin gene LmjF22.0230. Blue represents residues A, C, F, I, M, V and W; Dark Blue represents residues H and Y; Pink represents residues E and D; Green represents residues N, Q, S and T; Yellow represents residue P; Orange represents residue G; Coral represents residues K and R; ‘*’ indicates position which have a single fully conserved residue; ‘:’indicates a strong group of conserved amino acids; ‘.’ indicates a weaker group of conserved amino acids.

https://doi.org/10.1371/journal.pntd.0010046.s002

(TIF)

S3 Fig. Phylogeny of AAT22 in Leishmania ssp.

The evolutionary history was inferred using the Neighbor-Joining method. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (500 replicates) are shown next to the branches. The tree is drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using the Poisson correction method and are in the units of the number of amino acid substitutions per site. All ambiguous positions were removed for each sequence pair (pairwise deletion option). Orthologs of AAT22 in Leptomonas seymouri, Crithidia fasciculata and Trypanosoma cruzi were used as outgroup.

https://doi.org/10.1371/journal.pntd.0010046.s003

(TIF)

S1 Table.

Genes with an increased (spreadsheet A) or decreased (spreadsheet B) copy number in L. tarentolae A150.1 as measured by read counts ratios in comparison to L. tarentolae WT.

https://doi.org/10.1371/journal.pntd.0010046.s004

(XLSX)

S2 Table. List of SNVs detected in the coding regions of L. tarentolae A150.1.

https://doi.org/10.1371/journal.pntd.0010046.s005

(XLSX)

S3 Table. List of SNVs detected in the intergenic regions of L. tarentolae A150.1.

https://doi.org/10.1371/journal.pntd.0010046.s006

(XLSX)

Acknowledgments

This research was enabled in part by computing infrastructure provided by Calcul Québec (www.calculquebec.ca) and Compute Canada (www.computecanada.ca). Infrastructure and equipment were provided by the Canadian Foundation for Innovation. We thank Fereshteh Fani for help with EC₅₀s measurements.

References

1. Arya R, Dhembla C, Makde RD, Sundd M, Kundu S. An overview of the fatty acid biosynthesis in the protozoan parasite Leishmania and its relevance as a drug target against leishmaniasis. Mol Biochem Parasitol. 2021;246:111416. pmid:34555376
- View Article
- PubMed/NCBI
- Google Scholar
2. Jain V, Jain K. Molecular targets and pathways for the treatment of visceral leishmaniasis. Drug Discov Today. 2018;23(1):161–70. pmid:28919438
- View Article
- PubMed/NCBI
- Google Scholar
3. Kourbeli V, Chontzopoulou E, Moschovou K, Pavlos D, Mavromoustakos T, Papanastasiou IP. An Overview on Target-Based Drug Design against Kinetoplastid Protozoan Infections: Human African Trypanosomiasis, Chagas Disease and Leishmaniases. Molecules. 2021;26(15). pmid:34361781
- View Article
- PubMed/NCBI
- Google Scholar
4. Parab AR, McCall LI. Tryp-ing Up Metabolism: Role of Metabolic Adaptations in Kinetoplastid Disease Pathogenesis. Infect Immun. 2021;89(4). pmid:33526564
- View Article
- PubMed/NCBI
- Google Scholar
5. Boitz JM, Ullman B, Jardim A, Carter NS. Purine salvage in Leishmania: complex or simple by design? Trends Parasitol. 2012;28(8):345–52. pmid:22726696
- View Article
- PubMed/NCBI
- Google Scholar
6. Wilson ZN, Gilroy CA, Boitz JM, Ullman B, Yates PA. Genetic dissection of pyrimidine biosynthesis and salvage in Leishmania donovani. J Biol Chem. 2012;287(16):12759–70. pmid:22367196
- View Article
- PubMed/NCBI
- Google Scholar
7. Ritt JF, Raymond F, Leprohon P, Legare D, Corbeil J, Ouellette M. Gene amplification and point mutations in pyrimidine metabolic genes in 5-fluorouracil resistant Leishmania infantum. PLoS Negl Trop Dis. 2013;7(11):e2564. pmid:24278495
- View Article
- PubMed/NCBI
- Google Scholar
8. Ali JA, Tagoe DN, Munday JC, Donachie A, Morrison LJ, de Koning HP. Pyrimidine biosynthesis is not an essential function for Trypanosoma brucei bloodstream forms. PLoS One. 2013;8(3):e58034. pmid:23505454
- View Article
- PubMed/NCBI
- Google Scholar
9. Valente M, Timm J, Castillo-Acosta VM, Ruiz-Perez LM, Balzarini T, Nettleship JE, et al. Cell cycle regulation and novel structural features of thymidine kinase, an essential enzyme in Trypanosoma brucei. Mol Microbiol. 2016;102(3):365–85. pmid:27426054
- View Article
- PubMed/NCBI
- Google Scholar
10. Leija C, Rijo-Ferreira F, Kinch LN, Grishin NV, Nischan N, Kohler JJ, et al. Pyrimidine Salvage Enzymes Are Essential for De Novo Biosynthesis of Deoxypyrimidine Nucleotides in Trypanosoma brucei. PLoS Pathog. 2016;12(11):e1006010. pmid:27820863
- View Article
- PubMed/NCBI
- Google Scholar
11. Nare B, Luba J, Hardy LW, Beverley S. New approaches to Leishmania chemotherapy: pteridine reductase 1 (PTR1) as a target and modulator of antifolate sensitivity. Parasitology. 1997;114 Suppl:S101–10. pmid:9309772
- View Article
- PubMed/NCBI
- Google Scholar
12. Ouellette M, Drummelsmith J, El-Fadili A, Kundig C, Richard D, Roy G. Pterin transport and metabolism in Leishmania and related trypanosomatid parasites. Int J Parasitol. 2002;32(4):385–98. pmid:11849635
- View Article
- PubMed/NCBI
- Google Scholar
13. Bhattacharya A, Sharma M, Pakkinathan C, Rosen BP, Leprohon P, Ouellette M. Genomewide Analysis of Mode of Action of the S-Adenosylmethionine Analogue Sinefungin in Leishmania infantum. mSystems. 2019;4(5). pmid:31615876
- View Article
- PubMed/NCBI
- Google Scholar
14. Mukherjee T, Roy K, Bhaduri A. Acivicin: a highly active potential chemotherapeutic agent against visceral leishmaniasis. Biochem Biophys Res Commun. 1990;170(2):426–32. pmid:2383250
- View Article
- PubMed/NCBI
- Google Scholar
15. Fox BA, Bzik DJ. Organisation and sequence determination of glutamine-dependent carbamoyl phosphate synthetase II in Toxoplasma gondii. Int J Parasitol. 2003;33(1):89–96. pmid:12547350
- View Article
- PubMed/NCBI
- Google Scholar
16. Fijolek A, Hofer A, Thelander L. Expression, purification, characterization, and in vivo targeting of trypanosome CTP synthetase for treatment of African sleeping sickness. J Biol Chem. 2007;282(16):11858–65. pmid:17331943
- View Article
- PubMed/NCBI
- Google Scholar
17. Nakamura J, Straub K, Wu J, Lou L. The glutamine hydrolysis function of human GMP synthetase. Identification of an essential active site cysteine. J Biol Chem. 1995;270(40):23450–5. pmid:7559506
- View Article
- PubMed/NCBI
- Google Scholar
18. White TC, Fase-Fowler F, van Luenen H, Calafat J, Borst P. The H circles of Leishmania tarentolae are a unique amplifiable system of oligomeric DNAs associated with drug resistance. J Biol Chem. 1988;263(32):16977–83. pmid:3182827
- View Article
- PubMed/NCBI
- Google Scholar
19. Papadopoulou B, Roy G, Ouellette M. Autonomous replication of bacterial DNA plasmid oligomers in Leishmania. Mol Biochem Parasitol. 1994;65(1):39–49. pmid:7935627
- View Article
- PubMed/NCBI
- Google Scholar
20. Fernandez-Prada C, Sharma M, Plourde M, Bresson E, Roy G, Leprohon P, et al. High-throughput Cos-Seq screen with intracellular Leishmania infantum for the discovery of novel drug-resistance mechanisms. Int J Parasitol Drugs Drug Resist. 2018;8(2):165–73. pmid:29602064
- View Article
- PubMed/NCBI
- Google Scholar
21. Grondin K, Roy G, Ouellette M. Formation of extrachromosomal circular amplicons with direct or inverted duplications in drug-resistant Leishmania tarentolae. Mol Cell Biol. 1996;16(7):3587–95. pmid:8668175
- View Article
- PubMed/NCBI
- Google Scholar
22. Goto Y, Kuroki A, Suzuki K, Yamagishi J. Draft Genome Sequence of Leishmania tarentolae Parrot Tar II, Obtained by Single-Molecule Real-Time Sequencing. Microbiol Resour Announc. 2020;9(21). pmid:32439660
- View Article
- PubMed/NCBI
- Google Scholar
23. Li H, Durbin R. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics. 2010;26(5):589–95. pmid:20080505
- View Article
- PubMed/NCBI
- Google Scholar
24. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303. pmid:20644199
- View Article
- PubMed/NCBI
- Google Scholar
25. Fernandez-Prada C, Vincent IM, Brotherton MC, Roberts M, Roy G, Rivas L, et al. Different Mutations in a P-type ATPase Transporter in Leishmania Parasites are Associated with Cross-resistance to Two Leading Drugs by Distinct Mechanisms. PLoS Negl Trop Dis. 2016;10(12):e0005171. pmid:27911896
- View Article
- PubMed/NCBI
- Google Scholar
26. Papadopoulou B, Roy G, Ouellette M. Frequent amplification of a short chain dehydrogenase gene as part of circular and linear amplicons in methotrexate resistant Leishmania. Nucleic Acids Res. 1993;21(18):4305–12. pmid:8414986
- View Article
- PubMed/NCBI
- Google Scholar
27. El Khoury JY, Boucher N, Bergeron MG, Leprohon P, Ouellette M. Penicillin induces alterations in glutamine metabolism in Streptococcus pneumoniae. Sci Rep. 2017;7(1):14587. pmid:29109543
- View Article
- PubMed/NCBI
- Google Scholar
28. Raymond F, Boisvert S, Roy G, Ritt JF, Legare D, Isnard A, et al. Genome sequencing of the lizard parasite Leishmania tarentolae reveals loss of genes associated to the intracellular stage of human pathogenic species. Nucleic Acids Res. 2012;40(3):1131–47. pmid:21998295
- View Article
- PubMed/NCBI
- Google Scholar
29. Laffitte MC, Genois MM, Mukherjee A, Legare D, Masson JY, Ouellette M. Formation of linear amplicons with inverted duplications in Leishmania requires the MRE11 nuclease. PLoS Genet. 2014;10(12):e1004805. pmid:25474106
- View Article
- PubMed/NCBI
- Google Scholar
30. Jackson AP. Origins of amino acid transporter loci in trypanosomatid parasites. BMC Evol Biol. 2007;7:26. pmid:17319943
- View Article
- PubMed/NCBI
- Google Scholar
31. Bhattacharya A, Leprohon P, Bigot S, Padmanabhan PK, Mukherjee A, Roy G, et al. Coupling chemical mutagenesis to next generation sequencing for the identification of drug resistance mutations in Leishmania. Nat Commun. 2019;10(1):5627. pmid:31819054
- View Article
- PubMed/NCBI
- Google Scholar
32. Paes LS, Galvez Rojas RL, Daliry A, Floeter-Winter LM, Ramirez MI, Silber AM. Active transport of glutamate in Leishmania (Leishmania) amazonensis. J Eukaryot Microbiol. 2008;55(5):382–7. pmid:19017058
- View Article
- PubMed/NCBI
- Google Scholar
33. Silber AM, Rojas RL, Urias U, Colli W, Alves MJ. Biochemical characterization of the glutamate transport in Trypanosoma cruzi. Int J Parasitol. 2006;36(2):157–63. pmid:16373069
- View Article
- PubMed/NCBI
- Google Scholar
34. Inbar E, Schlisselberg D, Suter Grotemeyer M, Rentsch D, Zilberstein D. A versatile proline/alanine transporter in the unicellular pathogen Leishmania donovani regulates amino acid homoeostasis and osmotic stress responses. Biochem J. 2013;449(2):555–66. pmid:22994895
- View Article
- PubMed/NCBI
- Google Scholar
35. Williams K, Cullati S, Sand A, Biterova EI, Barycki JJ. Crystal structure of acivicin-inhibited gamma-glutamyltranspeptidase reveals critical roles for its C-terminus in autoprocessing and catalysis. Biochemistry. 2009;48(11):2459–67. pmid:19256527
- View Article
- PubMed/NCBI
- Google Scholar
36. Chittur SV, Klem TJ, Shafer CM, Davisson VJ. Mechanism for acivicin inactivation of triad glutamine amidotransferases. Biochemistry. 2001;40(4):876–87. pmid:11170408
- View Article
- PubMed/NCBI
- Google Scholar
37. Miles BW, Thoden JB, Holden HM, Raushel FM. Inactivation of the amidotransferase activity of carbamoyl phosphate synthetase by the antibiotic acivicin. J Biol Chem. 2002;277(6):4368–73. pmid:11729189
- View Article
- PubMed/NCBI
- Google Scholar
38. Pinto A, Tamborini L, Cullia G, Conti P, De Micheli C. Inspired by Nature: The 3-Halo-4,5-dihydroisoxazole Moiety as a Novel Molecular Warhead for the Design of Covalent Inhibitors. ChemMedChem. 2016;11(1):10–4. pmid:26607551
- View Article
- PubMed/NCBI
- Google Scholar
39. Laffitte MN, Leprohon P, Papadopoulou B, Ouellette M. Plasticity of the Leishmania genome leading to gene copy number variations and drug resistance. F1000Res. 2016;5:2350. pmid:27703673
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. Arya R, Dhembla C, Makde RD, Sundd M, Kundu S. An overview of the fatty acid biosynthesis in the protozoan parasite Leishmania and its relevance as a drug target against leishmaniasis. Mol Biochem Parasitol. 2021;246:111416. pmid:34555376
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Jain V, Jain K. Molecular targets and pathways for the treatment of visceral leishmaniasis. Drug Discov Today. 2018;23(1):161–70. pmid:28919438
View Article
PubMed/NCBI
Google Scholar

[6] View Article

[7] PubMed/NCBI

[8] Google Scholar

[ref3] 3. Kourbeli V, Chontzopoulou E, Moschovou K, Pavlos D, Mavromoustakos T, Papanastasiou IP. An Overview on Target-Based Drug Design against Kinetoplastid Protozoan Infections: Human African Trypanosomiasis, Chagas Disease and Leishmaniases. Molecules. 2021;26(15). pmid:34361781
View Article
PubMed/NCBI
Google Scholar

[10] View Article

[11] PubMed/NCBI

[12] Google Scholar

[ref4] 4. Parab AR, McCall LI. Tryp-ing Up Metabolism: Role of Metabolic Adaptations in Kinetoplastid Disease Pathogenesis. Infect Immun. 2021;89(4). pmid:33526564
View Article
PubMed/NCBI
Google Scholar

[14] View Article

[15] PubMed/NCBI

[16] Google Scholar

[ref5] 5. Boitz JM, Ullman B, Jardim A, Carter NS. Purine salvage in Leishmania: complex or simple by design? Trends Parasitol. 2012;28(8):345–52. pmid:22726696
View Article
PubMed/NCBI
Google Scholar

[18] View Article

[19] PubMed/NCBI

[20] Google Scholar

[ref6] 6. Wilson ZN, Gilroy CA, Boitz JM, Ullman B, Yates PA. Genetic dissection of pyrimidine biosynthesis and salvage in Leishmania donovani. J Biol Chem. 2012;287(16):12759–70. pmid:22367196
View Article
PubMed/NCBI
Google Scholar

[22] View Article

[23] PubMed/NCBI

[24] Google Scholar

[ref7] 7. Ritt JF, Raymond F, Leprohon P, Legare D, Corbeil J, Ouellette M. Gene amplification and point mutations in pyrimidine metabolic genes in 5-fluorouracil resistant Leishmania infantum. PLoS Negl Trop Dis. 2013;7(11):e2564. pmid:24278495
View Article
PubMed/NCBI
Google Scholar

[26] View Article

[27] PubMed/NCBI

[28] Google Scholar

[ref8] 8. Ali JA, Tagoe DN, Munday JC, Donachie A, Morrison LJ, de Koning HP. Pyrimidine biosynthesis is not an essential function for Trypanosoma brucei bloodstream forms. PLoS One. 2013;8(3):e58034. pmid:23505454
View Article
PubMed/NCBI
Google Scholar

[30] View Article

[31] PubMed/NCBI

[32] Google Scholar

[ref9] 9. Valente M, Timm J, Castillo-Acosta VM, Ruiz-Perez LM, Balzarini T, Nettleship JE, et al. Cell cycle regulation and novel structural features of thymidine kinase, an essential enzyme in Trypanosoma brucei. Mol Microbiol. 2016;102(3):365–85. pmid:27426054
View Article
PubMed/NCBI
Google Scholar

[34] View Article

[35] PubMed/NCBI

[36] Google Scholar

[ref10] 10. Leija C, Rijo-Ferreira F, Kinch LN, Grishin NV, Nischan N, Kohler JJ, et al. Pyrimidine Salvage Enzymes Are Essential for De Novo Biosynthesis of Deoxypyrimidine Nucleotides in Trypanosoma brucei. PLoS Pathog. 2016;12(11):e1006010. pmid:27820863
View Article
PubMed/NCBI
Google Scholar

[38] View Article

[39] PubMed/NCBI

[40] Google Scholar

[ref11] 11. Nare B, Luba J, Hardy LW, Beverley S. New approaches to Leishmania chemotherapy: pteridine reductase 1 (PTR1) as a target and modulator of antifolate sensitivity. Parasitology. 1997;114 Suppl:S101–10. pmid:9309772
View Article
PubMed/NCBI
Google Scholar

[42] View Article

[43] PubMed/NCBI

[44] Google Scholar

[ref12] 12. Ouellette M, Drummelsmith J, El-Fadili A, Kundig C, Richard D, Roy G. Pterin transport and metabolism in Leishmania and related trypanosomatid parasites. Int J Parasitol. 2002;32(4):385–98. pmid:11849635
View Article
PubMed/NCBI
Google Scholar

[46] View Article

[47] PubMed/NCBI

[48] Google Scholar

[ref13] 13. Bhattacharya A, Sharma M, Pakkinathan C, Rosen BP, Leprohon P, Ouellette M. Genomewide Analysis of Mode of Action of the S-Adenosylmethionine Analogue Sinefungin in Leishmania infantum. mSystems. 2019;4(5). pmid:31615876
View Article
PubMed/NCBI
Google Scholar

[50] View Article

[51] PubMed/NCBI

[52] Google Scholar

[ref14] 14. Mukherjee T, Roy K, Bhaduri A. Acivicin: a highly active potential chemotherapeutic agent against visceral leishmaniasis. Biochem Biophys Res Commun. 1990;170(2):426–32. pmid:2383250
View Article
PubMed/NCBI
Google Scholar

[54] View Article

[55] PubMed/NCBI

[56] Google Scholar

[ref15] 15. Fox BA, Bzik DJ. Organisation and sequence determination of glutamine-dependent carbamoyl phosphate synthetase II in Toxoplasma gondii. Int J Parasitol. 2003;33(1):89–96. pmid:12547350
View Article
PubMed/NCBI
Google Scholar

[58] View Article

[59] PubMed/NCBI

[60] Google Scholar

[ref16] 16. Fijolek A, Hofer A, Thelander L. Expression, purification, characterization, and in vivo targeting of trypanosome CTP synthetase for treatment of African sleeping sickness. J Biol Chem. 2007;282(16):11858–65. pmid:17331943
View Article
PubMed/NCBI
Google Scholar

[62] View Article

[63] PubMed/NCBI

[64] Google Scholar

[ref17] 17. Nakamura J, Straub K, Wu J, Lou L. The glutamine hydrolysis function of human GMP synthetase. Identification of an essential active site cysteine. J Biol Chem. 1995;270(40):23450–5. pmid:7559506
View Article
PubMed/NCBI
Google Scholar

[66] View Article

[67] PubMed/NCBI

[68] Google Scholar

[ref18] 18. White TC, Fase-Fowler F, van Luenen H, Calafat J, Borst P. The H circles of Leishmania tarentolae are a unique amplifiable system of oligomeric DNAs associated with drug resistance. J Biol Chem. 1988;263(32):16977–83. pmid:3182827
View Article
PubMed/NCBI
Google Scholar

[70] View Article

[71] PubMed/NCBI

[72] Google Scholar

[ref19] 19. Papadopoulou B, Roy G, Ouellette M. Autonomous replication of bacterial DNA plasmid oligomers in Leishmania. Mol Biochem Parasitol. 1994;65(1):39–49. pmid:7935627
View Article
PubMed/NCBI
Google Scholar

[74] View Article

[75] PubMed/NCBI

[76] Google Scholar

[ref20] 20. Fernandez-Prada C, Sharma M, Plourde M, Bresson E, Roy G, Leprohon P, et al. High-throughput Cos-Seq screen with intracellular Leishmania infantum for the discovery of novel drug-resistance mechanisms. Int J Parasitol Drugs Drug Resist. 2018;8(2):165–73. pmid:29602064
View Article
PubMed/NCBI
Google Scholar

[78] View Article

[79] PubMed/NCBI

[80] Google Scholar

[ref21] 21. Grondin K, Roy G, Ouellette M. Formation of extrachromosomal circular amplicons with direct or inverted duplications in drug-resistant Leishmania tarentolae. Mol Cell Biol. 1996;16(7):3587–95. pmid:8668175
View Article
PubMed/NCBI
Google Scholar

[82] View Article

[83] PubMed/NCBI

[84] Google Scholar

[ref22] 22. Goto Y, Kuroki A, Suzuki K, Yamagishi J. Draft Genome Sequence of Leishmania tarentolae Parrot Tar II, Obtained by Single-Molecule Real-Time Sequencing. Microbiol Resour Announc. 2020;9(21). pmid:32439660
View Article
PubMed/NCBI
Google Scholar

[86] View Article

[87] PubMed/NCBI

[88] Google Scholar

[ref23] 23. Li H, Durbin R. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics. 2010;26(5):589–95. pmid:20080505
View Article
PubMed/NCBI
Google Scholar

[90] View Article

[91] PubMed/NCBI

[92] Google Scholar

[ref24] 24. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303. pmid:20644199
View Article
PubMed/NCBI
Google Scholar

[94] View Article

[95] PubMed/NCBI

[96] Google Scholar

[ref25] 25. Fernandez-Prada C, Vincent IM, Brotherton MC, Roberts M, Roy G, Rivas L, et al. Different Mutations in a P-type ATPase Transporter in Leishmania Parasites are Associated with Cross-resistance to Two Leading Drugs by Distinct Mechanisms. PLoS Negl Trop Dis. 2016;10(12):e0005171. pmid:27911896
View Article
PubMed/NCBI
Google Scholar

[98] View Article

[99] PubMed/NCBI

[100] Google Scholar

[ref26] 26. Papadopoulou B, Roy G, Ouellette M. Frequent amplification of a short chain dehydrogenase gene as part of circular and linear amplicons in methotrexate resistant Leishmania. Nucleic Acids Res. 1993;21(18):4305–12. pmid:8414986
View Article
PubMed/NCBI
Google Scholar

[102] View Article

[103] PubMed/NCBI

[104] Google Scholar

[ref27] 27. El Khoury JY, Boucher N, Bergeron MG, Leprohon P, Ouellette M. Penicillin induces alterations in glutamine metabolism in Streptococcus pneumoniae. Sci Rep. 2017;7(1):14587. pmid:29109543
View Article
PubMed/NCBI
Google Scholar

[106] View Article

[107] PubMed/NCBI

[108] Google Scholar

[ref28] 28. Raymond F, Boisvert S, Roy G, Ritt JF, Legare D, Isnard A, et al. Genome sequencing of the lizard parasite Leishmania tarentolae reveals loss of genes associated to the intracellular stage of human pathogenic species. Nucleic Acids Res. 2012;40(3):1131–47. pmid:21998295
View Article
PubMed/NCBI
Google Scholar

[110] View Article

[111] PubMed/NCBI

[112] Google Scholar

[ref29] 29. Laffitte MC, Genois MM, Mukherjee A, Legare D, Masson JY, Ouellette M. Formation of linear amplicons with inverted duplications in Leishmania requires the MRE11 nuclease. PLoS Genet. 2014;10(12):e1004805. pmid:25474106
View Article
PubMed/NCBI
Google Scholar

[114] View Article

[115] PubMed/NCBI

[116] Google Scholar

[ref30] 30. Jackson AP. Origins of amino acid transporter loci in trypanosomatid parasites. BMC Evol Biol. 2007;7:26. pmid:17319943
View Article
PubMed/NCBI
Google Scholar

[118] View Article

[119] PubMed/NCBI

[120] Google Scholar

[ref31] 31. Bhattacharya A, Leprohon P, Bigot S, Padmanabhan PK, Mukherjee A, Roy G, et al. Coupling chemical mutagenesis to next generation sequencing for the identification of drug resistance mutations in Leishmania. Nat Commun. 2019;10(1):5627. pmid:31819054
View Article
PubMed/NCBI
Google Scholar

[122] View Article

[123] PubMed/NCBI

[124] Google Scholar

[ref32] 32. Paes LS, Galvez Rojas RL, Daliry A, Floeter-Winter LM, Ramirez MI, Silber AM. Active transport of glutamate in Leishmania (Leishmania) amazonensis. J Eukaryot Microbiol. 2008;55(5):382–7. pmid:19017058
View Article
PubMed/NCBI
Google Scholar

[126] View Article

[127] PubMed/NCBI

[128] Google Scholar

[ref33] 33. Silber AM, Rojas RL, Urias U, Colli W, Alves MJ. Biochemical characterization of the glutamate transport in Trypanosoma cruzi. Int J Parasitol. 2006;36(2):157–63. pmid:16373069
View Article
PubMed/NCBI
Google Scholar

[130] View Article

[131] PubMed/NCBI

[132] Google Scholar

[ref34] 34. Inbar E, Schlisselberg D, Suter Grotemeyer M, Rentsch D, Zilberstein D. A versatile proline/alanine transporter in the unicellular pathogen Leishmania donovani regulates amino acid homoeostasis and osmotic stress responses. Biochem J. 2013;449(2):555–66. pmid:22994895
View Article
PubMed/NCBI
Google Scholar

[134] View Article

[135] PubMed/NCBI

[136] Google Scholar

[ref35] 35. Williams K, Cullati S, Sand A, Biterova EI, Barycki JJ. Crystal structure of acivicin-inhibited gamma-glutamyltranspeptidase reveals critical roles for its C-terminus in autoprocessing and catalysis. Biochemistry. 2009;48(11):2459–67. pmid:19256527
View Article
PubMed/NCBI
Google Scholar

[138] View Article

[139] PubMed/NCBI

[140] Google Scholar

[ref36] 36. Chittur SV, Klem TJ, Shafer CM, Davisson VJ. Mechanism for acivicin inactivation of triad glutamine amidotransferases. Biochemistry. 2001;40(4):876–87. pmid:11170408
View Article
PubMed/NCBI
Google Scholar

[142] View Article

[143] PubMed/NCBI

[144] Google Scholar

[ref37] 37. Miles BW, Thoden JB, Holden HM, Raushel FM. Inactivation of the amidotransferase activity of carbamoyl phosphate synthetase by the antibiotic acivicin. J Biol Chem. 2002;277(6):4368–73. pmid:11729189
View Article
PubMed/NCBI
Google Scholar

[146] View Article

[147] PubMed/NCBI

[148] Google Scholar

[ref38] 38. Pinto A, Tamborini L, Cullia G, Conti P, De Micheli C. Inspired by Nature: The 3-Halo-4,5-dihydroisoxazole Moiety as a Novel Molecular Warhead for the Design of Covalent Inhibitors. ChemMedChem. 2016;11(1):10–4. pmid:26607551
View Article
PubMed/NCBI
Google Scholar

[150] View Article

[151] PubMed/NCBI

[152] Google Scholar

[ref39] 39. Laffitte MN, Leprohon P, Papadopoulou B, Ouellette M. Plasticity of the Leishmania genome leading to gene copy number variations and drug resistance. F1000Res. 2016;5:2350. pmid:27703673
View Article
PubMed/NCBI
Google Scholar

[154] View Article

[155] PubMed/NCBI

[156] Google Scholar

Figures

Abstract

Author summary

Introduction

Methods

Cells lines, culture conditions

DNA manipulation

Pulse-field gel electrophoresis (PFGE)

Whole genome sequencing and analysis

Cellular transport assays

Glutamate and glutamine measurements in Leishmania by LC-MS

Statistical analysis

Results

Leishmania tarentolae acivicin resistant mutant and gene copy number variation

Single nucleotide polymorphisms in L. tarentolae acivicin resistant A150.1 mutant

Genes involved in acivicin resistance

Discussion

Supporting information

S1 Fig. Chromosomal amplification and deletion in L. tarentolae A150.1.

S2 Fig. Alignment of SNAT6 homologs.

S3 Fig. Phylogeny of AAT22 in Leishmania ssp.

S1 Table.

S2 Table. List of SNVs detected in the coding regions of L. tarentolae A150.1.

S3 Table. List of SNVs detected in the intergenic regions of L. tarentolae A150.1.

Acknowledgments

References