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

Genomic and Transcriptomic Analysis for Identification of Genes and Interlinked Pathways Mediating Artemisinin Resistance in Leishmania donovani

Current therapy for visceral leishmaniasis (VL), compromised by drug resistance, toxicity, and high cost, demands for more effective, safer, and low-cost drugs. Artemisinin has been found to be an effectual drug alternative in experimental models of leishmaniasis. Comparative genome and transcriptome analysis of in vitro-adapted artesunate-resistant (K133AS-R) and -sensitive wild-type (K133WT) Leishmania donovani parasites was carried out using next-generation sequencing and single-color DNA microarray technology, respectively, to identify genes and interlinked pathways contributing to drug resistance. Whole-genome sequence analysis of K133WT vs. K133AS-R parasites revealed substantial variation among the two and identified 240 single nucleotide polymorphisms (SNPs), 237 insertion deletions (InDels), 616 copy number variations (CNVs) (377 deletions and 239 duplications), and trisomy of chromosome 12 in K133AS-R parasites. Transcriptome analysis revealed differential expression of 208 genes (fold change ≥ 2) in K133AS-R parasites. Functional categorization and analysis of modulated genes of interlinked pathways pointed out plausible adaptations in K133AS-R parasites, such as (i) a dependency on lipid and amino acid metabolism for generating energy, (ii) reduced DNA and protein synthesis leading to parasites in the quiescence state, and (iii) active drug efflux. The upregulated expression of cathepsin-L like protease, amastin-like surface protein, and amino acid transporter and downregulated expression of the gene encoding ABCG2, pteridine receptor, adenylatecyclase-type receptor, phosphoaceylglucosamine mutase, and certain hypothetical proteins are concordant with genomic alterations suggesting their potential role in drug resistance. The study provided an understanding of the molecular basis linked to artemisinin resistance in Leishmania parasites, which may be advantageous for safeguarding this drug for future use.

Ascorbate-Dependent Peroxidase (APX) from Leishmania amazonensis Is a Reactive Oxygen Species-Induced Essential Enzyme That Regulates Virulence

The molecular mechanisms underlying biological differences between two Leishmania species that cause cutaneous disease, L. major and L. amazonensis, are poorly understood. In L. amazonensis, reactive oxygen species (ROS) signaling drives differentiation of nonvirulent promastigotes into forms capable of infecting host macrophages. Tight spatial and temporal regulation of H₂O₂ is key to this signaling mechanism, suggesting a role for ascorbate-dependent peroxidase (APX), which degrades mitochondrial H₂O₂ Earlier studies showed that APX-null L. major parasites are viable, accumulate higher levels of H₂O₂, generate a greater yield of infective metacyclic promastigotes, and have increased virulence. In contrast, we found that in L. amazonensis, the ROS-inducible APX is essential for survival of all life cycle stages. APX-null promastigotes could not be generated, and parasites carrying a single APX allele were impaired in their ability to infect macrophages and induce cutaneous lesions in mice. Similar to what was reported for L. major, APX depletion in L. amazonensis enhanced differentiation of metacyclic promastigotes and amastigotes, but the parasites failed to replicate after infecting macrophages. APX expression restored APX single-knockout infectivity, while expression of catalytically inactive APX drastically reduced virulence. APX overexpression in wild-type promastigotes reduced metacyclogenesis, but enhanced intracellular survival following macrophage infection or inoculation into mice. Collectively, our data support a role for APX-regulated mitochondrial H₂O₂ in promoting differentiation of virulent forms in both L. major and L. amazonensis Our results also uncover a unique requirement for APX-mediated control of ROS levels for survival and successful intracellular replication of L. amazonensis.

Identification of divergent Leishmania (Viannia) braziliensis ecotypes derived from a geographically restricted area through whole genome analysis

