Amplification of a new region of DNA in an unselected laboratory stock of L. tarentolae: the T region

A new DNA amplification is described from an isolate of the lizard parasite Leishmania tarentolae. This DNA is present in up to 50 copies in the Trager line of this species and present but not amplified in all other lines tested. This amplification has been named the T amplification (for Tarentolae/Trager). Restriction enzyme digestion and electrophoresis of total DNA reveal amplified fragments totalling 19 kb following staining with ethidium bromide, a finding confirmed by the use of specific hybridization probes. Much of the amplified T DNA occurs as extra-chromosomal circular molecules. No cross-hybridization was observed between the T region and other amplified DNA of Leishmania, or the maxicircle of L. tarentolae, nor was resistance to methotrexate, chloroquine or primaquine detected in the T-amplified line. Combined with our previous results showing H region amplification in 2 other unselected lab stocks, these data demonstrate the prevalence of apparently spontaneous gene amplifications in L. tarentolae.     

Amplified DNAs in laboratory stocks of Leishmania tarentolae: extrachromosomal circles structurally and functionally similar to the inverted-H-region amplification of methotrexate-resistant Leishmania major.

We describe the structure of amplified DNA that was discovered in two laboratory stocks of the protozoan parasite Leishmania tarentolae. Restriction mapping and molecular cloning revealed that a region of 42 kilobases was amplified 8- to 30-fold in these lines. Southern blot analyses of digested DNAs or chromosomes separated by pulsed-field electrophoresis showed that the amplified DNA corresponded to the H region, a locus defined originally by its amplification in methotrexate-resistant Leishmania major (S. M. Beverley, J. A. Coderre, D. V. Santi, and R. T. Schimke, Cell 38:431-439, 1984). Similarities between the amplified DNA of the two species included (i) extensive cross-hybridization; (ii) approximate conservation of sequence order; (iii) extrachromosomal localization; (iv) an overall inverted, head-to-head configuration as a circular 140-kilobase tetrameric molecule; (v) two regions of DNA sequence rearrangement, each of which was closely associated with the two centers of the inverted repeats; (vi) association with methotrexate resistance; and (vii) phenotypically conservative amplification, in which the wild-type chromosomal arrangement was retained without apparent modification. Our data showed that amplified DNA mediating drug resistance arose in unselected L. tarentolae, although the pressures leading to apparently spontaneous amplification and maintenance of the H region are not known. The simple structure and limited extent of DNA amplified in these and other Leishmania lines suggests that the study of gene amplification in Leishmania spp. offers an attractive model system for the study of amplification in cultured mammalian cells and tumors. We also introduced a method for measuring the size of large circular DNAs, using gamma-irradiation to introduce limited double-strand breaks followed by sizing of the linear DNAs by pulsed-field electrophoresis.     

Multiple drug resistance and conservative amplification of the H region in Leishmania major

Amplification of the H region has been previously observed in methotrexate (MTX)-resistant strains of Leishmania major and in unselected laboratory stocks of L. tarentolae. We now show that selection of L. major with the structurally unrelated drugs primaquine or terbinafine generated resistant lines exhibiting H region amplification and 23- and 12-fold cross-resistance to MTX, respectively. These and other drug-resistant lines bearing H region amplification also exhibited weak cross-resistance to primaquine and terbinafine, associating the amplified H region with pleiotropic resistance to MTX and other drugs. In contrast, lines selected for chloroquine or pentamidine resistance did not show H region amplification or this pattern of drug cross-resistance. The primaquine- and terbinafine-selected lines exhibited wild-type levels of dihydrofolate reductase-thymidylate synthase and normal uptake and accumulation of MTX, and the MTX resistance of these lines was not reversed by verapamil. These data suggest that the mechanism of MTX cross-resistance associated with H region amplification is novel and distinct from that mediated by overexpression of MDR genes in multidrug-resistant mammalian cells. Structural studies indicated that the amplified H region DNA in these L. major lines was largely (possibly exclusively) extra-chromosomal and consisted of circular inverted repeats joined at two DNA rearrangement junctions. Southern blot analyses showed that these rearrangement junctions were identical in four independent cell lines, suggesting that these sites are "hotspots" for DNA rearrangement. H region amplification in all of these lines was conservative, defined as retention of the chromosomal H region locus without structural alteration or reduction in copy number. This finding is consistent with an over-replication/recombination model for amplification of the H region.     

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.     

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

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

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

