Formation of extrachromosomal circular amplicons with direct or inverted duplications in drug-resistant Leishmania tarentolae

Circular extrachromosomal DNA in Euglena gracilis under normal and stress conditions

Extrachromosomal circular DNA (eccDNA) enhances genomic plasticity, augmenting its coding and regulatory potential. Advances in high-throughput sequencing have enabled the investigation of these structural variants. Although eccDNAs have been investigated in numerous taxa, they remained understudied in euglenids. Therefore, we examined eccDNAs predicted from Illumina sequencing data of Euglena gracilis Z SAG 1224-5/25, grown under optimal photoperiod and exposed to UV irradiation. We identified approximately 1000 unique eccDNA candidates, about 20% of which were shared across conditions. We also observed a significant enrichment of mitochondrially encoded eccDNA in the UV-irradiated sample. Furthermore, we found that the heterogeneity of eccDNA was reduced in UV-exposed samples compared to cells that were grown in optimal conditions. Hence, eccDNA appears to play a role in the response to oxidative stress in Euglena, as it does in other studied organisms. In addition to contributing to the understanding of Euglena genomes, our results contribute to the validation of bioinformatics pipelines on a large, non-model genome.

Antigenic diversity in malaria parasites is maintained on extrachromosomal DNA

Sequence variation among antigenic var genes enables Plasmodium falciparum malaria parasites to evade host immunity. Using long sequence reads from haploid clones from a mutation accumulation experiment, we detect var diversity inconsistent with simple chromosomal inheritance. We discover putatively circular DNA that is strongly enriched for var genes, which exist in multiple alleles per locus separated by recombination and indel events. Extrachromosomal DNA likely contributes to rapid antigenic diversification in P. falciparum.

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.

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.

Application of CRISPR/Cas9-Mediated Genome Editing in Leishmania

CRISPR-Cas9 is an RNA guided endonuclease derived from the bacterium Streptococcus pyogenes. Due to its simplicity, versatility, and high efficiency, it has been widely used for genome editing in a variety of organisms including the protozoan parasite Leishmania, the causative agent of human leishmaniasis. Compared to the traditional homologous recombination gene targeting method, CRISPR-Cas9 has been shown to be a more efficient method to delete or disrupt Leishmania genes, generate point mutations, and add tags to endogenous genes. Notably, the stable CRISPR expression systems were shown to delete multicopy family Leishmania genes and genes present in multiploid chromosomes, identify essential Leishmania genes, and create specific chromosome translocations. In this chapter, we describe detailed procedures on using the stable CRISPR expression system for genome editing in Leishmania. These procedures include CRISPR targeting site selection, gRNA design, cloning single and double gRNA coding sequences into the Leishmania CRISPR vector pLdCN, oligonucleotide donor and drug resistance selection donor design, Leishmania cell transfection, screening, and isolation of CRISPR-edited mutants. As the principles of gene editing are generally similar, many of these procedures could also apply to the transient Leishmania CRISPR systems described by other labs.

Single-Strand Annealing Plays a Major Role in Double-Strand DNA Break Repair following CRISPR-Cas9 Cleavage in Leishmania

CRISPR-Cas9 genome editing relies on an efficient double-strand DNA break (DSB) and repair. Contrary to mammalian cells, the protozoan parasite Leishmania lacks the most efficient nonhomologous end-joining pathway and uses microhomology-mediated end joining (MMEJ) and, occasionally, homology-directed repair to repair DSBs. Here, we reveal that Leishmania predominantly uses single-strand annealing (SSA) (>90%) instead of MMEJ (<10%) for DSB repair (DSBR) following CRISPR targeting of the miltefosine transporter gene, resulting in 9-, 18-, 20-, and 29-kb sequence deletions and multiple gene codeletions. Strikingly, when targeting the Leishmania donovani LdBPK_241510 gene, SSA even occurred by using direct repeats 77 kb apart, resulting in the codeletion of 15 Leishmania genes, though with a reduced frequency. These data strongly indicate that DSBR is not efficient in Leishmania, which explains why more than half of DSBs led to cell death and why the CRISPR gene-targeting efficiency is low compared with that in other organisms. Since direct repeat sequences are widely distributed in the Leishmania genome, we predict that many DSBs created by CRISPR are repaired by SSA. It is also revealed that DNA polymerase theta is involved in both MMEJ and SSA in Leishmania Collectively, this study establishes that DSBR mechanisms and their competence in an organism play an important role in determining the outcome and efficacy of CRISPR gene targeting. These observations emphasize the use of donor DNA templates to improve gene editing specificity and efficiency in Leishmania In addition, we developed a novel Staphylococcus aureus Cas9 constitutive expression vector (pLdSaCN) for gene targeting in Leishmania IMPORTANCE Due to differences in double-strand DNA break (DSB) repair mechanisms, CRISPR-Cas9 gene editing efficiency can vary greatly in different organisms. In contrast to mammalian cells, the protozoan parasite Leishmania uses microhomology-mediated end joining (MMEJ) and, occasionally, homology-directed repair (HDR) to repair DSBs but lacks the nonhomologous end-joining pathway. Here, we show that Leishmania predominantly uses single-strand annealing (SSA) instead of MMEJ for DSB repairs (DSBR), resulting in large deletions that can include multiple genes. This strongly indicates that the overall DSBR in Leishmania is inefficient and therefore can influence the outcome of CRISPR-Cas9 gene editing, highlighting the importance of using a donor DNA to improve gene editing fidelity and efficiency in Leishmania.

Predominance of Leishmania major and rare occurrence of Leishmania tropica with haplotype variability at the center of Iran

BACKGROUND

