FISH analysis reveals aneuploidy and continual generation of chromosomal mosaicism in Leishmania major

Role of PIP39 in oxidative stress response appears conserved in kinetoplastids

Kinetoplastids, a group of flagellated protists that are often insect intestinal parasites, encounter various sources of oxidative stress. Such stressors include reactive oxygen species, both internally produced within the protist, and induced externally by host immune responses. This investigation focuses on the role of a highly conserved aspartate-based protein phosphatase, PTP-Interacting protein (PIP39) in managing oxidative stress. In addition to its well accepted role in a Trypanosoma brucei life stage transition, there is evidence of PIP39 participation in the T. brucei oxidative stress response. To examine whether this latter PIP39 role may exist more broadly, we aimed to elucidate PIP39's contribution to redox homeostasis in the monoxenous parasite Leptomonas seymouri. Utilizing CRISPR-Cas9-mediated elimination of PIP39 in conjunction with oxidative stress assays, we demonstrate that PIP39 is required for cellular tolerance to oxidative stress in L. seymouri, positing it as a putative regulatory node for adaptive stress responses. We propose that future analysis of L. seymouri PIP39 enzymatic activity, regulation, and potential localization to a specialized organelle termed a glycosome will contribute to a deeper understanding of the molecular mechanisms by which protozoan parasites adapt to oxidative environments. Our study also demonstrates success at using gene editing tools developed for Leishmania for the related L. seymouri.

The intersection of host in vivo metabolism and immune responses to infection with kinetoplastid and apicomplexan parasites

SUMMARYProtozoan parasite infection dramatically alters host metabolism, driven by immunological demand and parasite manipulation strategies. Immunometabolic checkpoints are often exploited by kinetoplastid and protozoan parasites to establish chronic infection, which can significantly impair host metabolic homeostasis. The recent growth of tools to analyze metabolism is expanding our understanding of these questions. Here, we review and contrast host metabolic alterations that occur in vivo during infection with Leishmania, trypanosomes, Toxoplasma, Plasmodium, and Cryptosporidium. Although genetically divergent, there are commonalities among these pathogens in terms of metabolic needs, induction of the type I immune responses required for clearance, and the potential for sustained host metabolic dysbiosis. Comparing these pathogens provides an opportunity to explore how transmission strategy, nutritional demand, and host cell and tissue tropism drive similarities and unique aspects in host response and infection outcome and to design new strategies to treat disease.

Recent Advances in CRISPR/Cas9-Mediated Genome Editing in Leishmania Strains

BACKGROUND

Genome manipulation of Leishmania species and the creation of modified strains are widely employed strategies for various purposes, including gene function studies, the development of live attenuated vaccines, and the engineering of host cells for protein production.

OBJECTIVE

Despite the introduction of novel manipulation approaches like CRISPR/Cas9 technology with significant advancements in recent years, the development of a reliable protocol for efficiently and precisely altering the genes of Leishmania strains remains a challenging endeavor. Following the successful adaptation of the CRISPR/Cas9 system for higher eukaryotic cells, several research groups have endeavored to apply this system to manipulate the genome of Leishmania.

RESULTS

Despite the substantial differences between Leishmania and higher eukaryotes, the CRISPR/Cas9 system has been effectively tested and applied in Leishmania.  CONCLUSION: This comprehensive review summarizes all the CRISPR/Cas9 systems that have been employed in Leishmania, providing details on their methods and the expression systems for Cas9 and gRNA. The review also explores the various applications of the CRISPR system in Leishmania, including the deletion of multicopy gene families, the development of the Leishmania vaccine, complete gene deletions, investigations into chromosomal translocations, protein tagging, gene replacement, large-scale gene knockout, genome editing through cytosine base replacement, and its innovative use in the detection of Leishmania. In addition, the review offers an up-to-date overview of all double-strand break repair mechanisms in Leishmania.

The adaptive roles of aneuploidy and polyclonality in Leishmania in response to environmental stress

Aneuploidy is generally considered harmful, but in some microorganisms, it can act as an adaptive mechanism against environmental stress. Here, we use Leishmania-a protozoan parasite with remarkable genome plasticity-to study the early steps of aneuploidy evolution under high drug pressure (using antimony or miltefosine as stressors). By combining single-cell genomics, lineage tracing with cellular barcodes, and longitudinal genome characterization, we reveal that aneuploidy changes under antimony pressure result from polyclonal selection of pre-existing karyotypes, complemented by further and rapid de novo alterations in chromosome copy number along evolution. In the case of miltefosine, early parasite adaptation is associated with independent point mutations in a miltefosine transporter gene, while aneuploidy changes only emerge later, upon exposure to increased drug levels. Therefore, polyclonality and genome plasticity are hallmarks of parasite adaptation, but the scenario of aneuploidy dynamics depends on the nature and strength of the environmental stress as well as on the existence of other pre-adaptive mechanisms.

Visceral Leishmaniasis Caused by Leishmania Tropica

PURPOSE

In Turkey, the main causative agent of visceral leishmaniasis (VL) is Leishmania. infantum and the main causative agent of cutaneous leishmaniasis (CL) is Leishmania tropica. In this study, we aimed to discuss the possible mechanisms, clinical aspects, and threat of visceralizing L. tropica.

METHODS

This study includes seven cases of VL caused by L. tropica.Five patients were male (71%) and four were adults (57%).

RESULTS

All the VL patients complained of fever and splenomegaly. Fatigue, pancytopenia, and hepatomegaly were present in six patients each (86%), while weight loss and gastrointestinal system (GIS) symptoms were present in 5 patients (71%).

CONCLUSIONS

In this study, we have evaluated seven cases of visceralized L. tropica (VLT) in the context of the changing leishmaniasis epidemiology in Turkey. We have evaluated the possible mechanisms of visceralization; inter- and intraspecies genetic exchange with all the old world leishmaniasis agents present in the region, stress induced by inappropriate use of drugs, and possible ongoing adaptation mechanisms of Leishmania spp. The threat posed by VLT is significant as L. tropica is the most widespread and most common cause of leishmaniasis in Turkey. We do not know the vectorial capacity of the sand flies for the transmission of VLT strains or if these strains are in circulation in Turkey. Future studies should be carried out to investigate these issues as the transition of L. tropica from a mild disease-causing agent to a mortal one poses a significant public health concern for Turkey and Europe.

Drug resistance in Leishmania: does it really matter?

Treatment failure (TF) jeopardizes the management of parasitic diseases, including leishmaniasis. From the parasite's point of view, drug resistance (DR) is generally considered as central to TF. However, the link between TF and DR, as measured by in vitro drug susceptibility assays, is unclear, some studies revealing an association between treatment outcome and drug susceptibility, others not. Here we address three fundamental questions aiming to shed light on these ambiguities. First, are the right assays being used to measure DR? Second, are the parasites studied, which are generally those that adapt to in vitro culture, actually appropriate? Finally, are other parasite factors - such as the development of quiescent forms that are recalcitrant to drugs - responsible for TF without DR?

The paradigm of intracellular parasite survival and drug resistance in leishmanial parasite through genome plasticity and epigenetics: Perception and future perspective

Leishmania is an intracellular, zoonotic, kinetoplastid eukaryote with more than 1.2 million cases all over the world. The leishmanial chromosomes are divided into polymorphic chromosomal ends, conserved central domains, and antigen-encoding genes found in telomere-proximal regions. The genome flexibility of chromosomal ends of the leishmanial parasite is known to cause drug resistance and intracellular survival through the evasion of host defense mechanisms. Therefore, in this review, we discuss the plasticity of Leishmania genome organization which is the primary cause of drug resistance and parasite survival. Moreover, we have not only elucidated the causes of such genome plasticity which includes aneuploidy, epigenetic factors, copy number variation (CNV), and post-translation modification (PTM) but also highlighted their impact on drug resistance and parasite survival.