Selection for methotrexate resistance in Leishmania spp. is often associated with amplification of the H locus short-chain dehydrogenase-reductase gene ptr1 as part of extrachromosomal elements. Extensive sequences are always coamplified and often contain inverted duplications, most likely formed by the annealing of inverted repeats present at the H locus. By gene targeting mediated by homologous recombination, several repeated sequences were introduced in the vicinity of ptr1. Selection for methotrexate resistance in these transfectants led to ptr1 amplification as part of small circles with direct or inverted duplications whether the integrated sequences consisted of direct or inverted repeats. Hence, for a region to he amplified in L. tarentolae during drug selection, a drug resistance gene is required and must be flanked by (any) homologous repeated sequences. The distance between these repeats and their orientation will determine the length of the amplicon and whether it contains direct or inverted duplications.     

Heavy metal resistance: a new role for P-glycoproteins in Leishmania.

P-glycoproteins are responsible for multidrug resistance in tumor cell lines and are thought to have a physiologic role in exporting cellular metabolites. We now report that a P-glycoprotein gene in the H region of the trypanosomatid protozoan Leishmania confers resistance to heavy metals when present in multiple copies. The Leishmania H region is frequently amplified in drug-resistant lines and is associated with metal resistance. Leishmania expression vectors were used to introduce multiple copies of segments of the Leishmania major H region into wild-type L. major promastigotes. Only constructs bearing a segment of L. major DNA containing the P-glycoprotein lmpgpA conferred arsenite resistance. Deletional analysis of the arsenite-resistant construct mapped resistance to the lmpgpA protein coding region. Lines expressing lmpgpA showed resistance to arsenite and trivalent antimonials, but not to pentavalent antimonials, zinc, cadmium, or the typical multidrug-resistant P-glycoprotein substrates vinblastine and puromycin. Transfection of the Leishmania tarentolae P-glycoprotein homologue ltpgpA resulted in a similar resistance profile. Thus, these pgpAs represent a functionally distinct group of P-glycoproteins which exhibit a substrate specificity similar to prokaryotic heavy metal pumps. Additionally, several arguments suggest that pgpAs may play a role in the susceptibility of Leishmania to clinically utilized antimonials.     

Maxicircle (mitochondrial) genome sequence (partial) of Leishmania major: gene content, arrangement and composition compared with Leishmania tarentolae

We report 8420 bp of DNA sequence data from the maxicircle (mitochondrial) genome of Leishmania major (MHOM/SU/73/5ASKH), a much larger portion of this genome than has been reported previously from any Leishmania species infecting humans. This region contains 10 partial and complete genes: 5 protein-encoding genes (COII, COIII, ND1, ND7 and Cyt b); two ribosomal RNA subunits (12S and 9S) and three unidentified open reading frames (MURF1, MURF4 (ATPase6) and MURF5), as in the lizard-infecting species L. tarentolae. The genes from L. major exhibit 85-87% identity with those of L. tarentolae at the nucleotide level and 71-94% identity at the amino acid level. Most differences between sequences from the two species are transversions. The gene order and arrangement within the maxicircle of L. major are similar to those in L. tarentolae, but base composition and codon usage differ between the species. Codons assigned for initiation for protein-coding genes available for comparison are similar in five genes in the two species. Pre-editing was identified in some of the protein-coding genes. Short intergenic non-coding regions are also present in L. major as they are in L. tarentolae. Intergenic regions between 9S rRNA and MURF5, MURF1 and ND1 genes are G+C rich and considered to be extensive RNA editing regions. The RNA editing process is likely to be conserved in similar pattern in L. major as in L. tarentolae.     

