Effect of large targeted deletions on the mitotic stability of an extra chromosome mediating drug resistance in Leishmania

Metabolic selection of a homologous recombination-mediated gene loss protects Trypanosoma brucei from ROS production by glycosomal fumarate reductase

The genome of trypanosomatids rearranges by using repeated sequences as platforms for amplification or deletion of genomic segments. These stochastic recombination events have a direct impact on gene dosage and foster the selection of adaptive traits in response to environmental pressure. We provide here such an example by showing that the phosphoenolpyruvate carboxykinase (PEPCK) gene knockout (Δpepck) leads to the selection of a deletion event between two tandemly arranged fumarate reductase (FRDg and FRDm2) genes to produce a chimeric FRDg-m2 gene in the Δpepck∗ cell line. FRDg is expressed in peroxisome-related organelles, named glycosomes, expression of FRDm2 has not been detected to date, and FRDg-m2 is nonfunctional and cytosolic. Re-expression of FRDg significantly impaired growth of the Δpepck∗ cells, but FRD enzyme activity was not required for this negative effect. Instead, glycosomal localization as well as the covalent flavinylation motif of FRD is required to confer growth retardation and intracellular accumulation of reactive oxygen species (ROS). The data suggest that FRDg, similar to Escherichia coli FRD, can generate ROS in a flavin-dependent process by transfer of electrons from NADH to molecular oxygen instead of fumarate when the latter is unavailable, as in the Δpepck background. Hence, growth retardation is interpreted as a consequence of increased production of ROS, and rearrangement of the FRD locus liberates Δpepck∗ cells from this obstacle. Interestingly, intracellular production of ROS has been shown to be required to complete the parasitic cycle in the insect vector, suggesting that FRDg may play a role in this process.

Identification of the centromeres of Leishmania major: revealing the hidden pieces

Leishmania affects millions of people worldwide. Its genome undergoes constitutive mosaic aneuploidy, a type of genomic plasticity that may serve as an adaptive strategy to survive distinct host environments. We previously found high rates of asymmetric chromosome allotments during mitosis that lead to the generation of such ploidy. However, the underlying molecular events remain elusive. Centromeres and kinetochores most likely play a key role in this process, yet their identification has failed using classical methods. Our analysis of the unconventional kinetochore complex recently discovered in Trypanosoma brucei (KKTs) leads to the identification of a Leishmania KKT gene candidate (LmKKT1). The GFP-tagged LmKKT1 displays "kinetochore-like" dynamics of intranuclear localization throughout the cell cycle. By ChIP-Seq assay, one major peak per chromosome is revealed, covering a region of 4 ±2 kb. We find two largely conserved motifs mapping to 14 of 36 chromosomes while a higher density of retroposons are observed in 27 of 36 centromeres. The identification of centromeres and of a kinetochore component of Leishmania chromosomes opens avenues to explore their role in mosaic aneuploidy.

Genome-wide mapping reveals single-origin chromosome replication in Leishmania, a eukaryotic microbe

BACKGROUND

DNA replication initiates on defined genome sites, termed origins. Origin usage appears to follow common rules in the eukaryotic organisms examined to date: all chromosomes are replicated from multiple origins, which display variations in firing efficiency and are selected from a larger pool of potential origins. To ask if these features of DNA replication are true of all eukaryotes, we describe genome-wide origin mapping in the parasite Leishmania.

RESULTS

Origin mapping in Leishmania suggests a striking divergence in origin usage relative to characterized eukaryotes, since each chromosome appears to be replicated from a single origin. By comparing two species of Leishmania, we find evidence that such origin singularity is maintained in the face of chromosome fusion or fission events during evolution. Mapping Leishmania origins suggests that all origins fire with equal efficiency, and that the genomic sites occupied by origins differ from related non-origins sites. Finally, we provide evidence that origin location in Leishmania displays striking conservation with Trypanosoma brucei, despite the latter parasite replicating its chromosomes from multiple, variable strength origins.

CONCLUSIONS

The demonstration of chromosome replication for a single origin in Leishmania, a microbial eukaryote, has implications for the evolution of origin multiplicity and associated controls, and may explain the pervasive aneuploidy that characterizes Leishmania chromosome architecture.

Drug resistance analysis by next generation sequencing in Leishmania

•

WGS revealed the complexity of resistance which is often polyclonal in Leishmania.

•

An impressive variation in chromosome copy numbers exist between Leishmania species.

•

The genotype heterogeneity complicates the analysis of resistance in field isolates.