Life in plastic, it's fantastic! How Leishmania exploit genome instability to shape gene expression

Leishmania are kinetoplastid pathogens that cause leishmaniasis, a debilitating and potentially life-threatening infection if untreated. Unusually, Leishmania regulate their gene expression largely post-transcriptionally due to the arrangement of their coding genes into polycistronic transcription units that may contain 100s of functionally unrelated genes. Yet, Leishmania are capable of rapid and responsive changes in gene expression to challenging environments, often instead correlating with dynamic changes in their genome composition, ranging from chromosome and gene copy number variations to the generation of extrachromosomal DNA and the accumulation of point mutations. Typically, such events indicate genome instability in other eukaryotes, coinciding with genetic abnormalities, but for Leishmania, exploiting these products of genome instability can provide selectable substrates to catalyse necessary gene expression changes by modifying gene copy number. Unorthodox DNA replication, DNA repair, replication stress factors and DNA repeats are recognised in Leishmania as contributors to this intrinsic instability, but how Leishmania regulate genome plasticity to enhance fitness whilst limiting toxic under- or over-expression of co-amplified and co-transcribed genes is unclear. Herein, we focus on fresh, and detailed insights that improve our understanding of genome plasticity in Leishmania. Furthermore, we discuss emerging models and factors that potentially circumvent regulatory issues arising from polycistronic transcription. Lastly, we highlight key gaps in our understanding of Leishmania genome plasticity and discuss future studies to define, in higher resolution, these complex regulatory interactions.

Identification of Leishmania infantum Epidemiology, Drug Resistance and Pathogenicity Biomarkers with Nanopore Sequencing

The emergence of drug-resistant strains of the parasite Leishmania infantum infecting dogs and humans represents an increasing threat. L. infantum genomes are complex and unstable with extensive structural variations, ranging from aneuploidies to multiple copy number variations (CNVs). These CNVs have recently been validated as biomarkers of Leishmania concerning virulence, tissue tropism, and drug resistance. As a proof-of-concept to develop a novel diagnosis platform (LeishGenApp), four L. infantum samples from humans and dogs were nanopore sequenced. Samples were epidemiologically typed within the Mediterranean L. infantum group, identifying members of the JCP5 and non-JCP5 subgroups, using the conserved region (CR) of the maxicircle kinetoplast. Aneuploidies were frequent and heterogenous between samples, yet only chromosome 31 tetrasomy was common between all the samples. A high frequency of aneuploidies was observed for samples with long passage history (MHOM/TN/80/IPT-1), whereas fewer were detected for samples maintained in vivo (MCRI/ES/2006/CATB033). Twenty-two genes were studied to generate a genetic pharmacoresistance profile against miltefosine, allopurinol, trivalent antimonials, amphotericin, and paromomycin. MHOM/TN/80/IPT-1 and MCRI/ES/2006/CATB033 displayed a genetic profile with potential resistance against miltefosine and allopurinol. Meanwhile, MHOM/ES/2016/CATB101 and LCAN/ES/2020/CATB102 were identified as potentially resistant against paromomycin. All four samples displayed a genetic profile for resistance against trivalent antimonials. Overall, this proof-of-concept revealed the potential of nanopore sequencing and LeishGenApp for the determination of epidemiological, drug resistance, and pathogenicity biomarkers in L. infantum.

Sex in protists: A new perspective on the reproduction mechanisms of trypanosomatids

The Protist kingdom individuals are the most ancestral representatives of eukaryotes. They have inhabited Earth since ancient times and are currently found in the most diverse environments presenting a great heterogeneity of life forms. The unicellular and multicellular algae, photosynthetic and heterotrophic organisms, as well as free-living and pathogenic protozoa represents the protist group. The evolution of sex is directly associated with the origin of eukaryotes being protists the earliest protagonists of sexual reproduction on earth. In eukaryotes, the recombination through genetic exchange is a ubiquitous mechanism that can be stimulated by DNA damage. Scientific evidences support the hypothesis that reactive oxygen species (ROS) induced DNA damage can promote sexual recombination in eukaryotes which might have been a decisive factor for the origin of sex. The fact that some recombination enzymes also participate in meiotic sex in modern eukaryotes reinforces the idea that sexual reproduction emerged as consequence of specific mechanisms to cope with mutations and alterations in genetic material. In this review we will discuss about origin of sex and different strategies of evolve sexual reproduction in some protists such that cause human diseases like malaria, toxoplasmosis, sleeping sickness, Chagas disease, and leishmaniasis.

Four layer multi-omics reveals molecular responses to aneuploidy in Leishmania

Aneuploidy causes system-wide disruptions in the stochiometric balances of transcripts, proteins, and metabolites, often resulting in detrimental effects for the organism. The protozoan parasite Leishmania has an unusually high tolerance for aneuploidy, but the molecular and functional consequences for the pathogen remain poorly understood. Here, we addressed this question in vitro and present the first integrated analysis of the genome, transcriptome, proteome, and metabolome of highly aneuploid Leishmania donovani strains. Our analyses unambiguously establish that aneuploidy in Leishmania proportionally impacts the average transcript- and protein abundance levels of affected chromosomes, ultimately correlating with the degree of metabolic differences between closely related aneuploid strains. This proportionality was present in both proliferative and non-proliferative in vitro promastigotes. However, as in other Eukaryotes, we observed attenuation of dosage effects for protein complex subunits and in addition, non-cytoplasmic proteins. Differentially expressed transcripts and proteins between aneuploid Leishmania strains also originated from non-aneuploid chromosomes. At protein level, these were enriched for proteins involved in protein metabolism, such as chaperones and chaperonins, peptidases, and heat-shock proteins. In conclusion, our results further support the view that aneuploidy in Leishmania can be adaptive. Additionally, we believe that the high karyotype diversity in vitro and absence of classical transcriptional regulation make Leishmania an attractive model to study processes of protein homeostasis in the context of aneuploidy and beyond.

Genome deletions to overcome the directed loss of gene function in Leishmania

With the global reach of the Neglected Tropical Disease leishmaniasis increasing, coupled with a tiny armory of therapeutics which all have problems with resistance, cost, toxicity and/or administration, the validation of new drug targets in the causative insect vector borne protozoa Leishmania spp is more important than ever. Before the introduction of CRISPR Cas9 technology in 2015 genetic validation of new targets was carried out largely by targeted gene knockout through homologous recombination, with the majority of genes targeted (~70%) deemed non-essential. In this study we exploit the ready availability of whole genome sequencing technology to reanalyze one of these historic cell lines, a L. major knockout in the catalytic subunit of serine palmitoyltransferase (LCB2), which causes a complete loss of sphingolipid biosynthesis but remains viable and infective. This revealed a number of Single Nucleotide Polymorphisms, but also the complete loss of several coding regions including a gene encoding a putative ABC3A orthologue, a putative sterol transporter. Hypothesizing that the loss of such a transporter may have facilitated the directed knockout of the catalytic subunit of LCB2 and the complete loss of de novo sphingolipid biosynthesis, we re-examined LCB2 in a L. mexicana line engineered for straightforward CRISPR Cas9 directed manipulation. Strikingly, LCB2 could not be knocked out indicating essentiality. However, simultaneous deletion of LCB2 and the putative ABC3A was possible. This indicated that the loss of the putative ABC3A facilitated the loss of sphingolipid biosynthesis in Leishmania, and suggested that we should re-examine the many other Leishmania knockout lines where genes were deemed non-essential.