Leishmania braziliensis, the main etiological agent of cutaneous leishmaniasis (CL) in Latin America, is characterized by major differences in basic biology in comparison with better-known Leishmania species. It is also associated with a high phenotypic and possibly genetic diversity that need to be more adequately defined. Here we used whole genome sequences to evaluate the genetic diversity of ten L. braziliensis isolates from a CL endemic area from Northeastern Brazil, previously classified by Multi Locus Enzyme Electrophoresis (MLEE) into ten distinct zymodemes. These sequences were first mapped using the L. braziliensis M2904 reference genome followed by identification of Single Nucleotide Polymorphisms (SNPs). A substantial level of diversity was observed when compared with the reference genome, with SNP counts ranging from ~95,000 to ~131,000 for the different isolates. When the genome data was used to infer relationship between isolates, those belonging to zymodemes Z72/Z75, recovered from forested environments, were found to cluster separately from the others, generally associated with more urban environments. Among the remaining isolates, those from zymodemes Z74/Z106 were also found to form a separate group. Phylogenetic analyses were also performed using Multi-Locus Sequence Analysis from genes coding for four metabolic enzymes used for MLEE as well as the gene sequence coding for the Hsp70 heat shock protein. All 10 isolates were firmly identified as L. braziliensis, including the zymodeme Z26 isolate previously classified as Leishmania shawi, with the clustering into three groups confirmed. Aneuploidy was also investigated but found in general restricted to chromosome 31, with a single isolate, from zymodeme Z27, characterized by extra copies for other chromosomes. Noteworthy, both Z72 and Z75 isolates are characterized by a much reduced heterozygosity. Our data is consistent with the existence of distinct evolutionary groups in the restricted area sampled and a substantial genetic diversity within L. braziliensis.

Genomic insights into virulence mechanisms of Leishmania donovani: evidence from an atypical strain

Background

Leishmaniasis is a neglected tropical disease with diverse clinical phenotypes, determined by parasite, host and vector interactions. Despite the advances in molecular biology and the availability of more Leishmania genome references in recent years, the association between parasite species and distinct clinical phenotypes remains poorly understood. We present a genomic comparison of an atypical variant of Leishmania donovani from a South Asian focus, where it mostly causes cutaneous form of leishmaniasis.

Results

Clinical isolates from six cutaneous leishmaniasis patients (CL-SL); 2 of whom were poor responders to antimony (CL-PR), and two visceral leishmaniasis patients (VL-SL) were sequenced on an Illumina MiSeq platform. Chromosome aneuploidy was observed in both groups but was more frequent in CL-SL. 248 genes differed by 2 fold or more in copy number among the two groups. Genes involved in amino acid use (LdBPK_271940) and energy metabolism (LdBPK_271950), predominated the VL-SL group with the same distribution pattern reflected in gene tandem arrays. Genes encoding amastins were present in higher copy numbers in VL-SL and CL-PR as well as being among predicted pseudogenes in CL-SL. Both chromosome and SNP profiles showed CL-SL and VL-SL to form two distinct groups. While expected heterozygosity was much higher in VL-SL, SNP allele frequency patterns did not suggest potential recent recombination breakpoints. The SNP/indel profile obtained using the more recently generated PacBio sequence did not vary markedly from that based on the standard LdBPK282A1 reference. Several genes previously associated with resistance to antimonials were observed in higher copy numbers in the analysis of CL-PR. H-locus amplification was seen in one cutaneous isolate which however did not belong to the CL-PR group.

Conclusions

The data presented suggests that intra species variations at chromosome and gene level are more likely to influence differences in tropism as well as response to treatment, and contributes to greater understanding of parasite molecular mechanisms underpinning these differences. These findings should be substantiated with a larger sample number and expression/functional studies.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5271-z) contains supplementary material, which is available to authorized users.

Genomic Analysis of Colombian Leishmania panamensis strains with different level of virulence

The establishment of Leishmania infection in mammalian hosts and the subsequent manifestation of clinical symptoms require internalization into macrophages, immune evasion and parasite survival and replication. Although many of the genes involved in these processes have been described, the genetic and genomic variability associated to differences in virulence is largely unknown. Here we present the genomic variation of four Leishmania (Viannia) panamensis strains exhibiting different levels of virulence in BALB/c mice and its application to predict novel genes related to virulence. De novo DNA sequencing and assembly of the most virulent strain allowed comparative genomics analysis with sequenced L. (Viannia) panamensis and L. (Viannia) braziliensis strains, and showed important variations at intra and interspecific levels. Moreover, the mutation detection and a CNV search revealed both base and structural genomic variation within the species. Interestingly, we found differences in the copy number and protein diversity of some genes previously related to virulence. Several machine-learning approaches were applied to combine previous knowledge with features derived from genomic variation and predict a curated set of 66 novel genes related to virulence. These genes can be prioritized for validation experiments and could potentially become promising drug and immune targets for the development of novel prophylactic and therapeutic interventions.

Insights on a putative aminoacyl-tRNA-protein transferase of Leishmania major