Leishmania major is a causative agent of zoonotic cutaneous leishmaniasis in the center of Iran, Abarkouh district. Molecular characterization and precise incrimination of Leishmania species was carried out to perform controlling measurements and to design treatment programs for zoonotic cutaneous leishmaniasis.

METHODS

All smears isolated from ulcers of suspected patients were examined under a light microscope and graded for amastigotes frequency. Extraction of DNA, PCR, RFLP and sequencing of ITS-rDNA genotype were done to increase the efficacy of Leishmania parasites identification at their species-specific level and to detect any Leishmania infections within.

RESULTS

Humans were found to be infected with L. major with high infection frequency and also Leishmania tropica was identified with low occurrence for the first time as non-native species using molecular analyses. The rates of infections was considerable with microscopic observation (n = 65, 73%) out of 89 smears prepared from suspected patients. Molecular analyses showed that the density of L. major was significantly higher (n = 48, 53.93%) than L. tropica (n = 4, 4.49%) (Mann-Whitney U test: p < 0.05) and two samples (2.25%) remained ambiguous after several sequencing. L. major did not have diversity with two common haplotypes but L. tropica were found to exhibit high diversity with three novel haplotypes.

CONCLUSION

L. major was considered the causative agent of leishmaniasis in the region, but the identification of a non-native L. tropica revealed the importance of further isolation of Leishmania parasites following molecular analyses and confirmation, and also revealed the importance of further isolation of Leishmania parasites from patients of the field areas who do not have easily access to health care centers for specialized treatment strategies.

Molecular Preadaptation to Antimony Resistance in Leishmania donovani on the Indian Subcontinent

The “antibiotic resistance crisis” is a major challenge for scientists and medical professionals. This steady rise in drug-resistant pathogens also extends to parasitic diseases, with antimony being the first anti-Leishmania drug that fell in the Indian subcontinent (ISC). Leishmaniasis is a major but neglected infectious disease with limited therapeutic options. Therefore, understanding how parasites became resistant to antimonials is of commanding importance. In this study, we experimentally characterized the dynamics of this resistance acquisition and show for the first time that some Leishmania populations of the ISC were preadapted to antimony resistance, likely driven by environmental factors or by drugs used in the 19th century.

Antimonials (Sb) were used for decades for chemotherapy of visceral leishmaniasis (VL). Now abandoned in the Indian subcontinent (ISC) because of Leishmania donovani resistance, this drug offers a unique model for understanding drug resistance dynamics. In a previous phylogenomic study, we found two distinct populations of L. donovani: the core group (CG) in the Gangetic plains and ISC1 in the Nepalese highlands. Sb resistance was only encountered within the CG, and a series of potential markers were identified. Here, we analyzed the development of resistance to trivalent antimonials (Sb^III) upon experimental selection in ISC1 and CG strains. We observed that (i) baseline Sb^III susceptibility of parasites was higher in ISC1 than in the CG, (ii) time to Sb^III resistance was higher for ISC1 parasites than for CG strains, and (iii) untargeted genomic and metabolomic analyses revealed molecular changes along the selection process: these were more numerous in ISC1 than in the CG. Altogether these observations led to the hypothesis that CG parasites are preadapted to Sb^III resistance. This hypothesis was experimentally confirmed by showing that only wild-type CG strains could survive a direct exposure to the maximal concentration of Sb^III. The main driver of this preadaptation was shown to be MRPA, a gene involved in Sb^III sequestration and amplified in an intrachromosomal amplicon in all CG strains characterized so far. This amplicon emerged around 1850 in the CG, well before the implementation of antimonials for VL chemotherapy, and we discuss here several hypotheses of selective pressure that could have accompanied its emergence.

IMPORTANCE The “antibiotic resistance crisis” is a major challenge for scientists and medical professionals. This steady rise in drug-resistant pathogens also extends to parasitic diseases, with antimony being the first anti-Leishmania drug that fell in the Indian subcontinent (ISC). Leishmaniasis is a major but neglected infectious disease with limited therapeutic options. Therefore, understanding how parasites became resistant to antimonials is of commanding importance. In this study, we experimentally characterized the dynamics of this resistance acquisition and show for the first time that some Leishmania populations of the ISC were preadapted to antimony resistance, likely driven by environmental factors or by drugs used in the 19th century.

Intrachromosomal amplification, locus deletion and point mutation in the aquaglyceroporin AQP1 gene in antimony resistant Leishmania (Viannia) guyanensis

BACKGROUND

Antimony resistance complicates the treatment of infections caused by the parasite Leishmania.

METHODOLOGY/PRINCIPAL FINDINGS

Using next generation sequencing, we sequenced the genome of four independent Leishmania guyanensis antimony-resistant (SbR) mutants and found different chromosomal alterations including aneuploidy, intrachromosomal gene amplification and gene deletion. A segment covering 30 genes on chromosome 19 was amplified intrachromosomally in three of the four mutants. The gene coding for the multidrug resistance associated protein A involved in antimony resistance was also amplified in the four mutants, most likely through chromosomal translocation. All mutants also displayed a reduced accumulation of antimony mainly due to genomic alterations at the level of the subtelomeric region of chromosome 31 harboring the gene coding for the aquaglyceroporin 1 (LgAQP1). Resistance involved the loss of LgAQP1 through subtelomeric deletions in three mutants. Interestingly, the fourth mutant harbored a single G133D point mutation in LgAQP1 whose role in resistance was functionality confirmed through drug sensitivity and antimony accumulation assays. In contrast to the Leishmania subspecies that resort to extrachromosomal amplification, the Viannia strains studied here used intrachromosomal amplification and locus deletion.

CONCLUSIONS/SIGNIFICANCE

This is the first report of a naturally occurred point mutation in AQP1 in antimony resistant parasites.