Comparison of the maxicircle (mitochondrial) genomes of Leishmania tarentolae and Trypanosoma brucei at the level of nucleotide sequence

The entire 16.7-kilobase (kb) transcribed region of the Leishmania tarentolae maxicircle was compared to the entire 15-kb transcribed region of the Trypanosoma brucei maxicircle at the nucleotide sequence level by dot matrix analysis and by alignments of individual genes. The L. tarentolae NADH dehydrogenase subunit 1 (ND1) gene was identified in a newly obtained 2.9-kb sequence. All but two regions which flank the cytochrome b gene are highly conserved in both species. One 3.1-kb region in L. tarentolae that contains the cytochrome oxidase subunit III (COIII) gene and several open reading frames corresponds to a 2-kb sequence in T. brucei with limited sequence homology that lacks the COIII gene. Another 0.6-kb region that comprises an unidentified open reading frame (open reading frame 12) in L. tarentolae is substituted by a nonhomologous 0.4-kb open reading frame in T. brucei. A short intergenic region between the ND1 gene and the maxicircle unidentified reading frame 1 gene shows limited sequence homology, and the regions between the ND4 and ND5 genes and between the COI and ND4 genes are not conserved. All of the intergenic regions share G + C richness and a similar pattern of G versus C strand bias. 1.8 kb of the L. tarentolae divergent region (DV) and around 3 kb of the T. brucei DV were also obtained. The T. brucei DV sequences were not homologous to the L. tarentolae DV sequence but were organized in a similar fashion with tandem repeats of varying complexity.     

The divergent region of the Leishmania tarentolae kinetoplast maxicircle DNA contains a diverse set of repetitive sequences.

A 2.76 kb segment of the 12 kb divergent region of the Leishmania tarentolae kinetoplast maxicircle DNA consists almost entirely of repeated sequences. The repeats can be grouped into six families, some of which are present throughout the remainder of the divergent region. The repeats are oriented in a head-to-tail fashion with the three simplest repeats clustered into large arrays. A 47 bp palindrome and two copies of a "supercluster" of three different types of repeats are also present in the sequenced region. A sequence change in the divergent region is described for a clonal strain of L. tarentolae which was passaged continuously for several years. The repetitive sequences found in the divergent region appear to be appropriate substrates for the presumed deletion/insertion/recombination events occurring in this rapidly evolving portion of the maxicircle.     

Direct and inverted DNA repeats associated with P-glycoprotein gene amplification in drug resistant Leishmania.

The H circle of Leishmania species contains a 30 kb inverted duplication separated by two unique DNA segments, a and b. The corresponding H region of chromosomal DNA has only one copy of the duplicated DNA. We show here that the chromosomal segments a and b are flanked by inverted repeats (198 and 1241 bp) and we discuss how these repeats could lead to formation of H circles from chromosomal DNA. Selection of Leishmania tarentolae for methotrexate resistance indeed resulted in the de novo formation of circles with long inverted duplication, but two mutants selected for arsenite resistance contained new H region plasmids without such duplications. One of these plasmids appears due to a homologous recombination between two P-glycoprotein genes with a high degree of sequence homology. Our results show how the same DNA region in Leishmania may be amplified to give plasmids with or without long inverted duplications and apparently by different mechanisms.     

Sequences of six genes and several open reading frames in the kinetoplast maxicircle DNA of Leishmania tarentolae