The N-end rule pathway leads to regulated proteolysis as an adaptive response to external stress and is ubiquitous from bacteria to mammals. In this study, we investigated a gene coding for a putative core enzyme of this post-translational regulatory pathway in Leishmania major, which may be crucial during cytodifferentiation and the environment adaptive responses of the parasite. Leucyl, phenylalanyl-tRNA protein transferase and arginyl-tRNA protein transferase are key components of this pathway in E. coli and eukaryotes, respectively. They catalyze the specific conjugation of leucine, phenylalanine or arginine to proteins containing exposed N-terminal amino acid residues, which are recognized by the machinery for the targeted proteolysis. Here, we characterized a conserved hypothetical protein coded by the LmjF.21.0725 gene in L. major. In silico analysis suggests that the LmjF.21.0725 protein is highly conserved among species of Leishmania and might belong to the Acyl CoA-N-acyltransferases (NAT) superfamily of proteins. Immunofluorescence cell imaging indicates that the cytosolic localization of the studied protein and the endogenous levels of the protein in promastigotes are barely detectable by western blotting assay. The knockout of the two alleles of LmjF.21.0725 by homologous recombination was only possible in the heterozygous transfectant expressing LmjF.21.0725 as a transgene from a plasmid. Moreover, the kinetics of loss of the plasmid in the absence of drug pressure suggests that maintenance of the gene is essential for promastigote survival. Here, evidence is provided that this putative aminoacyl tRNA-protein transferase is essential for parasite survival. The enzyme activity and corresponding post-translational regulatory pathway are yet to be investigated.

Whole genome sequencing of live attenuated Leishmania donovani parasites reveals novel biomarkers of attenuation and enables product characterization

No licensed human vaccines are currently available against leishmaniasis. Several anti-leishmanial vaccines are currently undergoing testing, including genetically modified live-attenuated parasite vaccines. Studies with live attenuated Leishmania vaccines such as centrin deleted Leishmania donovani parasites (LdCen^−/−) showed protective immunity in animal models. Such studies typically examined the biomarkers of protective immunity however the biomarkers of attenuation in the parasite preparations have not received adequate attention. As several candidate vaccines enter clinical trials, a more complete product characterization to enable maintenance of product quality will help meet regulatory requirements. Towards this goal, we have determined the complete genome sequence of LdCen^−/− and its parent strain Ld1S-2D (LdWT) and characterized the LdCen^−/− vaccine strain using bioinformatics tools. Results showed that the LdCen^−/− parasites, in addition to loss of the centrin gene, have additional deletions ranging from 350 bp to 6900 bp in non-contiguous loci on several chromosomes, most commonly in untranslated regions. We have experimentally verified a subset of these adventitious deletions that had no impact on the attenuation of the LdCen^−/− parasites. Our results identified hitherto unknown features of attenuation of virulence that could be used as markers of product quality in production lots and highlight the importance of product characterization in parasitic vaccines.

Chromosomal Translocations in the Parasite Leishmania by a MRE11/RAD50-Independent Microhomology-Mediated End Joining Mechanism

The parasite Leishmania often relies on gene rearrangements to survive stressful environments. However, safeguarding a minimum level of genome integrity is important for cell survival. We hypothesized that maintenance of genomic integrity in Leishmania would imply a leading role of the MRE11 and RAD50 proteins considering their role in DNA repair, chromosomal organization and protection of chromosomes ends in other organisms. Attempts to generate RAD50 null mutants in a wild-type background failed and we provide evidence that this gene is essential. Remarkably, inactivation of RAD50 was possible in a MRE11 null mutant that we had previously generated, providing good evidence that RAD50 may be dispensable in the absence of MRE11. Inactivation of the MRE11 and RAD50 genes led to a decreased frequency of homologous recombination and analysis of the null mutants by whole genome sequencing revealed several chromosomal translocations. Sequencing of the junction between translocated chromosomes highlighted microhomology sequences at the level of breakpoint regions. Sequencing data also showed a decreased coverage at subtelomeric locations in many chromosomes in the MRE11^-/-RAD50^-/- parasites. This study demonstrates an MRE11-independent microhomology-mediated end-joining mechanism and a prominent role for MRE11 and RAD50 in the maintenance of genomic integrity. Moreover, we suggest the possible involvement of RAD50 in subtelomeric regions stability.

The parasite Leishmania relies on gene rearrangements to survive stressful conditions. However, maintaining a minimum level of genomic integrity is crucial for cell survival. Studies in other organisms have provided evidence that the DNA repair proteins MRE11 and RAD50 are involved in chromosomes organization, protection of chromosomes ends and therefore in the maintenance of genomic integrity. In this manuscript, we present the conditional inactivation of the Leishmania infantum RAD50 gene that was only possible in MRE11 deficient cells and suggest the genetic background is crucial for RAD50 inactivation. We demonstrate the occurrence of chromosomal translocations in the MRE11 and RAD50 deficient cells and described a MRE11-independent microhomology-mediated end-joining mechanism at the level of translocation breakpoints. We also suggest a possible involvement of RAD50 in subtelomeric regions stability. Our results highlight that both MRE11 and RAD50 are important for the maintenance of genomic integrity in Leishmania.