The use of next generation sequencing has the power to expedite the identification of drug resistance determinants and biomarkers and was applied successfully to drug resistance studies in Leishmania. This allowed the identification of modulation in gene expression, gene dosage alterations, changes in chromosome copy numbers and single nucleotide polymorphisms that correlated with resistance in Leishmania strains derived from the laboratory and from the field. An impressive heterogeneity at the population level was also observed, individual clones within populations often differing in both genotypes and phenotypes, hence complicating the elucidation of resistance mechanisms. This review summarizes the most recent highlights that whole genome sequencing brought to our understanding of Leishmania drug resistance and likely new directions.

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.

Novel insights into genome plasticity in Eukaryotes: mosaic aneuploidy in Leishmania

Leishmania are unicellular eukaryotes that have many markedly original molecular features compared with other uni- or multicellular eukaryotes like yeasts or mammals. Genome plasticity in this parasite has been the subject of many publications, and has been associated with drug resistance or adaptability. Aneuploidy has been suspected by several authors and it is now confirmed using state-of-the-art technologies such as high-throughput DNA sequencing. The analysis of genome contents at the single cell level using fluorescence in situ hybridization (FISH) has brought a new light on the genome organization: within a cell population, every chromosome, in every cell, may be present in at least two ploidy states (being either monosomic, disomic or trisomic), and the chromosomal content varies greatly from cell to cell, thus generating a constitutive intra-strain genomic heterogeneity, here termed 'mosaic aneuploidy'. Mosaic aneuploidy deeply affects the genetics of these organisms, leading, for example, to an extreme degree of intra-strain genomic diversity, as well as to a clearance of heterozygous cells in the population without however affecting genetic heterogeneity. Second, mosaic aneuploidy might be considered as a powerful strategy evolved by the parasite for adapting to modifications of environment conditions as well as for the emergence of drug resistance. On the whole, mosaic aneuploidy may be considered as a novel mechanism for generating phenotypic diversity driven by genomic plasticity.

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.

The SIDER2 elements, interspersed repeated sequences that populate the Leishmania genomes, constitute subfamilies showing chromosomal proximity relationship

Background

Protozoan parasites of the genus Leishmania are causative agents of a diverse spectrum of human diseases collectively known as leishmaniasis. These eukaryotic pathogens that diverged early from the main eukaryotic lineage possess a number of unusual genomic, molecular and biochemical features. The completion of the genome projects for three Leishmania species has generated invaluable information enabling a direct analysis of genome structure and organization.

Results

By using DNA macroarrays, made with Leishmania infantum genomic clones and hybridized with total DNA from the parasite, we identified a clone containing a repeated sequence. An analysis of the recently completed genome sequence of L. infantum, using this repeated sequence as bait, led to the identification of a new class of repeated elements that are interspersed along the different L. infantum chromosomes. These elements turned out to be homologues of SIDER2 sequences, which were recently identified in the Leishmania major genome; thus, we adopted this nomenclature for the Leishmania elements described herein. Since SIDER2 elements are very heterogeneous in sequence, their precise identification is rather laborious. We have characterized 54 LiSIDER2 elements in chromosome 32 and 27 ones in chromosome 20. The mean size for these elements is 550 bp and their sequence is G+C rich (mean value of 66.5%). On the basis of sequence similarity, these elements can be grouped in subfamilies that show a remarkable relationship of proximity, i.e. SIDER2s of a given subfamily locate close in a chromosomal region without intercalating elements. For comparative purposes, we have identified the SIDER2 elements existing in L. major and Leishmania braziliensis chromosomes 32. While SIDER2 elements are highly conserved both in number and location between L. infantum and L. major, no such conservation exists when comparing with SIDER2s in L. braziliensis chromosome 32.

Conclusion

SIDER2 elements constitute a relevant piece in the Leishmania genome organization. Sequence characteristics, genomic distribution and evolutionarily conservation of SIDER2s are suggestive of relevant functions for these elements in Leishmania. Apart from a proved involvement in post-trancriptional mechanisms of gene regulation, SIDER2 elements could be involved in DNA amplification processes and, perhaps, in chromosome segregation as centromeric sequences.

Clonal propagation and the fast generation of karyotype diversity: An in vitro Leishmania model

In the present work we studied the karyotype stability during long-term in vitro maintenance in 3 cloned strains of Leishmania (Viannia) peruviana, Leishmania (Viannia) braziliensis and a hybrid between both species. Only the L. (V.) peruviana strain showed an unstable karyotype, even after subcloning. Four chromosomes were studied in detail, each of them characterized by homologous chromosomes of different size (heteromorphy). Variations in chromosome patterns during in vitro maintenance were rapid and discrete, involving loss of heteromorphy or appearance of additional chromosome size variants. The resulting pattern was not the same according to experimental conditions (subinoculation rate or incubation temperature), and interestingly, this was associated with differences in growth behaviour of the respective parasites. No change in total ploidy of the cells was observed by flow cytometry. We discuss several mechanisms that might account for this variation of chromosome patterns, but we favour the occurrence of aneuploidy, caused by aberrant chromosome segregation during mitosis. Our results provide insight into the generation of karyotype diversity in natural conditions and highlight the relativity of the clone concept in parasitology.