The DNA sequence of approximately 80% of the transcribed region of the kinetoplast maxicircle DNA of Leishmania tarentolae was obtained, and structural genes were localized by comparison of the translated amino acid sequences with those of known mitochondrial genes from other organisms. By this method, the genes for cytochrome oxidase subunits I, II, and III, cytochrome b, and human mitochondrial unidentified reading frames 4 and 5 were identified. By comparing the amino acid sequences of the putative L. tarentolae genes with those of known genes, we conclude that TGA codes for tryptophan, as in most other mitochondrial systems. This is the only apparent change from the universal genetic code. The six identified structural genes show various degrees of divergence from the homologous genes in other species, with cytochrome oxidase subunit I being the most conserved and cytochrome oxidase subunit III being the least conserved. A comparison of the cytochrome b genes from L. tarentolae and Trypanosoma brucei showed that the ratio of transversions to transitions is 1:1, suggesting that these species diverged from each other more than 80 X 10(6) years ago. Several as yet unidentified open reading frames were also present in the maxicircle sequence. These data confirm that maxicircle DNA has a coding potential which typifies other mitochondrial systems.     

Gene amplification in amphotericin B-resistant Leishmania tarentolae

Two Leishmania tarentolae cells were selected step by step for resistance to the polyene antibiotic amphotericin B, a second-line drug against the parasite Leishmania. One of the mutants was cross-resistant to ketoconazole. DNA amplification was observed in both mutants. The amplicons were extrachromosomal circles and were derived from different chromosomes. In one mutant the circle was unusually stable as it remained within the cell despite numerous passages in the absence of the drug. A circumstantial link between the copy number of amplicons and the resistance levels was established. Gene transfection experiments indicated that the link between the locus amplified and the resistance levels was not straightforward and possibly several mutations act together to lead to amphotericin B resistance.     

Drug uptake and modulation of drug resistance in Leishmania by an aquaglyceroporin

Leishmaniasis is a protozoan parasitic disease that affects 12 million people worldwide. The first line choice for the treatment of this disease is antimonial drugs. In the endemic regions, resistance to this class of drugs is a major impediment to treatment. Microbes often become resistant to drugs by mutation or down-regulation of uptake systems, but the uptake system for the antimonial drugs in Leishmania is unknown. In other organisms, aquaglyceroporins have been shown to facilitate uptake of trivalent metalloids. In this study, we report the identification and characterization of aquaglyceroporins from Leishmania major (LmAQP1) and Leishmania tarentolae (LtAQP1), respectively. These Leishmania proteins have the conserved signature motifs of aquaglyceroporins. Transfection of LmAQP1 into three species of Leishmania, L. tarentolae, Leishmania infantum, and L. major, produced hypersensitivity to both As(III) and Sb(III) in all three strains. Increased production of LmAQP1 was detected by immunoblotting. Drug-resistant parasites with various mutations leading to resistance mechanisms became hypersensitive to both metalloids after expression of LmAQP1. Increased rates of uptake of As(III) or Sb(III) correlated with metalloid sensitivity of the wild type and drug-resistant transfectants. Transfection of LmAQP1 in a Pentostam-resistant field isolate also sensitized the parasite in the macrophage-associated amastigote form. One allele of LmAQP1 was disrupted in L. major, and the resulting cells became 10-fold more resistant to Sb(III). This is the first report of the uptake of a metalloid drug by an aquaglyceroporin in Leishmania, suggesting a strategy to reverse resistance in the field.     

Leishmania tarentolae: purification and characterization of tubulin and its suitability for antileishmanial drug screening

Previously, tubulin has been purified from Leishmania amazonensis and used to identify novel molecules with selective antimitotic activity. However, L. amazonensis is pathogenic and requires a relatively expensive medium for large-scale cultivation. Herein, the purification and characterization of tubulin from the non-pathogenic Leishmania tarentolae is reported, together with the sequence of alpha- and beta-tubulin from this organism. This protein was purified by sonication, diethylaminoethyl-Sepharose chromatography, and one assembly disassembly cycle in 1% overall recovery based on total cellular protein. Leishmania tarentolae tubulin was indistinguishable from the corresponding L. amazonensis protein in terms of binding affinity for dinitroaniline sulfanilamides and sensitivity to assembly inhibition by these compounds. The amino acid sequences derived from the L. tarentolae alpha- and beta-tubulin genes were 99.6 and 99.4% identical to the corresponding amino acid sequences from the Leishmania major Friedlin strain. These results indicate that tubulin from L. tarentolae is suitable for use in drug screening.     