Disclosing the essentiality of ribose-5-phosphate isomerase B in Trypanosomatids

Ribose-5-phosphate isomerase (RPI) belongs to the non-oxidative branch of the pentose phosphate pathway, catalysing the inter-conversion of D-ribose-5-phosphate and D-ribulose-5-phosphate. Trypanosomatids encode a type B RPI, whereas humans have a structurally unrelated type A, making RPIB worthy of exploration as a potential drug target. Null mutant generation in Leishmania infantum was only possible when an episomal copy of RPIB gene was provided, and the latter was retained both in vitro and in vivo in the absence of drug pressure. This suggests the gene is essential for parasite survival. Importantly, the inability to remove the second allele of RPIB gene in sKO mutants complemented with an episomal copy of RPIB carrying a mutation that abolishes isomerase activity suggests the essentiality is due to its metabolic function. In vitro, sKO promastigotes exhibited no defect in growth, metacyclogenesis or macrophage infection, however, an impairment in intracellular amastigotes' replication was observed. Additionally, mice infected with sKO mutants rescued by RPIB complementation had a reduced parasite burden in the liver. Likewise, Trypanosoma brucei is resistant to complete RPIB gene removal and mice infected with sKO mutants showed prolonged survival upon infection. Taken together our results genetically validate RPIB as a potential drug target in trypanosomatids.

A Trypanosomatid Iron Transporter that Regulates Mitochondrial Function Is Required for Leishmania amazonensis Virulence

Iron, an essential co-factor of respiratory chain proteins, is critical for mitochondrial function and maintenance of its redox balance. We previously reported a role for iron uptake in differentiation of Leishmania amazonensis into virulent amastigotes, by a mechanism that involves reactive oxygen species (ROS) production and is independent of the classical pH and temperature cues. Iron import into mitochondria was proposed to be essential for this process, but evidence supporting this hypothesis was lacking because the Leishmania mitochondrial iron transporter was unknown. Here we describe MIT1, a homolog of the mitochondrial iron importer genes mrs3 (yeast) and mitoferrin-1 (human) that is highly conserved among trypanosomatids. MIT1 expression was essential for the survival of Trypanosoma brucei procyclic but not bloodstream forms, which lack functional respiratory complexes. L. amazonensis LMIT1 null mutants could not be generated, suggesting that this mitochondrial iron importer is essential for promastigote viability. Promastigotes lacking one LMIT1 allele (LMIT1/Δlmit1) showed growth defects and were more susceptible to ROS toxicity, consistent with the role of iron as the essential co-factor of trypanosomatid mitochondrial superoxide dismutases. LMIT1/Δlmit1 metacyclic promastigotes were unable to replicate as intracellular amastigotes after infecting macrophages or cause cutaneous lesions in mice. When induced to differentiate axenically into amastigotes, LMIT1/Δlmit1 showed strong defects in iron content and function of mitochondria, were unable to upregulate the ROS-regulatory enzyme FeSOD, and showed mitochondrial changes suggestive of redox imbalance. Our results demonstrate the importance of mitochondrial iron uptake in trypanosomatid parasites, and highlight the role of LMIT1 in the iron-regulated process that orchestrates differentiation of L. amazonensis into infective amastigotes.

Leishmaniasis is a serious parasitic disease that affects 12 million people worldwide, with clinical manifestations ranging from self-healing cutaneous lesions to deadly visceralizing disease. A vaccine is not available, and new and less toxic drugs against this protozoan parasite are urgently needed. Following introduction into vertebrate hosts during a sand fly blood meal, Leishmania parasites undergo extensive changes in morphology and metabolism that are critical for adaptation to life inside host macrophages and replication as amastigotes. Earlier studies identified major events that occur during amastigote differentiation, but the signaling mechanism initiating this process remained poorly understood. Previously we demonstrated a novel role for the reactive oxygen species (ROS) H₂O₂ in initiating amastigote differentiation, a process proposed to be dependent on iron availability inside the parasite’s mitochondria. In this study we identify LMIT1, a Leishmania transmembrane protein that functions as a mitochondrial iron transporter and is conserved in other trypanosomatid protozoan parasites. Reduced LMIT1 expression impairs mitochondrial function in the infective amastigote stage, abolishing parasite virulence. Our findings identify LMIT1 as a promising new drug target, and support the conclusion that iron-dependent ROS signals generated in the mitochondria regulate differentiation of virulent Leishmania amastigotes.