Formation of linear amplicons with inverted duplications in Leishmania requires the MRE11 nuclease

Extrachromosomal DNA amplification is frequent in the protozoan parasite Leishmania selected for drug resistance. The extrachromosomal amplified DNA is either circular or linear, and is formed at the level of direct or inverted homologous repeated sequences that abound in the Leishmania genome. The RAD51 recombinase plays an important role in circular amplicons formation, but the mechanism by which linear amplicons are formed is unknown. We hypothesized that the Leishmania infantum DNA repair protein MRE11 is required for linear amplicons following rearrangements at the level of inverted repeats. The purified LiMRE11 protein showed both DNA binding and exonuclease activities. Inactivation of the LiMRE11 gene led to parasites with enhanced sensitivity to DNA damaging agents. The MRE11^−/− parasites had a reduced capacity to form linear amplicons after drug selection, and the reintroduction of an MRE11 allele led to parasites regaining their capacity to generate linear amplicons, but only when MRE11 had an active nuclease activity. These results highlight a novel MRE11-dependent pathway used by Leishmania to amplify portions of its genome to respond to a changing environment.

Extrachromosomal DNA amplification is frequent in the human protozoan parasite Leishmania when challenged with drug or other stressful conditions. DNA amplicons, either circular or linear, are formed by recombination between direct or inverted repeats spread throughout the genome of the parasite. The recombinase RAD51 is involved in the formation of circular amplicons, but the mechanism by which linear amplicons are formed is still unknown in this parasite. Studies in other organisms have provided some evidence that a DNA break is required for linear amplifications, and that the DNA repair protein MRE11 can be involved in this process. In this work, we present our biochemical, cellular and molecular characterization of the Leishmania infantum MRE11 orthologue and provide evidence that this nuclease is involved in the formation of linear amplicons in Leishmania. Our results highlight a novel MRE11-dependent pathway used by Leishmania to amplify portions of its genome to respond to a changing environment.

DNA repair pathways in trypanosomatids: from DNA repair to drug resistance

All living organisms are continuously faced with endogenous or exogenous stress conditions affecting genome stability. DNA repair pathways act as a defense mechanism, which is essential to maintain DNA integrity. There is much to learn about the regulation and functions of these mechanisms, not only in human cells but also equally in divergent organisms. In trypanosomatids, DNA repair pathways protect the genome against mutations but also act as an adaptive mechanism to promote drug resistance. In this review, we scrutinize the molecular mechanisms and DNA repair pathways which are conserved in trypanosomatids. The recent advances made by the genome consortiums reveal the complete genomic sequences of several pathogens. Therefore, using bioinformatics and genomic sequences, we analyze the conservation of DNA repair proteins and their key protein motifs in trypanosomatids. We thus present a comprehensive view of DNA repair processes in trypanosomatids at the crossroads of DNA repair and drug resistance.

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

The human parasite Leishmania uses adaptive gene rearrangements and amplification involving repeated sequences on a genome-wide scale as one strategy to adapt to a changing environment.

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.

Variations in the copy number of DNA segments account for a substantial amount of genome diversity of most organisms. DNA amplification, a contributor to copy number variation, can occur in response to various stresses or after altered growth conditions, leading to extensive and often reversible genetic variation. DNA amplification in the parasite Leishmania occurs outside the normal chromosomes and arises by DNA rearrangements involving homologous repeated sequences. We show here that such repeated sequences are widespread in the Leishmania genome and that most of the Leishmania genome is subject to stochastic gene rearrangements mediated by these low-copy repeat sequences. Thus, although cells in the population have a common core genome, many individual cells will differ from the rest of the population by carrying one or more distinct extrachromosomal amplicon. Upon selection with either drugs or culture conditions, a subpopulation can emerge where the amplicon copy number per cell increases, and this clone of cells can then expand to dominate the population. We propose that Leishmania uses adaptive gene amplification at a genome-wide scale as one strategy to adapt to a changing environment.

Antimony transport mechanisms in resistant leishmania parasites

Antimonial compounds have been used for more than a century in the treatment of the parasitic disease leishmaniasis. Although pentavalent antimonials are still first-line drugs in several developing countries, this class of drugs is no longer recommended in the Indian sub-continent because of the emergence of drug resistance. The precise mechanisms involved in the resistance of leishmania parasites to antimony are still subject to debate. It is now well documented that drug resistance in leishmania parasites is a multifactorial phenomenon involving multiple genes whose expression pattern synergistically leads to the resistance phenotype. The reduction of intracellular antimony accumulation is a frequent change observed in resistant leishmania cells; however, no comprehensive transport model has been presented so far to explain this change and its contribution to Leishmania resistance. The present review firstly covers the actual knowledge on the metabolism of antimonial drugs, the mechanisms of their transmembrane transport and intracellular processing in Leishmania. It further describes both the functional and molecular changes associated with Sb resistance in this organism. Possible transport models based on the actual knowledge are then presented, as well as their functional implications. Biophysical and pharmacological strategies are finally proposed for the precise identification of the transport pathways.

Telomeric gene deletion and intrachromosomal amplification in antimony-resistant Leishmania

Antimonials are still the mainstay of treatment against leishmaniasis but drug resistance is increasing. We carried out short read next-generation sequencing (NGS) and comparative genomic hybridization (CGH) of three independent Leishmania major antimony-resistant mutants. Copy number variations were consistently detected with both NGS and CGH. A major attribute of antimony resistance was a novel terminal deletion of variable length (67 kb to 204 kb) of the polyploid chromosome 31 in the three mutants. Terminal deletions in two mutants occurred at the level of inverted repeated sequences. The AQP1 gene coding for an aquaglyceroporin was part of the deleted region and its transfection into resistant mutants reverted resistance to SbIII. We also highlighted an intrachromosomal amplification of a subtelomeric locus on chromosome 34 in one mutant. This region encoded for ascorbate-dependent peroxidase (APX) and glucose-6-phosphate dehydrogenase (G6PDH). Overexpression of these genes in revertant backgrounds demonstrated resistance to SbIII and protection from reactive oxygen species (ROS). Generation of a G6PDH null mutant in one revertant exhibited SbIII sensitivity and a decreased protection of ROS. Our genomic analyses and functional validation highlighted novel genomic rearrangements, functionally important resistant loci and the implication of new genes in antimony resistance in Leishmania.