Genomic diversity and genetic variation of Leishmania panamensis within its endemic range

Species belonging to the Leishmania (Viannia) subgenus are important causative agents of cutaneous and mucocutaneous leishmaniasis in Central and South America. These parasites possess several distinctive biological features that are influenced by their genetics, population structure, and genome instability. To date, several studies have revealed varying degrees of genetic diversity within Leishmania species. Particularly, in species of the L. (Viannia) subgenus, a generalized high intraspecific genetic diversity has been reported, although, conflicting conclusions have been drawn using different molecular techniques. Despite being the most common Leishmania species circulating in Panama and Colombia, few studies have analyzed clinical samples of Leishmania panamensis using whole-genome sequencing, and their restricted number of samples has limited the information they can provide to understand the population structure of L. panamensis. Here, we used next generation sequencing (NGS) to explore the genetic diversity of L. panamensis within its endemic range, analyzing data from 43 isolates of Colombian and Panamanian origin. Our results show the occurrence of three well-defined geographically correlated groups, and suggests the possible occurrence of additional phylogeographic groups. Furthermore, these results support the existence of a mixed mode of reproduction in L. panamensis, with varying frequencies of events of genetic recombination occurring primarily within subpopulations of closely related strains. This study offers important insights into the population genetics and reproduction mode of L. panamensis, paving the way to better understand their population structure and the emergence and maintenance of key eco-epidemiological traits.

Replacement of Leishmania (Leishmania) infantum Populations in an Endemic Focus of Visceral Leishmaniasis in Brazil

Visceral leishmaniasis is an important global health problem with an estimated of 50,000 to 90,000 new cases per year. VL is the most serious form of leishmaniasis as it can be fatal in 95% of the cases if it remains untreated. VL is a particularly acute problem in Brazil which contributed with 97% of all cases reported in 2020 in the Americas. In this country, VL affects mainly the poorest people in both urban and rural areas and continues to have a high mortality rate estimated around 8.15%. Here, we performed a temporal parasite population study using whole genome sequence data from a set of 34 canine isolates sampled in 2008, 2012 and 2015 from a re-emergent focus in Southeastern Brazil. Our study found the presence of two distinct sexual subpopulations that corresponded to two isolation periods. These subpopulations diverged hundreds of years ago with no apparent gene flow between them suggesting a process of rapid replacement during a two-year period. Sequence comparisons and analysis of nucleotide diversity also showed evidence of balancing selection acting on transport-related genes and antigenic families. To our knowledge this is the first population genomic study showing a turn-over of parasite populations in an endemic region for leishmaniasis. The complexity and rapid adaptability of these parasites pose new challenges to control activities and demand more integrated approaches to understand this disease in New World foci.

Genome plasticity driven by aneuploidy and loss of heterozygosity in Trypanosoma cruzi

Trypanosoma cruzi the causative agent of Chagas disease shows a marked genetic diversity and divided into at least six Discrete Typing Units (DTUs). High intra genetic variability has been observed in the TcI DTU, the most widely distributed DTU, where patterns of genomic diversity can provide information on ecological and evolutionary processes driving parasite population structure and genome organization. Chromosomal aneuploidies and rearrangements across multigene families represent an evidence of T. cruzi genome plasticity. We explored genomic diversity among 18 Colombian T. cruzi I clones and 15 T. cruzi I South American strains. Our results confirm high genomic variability, heterozygosity and presence of a clade compatible with the TcI_dom genotype, described for strains from humans in Colombia and Venezuela. TcI showed high structural plasticity across the geographical region studied. Differential events of whole and segmental aneuploidy (SA) along chromosomes even between clones from the same strain were found and corroborated by the depth and allelic frequency. We detected loss of heterozygosity (LOH) events in different chromosomes, however, the size and location of segments under LOH varied between clones. Genes adjacent to breakpoints were evaluated, and retrotransposon hot spot genes flanked the beginning of segmental aneuploidies. Our results suggest that T. cruzi genomes, like those of Leishmania, may have a highly unstable structure and there is now an urgent need to design experiments to explore any potential adaptive role for the plasticity observed.

Trypanosoma cruzi Genomic Variability: Array Comparative Genomic Hybridization Analysis of Clone and Parental Strain

Trypanosoma cruzi, the etiological agent of Chagas disease, exhibits extensive inter- and intrastrain genetic diversity. As we have previously described, there are some genetic differences between the parental G strain and its clone D11, which was isolated by the limiting dilution method and infection of cultured mammalian cells. Electrophoretic karyotyping and Southern blot hybridization of chromosomal bands with specific markers revealed chromosome length polymorphisms of small size with additional chromosomal bands in clone D11 and the maintenance of large syntenic groups. Both G strain and clone D11 belong to the T. cruzi lineage TcI. Here, we designed intraspecific array-based comparative genomic hybridization (aCGH) to identify chromosomal regions harboring copy-number variations between clone D11 and the G strain. DNA losses were more extensive than DNA gains in clone D11. Most alterations were flanked by repeated sequences from multigene families that could be involved in the duplication and deletion events. Several rearrangements were detected by chromoblot hybridization and confirmed by aCGH. We have integrated the information of genomic sequence data obtained by aCGH to the electrophoretic karyotype, allowing the reconstruction of possible recombination events that could have generated the karyotype of clone D11. These rearrangements may be explained by unequal crossing over between sister or homologous chromatids mediated by flanking repeated sequences and unequal homologous recombination via break-induced replication. The genomic changes detected by aCGH suggest the presence of a dynamic genome that responds to environmental stress by varying the number of gene copies and generating segmental aneuploidy.

Axenic interspecies and intraclonal hybrid formation in Leishmania: Successful crossings between visceral and cutaneous strains

Diseases caused by trypanosomatids are serious public health concerns in low-income endemic countries. Leishmaniasis is presented in two main clinical forms, visceral leishmaniasis-caused by L. infantum and L. donovani-and cutaneous leishmaniasis-caused by many species, including L. major, L. tropica and L. braziliensis. As for certain other trypanosomatids, sexual reproduction has been confirmed in these parasites, and formation of hybrids can contribute to virulence, drug resistance or adaptation to the host immune system. In the present work, the capability of intraclonal and interspecies genetic exchange has been investigated using three parental strains: L. donovani, L. tropica and L. major, which have been engineered to express different fluorescent proteins and antibiotic resistance markers in order to facilitate the phenotypic selection of hybrid parasites after mating events. Stationary and exponential-phase promastigotes of each species were used, in in vitro experiments, some of them containing LULO cells (an embryonic cell line derived from Lutzomyia longipalpis). Several intraclonal hybrids were obtained with L. tropica as crossing progenitor, but not with L. donovani or L. major. In interspecies crossings, three L. donovani x L. major hybrids and two L. donovani x L. tropica hybrids were isolated, thereby demonstrating the feasibility to obtain in vitro hybrids of parental lines causing different tropism of leishmaniasis. Ploidy analysis revealed an increase in DNA content in all hybrids compared to the parental strains, and nuclear analysis showed that interspecies hybrids are complete hybrids, i.e. each of them showing at least one chromosomal set from each parental. Regarding kDNA inheritance, discrepancies were observed between maxi and minicircle heritage. Finally, phenotypic studies showed either intermediate phenotypes in terms of growth profiles, or a decreased in vitro infection capacity compared to the parental cells. To the best of our knowledge, this is the first time that in vitro interspecies outcrossing has been demonstrated between Leishmania species with different tropism, thus contributing to shed light on the mechanisms underlying sexual reproduction in these parasites.

Nitric Oxide Resistance in Leishmania (Viannia) braziliensis Involves Regulation of Glucose Consumption, Glutathione Metabolism and Abundance of Pentose Phosphate Pathway Enzymes