Identification of a DNA fragment that increases mitotic stability of episomal linear DNAs in Leishmania major

The centromere is a specialized region of eukaryotic chromosomes, the site of kinetochore formation, spindle attachment and regulation of chromosome segregation during mitotic and meiotic cell divisions. To identify sequences which increase mitotic stability and/or act as potential centromeres in Leishmania major, we first generated libraries of Leishmania linear artificial chromosomes (LACs) bearing 30 kb inserts of randomly selected genomic DNAs. These were introduced into parasites, and then their stability was assessed following a period of 10 passages of growth in the absence of selective pressure. Approximately 80% of the 108 transfectants tested lost their LACs promptly and only 20% of the recombinants were retained; of these six showed strong but partial stability (maintained in 30-46% of cells). Mapping and sequencing of one clone (cSC10), which confers the highest degree of maintenance, revealed the presence of a sequence that was found within another stable episome, and which is dispersed in the genome of L. major. The implications of these data to the possible mechanisms of chromosomal maintenance are discussed.

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.

Functional mapping of a trypanosome centromere by chromosome fragmentation identifies a 16-kb GC-rich transcriptional "strand-switch" domain as a major feature

Trypanosomatids are an ancient family that diverged from the main eukaryotic lineage early in evolution, which display several unique features of gene organization and expression. Although genome sequencing is now complete, the nature of centromeres in these and other parasitic protozoa has not been resolved. Here, we report the functional mapping of a centromere in the American trypanosome, Trypanosoma cruzi, a parasite with an unusual mechanism of genetic exchange that involves the generation of aneuploidy by nuclear hybridization. Using a telomere-associated chromosome fragmentation approach, we show that the region required for the mitotic stability of chromosome 3 encompasses a transcriptional "strand-switch" domain constituted by a 16-kb GC-rich island. The domain contains several degenerate retrotransposon-like insertions, but atypically, lacks the arrays of satellite repeats normally associated with centromeric regions. This unusual type of organization may represent a paradigm for centromeres in T. cruzi and other primitive eukaryotes.

The small chromosomes of Trypanosoma brucei involved in antigenic variation are constructed around repetitive palindromes

Most eukaryotic genomes contain large regions of satellite DNA. These arrays are often associated with essential chromosomal functions, but remain largely absent from genome projects because of difficulties in cloning and sequence assembly. The numerous small chromosomes of the parasite Trypanosoma brucei fall into this category, yet are critical to understanding the genome because of their role in antigenic variation. Their relatively small size, however, makes them particularly amenable to physical mapping. We have produced fine-resolution maps of 17 complete minichromosomes and partial maps of two larger intermediate-sized chromosomes. This revealed a canonical structure shared by both chromosomal classes based around a large central core of 177-bp repeats. Around the core are variable-length genic regions, the lengths of which define chromosomal class. We show the core region to be a repetitive palindrome with a single inversion point common to all the chromosomes of both classes, suggesting a mechanism of genesis for these chromosomes. Moreover, palindromy appears to be a feature of (peri)centromeres in other species that can be easily overlooked. We propose that sequence inversion is one of the higher-order sequence motifs that confer chromosomal stability.

Physiological consequences of drug resistance in Leishmania and their relevance for chemotherapy

In the early twentieth century, infectious diseases were a leading cause of death worldwide. Through the following years, morbidity and mortality caused by infectious diseases decreased considerably in the developed world, but not in the developing world, where infectious diseases remain an important reason for concern. For example, leishmaniosis has become into a serious Third World problem. This is mainly due to an increasing frequency of drug-resistance in Leishmania and an enhanced risk of co-infection with HIV. Drug-resistance is usually associated with an increased expression of specific P-glycoproteins involved in membrane transport. The present review summarizes information which shows that drug-resistance is also associated with changes in physiological events such as parasite infectivity, incorporation of metabolites, xenobiotics conjugation and traffic, intracellular metabolism, host-parasite interaction, parasite cell shape and promastigote-amastigote differentiation. Furthermore, these events may change in a coordinated manner. An understanding of these physiological events may be helpful for designing chemotherapeutic approaches to multiple cellular targets, identifying strategies to circumvent Leishmania drug-resistance and succesfully treating leishmaniosis.

Leishmania major chromosome 3 contains two long convergent polycistronic gene clusters separated by a tRNA gene