Whole genome sequencing of multiple Leishmania donovani clinical isolates provides insights into population structure and mechanisms of drug resistance

Visceral leishmaniasis is a potentially fatal disease endemic to large parts of Asia and Africa, primarily caused by the protozoan parasite Leishmania donovani. Here, we report a high-quality reference genome sequence for a strain of L. donovani from Nepal, and use this sequence to study variation in a set of 16 related clinical lines, isolated from visceral leishmaniasis patients from the same region, which also differ in their response to in vitro drug susceptibility. We show that whole-genome sequence data reveals genetic structure within these lines not shown by multilocus typing, and suggests that drug resistance has emerged multiple times in this closely related set of lines. Sequence comparisons with other Leishmania species and analysis of single-nucleotide diversity within our sample showed evidence of selection acting in a range of surface- and transport-related genes, including genes associated with drug resistance. Against a background of relative genetic homogeneity, we found extensive variation in chromosome copy number between our lines. Other forms of structural variation were significantly associated with drug resistance, notably including gene dosage and the copy number of an experimentally verified circular episome present in all lines and described here for the first time. This study provides a basis for more powerful molecular profiling of visceral leishmaniasis, providing additional power to track the drug resistance and epidemiology of an important human pathogen.

Antimony resistance in leishmania, focusing on experimental research

Leishmaniases are parasitic diseases that spread in many countries with a prevalence of 12 million cases. There are few available treatments and antimonials are still of major importance in the therapeutic strategies used in most endemic regions. However, resistance toward these compounds has recently emerged in areas where the replacement of these drugs is mainly limited by the cost of alternative molecules. In this paper, we reviewed the studies carried out on antimonial resistance in Leishmania. Several common limitations of these works are presented before prevalent approaches to evidence antimonial resistance are related. Afterwards, phenotypic determination of resistance is described, then confronted to clinical outcome. Finally, we detail molecular mechanisms and targets involved in resistance and already identified in vitro within selected mutant strains or in clinical isolates.

Intrachromosomal tandem duplication and repeat expansion during attempts to inactivate the subtelomeric essential gene GSH1 in Leishmania

Gamma-glutamylcysteine synthetase encoded by GSH1 is the rate-limiting enzyme in the biosynthesis of glutathione and trypanothione in Leishmania. Attempts to generate GSH1 null mutants by gene disruption failed in Leishmania infantum. Removal of even a single allele invariably led to the generation of an extra copy of GSH1, maintaining two intact wild-type alleles. In the second and even third round of inactivation, the markers integrated at the homologous locus but always preserved two intact copies of GSH1. We probed into the mechanism of GSH1 duplication. GSH1 is subtelomeric on chromosome 18 and Southern blot analysis indicated that a 10-kb fragment flanked by 466-bp direct repeated sequences was duplicated in tandem on the same chromosomal allele each time GSH1 was targeted. Polymerase chain reaction analysis and sequencing confirmed the generation of novel junctions created at the level of the 466-bp repeats consequent to locus duplication. In loss of heterozygosity attempts, the same repeated sequences were utilized for generating extrachromosomal circular amplicons. Our results are consistent with break-induced replication as a mechanism for the generation of this regional polyploidy to compensate for the inactivation of an essential gene. This chromosomal repeat expansion through repeated sequences could be implicated in locus duplication in Leishmania.

Gene expression profiling and molecular characterization of antimony resistance in Leishmania amazonensis

Background

Drug resistance is a major problem in leishmaniasis chemotherapy. RNA expression profiling using DNA microarrays is a suitable approach to study simultaneous events leading to a drug-resistance phenotype. Genomic analysis has been performed primarily with Old World Leishmania species and here we investigate molecular alterations in antimony resistance in the New World species L. amazonensis.

Methods/Principal Findings

We selected populations of L. amazonensis promastigotes for resistance to antimony by step-wise drug pressure. Gene expression of highly resistant mutants was studied using DNA microarrays. RNA expression profiling of antimony-resistant L. amazonensis revealed the overexpression of genes involved in drug resistance including the ABC transporter MRPA and several genes related to thiol metabolism. The MRPA overexpression was validated by quantitative real-time RT-PCR and further analysis revealed that this increased expression was correlated to gene amplification as part of extrachromosomal linear amplicons in some mutants and as part of supernumerary chromosomes in other mutants. The expression of several other genes encoding hypothetical proteins but also nucleobase and glucose transporter encoding genes were found to be modulated.

Conclusions/Significance

Mechanisms classically found in Old World antimony resistant Leishmania were also highlighted in New World antimony-resistant L. amazonensis. These studies were useful to the identification of resistance molecular markers.

Leishmania are unicellular microorganisms that can be transmitted to humans by the bite of sandflies. They cause a spectrum of diseases called leishmaniasis, which are classified as neglected tropical diseases by the World Health Organization. The treatment of leishmaniasis is based on the administration of antimony-containing drugs. These drugs have been used since 1947 and still constitute the mainstay for leishmaniasis treatment in several countries. One of the problems with these compounds is the emergence of resistance. Our work seeks to understand how these parasites become resistant to the drug. We studied antimony-resistant Leishmania amazonensis mutants. We analyzed gene expression at the whole genome level in antimony-resistant parasites and identified mechanisms used by Leishmania for resistance. This work could help us in developing new strategies for treatment in endemic countries where people are unresponsive to antimony-based chemotherapy. The identification of common mechanisms among different species of resistant parasites may also contribute to the development of diagnostic kits to identify and monitor the spread of resistance.