In American Tegumentary Leishmaniasis production of cytokines, reactive oxygen species and nitric oxide (NO) by host macrophages normally lead to parasite death. However, some Leishmania braziliensis strains exhibit natural NO resistance. NO-resistant strains cause more lesions and are frequently more resistant to antimonial treatment than NO-susceptible ones, suggesting that NO-resistant parasites are endowed with specific mechanisms of survival and persistence. To tests this, we analyzed the effect of pro- and antioxidant molecules on the infectivity in vitro of L. braziliensis strains exhibiting polar phenotypes of resistance or susceptibility to NO. In addition, we conducted a comprehensive quantitative mass spectrometry-based proteomics analysis of those parasites. NO-resistant parasites were more infective to peritoneal macrophages, even in the presence of high levels of reactive species. Principal component analysis of protein concentration values clearly differentiated NO-resistant from NO-susceptible parasites, suggesting that there are natural intrinsic differences at molecular level among those strains. Upon NO exposure, NO-resistant parasites rapidly modulated their proteome, increasing their total protein content and glutathione (GSH) metabolism. Furthermore, NO-resistant parasites showed increased glucose analogue uptake, and increased abundance of phosphotransferase and G6PDH after nitrosative challenge, which can contribute to NADPH pool maintenance and fuel the reducing conditions for the recovery of GSH upon NO exposure. Thus, increased glucose consumption and GSH-mediated redox capability may explain the natural resistance of L. braziliensis against NO.

Stress conditions promote Leishmania hybridization in vitro marked by expression of the ancestral gamete fusogen HAP2 as revealed by single-cell RNA-seq

Leishmania are protozoan parasites transmitted by the bite of sand fly vectors producing a wide spectrum of diseases in their mammalian hosts. These diverse clinical outcomes are directly associated with parasite strain and species diversity. Although Leishmania reproduction is mainly clonal, a cryptic sexual cycle capable of producing hybrid genotypes has been inferred from population genetic studies and directly demonstrated by laboratory crosses. Experimentally, mating competence has been largely confined to promastigotes developing in the sand fly midgut. The ability to hybridize culture promastigotes in vitro has been limited so far to low-efficiency crosses between two Leishmania tropica strains, L747 and MA37, that mate with high efficiency in flies. Here, we show that exposure of promastigote cultures to DNA damage stress produces a remarkably enhanced efficiency of in vitro hybridization of the L. tropica strains and extends to other species, including Leishmania donovani, Leishmania infantum, and Leishmania braziliensis, a capacity to generate intra- and interspecific hybrids. Whole-genome sequencing and total DNA content analyses indicate that the hybrids are in each case full genome, mostly tetraploid hybrids. Single-cell RNA sequencing of the L747 and MA37 parental lines highlights the transcriptome heterogeneity of culture promastigotes and reveals discrete clusters that emerge post-irradiation in which genes potentially involved in genetic exchange are expressed, including the ancestral gamete fusogen HAP2. By generating reporter constructs for HAP2, we could select for promastigotes that could either hybridize or not in vitro. Overall, this work reveals that there are specific populations involved in Leishmania hybridization associated with a discernible transcriptomic signature, and that stress facilitated in vitro hybridization can be a transformative approach to generate large numbers of hybrid genotypes between diverse species and strains.

High genome plasticity and frequent genetic exchange in Leishmania tropica isolates from Afghanistan, Iran and Syria

Background

The kinetoplastid protozoan Leishmania tropica mainly causes cutaneous leishmaniasis in humans in the Middle East, and relapse or treatment failure after treatment are common in this area. L. tropica’s digenic life cycle includes distinct stages in the vector sandfly and the mammalian host. Sexual reproduction and genetic exchange appear to occur more frequently than in other Leishmania species. Understanding these processes is complicated by chromosome instability during cell division that yields aneuploidy, recombination and heterozygosity. This combination of rare recombination and aneuploid permits may reveal signs of hypothetical parasexual mating, where diploid cells fuse to form a transient tetraploid that undergoes chromosomal recombination and gradual chromosomal loss.

Methodology/principal findings

The genome-wide SNP diversity from 22 L. tropica isolates showed chromosome-specific runs of patchy heterozygosity and extensive chromosome copy number variation. All these isolates were collected during 2007–2017 in Sweden from patients infected in the Middle East and included isolates from a patient possessing two genetically distinct leishmaniasis infections three years apart with no evidence of re-infection. We found differing ancestries on the same chromosome (chr36) across multiple samples: matching the reference genome with few derived alleles, followed by blocks of heterozygous SNPs, and then by clusters of homozygous SNPs with specific recombination breakpoints at an inferred origin of replication. Other chromosomes had similar marked changes in heterozygosity at strand-switch regions separating polycistronic transcriptional units.

Conclusion/significance

These large-scale intra- and inter-chromosomal changes in diversity driven by recombination and aneuploidy suggest multiple mechanisms of cell reproduction and diversification in L. tropica, including mitotic, meiotic and parasexual processes. It underpins the need for more genomic surveillance of Leishmania, to detect emerging hybrids that could spread more widely and to better understand the association between genetic variation and treatment outcome. Furthering our understanding of Leishmania genome evolution and ancestry will aid better diagnostics and treatment for cutaneous leishmaniasis caused by L.tropica in the Middle East.

Cutaneous leishmaniasis is mainly caused by Leishmania tropica in the Middle East, where it is known for treatment failure and a need for prolonged and/or multiple treatments. Several factors affect the clinical presentation and treatment outcome, such as host genetic variability and specific immune response, as well as environmental factors and the vector species. Little is known about the parasite genome and its influence on treatment response. By analysing the genome of 22 isolates of L. tropica, we have revealed extensive genomic variation and a complex population structure with evidence of genetic exchange within and among the isolates, indicating a possible presence of sexual or parasexual mechanisms. Understanding the Leishmania genome better may improve future treatment and better understanding of treatment failure and relapse.

High throughput single-cell genome sequencing gives insights into the generation and evolution of mosaic aneuploidy in Leishmania donovani

Leishmania, a unicellular eukaryotic parasite, is a unique model for aneuploidy and cellular heterogeneity, along with their potential role in adaptation to environmental stresses. Somy variation within clonal populations was previously explored in a small subset of chromosomes using fluorescence hybridization methods. This phenomenon, termed mosaic aneuploidy (MA), might have important evolutionary and functional implications but remains under-explored due to technological limitations. Here, we applied and validated a high throughput single-cell genome sequencing method to study for the first time the extent and dynamics of whole karyotype heterogeneity in two clonal populations of Leishmania promastigotes representing different stages of MA evolution in vitro. We found that drastic changes in karyotypes quickly emerge in a population stemming from an almost euploid founder cell. This possibly involves polyploidization/hybridization at an early stage of population expansion, followed by assorted ploidy reduction. During further stages of expansion, MA increases by moderate and gradual karyotypic alterations, affecting a defined subset of chromosomes. Our data provide the first complete characterization of MA in Leishmania and pave the way for further functional studies.

Leishmania and the Model of Predominant Clonal Evolution

As it is the case for other pathogenic microorganisms, the respective impact of clonality and genetic exchange on Leishmania natural populations has been the object of lively debates since the early 1980s. The predominant clonal evolution (PCE) model states that genetic exchange in these parasites’ natural populations may have a high relevance on an evolutionary scale, but is not sufficient to erase a persistent phylogenetic signal and the existence of bifurcating trees. Recent data based on high-resolution markers and genomic polymorphisms fully confirm the PCE model down to a microevolutionary level.