Roles of Rad51 paralogs for promoting homologous recombination in Leishmania infantum

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

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.

Gene-deleted live-attenuated Trypanosoma cruzi parasites as vaccines to protect against Chagas disease

Chagas disease is a neglected tropical disease caused by the protozoan parasite Trypanosoma cruzi. This illness is now becoming global, mainly due to congenital transmission, and so far, there are no prophylactic or therapeutic vaccines available to either prevent or treat Chagas disease. Therefore, different approaches aimed at identifying new protective immunogens are urgently needed. Live vaccines are likely to be more efficient in inducing protection, but safety issues linked with their use have been raised. The development of improved protozoan genetic manipulation tools and genomic and biological information has helped to increase the safety of live vaccines. These advances have generated a renewed interest in the use of genetically attenuated parasites as vaccines against Chagas disease. This review discusses the protective capacity of genetically attenuated parasite vaccines and the challenges and perspectives for the development of an effective whole-parasite Chagas disease vaccine.

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.

LeishMicrosatDB: open source database of repeat sequences detected in six fully sequenced Leishmania genomes

A Leishmania Microsatellite Database (LeishMicrosatDB) is reported for genome wise mining of microsatellites in six Leishmania species, using in silico techniques. This was created to provide parasitologists a platform to understand the genome characterization, mapping, phylogeny and evolutionary analysis. The present version of the database contains 1,738,669 simple sequence repeats of which 181 s756 repeats are present in compound form. The repeats can be sought in a chromosome using input parameters such as repeat type (mono- hexa), coding status, repeat unit length and repeat sequence motif. The genic repeats have been further hyperlinked with their corresponding locus id, and the database is appended with primer3 plus for primer designing of selected repeats with left and right flanking sequences up to 250 bp. Information on clustering and polymorphic repeats can also be retrieved. This database may also be adopted as a tool to study the relative occurrence and distribution of microsatellites across the parasitic genome. The database can enable a biologist to select markers at desired intervals over the chromosomes, and can be accessed as an open source repository at http://biomedinformri.com/leishmicrosat.

DATABASE URL

http://biomedinformri.com/leishmicrosat.

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.

Genomic analysis of sequence-dependent DNA curvature in Leishmania

Leishmania major is a flagellated protozoan parasite of medical importance. Like other members of the Trypanosomatidae family, it possesses unique mechanisms of gene expression such as constitutive polycistronic transcription of directional gene clusters, gene amplification, mRNA trans-splicing, and extensive editing of mitochondrial transcripts. The molecular signals underlying most of these processes remain under investigation. In order to investigate the role of DNA secondary structure signals in gene expression, we carried out a genome-wide in silico analysis of the intrinsic DNA curvature. The L. major genome revealed a lower frequency of high intrinsic curvature regions as well as inter- and intra- chromosomal distribution heterogeneity, when compared to prokaryotic and eukaryotic organisms. Using a novel method aimed at detecting region-integrated intrinsic curvature (RIIC), high DNA curvature was found to be associated with regions implicated in transcription initiation. Those include divergent strand-switch regions between directional gene clusters and regions linked to markers of active transcription initiation such as acetylated H3 histone, TRF4 and SNAP50. These findings suggest a role for DNA curvature in transcription initiation in Leishmania supporting the relevance of DNA secondary structures signals.

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.

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.

Genetically attenuated Trypanosoma cruzi parasites as a potential vaccination tool

Chagas disease is the clinical manifestation of the infection produced by the parasite Trypanosoma cruzi. Currently there is no vaccine to prevent this disease and the protection attained with vaccines containing non-replicating parasites is limited. Genetically attenuated trypanosomatid parasites can be obtained by deletion of selected genes. Gene deletion takes advantage of the fact that this parasite can undergo homologous recombination between endogenous and foreign DNA sequences artificially introduced in the cells. This approach facilitated the discovery of several unknown gene functions, as well as allowing us to speculate about the potential for genetically attenuated live organisms as experimental immunogens. Vaccination with live attenuated parasites has been used effectively in mice to reduce parasitemia and histological damage, and in dogs, to prevent vector-delivered infection in the field. However, the use of live parasites as immunogens is controversial due to the risk of reversion to a virulent phenotype. Herein, we present our results from experiments on genetic manipulation of two T. cruzi strains to produce parasites with impaired replication and infectivity, and using the mutation of the dhfr-ts gene as a safety device against reversion to virulence.