Gene expression modulation is associated with gene amplification, supernumerary chromosomes and chromosome loss in antimony-resistant Leishmania infantum

Antimonials remain the first line drug against the protozoan parasite Leishmania but their efficacy is threatened by resistance. We carried out a RNA expression profiling analysis comparing an antimony-sensitive and -resistant (Sb2000.1) strain of Leishmania infantum using whole-genome 70-mer oligonucleotide microarrays. Several genes were differentially expressed between the two strains, several of which were found to be physically linked in the genome. MRPA, an ATP-binding cassette (ABC) gene known to be involved in antimony resistance, was overexpressed in the antimony-resistant mutant along with three other tandemly linked genes on chromosome 23. This four gene locus was flanked by 1.4 kb repeated sequences from which an extrachromosomal circular amplicon was generated in the resistant cells. Interestingly, gene expression modulation of entire chromosomes occurred in the antimony-resistant mutant. Southern blots analyses and comparative genomic hybridizations revealed that this was either due to the presence of supernumerary chromosomes or to the loss of one chromosome. Leishmania parasites with haploid chromosomes were viable. Changes in copy number for some of these chromosomes were confirmed in another antimony-resistant strain. Selection of a partial revertant line correlated antimomy resistance levels and the copy number of aneuploid chromosomes, suggesting a putative link between aneuploidy and drug resistance in Leishmania.

Modulation of gene expression in drug resistant Leishmania is associated with gene amplification, gene deletion and chromosome aneuploidy

Gene expression and DNA copy number analyses using full genome oligonucleotide microarrays of Leishmania reveal molecular mechanisms of methotrexate resistance.

Background

Drug resistance can be complex, and several mutations responsible for it can co-exist in a resistant cell. Transcriptional profiling is ideally suited for studying complex resistance genotypes and has the potential to lead to novel discoveries. We generated full genome 70-mer oligonucleotide microarrays for all protein coding genes of the human protozoan parasites Leishmania major and Leishmania infantum. These arrays were used to monitor gene expression in methotrexate resistant parasites.

Results

Leishmania is a eukaryotic organism with minimal control at the level of transcription initiation and few genes were differentially expressed without concomitant changes in DNA copy number. One exception was found in Leishmania major, where the expression of whole chromosomes was down-regulated. The microarrays highlighted several mechanisms by which the copy number of genes involved in resistance was altered; these include gene deletion, formation of extrachromosomal circular or linear amplicons, and the presence of supernumerary chromosomes. In the case of gene deletion or gene amplification, the rearrangements have occurred at the sites of repeated (direct or inverted) sequences. These repeats appear highly conserved in both species to facilitate the amplification of key genes during environmental changes. When direct or inverted repeats are absent in the vicinity of a gene conferring a selective advantage, Leishmania will resort to supernumerary chromosomes to increase the levels of a gene product.

Conclusion

Aneuploidy has been suggested as an important cause of drug resistance in several organisms and additional studies should reveal the potential importance of this phenomenon in drug resistance in Leishmania.

Genome position and gene amplification

Genomic analyses of human cells expressing dihydrofolate reductase provide insight into the effects of genome position on the propensity for a drug-resistance gene to amplify in human cells.

Background

Amplifications, regions of focal high-level copy number change, lead to overexpression of oncogenes or drug resistance genes in tumors. Their presence is often associated with poor prognosis; however, the use of amplification as a mechanism for overexpression of a particular gene in tumors varies. To investigate the influence of genome position on propensity to amplify, we integrated a mutant form of the gene encoding dihydrofolate reductase into different positions in the human genome, challenged cells with methotrexate and then studied the genomic alterations arising in drug resistant cells.

Results

We observed site-specific differences in methotrexate sensitivity, amplicon organization and amplification frequency. One site was uniquely associated with a significantly enhanced propensity to amplify and recurrent amplicon boundaries, possibly implicating a rare folate-sensitive fragile site in initiating amplification. Hierarchical clustering of gene expression patterns and subsequent gene enrichment analysis revealed two clusters differing significantly in expression of MYC target genes independent of integration site.

Conclusion

These studies suggest that genome context together with the particular challenges to genome stability experienced during the progression to cancer contribute to the propensity to amplify a specific oncogene or drug resistance gene, whereas the overall functional response to drug (or other) challenge may be independent of the genomic location of an oncogene.

Drug resistance in Leishmania: similarities and differences to other organisms

The main line of defense available against parasitic protozoa is chemotherapy. Drug resistance has emerged however, as a primary obstacle to the successful treatment and control of parasitic diseases. Leishmania spp., the causative agents of leishmaniasis, have served as a useful model for studying mechanisms of drug resistance in vitro. Antimonials and amphotericin B are the first line drugs to treat Leishmania followed by pentamidine and a number of other drugs. Parasites resistant against all these classes of drugs have been selected under laboratory conditions. A multiplicity of resistance mechanisms has been detected, the most prevalent being gene amplification and transport mutations. With the tools now available, it should be possible to elucidate the mechanisms that govern drug resistance in field isolates and develop more effective chemotherapeutic agents.

Organization of H locus conserved repeats in Leishmania (Viannia) braziliensis correlates with lack of gene amplification and drug resistance