Unpicking the Roles of DNA Damage Protein Kinases in Trypanosomatids

To preserve genome integrity when faced with DNA lesions, cells activate and coordinate a multitude of DNA repair pathways to ensure timely error correction or tolerance, collectively called the DNA damage response (DDR). These interconnecting damage response pathways are molecular signal relays, with protein kinases (PKs) at the pinnacle. Focused efforts in model eukaryotes have revealed intricate aspects of DNA repair PK function, including how they direct DDR pathways and how repair reactions connect to wider cellular processes, including DNA replication and transcription. The Kinetoplastidae, including many parasites like Trypanosoma spp. and Leishmania spp. (causative agents of debilitating, neglected tropical infections), exhibit peculiarities in several core biological processes, including the predominance of multigenic transcription and the streamlining or repurposing of DNA repair pathways, such as the loss of non-homologous end joining and novel operation of nucleotide excision repair (NER). Very recent studies have implicated ATR and ATM kinases in the DDR of kinetoplastid parasites, whereas DNA-dependent protein kinase (DNA-PKcs) displays uncertain conservation, questioning what functions it fulfills. The wide range of genetic manipulation approaches in these organisms presents an opportunity to investigate DNA repair kinase roles in kinetoplastids and to ask if further kinases are involved. Furthermore, the availability of kinase inhibitory compounds, targeting numerous eukaryotic PKs, could allow us to test the suitability of DNA repair PKs as novel chemotherapeutic targets. Here, we will review recent advances in the study of trypanosomatid DNA repair kinases.

Comparative analysis of the transcriptional responses of five Leishmania species to trivalent antimony

BACKGROUND

Leishmaniasis is a neglected tropical disease caused by several species of Leishmania. The resistance phenotype of these parasites depends on the characteristics of each species, which contributes to increased therapeutic failures. Understanding the mechanism used by the parasite to survive under treatment pressure in order to identify potential common and specific therapeutic targets is essential for the control of leishmaniasis. The aim of this study was to investigate the expression profiles and potential shared and specific resistance markers of the main Leishmania species of medical importance [subgenus L. (Leishmania): L. donovani, L. infantum and L. amazonensis; subgenus L. (Viannia): L. panamensis and L. braziliensis)] resistant and sensitive to trivalent stibogluconate (SbIII).

METHODS

We conducted comparative analysis of the transcriptomic profiles (only coding sequences) of lines with experimentally induced resistance to SbIII from biological replicates of five Leishmania species available in the databases of four articles based on ortholog attribution. Simultaneously, we carried out functional analysis of ontology and reconstruction of metabolic pathways of the resulting differentially expressed genes (DEGs).

RESULTS

Resistant lines for each species had differential responses in metabolic processes, compound binding, and membrane components concerning their sensitive counterpart. One hundred and thirty-nine metabolic pathways were found, with the three main pathways comprising cysteine and methionine metabolism, glycolysis, and the ribosome. Differentially expressed orthologous genes assigned to species-specific responses predominated, with 899 self-genes. No differentially expressed genes were found in common among the five species. Two common upregulated orthologous genes were found among four species (L. donovani, L. braziliensis, L. amazonensis, and L. panamensis) related to an RNA-binding protein and the NAD(P)H cytochrome-B5-oxidoreductase complex, associated with transcriptional control and de novo synthesis of linoleic acid, critical mechanisms in resistance to antimonials.

CONCLUSION

Herein, we identified potential species-specific genes related to resistance to SbIII. Therefore, we suggest that future studies consider a treatment scheme that is species-specific. Despite the limitations of our study, this is the first approach toward unraveling the pan-genus genetic mechanisms of resistance in leishmaniasis.

Molecular tools confirm natural Leishmania (Viannia) guyanensis/L. (V.) shawi hybrids causing cutaneous leishmaniasis in the Amazon region of Brazil

Seven isolates from patients with American cutaneous leishmaniasis in the Amazon region of Brazil were phenotypically suggestive of Leishmania (Viannia) guyanensis/L. (V.) shawi hybrids. In this work, two molecular targets were employed to check the hybrid identity of the putative hybrids. Heat shock protein 70 (hsp70) gene sequences were analyzed by three different polymerase chain reaction (PCR) approaches, and two different patterns of inherited hsp70 alleles were found. Three isolates presented heterozygous L. (V.) guyanensis/L. (V.) shawi patterns, and four presented homozygous hsp70 patterns involving only L. (V.) shawi alleles. The amplicon sequences confirmed the RFLP patterns. The high-resolution melting method detected variant heterozygous and homozygous profiles. Single-nucleotide polymorphism genotyping/cleaved amplified polymorphic site analysis suggested a higher contribution from L. (V.) guyanensis in hsp70 heterozygous hybrids. Additionally, PCR-RFLP analysis targeting the enzyme mannose phosphate isomerase (mpi) gene indicated heterozygous and homozygous cleavage patterns for L. (V.) shawi and L. (V.) guyanensis, corroborating the hsp70 findings. In this communication, we present molecular findings based on partial informative regions of the coding sequences of hsp70 and mpi as markers confirming that some of the parasite strains from the Brazilian Amazon region are indeed hybrids between L. (V.) guyanensis and L. (V.) shawi.

Genomics and functional genomics in Leishmania and Trypanosoma cruzi: statuses, challenges and perspectives

The availability of Trypanosomatid genomic data in public databases has opened myriad experimental possibilities that have contributed to a more comprehensive understanding of the biology of these parasites and their interactions with hosts. In this review, after brief remarks on the history of the Trypanosoma cruzi and Leishmania genome initiatives, we present an overview of the relevant contributions of genomics, transcriptomics and functional genomics, discussing the primary obstacles, challenges, relevant achievements and future perspectives of these technologies.

Leishmania Sexual Reproductive Strategies as Resolved through Computational Methods Designed for Aneuploid Genomes

A cryptic sexual reproductive cycle in Leishmania has been inferred through population genetic studies revealing the presence of hybrid genotypes in natural isolates, with attempts made to decipher sexual strategies by studying complex chromosomal inheritance patterns. A more informative approach is to study the products of controlled, laboratory-based experiments where known strains or species are crossed in the sand fly vector to generate hybrid progeny. These hybrids can be subsequently studied through high resolution sequencing technologies and software suites such as PAINT that disclose inheritance patterns including ploidies, parental chromosome contributions and recombinations, all of which can inform the sexual strategy. In this work, we discuss the computational methods in PAINT that can be used to interpret the sexual strategies adopted specifically by aneuploid organisms and summarize how PAINT has been applied to the analysis of experimental hybrids to reveal meiosis-like sexual recombination in Leishmania.

Leishmania: Responding to environmental signals and challenges without regulated transcription

Here we describe the non-canonical control of gene expression in Leishmania, a single-cell parasite that is responsible for one of the major neglected tropical diseases. We discuss the lack of regulated RNA synthesis, the post-transcriptional gene regulation including RNA stability and regulated translation. We also show that genetic adaptations such as mosaic aneuploidy, gene copy number variations and DNA sequence polymorphisms are important means for overcoming drug challenge and environmental diversity. These mechanisms are discussed in the context of the unique flow of genetic information found in Leishmania and related protists.

Nuclear and mitochondrial genome sequencing of North-African Leishmania infantum isolates from cured and relapsed visceral leishmaniasis patients reveals variations correlating with geography and phenotype