Leishmania parasites (order Kinetoplastida, family Trypanosomatidae) cause a spectrum of human diseases ranging from asymptomatic to lethal. The approximately 33.6 Mb genome is distributed among 36 chromosome pairs that range in size from approximately 0.3 to 2.8 Mb. The complete nucleotide sequence of Leishmania major Friedlin chromosome 1 revealed 79 protein-coding genes organized into two divergent polycistronic gene clusters with the mRNAs transcribed towards the telomeres. We report here the complete nucleotide sequence of chromosome 3 (384 518 bp) and an analysis revealing 95 putative protein-coding ORFs. The ORFs are primarily organized into two large convergent polycistronic gene clusters (i.e. transcribed from the telomeres). In addition, a single gene at the left end is transcribed divergently towards the telomere, and a tRNA gene separates the two convergent gene clusters. Numerous genes have been identified, including those for metabolic enzymes, kinases, transporters, ribosomal proteins, spliceosome components, helicases, an RNA-binding protein and a DNA primase subunit.

The identification of circular extrachromosomal DNA in the nuclear genome of Trypanosoma brucei

Nuclear extrachromosomal DNA elements have been identified in several kinetoplastids such as Leishmania and Trypanosoma cruzi, but never in Trypanosoma brucei. They can occur naturally or arise spontaneously as the result of sublethal drug exposure of parasites. In most cases, they are represented as circular elements and are mitotically unstable. In this study we describe the presence of circular DNA in the nucleus of Trypanosoma brucei. This novel type of DNA was termed NR-element (NlaIII repeat element). In contrast to drug-induced episomes in other kinetoplastids, the T. brucei extrachromosomal NR-element is not generated by drug selection. Furthermore, the element is stable during mitosis over many generations. Restriction analysis of tagged NR-element DNA, unusual migration patterns during pulsed field gel electrophoresis (PFGE) and CsCl/ethidium bromide equilibrium centrifugation demonstrates that the NR-element represents circular DNA. Whereas it has been found in all field isolates of the parasites we analysed, it is not detectable in some laboratory strains notably the genome reference strain 927. The DNA sequence of this element is related to a 29 bp repeat present in the subtelomeric region of VSG-bearing chromosomes of T. brucei. It has been suggested that this subtelomeric region is part of a transition zone on chromosomes separating the relatively stable telomeric repeats from the recombinationaly active region downstream of VSG genes. Therefore, we discuss a functional connection between the occurrence of this circular DNA and subtelomeric recombination events in T. brucei.

The switch region on Leishmania major chromosome 1 is not required for mitotic stability or gene expression, but appears to be essential

The Leishmania genome project reference strain, Leishmania major Friedlin, is trisomic for chromosome 1. The complete sequence of this chromosome has revealed that the genes are grouped into two large clusters of the polycistronic type, each borne by one DNA strand and located on each side of a 1.6-kb sequence often termed the switch region. Several hypotheses concerning the role of this switch region have been put forward (region of initiation of transcription for both gene clusters, origin of replication or centromeric sequence). In the present study, we have deleted this region on the three copies of chromosome 1 by sequential targeted replacements. The absence of the switch region did not alter the mitotic stability of the three deleted chromosomes. This region therefore does not appear necessary for chromosomal replication or segregation. However, during the third targeting round, which aimed at knocking out the last switch region, a fourth copy of chromosome 1 that retained this region appeared in all clones analysed. This suggests that the persistence of this switch region is necessary for parasite survival. We then showed that the presence/absence of the switch region did not act upon the expression of a resistance marker gene inserted beforehand into the left gene cluster of the same chromosomal molecule. This result suggests that the presence of this 1.6-kb sequence is not necessary for the expression of all genes on chromosome 1.

Mitotic stability of a coding DNA sequence-free version of Leishmania major chromosome 1 generated by targeted chromosome fragmentation

The deletion of a 260-kb segment containing all the coding DNA sequences (CDS) of chromosome 1 of Leishmania major Friedlin strain was performed through homologous recombination during a transfection experiment. This allowed the selection of a mutant clone containing a linear extra chromosome sizing 155 kb (XC155). The structure of XC155 was determined by restriction analysis and DNA cloning and sequencing of the gel-purified chromosome: it is made of a 'mirror' inverted duplication of the 'right' end of chromosome 1a (approximately 25 kb at each end), and in its central part of a complex tandem amplification of the linearized transfection vector containing the hygromycin resistance gene (over approximately 105 kb). No sequence of the coding region of chromosome 1 (including the 1.6-kb 'switch' region) was found. By contrast, XC155 contains two large (approximately 13 kb) clusters of tandemly repeated subtelomeric sequences (272-bp 'satellite' DNA) as well as telomeric hexamer repeats. This extra chromosome was found to be mitotically stable after >150 generations without selective pressure in vitro. Two sequence elements are considered which may have an effect on mitotic stability and participate to centromeric function in this extra chromosome: the amplification of the input vector and the 272-bp 'satellite' DNA bound by telomeric repeats.