Resistance to antimonials is a major problem when treating visceral leishmaniasis in India and has already been described for New World parasites. Clinical response to meglumine antimoniate in patients infected with parasites of the Viannia sub-genus can be widely variable, suggesting the presence of mechanisms of drug resistance. In this work, we have compared L. major and L. braziliensis mutants selected in different drugs. The cross-resistance profiles of some cell lines resembled those of mutants bearing H locus amplicons. However, amplified episomal molecules were exclusively detected in L. major mutants. The analysis of the L. braziliensis H region revealed a strong conservation of gene synteny. The typical intergenic repeats that are believed to mediate the amplification of the H locus in species of the Leishmania sub-genus are partially conserved in the Viannia species. The conservation of these non-coding elements in equivalent positions in both species is indicative of their relevance within this locus. The absence of amplicons in L. braziliensis suggests that this species may not favour extra-chromosomal gene amplification as a source of phenotypic heterogeneity and fitness maintenance in changing environments.

Intrastrand annealing leads to the formation of a large DNA palindrome and determines the boundaries of genomic amplification in human cancer

Amplification of large chromosomal regions (gene amplification) is a common somatic alteration in human cancer cells and often is associated with advanced disease. A critical event initiating gene amplification is a DNA double-strand break (DSB), which is immediately followed by the formation of a large DNA palindrome. Large DNA palindromes are frequent and nonrandomly distributed in the genomes of cancer cells and facilitate a further increase in copy number. Although the importance of the formation of large DNA palindromes as a very early event in gene amplification is widely recognized, it is not known how a DSB is resolved to form a large DNA palindrome and whether any local DNA structure determines the location of large DNA palindromes. We show here that intrastrand annealing following a DNA double-strand break leads to the formation of large DNA palindromes and that DNA inverted repeats in the genome determine the efficiency of this event. Furthermore, in human Colo320DM cancer cells, a DNA inverted repeat in the genome marks the border between amplified and nonamplified DNA. Therefore, an early step of gene amplification is a regulated process that is facilitated by DNA inverted repeats in the genome.

LD1 amplifications in Leishmania

The genome of Leishmania is quite plastic. Chromosomal rearrangements and DNA amplifications are common events in all the species of the genus. Gene amplification occurs both as a mechanism of drug resistance and in the absence of drug pressure. The best known spontaneous amplification in Leishmania is the so-called LD1 family of amplicons. In the past few years there have been great advances in our knowledge of LD1 elements; here, Manuel Segovia and Ginés Ortiz review all the available data.

Formation of linear inverted repeat amplicons following targeting of an essential gene in Leishmania

Attempts to inactivate an essential gene in the protozoan parasite Leishmania have often led to the generation of extra copies of the wild-type alleles of the gene. In experiments with Leishmania tarentolae set up to disrupt the gene encoding the J-binding protein 1 (JBP1), a protein binding to the unusual base beta-D-glucosyl-hydroxymethyluracil (J) of Leishmania, we obtained JBP1 mutants containing linear DNA elements (amplicons) of approximately 100 kb. These amplicons consist of a long inverted repeat with telomeric repeats at both ends and contain either the two different targeting cassettes used to inactivate JBP1, or one cassette and one JBP1 gene. Each long repeat within the linear amplicons corresponds to sequences covering the JBP1 locus, starting at the telomeres upstream of JBP1 and ending in a approximately 220 bp sequence repeated in an inverted (palindromic) orientation downstream of the JBP1 locus. We propose that these amplicons have arisen by a template switch inside a DNA replication fork involving the inverted DNA repeats and helped by the gene targeting.

Widespread and nonrandom distribution of DNA palindromes in cancer cells provides a structural platform for subsequent gene amplification

Breakage-fusion-bridge cycles contribute to chromosome instability and generate large DNA palindromes that facilitate gene amplification in human cancers. The prevalence of large DNA palindromes in cancer is not known. Here, by using a new microarray-based approach called genome-wide analysis of palindrome formation, we show that palindromes occur frequently and are widespread in human cancers. Individual tumors seem to have a nonrandom distribution of palindromes in their genomes, and a subset of palindromic loci is associated with gene amplification. This indicates that the location of palindromes in the cancer genome can serve as a structural platform that supports subsequent gene amplification. Genome-wide analysis of palindrome formation is a new approach to identify structural chromosome aberrations associated with cancer.

The overexpression of a new ABC transporter in Leishmania is related to phospholipid trafficking and reduced infectivity

This paper reports the characterization of a new ABC transporter (LtrABC1.1), related to the human ABCA subfamily, in the protozoan parasite Leishmania tropica. LtrABC1.1 is a tandem duplicated gene flanked by inverted repeats. LtrABC1.1 is expressed mainly in the flagellar pocket of the parasite. Drug resistance studies in Leishmania overexpressing LtrABC1.1 showed the transporter not to confer resistance to a range of unrelated drugs. LtrABC1.1 appears to be involved in lipid movements across the plasma membrane of the parasite since overexpression reduces the accumulation of fluorescent phospholipid analogues. The activity of this protein may also affect membrane movement processes since secreted acid phosphatase (SAP) activity was significantly lower in promastigotes overexpressing LtrABC1.1. In vitro infection experiments with macrophages indicated LtrABC1.1-transfected parasites to be significantly less infective. Together, these results suggest that this new ABC transporter could play a role in lipid movements across the plasma membrane, and that its activity might influence vesicle trafficking. This is the first ABCA-like transporter described in unicellular eukaryotes.

Structural and functional analysis of an amplification containing a PGPA gene in a glucantime-resistant Leishmania (Viannia) guyanensis cell line