Although several studies have investigated genetic diversity of Leishmania infantum in North Africa, genome-wide analyses are lacking. Here, we conducted comparative analyses of nuclear and mitochondrial genomes of seven L. infantum isolates from Tunisia with the aim to gain insight into factors that drive genomic and phenotypic adaptation. Isolates were from cured (n=4) and recurrent (n=3) visceral leishmaniasis (VL) cases, originating from northern (n=2) and central (n=5) Tunisia, where respectively stable and emerging VL foci are observed. All isolates from relapsed patients were from Kairouan governorate (Centre); one showing resistance to the anti-leishmanial drug Meglumine antimoniate. Nuclear genome diversity of the isolates was analysed by comparison to the L. infantum JPCM5 reference genome. Kinetoplast maxi and minicircle sequences (1 and 59, respectively) were extracted from unmapped reads and identified by blast analysis against public data sets. The genome variation analysis grouped together isolates from the same geographical origins. Strains from the North were very different from the reference showing more than 34 587 specific single nucleotide variants, with one isolate representing a full genetic hybrid as judged by variant frequency. Composition of minicircle classes within isolates corroborated this geographical population structure. Read depth analysis revealed several significant gene copy number variations correlating with either geographical origin (amastin and Hsp33 genes) or relapse (CLN3 gene). However, no specific gene copy number variation was found in the drug-resistant isolate. In contrast, resistance was associated with a specific minicircle pattern suggesting Leishmania mitochondrial DNA as a potential novel source for biomarker discovery.

Expanded genome-wide comparisons give novel insights into population structure and genetic heterogeneity of Leishmania tropica complex

Leishmania tropica is one of the main causative agents of cutaneous leishmaniasis (CL). Population structures of L. tropica appear to be genetically highly diverse. However, the relationship between L. tropica strains genomic diversity, protein coding gene evolution and biogeography are still poorly understood. In this study, we sequenced the genomes of three new clinical L. tropica isolates, two derived from a recent outbreak of CL in camps hosting Syrian refugees in Lebanon and one historical isolate from Azerbaijan to further refine comparative genome analyses. In silico multilocus microsatellite typing (MLMT) was performed to integrate the current diversity of genome sequence data in the wider available MLMT genetic population framework. Single nucleotide polymorphism (SNPs), gene copy number variations (CNVs) and chromosome ploidy were investigated across the available 18 L. tropica genomes with a main focus on protein coding genes. MLMT divided the strains in three populations that broadly correlated with their geographical distribution but not populations defined by SNPs. Unique SNPs profiles divided the 18 strains into five populations based on principal component analysis. Gene ontology enrichment analysis of the protein coding genes with population specific SNPs profiles revealed various biological processes, including iron acquisition, sterols synthesis and drug resistance. This study further highlights the complex links between L. tropica important genomic heterogeneity and the parasite broad geographic distribution. Unique sequence features in protein coding genes identified in distinct populations reveal potential novel markers that could be exploited for the development of more accurate typing schemes to further improve our knowledge of the evolution and epidemiology of the parasite as well as highlighting protein variants of potential functional importance underlying L. tropica specific biology.

Genome duplication in Leishmania major relies on persistent subtelomeric DNA replication

DNA replication is needed to duplicate a cell’s genome in S phase and segregate it during cell division. Previous work in Leishmania detected DNA replication initiation at just a single region in each chromosome, an organisation predicted to be insufficient for complete genome duplication within S phase. Here, we show that acetylated histone H3 (AcH3), base J and a kinetochore factor co-localise in each chromosome at only a single locus, which corresponds with previously mapped DNA replication initiation regions and is demarcated by localised G/T skew and G4 patterns. In addition, we describe previously undetected subtelomeric DNA replication in G2/M and G1-phase-enriched cells. Finally, we show that subtelomeric DNA replication, unlike chromosome-internal DNA replication, is sensitive to hydroxyurea and dependent on 9-1-1 activity. These findings indicate that Leishmania’s genome duplication programme employs subtelomeric DNA replication initiation, possibly extending beyond S phase, to support predominantly chromosome-internal DNA replication initiation within S phase.

In-depth quantitative proteomics uncovers specie-specific metabolic programs in Leishmania (Viannia) species

Leishmania species are responsible for a broad spectrum of diseases, denominated Leishmaniasis, affecting over 12 million people worldwide. During the last decade, there have been impressive efforts for sequencing the genome of most of the pathogenic Leishmania spp. as well as hundreds of strains, but large-scale proteomics analyses did not follow these achievements and the Leishmania proteome remained mostly uncharacterized. Here, we report a comprehensive comparative study of the proteomes of strains representing L. braziliensis, L. panamensis and L. guyanensis species. Proteins extracted by SDS-mediated lysis were processed following the multi-enzyme digestion-filter aided sample preparation (FASP) procedure and analysed by high accuracy mass spectrometry. "Total Protein Approach" and "Proteomic Ruler" were applied for absolute quantification of proteins. Principal component analysis demonstrated very high reproducibility among biological replicates and a very clear differentiation of the three species. Our dataset comprises near 7000 proteins, representing the most complete Leishmania proteome yet known, and provides a comprehensive quantitative picture of the proteomes of the three species in terms of protein concentration and copy numbers. Analysis of the abundance of proteins from the major energy metabolic processes allow us to highlight remarkably differences among the species and suggest that these parasites depend on distinct energy substrates to obtain ATP. Whereas L. braziliensis relies the more on glycolysis, L. panamensis and L. guyanensis seem to depend mainly on mitochondrial respiration. These results were confirmed by biochemical assays showing opposite profiles for glucose uptake and O2 consumption in these species. In addition, we provide quantitative data about different membrane proteins, transporters, and lipids, all of which contribute for significant species-specific differences and provide rich substrate for explore new molecules for diagnosing purposes. Data are available via ProteomeXchange with identifier PXD017696.

Application of next generation sequencing (NGS) for descriptive analysis of 30 genomes of Leishmania infantum isolates in Middle-North Brazil

Visceral leishmaniasis (VL) is a life-threatening disease caused by the protozoa Leishmania donovani and L. infantum. Likely, L. infantum was introduced in the New World by the Iberic colonizers. Due to recent introduction, the genetic diversity is low. Access to genomic information through the sequencing of Leishmania isolates allows the characterization of populations through the identification and analysis of variations. Population structure information may reveal important data on disease dynamics. Aiming to describe the genetic diversity of L. infantum from the Middle-North, Brazil, next generation sequencing of 30 Leishmania isolates obtained in the city of Teresina, from where the disease dispersed, was performed. The variations were categorized accordingly to the genome region and impact and provided the basis for chromosomal ploidy and population structure analysis. The results showed low diversity between the isolates and the Iberic reference genome JPCM5. Most variations were seen in non-coding regions, with modifying impact. The ploidy number analysis showed aneuploid profile. The population structure analysis revealed the presence of two L. infantum populations identified in Teresina. Further population genetics studies with a larger number of isolates should be performed in order to identify the genetic background associated with virulence and parasite ecology.

Unravelling the proteomic signature of extracellular vesicles released by drug-resistant Leishmania infantum parasites

Leishmaniasis constitutes the 9^th largest disease burden among all infectious diseases. Control of this disease is based on a short list of chemotherapeutic agents headed by pentavalent antimonials, followed by miltefosine and amphotericin B; drugs that are far from ideal due to host toxicity, elevated cost, limited access, and high rates of drug resistance. Knowing that the composition of extracellular vesicles (EVs) can vary according to the state of their parental cell, we hypothesized that EVs released by drug-resistant Leishmania infantum parasites could contain unique and differently enriched proteins depending on the drug-resistance mechanisms involved in the survival of their parental cell line. To assess this possibility, we studied EV production, size, morphology, and protein content of three well-characterized drug-resistant L. infantum cell lines and a wild-type strain. Our results are the first to demonstrate that drug-resistance mechanisms can induce changes in the morphology, size, and distribution of L. infantum EVs. In addition, we identified L. infantum’s core EV proteome. This proteome is highly conserved among strains, with the exception of a handful of proteins that are enriched differently depending on the drug responsible for induction of antimicrobial resistance. Furthermore, we obtained the first snapshot of proteins enriched in EVs released by antimony-, miltefosine- and amphotericin-resistant parasites. These include several virulence factors, transcription factors, as well as proteins encoded by drug-resistance genes. This detailed study of L. infantum EVs sheds new light on the potential roles of EVs in Leishmania biology, particularly with respect to the parasite’s survival in stressful conditions. This work outlines a crucial first step towards the discovery of EV-based profiles capable of predicting response to antileishmanial agents.