Minicircle variable region probes for characterization of Leishmania (Viannia) species.

The minicircle molecules present in the kinetoplast DNA (kDNA) network constitute a particularly useful molecular tool because they are a multicopy target and present a variable region that differs among minicircle classes in the same network. Using the polymerase chain reaction (PCR) and a set of primers directed outwardly from the minicircle conserved region, it is possible to prepare molecular probes representing the pool of variable regions from the different minicircle classes in the kDNA. In order to examine the specificity of the minicircle variable region as hybridization probes in Leishmania (Viannia) species, such fragments were amplified from reference strains and from a panel of isolates representing the zymodeme diversity of Leishmania (Viannia) in Colombia. The size of the amplified products was conserved in Leishmania (Viannia) braziliensis, Leishmania (Viannia) guyanensis, and Leishmania (Viannia) panamensis (650 bp) and diverged in Leishmania (Viannia) equatorensis and Leishmania (Viannia) colombiensis (850 bp). The amplified products were further hybridized to variable region pools of Leishmania braziliensis, Leishmania panamensis, Leishmania guyanensis, and Leishmania equatorensis reference strains. The results obtained from the hybridization experiments support this approach as a means of defining relationships among strains. Hybridization allowed homologies to be perceived, whereas restriction fragment length analysis of the amplified products yielded strain-specific profiles. Apparently, L. (V.) equatorensis and L. (V.) colombiensis minicircle variable regions have no or only low homology with those of other Leishmania (Viannia) species, showing the divergence of those species within the subgenus.     

The ATPase subunit 6 gene sequence predicts that RNA editing is conserved between lizard- and human-infecting Leishmania.

Here we investigate the similarities in the kinetoplastid RNA editing process between human- and lizard-infecting Leishmania species. We present the sequence of the maxicircle-encoded ATPase subunit 6 gene from L. (V.) panamensis, L. (L.) mexicana and L. (L.) donovani species of human-infecting Leishmania. These represent the first available sequences of this gene from Leishmania species other than the lizard-infecting L. tarentolae. The gene sequences are highly conserved, both over the edited and unedited parts of the gene, implying that the RNA editing process is likely to be highly conserved between Leishmania species. Indeed, the first editing domain is absolutely conserved in all three Leishmania species studied and L. tarentolae. A phylogeny based on part of the ATPase subunit 6 gene placed the lizard-infecting Leishmania within the monophyletic Leishmania genus, supporting previous data which suggest that lizard- and human-infecting Leishmania species are closely related.     

The H circles of Leishmania tarentolae are a unique amplifiable system of oligomeric DNAs associated with drug resistance

We have induced drug resistance against methotrexate, an inhibitor of dihydrofolate reductase, and CB3717, an inhibitor of thymidylate synthetase in a strain of Leishmania tarentolae. The drug-resistant strains contain extrachromosomal DNA circles of 68 kilobases with a 30-kilobase inverted duplication flanked by 4- and 5 kilobase unique segments. We show that these circles are highly homologous to the drug-induced H circles of L. tropica (1). All three L. tarentolae strains analyzed contain a chromosomal copy of the H region without duplication, but two of the three strains contain extrachromosomal H circles as well, predominantly present as H circle dimers in one strain and as tetramers in the other. After induction of methotrexate resistance, monomeric circles, presumably derived from the oligomers, become the major type of circle. Our results indicate that the H region represents a genomic region that can be copied at very low frequency to yield circles by a precise, but unusual mechanism under natural conditions in wild-type cells. Although superficially analogous to the episomes of bacteria, the system is without precedent in nature.