Drug resistance is a complex phenomenon in Leishmania and commonly involves gene amplification. Active efflux and metal sequestration through a P-glycoprotein have been pointed to as the major mechanisms used by drug-resistant Leishmania. A gene amplification from a glucantime-resistant Leishmania (Viannia) guyanensis cell line was characterised in an attempt to understand the mechanism of metal resistance in this pathogenic species. We show that the amplification is present as an extrachromosomal amplicon of 30 kb and contains a PGPA gene ( LgPGPA), which is overexpressed in the resistant line as shown by Northern and Western blot analyses. In addition, we gathered evidence from transfection experiments for the role of the LgPGPA gene in oxyanion resistance in L. (V.) guyanensis. Our work indicates that, in this pathogenic New World Leishmania species, amplification of the PGPA gene is the major determinant in oxyanion resistance.

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.

Mechanisms of sod2 gene amplification in Schizosaccharomyces pombe

Gene amplification in eukaryotes plays an important role in drug resistance, tumorigenesis, and evolution. The Schizosaccharomyces pombe sod2 gene provides a useful model system to analyze this process. sod2 is near the telomere of chromosome I and encodes a plasma membrane Na(+)(Li(+))/H(+) antiporter. When sod2 is amplified, S. pombe survives otherwise lethal concentrations of LiCl, and >90% of the amplified sod2 genes are found in 180- and 225-kilobase (kb) linear amplicons. The sequence of the novel joint of the 180-kb amplicon indicates that it is formed by recombination between homologous regions near the telomeres of the long arm of chromosome I and the short arm of chromosome II. The 225-kb amplicon, isolated three times more frequently than the 180-kb amplicon, is a palindrome derived from a region near the telomere of chromosome I. The center of symmetry of this palindrome contains an inverted repeat consisting of two identical 134-base pair sequences separated by a 290-base pair spacer. LiCl-resistant mutants arise 200-600 times more frequently in strains deficient for topoisomerases or DNA ligase activity than in wild-type strains, but the mutant cells contain the same amplicons. These data suggest that amplicon formation may begin with DNA lesions such as breaks. In the case of the 225-kb amplicon, the breaks may lead to a hairpin structure, which is then replicated to form a double-stranded linear amplicon, or to a cruciform structure, which is then resolved to yield the same amplicon.

Role of the locus and of the resistance gene on gene amplification frequency in methotrexate resistant Leishmania tarentolae

The protozoan parasite Leishmania resists the antifolate methotrexate (MTX) by amplifying the R locus dihydrofolate reductase-thymidylate synthase ( dhfr-ts ) gene, the H locus ptr1 pterin reductase gene, and finally by mutation in a common folate/MTX transporter. Amplification of dhfr-ts has never been observed in Leishmania tarentolae MTX resistant mutants while ptr1 amplification is common. We have selected a L.tarentolae ptr1 null mutant for MTX resistance and observed dhfr-ts amplification in this mutant demonstrating that once a preferred resistance mechanism is unavailable, a second one will take over. By introducing the ptr1 gene at the R locus and the dhfr-ts gene at the H locus by gene targeting, we investigated the role of the resistance gene and the locus on the rate of gene amplification. Transfection studies indicated that ptr1 gave higher levels of MTX resistance than dhfr-ts. Consistent with this, when ptr1 was present as part of either the H locus or the R locus it was invariably amplified, while dhfr-ts was only amplified when ptr1 was inactivated. When dhfr-ts was present in a ptr1 null background on both the H locus and the R locus, amplification from the H locus was more frequent suggesting that both the gene and the locus are determining the frequency of gene amplification in Leishmania.

Increased transport of pteridines compensates for mutations in the high affinity folate transporter and contributes to methotrexate resistance in the protozoan parasite Leishmania tarentolae

Functional cloning led to the isolation of a novel methotrexate (MTX) resistance gene in the protozoan parasite Leishmania. The gene corresponds to orfG, an open reading frame (ORF) of the LD1/CD1 genomic locus that is frequently amplified in several Leishmania stocks. A functional ORF G-green fluorescence protein fusion was localized to the plasma membrane. Transport studies indicated that ORF G is a high affinity biopterin transporter. ORF G also transports folic acid, with a lower affinity, but does not transport the drug analog MTX. Disruption of both alleles of orfG led to a mutant strain that became hypersensitive to MTX and had no measurable biopterin transport. Leishmania tarentolae MTX-resistant cells without their high affinity folate transporters have a rearranged orfG gene and increased orfG RNA levels. Overexpression of orfG leads to increased biopterin uptake and, in folate-rich medium, to increased folate uptake. MTX-resistant cells compensate for mutations in their high affinity folate/MTX transporter by overexpressing ORF G, which increases the uptake of pterins and selectively increases the uptake of folic acid, but not MTX.

A dominant negative mutation in the GIM1 gene of Leishmania donovani is responsible for defects in glycosomal protein localization

Kinetoplastid protozoa contain a unique microbody organelle called the glycosome. Several important metabolic pathways are compartmentalized within the glycosome that are found in the cytoplasm of higher eukaryotes. We have previously reported the identification of a Leishmania donovani cell line called gim1-1, in which several normally glycosomal proteins are partially mislocalized to the cytoplasm. The GIM1 gene complements the defect and restores import of proteins to the glycosome. Here we demonstrate that GIM1 encodes an integral membrane protein of the glycosome. We also report that the mutant gim1-1 allele behaves as a dominant negative mutation. Introducing the gim1-1 allele extrachromasomally led to mislocalization of a glycosomal reporter protein even in wild-type cells. Gene disruption experiments in heterozygous GIM1/gim1-1 cells showed that when the mutant gim1-1 allele was lost, cells re-established normal glycosomal protein localization. Interestingly, no disruptions of the wild-type allele were obtained. These data indicate that a dominant negative mutation in the GIM1 gene is the sole genetic lesion responsible for the glycosomal defects in gim1-1, and suggest that GIM1 is an essential gene in Leishmania.