Visceral leishmaniasis is a life-threatening disease caused by Leishmania infantum parasites, which are transmitted by sand flies. In the absence of vaccines, current control of this disease is based on chemotherapy, which is comprised of a very limited arsenal threatened by the emergence and spread of drug-resistant strains. In the shadow of growing concern and treatment failure due to resistance, the characterization of extracellular vesicles (EVs) released by drug-resistant L. infantum parasites could shed some light on the complex nature of drug resistance in Leishmania and increase our understanding of the biology of the parasite. EVs are vesicles secreted by all eukaryotic cells whose contents (proteins, DNA/RNAs, lipids) vary as a function of their cellular origin. Our results demonstrate for the first time that EVs released by drug-resistant parasites are enriched in unique protein markers that reflect the drug-resistance mechanisms involved in the survival of parental cells. These unique proteins included several virulence and transcription factors, as well as drug-resistance genes; this offers a potential benefit for drug-resistant parasites in terms of parasite-to-parasite communication and host-parasite interactions. Collectively, our initial results could serve as a jumping-off point for the future development of novel EV-based diagnostic tools for the detection and appraisal of antimicrobial-resistant Leishmania populations.

MRPA-independent mechanisms of antimony resistance in Leishmania infantum

Control of both human and canine leishmaniasis is based on a very short list of chemotherapeutic agents, headed by antimonial derivatives (Sb). The utility of these molecules is severely threatened by high rates of drug resistance. The ABC transporter MRPA is one of the few key Sb resistance proteins described to date, whose role in detoxification has been thoroughly studied in Leishmania parasites. Nonetheless, its rapid amplification during drug selection complicates the discovery of other mechanisms potentially involved in Sb resistance. In this study, stepwise drug-resistance selection and next-generation sequencing were combined in the search for novel Sb-resistance mechanisms deployed by parasites when MRPA is abolished by targeted gene disruption. The gene mrpA is not essential in L. infantum, and its disruption leads to an Sb hypersensitive phenotype in both promastigotes and amastigotes. Five independent mrpA^-/- mutants were selected for antimony resistance. These mutants displayed major changes in their ploidy, as well as extrachromosomal linear amplifications of the subtelomeric region of chromosome 23, which includes the genes coding for ABCC1 and ABCC2. Overexpression of ABCC2, but not of ABCC1, resulted in increased Sb tolerance in the mrpA^-/- mutant. SNP analyses revealed three different heterozygous mutations in the gene coding for a serine acetyltransferase (SAT) involved in de novo cysteine synthesis in Leishmania. Overexpression of sat^Q390K, sat^G321R and sat^G325R variants led to a 2-3.2 -fold increase in Sb resistance in mrpA^-/- parasites. Only sat^G321R and sat^G325R induced increased Sb resistance in wild-type parasites. These results reinforce and expand knowledge on the complex nature of Sb resistance in Leishmania parasites.

Genomic analysis of natural intra-specific hybrids among Ethiopian isolates of Leishmania donovani

Parasites of the genus Leishmania (Kinetoplastida: Trypanosomatidae) cause widespread and devastating human diseases. Visceral leishmaniasis due to Leishmania donovani is endemic in Ethiopia where it has also been responsible for major epidemics. The presence of hybrid genotypes has been widely reported in surveys of natural populations, genetic variation reported in a number of Leishmania species, and the extant capacity for genetic exchange demonstrated in laboratory experiments. However, patterns of recombination and the evolutionary history of admixture that produced these hybrid populations remain unclear. Here, we use whole-genome sequence data to investigate Ethiopian L. donovani isolates previously characterized as hybrids by microsatellite and multi-locus sequencing. To date there is only one previous study on a natural population of Leishmania hybrids based on whole-genome sequences. We propose that these hybrids originate from recombination between two different lineages of Ethiopian L. donovani occurring in the same region. Patterns of inheritance are more complex than previously reported with multiple, apparently independent, origins from similar parents that include backcrossing with parental types. Analysis indicates that hybrids are representative of at least three different histories. Furthermore, isolates were highly polysomic at the level of chromosomes with differences between parasites recovered from a recrudescent infection from a previously treated individual. The results demonstrate that recombination is a significant feature of natural populations and contributes to the growing body of data that shows how recombination, and gene flow, shape natural populations of Leishmania.

Global genome diversity of the Leishmania donovani complex

Protozoan parasites of the Leishmania donovani complex - L. donovani and L. infantum - cause the fatal disease visceral leishmaniasis. We present the first comprehensive genome-wide global study, with 151 cultured field isolates representing most of the geographical distribution. L. donovani isolates separated into five groups that largely coincide with geographical origin but vary greatly in diversity. In contrast, the majority of L. infantum samples fell into one globally-distributed group with little diversity. This picture is complicated by several hybrid lineages. Identified genetic groups vary in heterozygosity and levels of linkage, suggesting different recombination histories. We characterise chromosome-specific patterns of aneuploidy and identified extensive structural variation, including known and suspected drug resistance loci. This study reveals greater genetic diversity than suggested by geographically-focused studies, provides a resource of genomic variation for future work and sets the scene for a new understanding of the evolution and genetics of the Leishmania donovani complex.

Trans-Atlantic Spill Over: Deconstructing the Ecological Adaptation of Leishmania infantum in the Americas

Pathogen fitness landscapes change when transmission cycles establish in non-native environments or spill over into new vectors and hosts. The introduction of Leishmania infantum in the Americas into the Neotropics during European colonization represents a unique case study to investigate the mechanisms of ecological adaptation of this important parasite. Defining the evolutionary trajectories that drive L. infantum fitness in this new environment are of great public health importance as they will allow unique insight into pathways of host/pathogen co-evolution and their consequences for region-specific changes in disease manifestation. This review summarizes current knowledge on L. infantum genetic and phenotypic diversity in the Americas and its possible role in the unique epidemiology of visceral leishmaniasis (VL) in the New World. We highlight the importance of appreciating adaptive molecular mechanisms in L. infantum to understand the parasites' successful establishment on the continent.

Meiotic sex in Chagas disease parasite Trypanosoma cruzi

Genetic exchange enables parasites to rapidly transform disease phenotypes and exploit new host populations. Trypanosoma cruzi, the parasitic agent of Chagas disease and a public health concern throughout Latin America, has for decades been presumed to exchange genetic material rarely and without classic meiotic sex. We present compelling evidence from 45 genomes sequenced from southern Ecuador that T. cruzi in fact maintains truly sexual, panmictic groups that can occur alongside others that remain highly clonal after past hybridization events. These groups with divergent reproductive strategies appear genetically isolated despite possible co-occurrence in vectors and hosts. We propose biological explanations for the fine-scale disconnectivity we observe and discuss the epidemiological consequences of flexible reproductive modes. Our study reinvigorates the hunt for the site of genetic exchange in the T. cruzi life cycle, provides tools to define the genetic determinants of parasite virulence, and reforms longstanding theory on clonality in trypanosomatid parasites.

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

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