Inversion/dimerization of plasmids mediated by inverted repeats

In contrast with earlier studies on the lambda and Escherichia coli genomes, recombination between inverted repeats on plasmids is highly efficient and shown to be recA-independent. In addition, the recombination product is exclusively a head-to-head inverted dimer. Here, we show that this recombination/rearrangement event can occur on different plasmid replicons and is not specific to the particular sequence within the inverted repeats. Transcription readthrough into the inverted repeats has little effect on this event. Genetic analysis has also indicated that most known recombination enzymes are not involved in this process. Specifically, single or double mutants defective in Holliday junction resolution systems (RuvABC and/or RecG/RusA) do not abolish this recombination/rearrangement event. This result does not support the previous models (i.e. the reciprocal-strand-switching and the cruciform-dumbbell models) in which intermediates containing Holliday junctions are proposed. Further analysis has demonstrated that the recombination/rearrangement frequency is dramatically (over 1000-fold) reduced if mismatches (2.8 %) are present within the inverted repeats. Mutations in dam, mutH and mutL genes partially or completely restored the recombination/rearrangement frequency to the level exhibited by the perfect inverted repeats, suggesting the formation of heteroduplexes during recombination/rearrangement. Sequencing analysis of the recombination/rearrangement products have indicated that the majority of the products do not involve crossing-over. We discuss a possible mechanism in which blockage of the lagging strand polymerase by a hairpin triggers recombination/rearrangement mediated by inverted repeats.

Linear amplicons as precursors of amplified circles in methotrexate-resistant Leishmania tarentolae

Gene amplification is frequently observed in Leishmania cells selected for drug resistance. By gene targeting we have tagged both alleles of the H locus of Leishmania tarentolae with the neomycin and hygromycin phosphotransferase genes ( neo and hyg ). Selection of these recombinant parasites for low level methotrexate resistance led to amplification of the H locus as part of linear amplicons. The availability of tags has permitted us to determine that both alleles can be amplified in the same cell and that chromosomal deletions are frequent. When methotrexate concentration was increased in subsequent selection steps, circles were observed in several mutants. We have introduced a hyg marker into linear amplicons to test whether the circles originated from linear amplicons. After selection with a high methotrexate concentration, circles with the hyg marker were observed, showing that circles can indeed be formed from linear amplicons. The tagging of H locus alleles permits appreciation of the extent of genetic rearrangements leading to amplicon formation in Leishmania cells selected for drug resistance.

Efflux systems and increased trypanothione levels in arsenite-resistant Leishmania

The mechanism of resistance to the metal arsenite has been studied and compared in L. mexicana, L. tropica, and L. tarentolae selected in a step by step manner for arsenite resistance. Amplification of the ABC transporter gene pgpA was found to be a frequent resistance mechanism in all species. Transfection of pgpA genes into different species indicated that both the origin of the pgpA gene and the recipient strain into which the gene is transfected seem important for resistance. An increase in the levels of trypanothione was also correlated with metal resistance in different Leishmania species. The mechanism used to increase the levels of trypanothione seems to differ, however, between the different species. This study points to a key role of transporters and thiol levels in metal resistance in Leishmania.

Organization, generation and replication of amphimeric genomes: a review

Genomes comprising a pair of separated inverted repeats and called 'amphimers' are reviewed. Amphimeric genomes are observed in a large variety of different organisms, ranging from archaebacteria to mammals. The widespread existence of amphimeric genomes in nature could be due to their particular dynamic structure. Amphimeric genomes containing long inverted segments may provide the only form in which a duplicated segment is stably retained in genomes. Amphimers are often found in amplified subgenomes, indicating that they could promote a special mechanism of DNA replication and amplification. The possible mechanisms of generation, isomerization and replication/amplification of different types of amphimeric genomes are discussed. The study of amphimeric mitochondrial petite genomes of yeast could be a good model system for the study of the role of inverted repeat sequences in genome dynamics.

A cruciform-dumbbell model for inverted dimer formation mediated by inverted repeats

Small inverted repeats (small palindromes) on plasmids have been shown to mediate a recombinational rearrangement event in Escherichia coli leading to the formation of inverted dimers (giant palindromes). This recombinational rearrangement event is efficient and independent of RecA and RecBCD. In this report, we propose a cruciform-dumbbell model to explain the inverted dimer formation mediated by inverted repeats. In this model, the inverted repeats promote the formation of a DNA cruciform which is processed by an endonuclease into a linear DNA with two hairpin loops at its ends. Upon DNA replication, this linear dumbbell-like DNA is then converted to the inverted dimer. In support of this model, linear dumbbell DNA molecules with unidirectional origin of DNA replication (ColE1 ori ) have been constructed and shown to transform E.coli efficiently resulting in the formation of the inverted dimer. The ability of linear dumbbell DNA to transform E.coli suggests that the terminal loops may be important in bypassing the requirement of DNA supercoiling for initiation of replication of the ColE1 ori.

Induction of large DNA palindrome formation in yeast: implications for gene amplification and genome stability in eukaryotes

Many amplified genes, including some oncogenes, are organized as large inverted repeats. How such giant palindromes are generated remains largely unknown. Recent studies of a palindrome in the ciliate Tetrahymena suggest a novel mechanism that requires chromosome breakage next to short inverted repeats. The prevalence of short inverted repeats in eukaryotic genomes raises the interesting possibility that this process may occur widely as a response to chromosome damage. Here we demonstrate that in Saccharomyces cerevisiae, large DNA palindromes are formed efficiently, probably by intramolecular recombination, when a double-strand break is introduced next to short inverted repeats. These results suggest a general mechanism for large palindromic DNA formation and reveal an important new source of genome instability resulting from chromosome breakage at selective sites.