Diagnosing Cutaneous leishmaniasis using Fluorescence in Situ Hybridization: the Sri Lankan Perspective

Cutaneous leishmaniasis (CL) caused by Leishmania donovani MON-37 is becoming a major public health problem in Sri Lanka, with 100 new cases per month being reported in endemic regions. Diagnosis of CL is challenging for several reasons. Due to relative specificity and rapidity we propose Fluorescence in Situ Hybridization as a diagnostic tool for CL. The objective was to evaluate the potential of Fluorescence in Situ Hybridization as a diagnostic method for Cutaneous leishmaniasis in Sri Lanka. Literature on current laboratory tests used to diagnose Cutaneous leishmaniasis in Sri Lanka and globally was reviewed. Sri Lankan data were reviewed systematically following the PRISMA guidelines. A narrative of the results is presented. There is currently no gold standard diagnostic method for Cutaneous leishmaniasis. Fluorescence in Situ Hybridization has been previously applied to detect dermal pathologies including those involving infectious agents, and its use to detect the Leishmania parasite in human cutaneous lesions reported in small number of studies, generally with limited numbers of subjects. Advantages of FISH has been specificity, cost and ease-of-use compared to the alternatives. Based on the available literature and our current work, FISH has potential for diagnosing CL and should now be evaluated in larger cohorts in endemic regions. FISH for CL diagnosis could find application in countries such as Sri Lanka, where laboratory facilities may be limited in rural areas where the disease burden is highest.

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.

Leishmania braziliensis is the main etiological agent of cutaneous leishmaniasis (CL) in Latin America. In the Pernambuco state, a highly endemic area for CL in Northeastern Brazil, ten zymodemes have been previously described, identified through the detection of variations in the mobility of selected enzymes on gel electrophoresis (MLEE). Here we used next-generation genome sequencing to reevaluate the genetic diversity and relatedness between isolates representing these zymodemes. The genetic analysis clustered these isolates into three distinct phylogenetic groups, with those circulating in forested environments more related to each other (Z72 and Z75) and separated from the other eight isolates from more urbanized environments, but nevertheless forming two further groups. All ten isolates were definitely identified as L. braziliensis, including one from zymodeme Z26, previously identified as L. shawi through MLEE. Chromosome copy number variation was observed in all isolates studied, but restricted mainly to chromosome 31, with a single isolate (from zymodeme Z27) showing more than two copies of other chromosomes. A low heterozygosity observed only for both Z72/Z75 isolates also confirm differences associated with them. Our findings confirm that MLEE, despite still being considered the gold standard for characterization of Leishmania spp., does not reflect relevant information on intra-specific variations. It also highlights the substantial diversity observed between L. braziliensis isolates.

Whole genome sequencing of experimental hybrids supports meiosis-like sexual recombination in Leishmania

Hybrid genotypes have been repeatedly described among natural isolates of Leishmania, and the recovery of experimental hybrids from sand flies co-infected with different strains or species of Leishmania has formally demonstrated that members of the genus possess the machinery for genetic exchange. As neither gamete stages nor cell fusion events have been directly observed during parasite development in the vector, we have relied on a classical genetic analysis to determine if Leishmania has a true sexual cycle. Here, we used whole genome sequencing to follow the chromosomal inheritance patterns of experimental hybrids generated within and between different strains of L. major and L. infantum. We also generated and sequenced the first experimental hybrids in L. tropica. We found that in each case the parental somy and allele contributions matched the inheritance patterns expected under meiosis 97–99% of the time. The hybrids were equivalent to F1 progeny, heterozygous throughout most of the genome for the markers that were homozygous and different between the parents. Rare, non-Mendelian patterns of chromosomal inheritance were observed, including a gain or loss of somy, and loss of heterozygosity, that likely arose during meiosis or during mitotic divisions of the progeny clones in the fly or culture. While the interspecies hybrids appeared to be sterile, the intraspecies hybrids were able to produce backcross and outcross progeny. Analysis of 5 backcross and outcross progeny clones generated from an L. major F1 hybrid, as well as 17 progeny clones generated from backcrosses involving a natural hybrid of L. tropica, revealed genome wide patterns of recombination, demonstrating that classical crossing over occurs at meiosis, and allowed us to construct the first physical and genetic maps in Leishmania. Altogether, the findings provide strong evidence for meiosis-like sexual recombination in Leishmania, presenting clear opportunities for forward genetic analysis and positional cloning of important genes.

Leishmania promastigotes are able to undergo genetic exchange during their growth and development in the sand fly vector, however, it is still not known if they have a true sexual cycle involving meiosis. Here, we used whole genome sequencing to follow the chromosomal inheritance patterns of 44 experimental hybrids generated between different strains of L. major, L. infantum, and L. tropica. In almost every case the number of chromosomes and the allele contributions from each parent matched the inheritance patterns expected under meiosis. Rare instances of hybrid chromosomes that did not fit Mendelian expectations were observed, including gain or loss of somy, and loss of heterozygosity. Strong evidence for a meiotic-like process was also obtained from the genome wide patterns of recombination observed in the offspring generated from backcrosses involving an experimental or natural hybrid, consistent with crossing over occurring between homologous chromosomes during meiosis. The frequency and position of the recombination breakpoints observed on each chromosome allowed us to construct the first physical and genetic maps in Leishmania. The results demonstrate that forward genetic approaches are possible for positional cloning of important genes.

TFIIIB Subunit Bdp1 Participates in RNA Polymerase III Transcription in the Protozoan Parasite Leishmania major

Leishmania major, a protozoan parasite that diverged early from the main eukaryotic lineage, exhibits unusual mechanisms of gene expression. Little is known in this organism about the transcription factors involved in the synthesis of tRNA, 5S rRNA, and snRNAs, transcribed by RNA Polymerase III (Pol III). Here we identify and characterize the TFIIIB subunit Bdp1 in L. major (LmBdp1). Bdp1 plays key roles in Pol III transcription initiation in other organisms, as it participates in Pol III recruitment and promoter opening. In silico analysis showed that LmBdp1 contains the typical extended SANT domain as well as other Bdp1 conserved regions. Nevertheless, LmBdp1 also displays distinctive features, including the presence of only one aromatic residue in the N-linker region. We were not able to produce null mutants of LmBdp1 by homologous recombination, as the obtained double replacement cell line contained an extra copy of LmBdp1, indicating that LmBdp1 is essential for the viability of L. major promastigotes. Notably, the mutant cell line showed reduced levels of the LmBdp1 protein, and its growth was significantly decreased in relation to wild-type cells. Nuclear run-on assays demonstrated that Pol III transcription was affected in the mutant cell line, and ChIP experiments showed that LmBdp1 binds to 5S rRNA, tRNA, and snRNA genes. Thus, our results indicate that LmBdp1 is an essential protein required for Pol III transcription in L. major.

Trypanosoma cruzi Genome Assemblies: Challenges and Milestones of Assembling a Highly Repetitive and Complex Genome

Trypanosoma cruzi present one of the most complex parasite genomes sequenced to date. Among its features are 600-kb-long repetitive multigene families' clusters, hybrid strains, and aneuploidies, which hampered genome assembly completeness and contiguity. Several approaches, such as Sanger sequencing in 2005, next-generation sequencing in 2011 and third-generation sequencing in 2018, were used to improve draft assemblies of different strains of this parasite. Hence, the study of T. cruzi genome assemblies' history is an excellent way to describe the evolution of genome sequencing methodologies and compare their efficiency and limitations to assembly complex genomes. In this book chapter, we summarize the principal findings and methodologies of T. cruzi genome assembly projects to date, highlighting the improvements and limitations of each approach.

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.