Chromosome and gene copy number variation allow major structural change between species and strains of Leishmania

LeishGEM: genome-wide deletion mutant fitness and protein localisations in Leishmania

LeishGEM is a genome-wide functional annotation community resource for Leishmania mexicana, where deletion mutant growth in vitro and in vivo is measured and protein localisation is determined by endogenous tagging and LOPIT-DC (localisation of organelle proteins by isotope tagging with differential centrifugation) spatial proteomics. Data are being made available pre-publication via http://leishgem.org which allows data-driven identification of the mechanisms for Leishmania parasitism.

Evaluation of the Leishmania Inositol Phosphorylceramide Synthase as a Drug Target Using a Chemical and Genetic Approach

The lack of effective vaccines and the development of resistance to the current treatments highlight the urgent need for new anti-leishmanials. Sphingolipid metabolism has been proposed as a promising source of Leishmania-specific targets as these lipids are key structural components of the eukaryotic plasma membrane and are involved in distinct cellular events. Inositol phosphorylceramide (IPC) is the primary sphingolipid in the Leishmania species and is the product of a reaction mediated by IPC synthase (IPCS). The antihistamine clemastine fumarate has been identified as an inhibitor of IPCS in L. major and a potent anti-leishmanial in vivo. Here we sought to further examine the target of this compound in the more tractable species L. mexicana, using an approach combining genomic, proteomic, metabolomic and lipidomic technologies, with molecular and biochemical studies. While the data demonstrated that the response to clemastine fumarate was largely conserved, unexpected disturbances beyond sphingolipid metabolism were identified. Furthermore, while deletion of the gene encoding LmxIPCS had little impact in vitro, it did influence clemastine fumarate efficacy and, importantly, in vivo pathogenicity. Together, these data demonstrate that clemastine does inhibit LmxIPCS and cause associated metabolic disturbances, but its primary target may lie elsewhere.

Genomic and epidemiological evidence for the emergence of a L. infantum/L. donovani hybrid with unusual epidemiology in northern Italy

Leishmania (L.) infantum is one of the main causative agents of animal and human leishmaniasis across many endemic areas in South America, Europe, North Africa, and Asia. Despite its clinical significance, little is known about the genetic diversity of L. infantum circulating in a given endemic area. Here, we investigate this important open question by applying a comparative genomics approach to seven L. infantum isolates from different hosts and Italian regions, including the northern part of the country (Emilia-Romagna, RER), Sicily, and Sardinia, as an initial attempt to explore the breadth of parasite genetic heterogeneity in Italy. Additionally, microsatellite analysis was carried out to compare the isolates from RER with other 70 L. infantum strains from the same region as well as 65 strains belonging to the L. donovani complex from other countries. We revealed important karyotypic instability and identified strain-specific changes in gene dosage, which affected important virulence factors such as amastins and surface antigen-like proteins. Single nucleotide polymorphism-based clustering analysis of these genomes together with over 80 publicly available L. infantum and L. donovani genomes placed the Italian isolates into three geographically distinct clusters within the Mediterranean basin and uncovered three isolates clustering with putative L. infantum/L. donovani hybrids isolated in Cyprus. As judged by microsatellite profiling, these hybrid isolates are representative of a sub-population of parasites circulating in northern Italy that preferentially infect humans but not dogs. Our results place Italy at the crossroads of L. infantum infection in the Mediterranean and call attention to the public health risk represented by the introduction of non-European Leishmania species.IMPORTANCEThis study closes important knowledge gaps with respect to Leishmania (L.) infantum genetic heterogeneity in a given endemic country, as exemplified here for Italy, and reveals genetic hybridization as a main cause for re-emerging human leishmaniasis in northern Italy. The observed high diversity of Leishmania parasites on the Italian peninsula suggests different geographical origins, with genomic adaptation to various ecologies affecting both pathogenicity and transmission potential. This is documented by the discovery of a putative L. infantum/L. donovani hybrid strain, which has been shown to preferentially infect humans but not dogs. Our results provide important information to health authorities, which need to consider the public health risk represented by the introduction of new Leishmania species into EU countries due to population displacement or travel from countries where exotic/allochthonous parasite species are endemic.

Streptococcus pyogenes Cas9 ribonucleoprotein delivery for efficient, rapid and marker-free gene editing in Trypanosoma and Leishmania

Kinetoplastids are unicellular eukaryotic flagellated parasites found in a wide range of hosts within the animal and plant kingdoms. They are known to be responsible in humans for African sleeping sickness (Trypanosoma brucei), Chagas disease (Trypanosoma cruzi), and various forms of leishmaniasis (Leishmania spp.), as well as several animal diseases with important economic impact (African trypanosomes, including Trypanosoma congolense). Understanding the biology of these parasites necessarily implies the ability to manipulate their genomes. In this study, we demonstrate that transfection of a ribonucleoprotein complex, composed of recombinant Streptococcus pyogenes Cas9 (SpCas9) and an in vitro-synthesized guide RNA, results in rapid and efficient genetic modifications of trypanosomatids, in marker-free conditions. This approach was successfully developed to inactivate, delete, and mutate candidate genes in various stages of the life cycle of T. brucei and T. congolense, and Leishmania promastigotes. The functionality of SpCas9 in these parasites now provides, to the research community working on these parasites, a rapid and efficient method of genome editing, without requiring plasmid construction and selection by antibiotics but requires only cloning and PCR screening of the clones. Importantly, this approach is adaptable to any wild-type parasite.

Ancestral aneuploidy and stable chromosomal duplication resulting in differential genome structure and gene expression control in trypanosomatid parasites

Aneuploidy is widely observed in both unicellular and multicellular eukaryotes, usually associated with adaptation to stress conditions. Chromosomal duplication stability is a tradeoff between the fitness cost of having unbalanced gene copies and the potential fitness gained from increased dosage of specific advantageous genes. Trypanosomatids, a family of protozoans that include species that cause neglected tropical diseases, are a relevant group to study aneuploidies. Their life cycle has several stressors that could select for different patterns of chromosomal duplications and/or losses, and their nearly universal use of polycistronic transcription increases their reliance on gene expansion/contraction, as well as post-transcriptional control as mechanisms for gene expression regulation. By evaluating the data from 866 isolates covering seven trypanosomatid genera, we have revealed that aneuploidy tolerance is an ancestral characteristic of trypanosomatids but has a reduced occurrence in a specific monophyletic clade that has undergone large genomic reorganization and chromosomal fusions. We have also identified an ancient chromosomal duplication that was maintained across these parasite's speciation, named collectively as the trypanosomatid ancestral supernumerary chromosome (TASC). TASC has most genes in the same coding strand, is expressed as a disomic chromosome (even having four copies), and has increased potential for functional variation, but it purges highly deleterious mutations more efficiently than other chromosomes. The evidence of stringent control over gene expression in this chromosome suggests that these parasites have adapted to mitigate the fitness cost associated with this ancient chromosomal duplication.

A comparative genomics approach reveals a local genetic signature of Leishmania tropica in Morocco

In Morocco, cutaneous leishmaniasis (CL) caused by Leishmania (L.) tropica is an important health problem. Despite the high incidence of CL in the country, the genomic heterogeneity of these parasites is still incompletely understood. In this study, we sequenced the genomes of 14 Moroccan isolates of L. tropica collected from confirmed cases of CL to investigate their genomic heterogeneity. Comparative genomics analyses were conducted by applying the recently established Genome Instability Pipeline (GIP), which allowed us to conduct phylogenomic and principal components analyses (PCA), and to assess genomic variations at the levels of the karyotype, gene copy number, single nucleotide polymorphisms (SNPs) and small insertions/deletions (INDELs) variants. Read-depth analyses revealed a mostly disomic karyotype, with the exception of the stable tetrasomy of chromosome 31. In contrast, we identified important gene copy number variations across all isolates, which affect known virulence genes and thus were probably selected in the field. SNP-based cluster analysis of the 14 isolates revealed a core group of 12 strains that formed a tight cluster and shared 45.1 % (87 751) of SNPs, as well as two strains (M3015, Ltr_16) that clustered separately from each other and the core group, suggesting the circulation of genetically highly diverse strains in Morocco. Phylogenetic analysis, which compared our 14 L. tropica isolates against 40 published genomes of L. tropica from a diverse array of locations, confirmed the genetic difference of our Moroccan isolates from all other isolates examined. In conclusion, our results indicate potential regional variations in SNP profiles that may differentiate Moroccan L. tropica from other L. tropica strains circulating in endemic countries in the Middle East. Our report paves the way for future research with a larger number of strains that will allow correlation of diverse phenotypes (resistance to treatments, virulence) and origins (geography, host species, year of isolation) to defined genomic signals such as gene copy number variations or SNP profiles that may represent interesting biomarker candidates.

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.

Comparative Analysis of Leishmania major Nucleoside Hydrolases Toward Selecting Multi-target Strategy

PURPOSE

Leishmania causes multiple types of leishmaniasis in different parts of the world. It has a lack of metabolic machine to produce purine bases. Therefore, the parasite produces purine bases through the breakdown of nutritional nucleotides and it makes the nucleoside hydrolases (NHs) good drug targets. They have different substrate-preferring (SP) types. Our objectives were modeling and comparative analysis of these protein structures for Leishmania major.

METHOD

In this work, available sequences for all SP types of L. major NH enzymes including inosine-uridine preferring NH (IUNH), inosine-guanosine preferring NH (IGNH), and inosine-adenosine-guanosine preferring NH (IAGNH) were used to make 24 structural models via SWISS-MODEL and LOMETS. After evaluating the structural models, three enzyme structures were finalized and used to analyze substrate-binding pockets.

RESULTS

The three SP types of L. major NH enzymes that can breakdown purine nucleosides were highly different in terms of sequence, structure, and profile of interacting residues within the substrate-binding pockets. In this study, new enzyme structures have been presented for three SP types and they have been compared in different aspects and it indicated that they were very different from each other.

CONCLUSION

Although, previously indicated that from these three SP types in genera other than Leishmania, the role of IGNH and IAGNH was greater than IUNH in supplying purine bases, till this work, just IUNH has been structurally studied and used in drug-design investigations for Leishmania. Therefore, we are offering to use all three SP types of NHs as multi-target strategy in anti-leishmaniosis drug-design studies.

Cutaneous Leishmania mexicana Infections in the United States: Defining Strains Through Endemic Human Pediatric Cases in Northern Texas

Over a six-month span, three patients under five years old with cutaneous leishmaniasis presented to the Pediatric Infectious Diseases Clinic at the University of Texas Southwestern Medical Center/Children's Health Dallas. None had traveled outside of the United States (US); all had confirmed L. mexicana infections by PCR. We provide case descriptions and images to increase the awareness of this disease among US physicians and scientists. Two patients responded to fluconazole, but one required topical paromomycin. Combining these cases with guidelines and our literature review, we suggest that: 1) higher doses (ten-twelve mg/kg/day) of fluconazole should be considered in young children to maximize likelihood and rapidity of response and 2) patients should transition to alternate agents if they do not respond to high-dose fluconazole within six weeks. Furthermore, and of particular interest to the broad microbiology community, we used samples from these cases as a proof-of-concept to propose a mechanism to strain-type US-endemic L. mexicana. For our analysis, we sequenced three housekeeping genes and the internal transcribed sequence 2 of the ribosomal RNA gene. We identified genetic changes that not only allow us to distinguish US-based L. mexicana strains from strains found in other areas of the Americas, but also establish polymorphisms that differ between US isolates. These techniques will allow documentation of genetic changes in this parasite as its range expands. Hence, our cases of cutaneous leishmaniasis provide significant evolutionary, treatment and public health implications as climate change increases exposure to formerly tropical diseases in previously non-endemic areas.

Diversity and dissemination of viruses in pathogenic protozoa

Viruses are the most abundant biological entities on Earth and play a significant role in the evolution of many organisms and ecosystems. In pathogenic protozoa, the presence of viruses has been linked to an increased risk of treatment failure and severe clinical outcome. Here, we studied the molecular epidemiology of the zoonotic disease cutaneous leishmaniasis in Peru and Bolivia through a joint evolutionary analysis of Leishmania braziliensis and their dsRNA Leishmania virus 1. We show that parasite populations circulate in tropical rainforests and are associated with single viral lineages that appear in low prevalence. In contrast, groups of hybrid parasites are geographically and ecologically more dispersed and associated with an increased prevalence, diversity and spread of viruses. Our results suggest that parasite gene flow and hybridization increased the frequency of parasite-virus symbioses, a process that may change the epidemiology of leishmaniasis in the region.

Drug Discovery for Cutaneous Leishmaniasis: A Review of Developments in the Past 15 Years

Leishmaniasis is a group of vector-borne, parasitic diseases caused by over 20 species of the protozoan Leishmania spp. The three major disease classifications, cutaneous, visceral, and mucocutaneous, have a range of clinical manifestations from self-healing skin lesions to hepatosplenomegaly and mucosal membrane damage to fatality. As a neglected tropical disease, leishmaniasis represents a major international health challenge, with nearly 350 million people living at risk of infection a year. The current chemotherapeutics used to treat leishmaniasis have harsh side effects, prolonged and costly treatment regimens, as well as emerging drug resistance, and are predominantly used for the treatment of visceral leishmaniasis. There is an undeniable need for the identification and development of novel chemotherapeutics targeting cutaneous leishmaniasis (CL), largely ignored by concerted drug development efforts. CL is mostly non-lethal and the most common presentation of this disease, with nearly 1 million new cases reported annually. Recognizing this unaddressed need, substantial yet fragmented progress in early drug discovery efforts for CL has occurred in the past 15 years and was outlined in this review. However, further work needs to be carried out to advance early discovery candidates towards the clinic. Importantly, there is a paucity of investment in the translation and development of therapies for CL, limiting the emergence of viable solutions to deal with this serious and complex international health problem.

Identification of a unique conserved region from a kinetoplastid genome of Leishmania orientalis (formerly named Leishmania siamensis) strain PCM2 in Thailand

Mitochondrial DNAs (mtDNAs) appear in almost all eukaryotic species and are useful molecular markers for phylogenetic studies and species identification. Kinetoplast DNAs (kDNAs) are structurally complex circular mtDNA networks in kinetoplastids, divided into maxicircles and minicircles. Despite several kDNAs of many Leishmania species being examined, the kDNAs of the new species, Leishmania orientalis (formerly named Leishmania siamensis) strain PCM2, have not been explored. This study aimed to investigate the maxicircle and minicircle DNAs of L. orientalis strain PCM2 using hybrid genome sequencing technologies and bioinformatic analyses. The kDNA sequences were isolated and assembled using the SPAdes hybrid assembler from the Illumina short-read and PacBio long-read data. Circular contigs of the maxicircle and minicircle DNAs were reconstructed and confirmed by BLASTn and rKOMICs programs. The kDNA genome was annotated by BLASTn before the genome comparison and phylogenetic analysis by progressiveMauve, MAFFT, and MEGA programs. The maxicircle of L. orientalis strain PCM2 (18,215 bp) showed 99.92% similarity and gene arrangement to Leishmania enriettii strain LEM3045 maxicircle with variation in the 12s rRNA gene and divergent region. Phylogenetics of the whole sequence, coding regions, divergent regions, and 12s rRNA gene also confirmed this relationship and subgenera separation. The identified 105 classes of minicircles (402-1177 bp) were clustered monophyletically and related to the Leishmania donovani minicircles. The kinetoplast maxicircle and minicircle DNAs of L. orientalis strain PCM2 contained a unique conserved region potentially useful for specific diagnosis of L. orientalis and further exploration of this parasite population genetics in Thailand and related regions.

Handheld Ultra-Fast Duplex Polymerase Chain Reaction Assays and Lateral Flow Detection and Identification of Leishmania Parasites for Cutaneous Leishmaniases Diagnosis

Early and accurate detection of infectious diseases is a key step for surveillance, epidemiology and control, notably timely disease diagnosis, patient management and follow-up. In this study, we aimed to develop handheld ultra-fast duplex PCR assays coupled to amplicon detection by lateral flow (LF) immunoassay to deliver a rapid and simple molecular diagnostic test for concomitant detection and identification of the main Leishmania parasites encountered in Tunisia. We selected two DNA targets to amplify L. major/L. tropica and L. infantum/L. tropica groups of species DNAs, respectively. We optimized the experimental conditions of a duplex ultra-fast PCR. The amplification is performed using a portable Palm convection PCR machine within 18 min, and the products are detected using an LF cassette within 10 min. The test allows the identification of the infecting species according to the position and number of test lines revealed. Tested on a selection of DNAs of representative Leishmania strains of the three studied species (N = 37), the ultra-fast duplex PCR-LF showed consistent, stable and reproducible results. The analytical limit of detection of the test was 0.4 pg for L. major, 4 pg for L. infantum and 40 pg for L. tropica.

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.

Complete assembly, annotation of virulence genes and CRISPR editing of the genome of Leishmania amazonensis PH8 strain

We report the sequencing and assembly of the PH8 strain of Leishmania amazonensis one of the etiological agents of leishmaniasis. After combining data from long Pacbio reads, short Illumina reads and synteny with the Leishmania mexicana genome, the sequence of 34 chromosomes with 8317 annotated genes was generated. Multigene families encoding three virulence factors, A2, amastins and the GP63 metalloproteases, were identified and compared to their annotation in other Leishmania species. As they have been recently recognized as virulence factors essential for disease establishment and progression of the infection, we also identified 14 genes encoding proteins involved in parasite iron and heme metabolism and compared to genes from other Trypanosomatids. To follow these studies with a genetic approach to address the role of virulence factors, we tested two CRISPR-Cas9 protocols to generate L. amazonensis knockout cell lines, using the Miltefosine transporter gene as a proof of concept.

Parasite Detection in Visceral Leishmaniasis Samples by Dye-Based qPCR Using New Gene Targets of Leishmania infantum and Crithidia

Visceral leishmaniasis (VL) is a neglected disease considered a serious public health problem, especially in endemic countries. Several studies have discovered monoxenous trypanosomatids (Leptomonas and Crithidia) in patients with VL. In different situations of leishmaniasis, investigations have examined cases of co-infection between Leishmania spp. and Crithidia spp. These coinfections have been observed in a wide range of vertebrate hosts, indicating that they are not rare. Diagnostic techniques require improvements and more robust tools to accurately detect the causative agent of VL. This study aimed to develop a real-time quantitative dye-based PCR (qPCR) assay capable of distinguishing Leishmania infantum from Crithidia-related species and to estimate the parasite load in samples of VL from humans and animals. The primer LinJ31_2420 targets an exclusive phosphatase of L. infantum; the primer Catalase_LVH60-12060_1F targets the catalase gene of Crithidia. Therefore, primers were designed to detect L. infantum and Crithidia sp. LVH60A (a novel trypanosomatid isolated from VL patients in Brazil), in samples related to VL. These primers were considered species-specific, based on sequence analysis using genome data retrieved from the TriTryp database and the genome assembling of Crithidia sp. LVH60A strain, in addition to experimental and clinical data presented herein. This novel qPCR assay was highly accurate in identifying and quantifying L. infantum and Crithidia sp. LVH60A in samples obtained experimentally (in vitro and in vivo) or collected from hosts (humans, dogs, cats, and vectors). Importantly, the screening of 62 cultured isolates from VL patients using these primers surprisingly revealed that 51 parasite cultures were PCR+ for Crithidia sp. In addition, qPCR assays identified the co-infection of L. infantum with Crithidia sp. LVH60A in two new VL cases in Brazil, confirming the suspicion of co-infection in a previously reported case of fatal VL. We believe that the species-specific genes targeted in this study can be helpful for the molecular diagnosis of VL, as well as for elucidating suspected co-infections with monoxenous-like trypanosomatids, which is a neglected fact of a neglected disease.

Kinetoplast Genome of Leishmania spp. Is under Strong Purifying Selection

Instability is an intriguing characteristic of many protist genomes, and trypanosomatids are not an exception in this respect. Some regions of trypanosomatid genomes evolve fast. For instance, the trypanosomatid mitochondrial (kinetoplast) genome consists of fairly conserved maxicircle and minicircle molecules that can, nevertheless, possess high nucleotide substitution rates between closely related strains. Recent experiments have demonstrated that rapid laboratory evolution can result in the non-functionality of multiple genes of kinetoplast genomes due to the accumulation of mutations or loss of critical genomic components. An example of a loss of critical components is the reported loss of entire minicircle classes in Leishmania tarentolae during laboratory cultivation, which results in an inability to generate some correctly encoded genes. In the current work, we estimated the evolutionary rates of mitochondrial and nuclear genome regions of multiple natural Leishmania spp. We analyzed synonymous and non-synonymous substitutions and, rather unexpectedly, found that the coding regions of kinetoplast maxicircles are among the most variable regions of both genomes. In addition, we demonstrate that synonymous substitutions greatly predominate among maxicircle coding regions and that most maxicircle genes show signs of purifying selection. These results imply that maxicircles in natural Leishmania populations remain functional despite their high mutation rate.

Translational reprogramming as a driver of antimony-drug resistance in Leishmania

Leishmania is a unicellular protozoan that has a limited transcriptional control and mostly uses post-transcriptional regulation of gene expression, although the molecular mechanisms of the process are still poorly understood. Treatments of leishmaniasis, pathologies associated with Leishmania infections, are limited due to drug resistance. Here, we report dramatic differences in mRNA translation in antimony drug-resistant and sensitive strains at the full translatome level. The major differences (2431 differentially translated transcripts) were demonstrated in the absence of the drug pressure supporting that complex preemptive adaptations are needed to efficiently compensate for the loss of biological fitness once they are exposed to the antimony. In contrast, drug-resistant parasites exposed to antimony activated a highly selective translation of only 156 transcripts. This selective mRNA translation is associated with surface protein rearrangement, optimized energy metabolism, amastins upregulation, and improved antioxidant response. We propose a novel model that establishes translational control as a major driver of antimony-resistant phenotypes in Leishmania.

Hit-to-lead optimization of a pyrazinylpiperazine series against Leishmania infantum and Leishmania braziliensis

An early hit-to-lead optimization of a novel pyrazinylpiperazine series against L. infantum and L. braziliensis has been performed after an extensive SAR focusing on the benzoyl fragment of hit (4). Deletion of the meta-Cl of (4) led to the obtention of the para-hydroxyl derivative (12), on which the design of most monosubstituted derivatives of the SAR was based. Further optimization of the series, involving disubstituted benzoyl fragments and the hydroxyl substituent of (12), allowed the obtention of a total of 15 compounds with increased antileishmanial potency (IC₅₀ < 10 μM), nine of which displayed activity in the low micromolar range (IC₅₀ < 5 μM). This optimization ultimately identified the ortho, meta-dihydroxyl derivative (46) as an early lead for this series (IC₅₀ (L. infantum) = 2.8 μM, IC₅₀ (L. braziliensis) = 0.2 μM). Additional assessment of some selected compounds against other trypanosomatid parasites revealed that this series is selective towards Leishmania parasites, and in silico ADMET predictions revealed satisfactory profiles for these compounds, allowing further lead optimization of the pyrazinylpiperazine class against Leishmania.

Leishmania infantum (JPCM5) Transcriptome, Gene Models and Resources for an Active Curation of Gene Annotations

Leishmania infantum is one of the causative agents of visceral leishmaniases, the most severe form of leishmaniasis. An improved assembly for the L. infantum genome was published five years ago, yet delineation of its transcriptome remained to be accomplished. In this work, the transcriptome annotation was attained by a combination of both short and long RNA-seq reads. The good agreement between the results derived from both methodologies confirmed that transcript assembly based on Illumina RNA-seq and further delimitation according to the positions of spliced leader (SAS) and poly-A (PAS) addition sites is an adequate strategy to annotate the transcriptomes of Leishmania, a procedure previously used for transcriptome annotation in other Leishmania species and related trypanosomatids. These analyses also confirmed that the Leishmania transcripts boundaries are relatively slippery, showing extensive heterogeneity at the 5′- and 3′-ends. However, the use of RNA-seq reads derived from the PacBio technology (referred to as Iso-Seq) allowed the authors to uncover some complex transcription patterns occurring at particular loci that would be unnoticed by the use of short RNA-seq reads alone. Thus, Iso-Seq analysis provided evidence that transcript processing at particular loci would be more dynamic than expected. Another noticeable finding was the observation of a case of allelic heterozygosity based on the existence of chimeric Iso-Seq reads that might be generated by an event of intrachromosomal recombination. In addition, we are providing the L. infantum gene models, including both UTRs and CDS regions, that would be helpful for undertaking whole-genome expression studies. Moreover, we have built the foundations of a communal database for the active curation of both gene/transcript models and functional annotations for genes and proteins.

Parasite hybridization promotes spreading of endosymbiotic viruses

Viruses are the most abundant biological entities on Earth and play a significant role in the evolution of many organisms and ecosystems. In pathogenic protozoa, the presence of endosymbiotic viruses has been linked to an increased risk of treatment failure and severe clinical outcome. Here, we studied the molecular epidemiology of the zoonotic disease cutaneous leishmaniasis in Peru and Bolivia through a joint evolutionary analysis of Leishmania braziliensis parasites and their endosymbiotic Leishmania RNA virus. We show that parasite populations circulate in isolated pockets of suitable habitat and are associated with single viral lineages that appear in low prevalence. In contrast, groups of hybrid parasites were geographically and ecologically dispersed, and commonly infected from a pool of genetically diverse viruses. Our results suggest that parasite hybridization, likely due to increased human migration and ecological perturbations, increased the frequency of endosymbiotic interactions known to play a key role in disease severity.

Leishmania allelic selection during experimental sand fly infection correlates with mutational signatures of oxidative DNA damage

Trypanosomatid pathogens are transmitted by blood-feeding insects, causing devastating human infections. These parasites show important phenotypic shifts that often impact parasite pathogenicity, tissue tropism, or drug susceptibility. The evolutionary mechanisms that allow for the selection of such adaptive phenotypes remain only poorly investigated. Here, we use Leishmania donovani as a trypanosomatid model pathogen to assess parasite evolutionary adaptation during experimental sand fly infection. Comparing the genome of the parasites before and after sand fly infection revealed a strong population bottleneck effect as judged by allele frequency analysis. Apart from random genetic drift caused by the bottleneck effect, our analyses revealed haplotype and allelic changes during sand fly infection that seem under natural selection given their convergence between independent biological replicates. Our analyses further uncovered signature mutations of oxidative DNA damage in the parasite genomes after sand fly infection, suggesting that Leishmania suffers from oxidative stress inside the insect digestive tract. Our results propose a model of Leishmania genomic adaptation during sand fly infection, with oxidative DNA damage and DNA repair processes likely driving haplotype and allelic selection. The experimental and computational framework presented here provides a useful blueprint to assess evolutionary adaptation of other eukaryotic pathogens inside their insect vectors, such as Plasmodium spp, Trypanosoma brucei, and Trypanosoma cruzi.

Selective whole-genome amplification reveals population genetics of Leishmania braziliensis directly from patient skin biopsies

In Brazil, Leishmania braziliensis is the main causative agent of the neglected tropical disease, cutaneous leishmaniasis (CL). CL presents on a spectrum of disease severity with a high rate of treatment failure. Yet the parasite factors that contribute to disease presentation and treatment outcome are not well understood, in part because successfully isolating and culturing parasites from patient lesions remains a major technical challenge. Here we describe the development of selective whole genome amplification (SWGA) for Leishmania and show that this method enables culture-independent analysis of parasite genomes obtained directly from primary patient skin samples, allowing us to circumvent artifacts associated with adaptation to culture. We show that SWGA can be applied to multiple Leishmania species residing in different host species, suggesting that this method is broadly useful in both experimental infection models and clinical studies. SWGA carried out directly on skin biopsies collected from patients in Corte de Pedra, Bahia, Brazil, showed extensive genomic diversity. Finally, as a proof-of-concept, we demonstrated that SWGA data can be integrated with published whole genome data from cultured parasite isolates to identify variants unique to specific geographic regions in Brazil where treatment failure rates are known to be high. SWGA provides a relatively simple method to generate Leishmania genomes directly from patient samples, unlocking the potential to link parasite genetics with host clinical phenotypes.

Genomic analysis of Leishmania turanica strains from different regions of Central Asia

The evolution in Leishmania is governed by the opposite forces of clonality and sexual reproduction, with vicariance being an important factor. As such, Leishmania spp. populations may be monospecific or mixed. Leishmania turanica in Central Asia is a good model to compare these two types. In most areas, populations of L. turanica are mixed with L. gerbilli and L. major. Notably, co-infection with L. turanica in great gerbils helps L. major to withstand a break in the transmission cycle. Conversely, the populations of L. turanica in Mongolia are monospecific and geographically isolated. In this work, we compare genomes of several well-characterized strains of L. turanica originated from monospecific and mixed populations in Central Asia in order to shed light on genetic factors, which may drive evolution of these parasites in different settings. Our results illustrate that evolutionary differences between mixed and monospecific populations of L. turanica are not dramatic. On the level of large-scale genomic rearrangements, we confirmed that different genomic loci and different types of rearrangements may differentiate strains originated from mixed and monospecific populations, with genome translocations being the most prominent example. Our data suggests that L. turanica has a significantly higher level of chromosomal copy number variation between the strains compared to its sister species L. major with only one supernumerary chromosome. This suggests that L. turanica (in contrast to L. major) is in the active phase of evolutionary adaptation.

Genome sequence of Leishmania mexicana MNYC/BZ/62/M379 expressing Cas9 and T7 RNA polymerase

We present the genome sequence of Leishmania mexicana MNYC/BZ/62/M379 modified to express Cas9 and T7 RNA-polymerase, revealing high similarity to the reference genome (MHOM/GT2001/U1103). Through RNAseq-based annotation of coding sequences and untranslated regions, we provide primer sequences for construct and sgRNA template generation for CRISPR-assisted gene deletion and endogenous tagging.

Immunoproteomics and phage display in the context of leishmaniasis complexity

Leishmaniasis is defined as a complex of diseases caused by protozoan parasites of the genus Leishmania, which comprises 20 parasite species pathogenic to mammalians, such as humans and dogs. From a clinical point of view, and considering the diversity and biological complexity of the parasites, vectors, and vertebrate hosts, leishmaniasis is classified according to the distinct clinical manifestations, such as tegumentary (involving the cutaneous, mucosal, and cutaneous-diffuse forms) and visceral leishmaniasis. Many issues and challenges remain unaddressed, which could be attributed to the complexity and diversity of the disease. The current demand for the identification of new Leishmania antigenic targets for the development of multicomponent-based vaccines, as well as for the production of specific diagnostic tests, is evident. In recent years, biotechnological tools have allowed the identification of several Leishmania biomarkers that might potentially be used for diagnosis and have an application in vaccine development. In this Mini Review, we discuss the different aspects of this complex disease that have been addressed by technologies such as immunoproteomics and phage display. It is extremely important to be aware of the potential applications of antigens selected in different screening context, so that they can be used appropriately, so understanding their performance, characteristics, and self-limitations.

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.

Parasite Genotype Is a Major Predictor of Mortality from Visceral Leishmaniasis

Visceral leishmaniasis (VL) is a potentially fatal disease caused mainly by Leishmania infantum in South America and Leishmania donovani in Asia and Africa. Disease outcomes have been associated with patient genotype, nutrition, age, sex, comorbidities, and coinfections. In this study, we examine the effects of parasite genetic variation on VL disease severity in Brazil. We collected and sequenced the genomes of 109 L. infantum isolates from patients in northeastern Brazil and retrieved matching patient clinical data from medical records, including mortality, sex, HIV coinfection, and laboratory data (creatinine, hemoglobin, and leukocyte and platelet counts). We identified genetic differences between parasite isolates, including single nucleotide polymorphisms (SNPs), small insertions/deletions (indels), and variations in genic, intergenic, and chromosome copy numbers (copy number variants [CNVs]). To describe associations between the parasite genotypes and clinical outcomes, we applied quantitative genetics methods of heritability and genome-wide association studies (GWAS), treating clinical outcomes as traits that may be influenced by parasite genotype. Multiple aspects of the genetic analysis indicate that parasite genotype affects clinical outcomes. We estimate that parasite genotype explains 83% chance of mortality (narrow-sense heritability [h²] = 0.83 ± 0.17) and has a significant relationship with patient sex (h² = 0.60 ± 0.27). Impacts of parasite genotype on other clinical traits are lower (h² ≤ 0.34). GWAS analysis identified multiple parasite genetic loci that were significantly associated with clinical outcomes; 17 CNVs were significantly associated with mortality, two with creatinine, and one with bacterial coinfection, jaundice, and HIV coinfection, and two SNPs/indels and six CNVs were associated with age, jaundice, HIV and bacterial coinfections, creatinine, and/or bleeding sites. Parasite genotype is an important factor in VL disease severity in Brazil. Our analysis indicates that specific genetic differences between parasites act as virulence factors, enhancing risks of severe disease and mortality. More detailed understanding of these virulence factors could be exploited for novel therapies. IMPORTANCE Multiple factors contribute to the risk of mortality from visceral leishmaniasis (VL), including, patient genotype, comorbidities, and nutrition. Many of these factors are influenced by socioeconomic biases. Our work suggests that the virulence of the infecting parasite is an important risk factor for mortality. We pinpoint some specific genomic markers that are associated with mortality, which can lead to a greater understanding of the molecular mechanisms that cause severe VL disease, to the identification of genetic markers for virulent parasites, and to the development of drug and vaccine therapies.

Self-Hybridization in Leishmania major

Genetic exchange between different Leishmania strains in the sand fly vector has been experimentally demonstrated and is supported by population genetic studies. In nature, opportunities for Leishmania interstrain mating are restricted to flies biting multiply infected hosts or through multiple bites of different hosts. In contrast, self-mating could occur in any infected sand fly. By crossing two recombinant lines derived from the same Leishmania major strain, each expressing a different drug-resistance marker, self-hybridization in L. major was confirmed in a natural sand fly vector, Phlebotomus duboscqi, and in frequencies comparable to interstrain crosses. We provide the first high resolution, whole-genome sequencing analysis of large numbers of selfing progeny, their parents, and parental subclones. Genetic exchange consistent with classical meiosis is supported by the biallelic inheritance of the rare homozygous single nucleotide polymorphisms (SNPs) that arose by mutation during the generation of the parental clones. In contrast, heterozygous SNPs largely failed to be transmitted in Mendelian ratios for reasons not understood. SNPs that were heterozygous in both parents, however, recombined to produce homozygous alleles in some hybrids. For trisomic chromosomes present in both parents, transmittal to the progeny was only altered by self-hybridization, involving a gain or loss of somy in frequencies predicted by a meiotic process. Whole-genome polyploidization was also observed in the selfing progeny. Thus, self-hybridization in Leishmania, with its potential to occur in any infected sand fly, may be an important source of karyotype variation, loss of heterozygosity, and functional diversity. IMPORTANCE Leishmania are parasitic protozoa that cause a wide spectrum of diseases collectively known as the leishmaniases. Sexual reproduction in Leishmania has been proposed as an important source of genetic diversity and has been formally demonstrated to occur inside the sand fly vector midgut. Nevertheless, in the wild, opportunities for genetic exchange between different Leishmania species or strains are restricted by the capacity of different Leishmania strains to colonize the same sand fly. In this work, we report the first high resolution, whole-genome sequence analysis of intraclonal genetic exchange as a type of self-mating in Leishmania. Our data reveal that self-hybridization can occur with comparable frequency as interstrain mating under experimental lab conditions, leading to important genomic alterations that can potentially take place within every naturally infected sand fly.

Genome sequence of Leishmania mexicana MNYC/BZ/62/M379 expressing Cas9 and T7 RNA polymerase

We present the genome sequence of Leishmania mexicana MNYC/BZ/62/M379 modified to express Cas9 and T7 RNA-polymerase, revealing high similarity to the reference genome (MHOM/GT2001/U1103). Through RNAseq-based annotation of coding sequences and untranslated regions, we provide primer sequences for construct and sgRNA template generation for CRISPR-assisted gene deletion and endogenous tagging.

ITS1 and cpb genetic polymorphisms in Algerian and Tunisian Leishmania infantum isolates from humans and dogs

Leishmania (L.) infantum strains, isolated from varying hosts and clinical manifestations (cutaneous, visceral and canine leishmaniasis), were investigated in order to understand the genetic polymorphisms within this species in Algeria and Tunisia. Two DNA-based typing methods were tested in order to evaluate their effectiveness against Multilocus enzyme electrophoresis (MLEE), widely considered as the reference method for Leishmania parasite typing. On the other hand, MLEE is cumbersome, high-cost, time consuming and frequently does not detect intra-species genetic polymorphisms. In this work, we used two molecular target regions to discriminate L. infantum strains, Internal transcribed spacer 1 (ITS1) and the cysteine proteinase B (cpb). The ITS1 region offers good resolution for Leishmania discrimination but does not spotlight intra-species polymorphisms. In contrast, cpbE and cpbF PCR-Sequencing demonstrated a certain variability within CL and VL Algerian and Tunisian L. infantum isolates. Following phylogenetic analyses of 44 L. infantum isolates, two main groups were identified, a group with 39 bp deletion in the cpb sequence, composed of cutaneous, visceral and canine isolates from both countries with no significant clinical or geographic distribution; these samples were typed as MON-1, MON-24, and MON-80 zymodemes. A second group which presents a clear clusterization of Tunisian cutaneous strains belonging to the L. infantum MON-24. This group, with no deletion in the mature domain of the cpb gene sequence, should be further explored with a higher number of samples.

Investigating the Leishmania donovani sacp Gene and Its Role in Macrophage Infection and Survival in Mice

The protozoan parasite Leishmania donovani is a causative agent of the neglected tropical disease known as visceral leishmaniasis, which can be lethal when untreated. Studying Leishmania viru-lence factors is crucial in determining how the parasite causes disease and identifying new targets for treatment. One potential virulence factor is L. donovani’s abundantly secreted protein: secreted acid phosphatase (SAcP). Whole-genome analysis revealed that the sacp gene was present in three copies in wild type L. donovani. Using CRISPR-Cas9 gene editing; we generated a sacp gene knockout termed LdΔSAcP, which demonstrated a loss of both the SAcP protein and an associated reduction in secreted acid phosphatase activity. Genome sequencing confirmed the precise dele-tion of the sacp gene in LdΔSAcP and identified several changes in the genome. LdΔSAcP demonstrated no significant changes in promastigote proliferation or its ability to infect and survive in macrophages compared to the wildtype strain. LdΔSAcP also demonstrated no change in murine liver infection; however, survival was impaired in the spleen. Taken together these results show that SAcP is not necessary for the survival of promastigotes in culture but may support long-term survival in the spleen. These observations also show that the use of CRISPR gene editing and WGS together are effective to investigate the function and phenotype of complex potential drug targets such as multicopy genes.

Three types of Leishmania mexicana amastigotes: Proteome comparison by quantitative proteomic analysis

Leishmania is the unicellular parasite transmitted by phlebotomine sand fly bite. It exists in two different forms; extracellular promastigotes, occurring in the gut of sand flies, and intracellular, round-shaped amastigotes residing mainly in vertebrate macrophages. As amastigotes originating from infected animals are often present in insufficient quality and quantity, two alternative types of amastigotes were introduced for laboratory experiments: axenic amastigotes and amastigotes from macrophages infected in vitro. Nevertheless, there is very little information about the degree of similarity/difference among these three types of amastigotes on proteomic level, whose comparison is crucial for assessing the suitability of using alternative types of amastigotes in experiments. In this study, L. mexicana amastigotes obtained from lesion of infected BALB/c mice were proteomically compared with alternatively cultivated amastigotes (axenic and macrophage-derived ones). Amastigotes of all three types were isolated, individually treated and analysed by LC-MS/MS proteomic analysis with quantification using TMT10-plex isobaric labeling. Significant differences were observed in the abundance of metabolic enzymes, virulence factors and proteins involved in translation and condensation of DNA. The most pronounced differences were observed between axenic amastigotes and lesion-derived amastigotes, macrophage-derived amastigotes were mostly intermediate between axenic and lesion-derived ones.

Laboratory diagnostics for human Leishmania infections: a polymerase chain reaction-focussed review of detection and identification methods

Leishmania infections span a range of clinical syndromes and impact humans from many geographic foci, but primarily the world's poorest regions. Transmitted by the bite of a female sand fly, Leishmania infections are increasing with human movement (due to international travel and war) as well as with shifts in vector habitat (due to climate change). Accurate diagnosis of the 20 or so species of Leishmania that infect humans can lead to the successful treatment of infections and, importantly, their prevention through modelling and intervention programs. A multitude of laboratory techniques for the detection of Leishmania have been developed over the past few decades, and although many have drawbacks, several of them show promise, particularly molecular methods like polymerase chain reaction. This review provides an overview of the methods available to diagnostic laboratories, from traditional techniques to the now-preferred molecular techniques, with an emphasis on polymerase chain reaction-based detection and typing methods.

Leishmania major Strain-Dependent Macrophage Activation Contributes to Pathogenicity in the Absence of Lymphocytes

Infection of C57BL/6 wild-type mice with Leishmania major 5-ASKH or Friedlin strains results in relatively similar pathogenicity with self-healing lesions within weeks. Parasite clearance depends on nitric oxide production by activated macrophages in response to cytokines produced mainly by CD4⁺ Th1 cells. In contrast, C57BL/6 Rag2 knockout mice, which lack T and B lymphocytes, show distinct pathologies during infection with these strains. Despite of the similar parasite number, the 5-ASKH infection induced severe inflammation rather than the Friedlin. To determine the immunological factors behind this phenomenon, we infected C57BL/6 Rag2 knockout mice with these two strains and compared immune cell kinetics and macrophage activation status. Compared with the Friedlin strain, the 5-ASKH strain elicited increased pathology associated with the accumulation of CD11b^high, Ly6G^high neutrophils by week four and increased the expression of macrophage activation markers. We then analyzed the differentially expressed transcripts in infected bone marrow-derived macrophages by RNA sequencing. It showed upregulation of multiple inflammatory transcripts, including Toll-like receptor 1/2 (TLR1/2), CD69, and CARD14, upon 5-ASKH infection. Our findings suggest that different L. major strains can trigger distinct macrophage activation, contributing to the disease outcome observed in the absence of lymphocytes but not in the presence of lymphocytes. IMPORTANCE Disease manifestations of cutaneous leishmaniasis (CL) range from self-healing cutaneous lesions to chronic forms of the disease, depending on the infecting Leishmania sp. and host immune protection. Previous works on mouse models of CL show the distinct pathogenicity of Leishmania major strains in the absence of lymphocytes. However, the mechanisms of this pathology remain uncovered. In the trial to understand the immunological process involved in lymphocyte-independent pathology, we have found a specific induction of macrophages by different L. major strains that affect their ability to mount innate responses leading to neutrophilic pathology when lymphocytes are ablated.

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.

Genome-wide analysis reveals allelic variation and chromosome copy number variation in paromomycin-resistant Leishmania donovani

In the absence of adequate diagnosis and treatment, leishmaniasis remains a major public health concern on a global scale. Drug resistance remains a key obstacle in controlling and eliminating visceral leishmaniasis. The therapeutic gap due to lack of target-specific medicine and vaccine can be minimized by obtaining parasite's genomic information. This study compared whole-genome sequence of paromomycin-resistant parasite (K133PMM) developed through in vitro adaptation and selection with sensitive Leishmania clinical isolate (K133WT). We found a large number of upstream and intergenic gene variations in K133PMM. There were 259 single nucleotide polymorphisms (SNPs), 187 insertion-deletion (InDels), and 546 copy number variations (CNVs) identified. Most of the genomic variations were found in the gene's upstream and non-coding regions. Ploidy estimation revealed chromosome 5 in tetrasomy and 6, 9, and 12 in trisomy, uniquely in K133PMM. These contain the genes for protein degradation, parasite motility, autophagy, cell cycle maintenance, and drug efflux membrane transporters. Furthermore, we also observed reduction in ploidy of chromosomes 15, 20, and 23, in the resistant parasite containing mostly the genes for hypothetical proteins and membrane transporters. We chronicled correlated genomic conversion and aneuploidy in parasites and hypothesize that this led to rapid evolutionary changes in response to drug induced pressure, which causes them to become resistant.

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.

Making the Most of Its Short Reads: A Bioinformatics Workflow for Analysing the Short-Read-Only Data of Leishmania orientalis (Formerly Named Leishmania siamensis) Isolate PCM2 in Thailand

Simple Summary

Leishmaniasis is a parasitic disease caused by flagellated protozoa of the genus Leishmania. Multiple genome sequencing platforms have been employed to complete Leishmania genomes at the expense of high cost. This study proposes an integrative bioinformatic workflow for assembling only the short-read data of Leishmania orientalis isolate PCM2 from Thailand and produce an acceptable-quality genome for further genomic analysis. This workflow gives extensive information required for identifying strain-specific markers and virulence-associated genes useful for drug and vaccine development before a more exhaustive and expensive investigation.

Abstract

Background: Leishmania orientalis (formerly named Leishmania siamensis) has been neglected for years in Thailand. The genomic study of L. orientalis has gained much attention recently after the release of the first high-quality reference genome of the isolate LSCM4. The integrative approach of multiple sequencing platforms for whole-genome sequencing has proven effective at the expense of considerably expensive costs. This study presents a preliminary bioinformatic workflow including the use of multi-step de novo assembly coupled with the reference-based assembly method to produce high-quality genomic drafts from the short-read Illumina sequence data of L. orientalis isolate PCM2. Results: The integrating multi-step de novo assembly by MEGAHIT and SPAdes with the reference-based method using the L. enriettii genome and salvaging the unmapped reads resulted in the 30.27 Mb genomic draft of L. orientalis isolate PCM2 with 3367 contigs and 8887 predicted genes. The results from the integrated approach showed the best integrity, coverage, and contig alignment when compared to the genome of L. orientalis isolate LSCM4 collected from the northern province of Thailand. Similar patterns of gene ratios and frequency were observed from the GO biological process annotation. Fifty GO terms were assigned to the assembled genomes, and 23 of these (accounting for 61.6% of the annotated genes) showed higher gene counts and ratios when results from our workflow were compared to those of the LSCM4 isolate. Conclusions: These results indicated that our proposed bioinformatic workflow produced an acceptable-quality genome of L. orientalis strain PCM2 for functional genomic analysis, maximising the usage of the short-read data. This workflow would give extensive information required for identifying strain-specific markers and virulence-associated genes useful for drug and vaccine development before a more exhaustive and expensive investigation.

Tissue tropism: Is it an intrinsic characteristic of Leishmania species?

The genus Leishmania comprises a wide range of species, some of which are pathogenic to humans and each of which has a different tissue preference, resulting in one of the three clinical forms of human leishmaniasis: visceral, cutaneous, or mucocutaneous. Although, all pathogenic species are deposited intradermally in the mammalian host upon an infectious sand fly bite, only the viscerotropic strains can leave the skin and reach the internal organs. We assume that Leishmania tissue tropism is not only the result of Leishmania genetic determinism but is also governed by the interaction of the parasite with different vectorial and human host elements. To shed light on these elements and key steps determining the course of the infection, we describe throughout this review the disease's progression from the early stages of infection taking place in the skin to the late stages succeeding in the parasite's visceral dissemination. Hence, we address the question of Leishmania tropism, through providing relevant hypotheses and answers gathered from the literature.

Bromodomain factor 5 is an essential regulator of transcription in Leishmania

Leishmania are unicellular parasites that cause human and animal diseases. Like other kinetoplastids, they possess large transcriptional start regions (TSRs) which are defined by histone variants and histone lysine acetylation. Cellular interpretation of these chromatin marks is not well understood. Eight bromodomain factors, the reader modules for acetyl-lysine, are found across Leishmania genomes. Using L. mexicana, Cas9-driven gene deletions indicate that BDF1–5 are essential for promastigotes. Dimerisable, split Cre recombinase (DiCre)-inducible gene deletion of BDF5 show it is essential for both promastigotes and murine infection. ChIP-seq identifies BDF5 as enriched at TSRs. XL-BioID proximity proteomics shows the BDF5 landscape is enriched for BDFs, HAT2, proteins involved in transcriptional activity, and RNA processing; revealing a Conserved Regulators of Kinetoplastid Transcription (CRKT) Complex. Inducible deletion of BDF5 causes global reduction in RNA polymerase II transcription. Our results indicate the requirement of Leishmania to interpret histone acetylation marks through the bromodomain-enriched CRKT complex for normal gene expression and cellular viability.

Leishmania use large (5–10 kb) transcriptional start regions, where the chromatin is highly enriched for acetylated histones, to drive the expression of polycistronic gene arrays. Here the authors show bromodomain-containing protein BDF5 is enriched at transcriptional start sites and its depletion leads to cell death in vitro and in murine infections, and they identify its interactors.

Stable Episomal Transfectant Leishmania infantum Promastigotes Over-Expressing the DEVH1 RNA Helicase Gene Down-Regulate Parasite Survival Genes

The compartmentalization of untranslated mRNA molecules in granules occurring in many eukaryotic organisms including trypanosomatids involves the formation of complexes between mRNA molecules and RNA-binding proteins (RBPs). The putative ATP-dependent DEAD/H RNA helicase (DEVH1) from Leishmania infantum (Kinetoplastida: Trypanosomatidae) is one such proteins. The objective of this research is finding differentially expressed genes in a stable episomal transfectant L. infantum promastigote line over-expressing DEVH1 in the stationary phase of growth in axenic culture to get insight into the biological roles of this RNA helicase in the parasite. Interestingly, genes related to parasite survival and virulence factors, such as the hydrophilic surface protein/small hydrophilic endoplasmic reticulum protein (HASP/SHERP) gene cluster, an amastin, and genes related to reactive oxygen species detoxification are down-regulated in DEVH1 transfectant promastigotes.

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.

Assembly of a Large Collection of Maxicircle Sequences and Their Usefulness for Leishmania Taxonomy and Strain Typing

Parasites of medical importance, such as Leishmania and Trypanosoma, are characterized by the presence of thousands of circular DNA molecules forming a structure known as kinetoplast, within the mitochondria. The maxicircles, which are equivalent to the mitochondrial genome in other eukaryotes, have been proposed as a promising phylogenetic marker. Using whole-DNA sequencing data, it is also possible to assemble maxicircle sequences as shown here and in previous works. In this study, based on data available in public databases and using a bioinformatics workflow previously reported by our group, we assembled the complete coding region of the maxicircles for 26 prototypical strains of trypanosomatid species. Phylogenetic analysis based on this dataset resulted in a robust tree showing an accurate taxonomy of kinetoplastids, which was also able to discern between closely related Leishmania species that are usually difficult to discriminate by classical methodologies. In addition, we provide a dataset of the maxicircle sequences of 60 Leishmania infantum field isolates from America, Western Europe, North Africa, and Eastern Europe. In agreement with previous studies, our data indicate that L. infantum parasites from Brazil are highly homogeneous and closely related to European strains, which were transferred there during the discovery of America. However, this study showed the existence of different L. infantum populations/clades within the Mediterranean region. A maxicircle signature for each clade has been established. Interestingly, two L. infantum clades were found coexisting in the same region of Spain, one similar to the American strains, represented by the Spanish JPCM5 reference strain, and the other, named “non-JPC like”, may be related to an important leishmaniasis outbreak that occurred in Madrid a few years ago. In conclusion, the maxicircle sequence emerges as a robust molecular marker for phylogenetic analysis and species typing within the kinetoplastids, which also has the potential to discriminate intraspecific variability.

Paving the Way: Contributions of Big Data to Apicomplexan and Kinetoplastid Research

In the age of big data an important question is how to ensure we make the most out of the resources we generate. In this review, we discuss the major methods used in Apicomplexan and Kinetoplastid research to produce big datasets and advance our understanding of Plasmodium, Toxoplasma, Cryptosporidium, Trypanosoma and Leishmania biology. We debate the benefits and limitations of the current technologies, and propose future advancements that may be key to improving our use of these techniques. Finally, we consider the difficulties the field faces when trying to make the most of the abundance of data that has already been, and will continue to be, generated.

Establish an allele-specific real-time PCR for Leishmania species identification

Background

Leishmaniasis is a serious neglected tropical disease that may lead to life-threatening outcome, which species are closely related to clinical diagnosis and patient management. The current Leishmania species determination method is not appropriate for clinical application. New Leishmania species identification tool is needed using clinical samples directly without isolation and cultivation of parasites.

Methods

A probe-based allele-specific real-time PCR assay was established for Leishmania species identification between Leishmania donovani and L. infantum for visceral leishmaniasis (VL) and among L. major, L. tropica and L. donovani/L. infantum for cutaneous leishmaniasis (CL), targeting hypoxanthine-guanine phosphoribosyl transferase (HGPRT) and spermidine synthase (SPDSYN) gene with their species-specific single nucleotide polymorphisms (SNPs). The limit of detection of this assay was evaluated based on 8 repeated tests with intra-assay standard deviation < 0.5 and inter-assay coefficients of variability < 5%. The specificity of this assay was tested with DNA samples obtained from Plasmodium falciparum, Toxoplasma gondii, Brucella melitensis and Orientia tsutsugamushi. Total 42 clinical specimens were used to evaluate the ability of this assay for Leishmania species identification. The phylogenetic tree was constructed using HGPRT and SPDSYN gene fragments to validate the performance of this assay.

Results

This new method was able to detect 3 and 12 parasites/reaction for VL and CL respectively, and exhibited no cross-reaction with P. falciparum, T. gondii, B. melitensis, O. tsutsugamushi and non-target species of Leishmania. Twenty-two samples from VL patients were identified as L. donovani (n = 3) and L. infantum (n = 19), and 20 specimens from CL patients were identified as L. major (n = 20), providing an agreement of 100% compared with sequencing results. For further validation, 29 sequences of HGPRT fragment from nine Leishmania species and 22 sequences from VL patients were used for phylogenetic analysis, which agreed with the results of this new method. Similar results were obtained with 43 sequences of SPDSYN fragment from 18 Leishmania species and 20 sequences from CL patients.

Conclusions

Our assay provides a rapid and accurate tool for Leishmania species identification which is applicable for species-adapted therapeutic schedule and patient management.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s40249-022-00992-y.

Leishmania diversity in bats from an endemic area for visceral and cutaneous leishmaniasis in Southeastern Brazil

This study aimed to determine the occurrence of Leishmania infection in bats in urban and wild areas in an endemic municipality for visceral and cutaneous leishmaniasis in Minas Gerais, Brazil. Between April 2014 to April 2015, 247 bats were captured and classified into 26 species belonging to Phyllostomidae (90.7%), Vespertilionidae (8.1%) and Molossidae (1.2%) families. Blood samples from 247 bats were collected and submitted to nested-PCR, targeting the variable V7-V8 region of the SSU rRNA gene, followed by sequencing of the PCR product. The overall infection rate of Leishmania spp. in bats was 4.4%. Of the eleven bats infected, ten were frugivorous bats: Artibeus planirostris (8/11), Artibeus lituratus (1/11) and Artibeus cinereus (1/11) and one a nectarivorous bat (Glossophaga soricina). None of the individuals exhibited macroscopic alterations in the skin, spleen or liver. Phylogenetic analysis separated Leishmania species in clades corresponding to the subgenera Viannia, Leishmania, and Mundinia, and supported that the isolates characterized in the present study clustered closely with Leishmania (Viannia) sp., Leishmania (Leishmania) infantum and Leishmania (Leishmania) amazonensis. Here we report for the first time the bat Artibeus cinereus as a host of Leishmania (Leishmania) amazonensis. In the study we found that the mean abundance of bats did not differ in wild habitats and urban areas and that bat-parasite interactions were similarly distributed in the two environments. On the other hand, further studies should be conducted in more recent times to verify whether there have been changes in these parameters.

Experimental evolution links post-transcriptional regulation to Leishmania fitness gain

The protozoan parasite Leishmania donovani causes fatal human visceral leishmaniasis in absence of treatment. Genome instability has been recognized as a driver in Leishmania fitness gain in response to environmental change or chemotherapy. How genome instability generates beneficial phenotypes despite potential deleterious gene dosage effects is unknown. Here we address this important open question applying experimental evolution and integrative systems approaches on parasites adapting to in vitro culture. Phenotypic analyses of parasites from early and late stages of culture adaptation revealed an important fitness tradeoff, with selection for accelerated growth in promastigote culture (fitness gain) impairing infectivity (fitness costs). Comparative genomics, transcriptomics and proteomics analyses revealed a complex regulatory network associated with parasite fitness gain, with genome instability causing highly reproducible, gene dosage-independent and -dependent changes. Reduction of flagellar transcripts and increase in coding and non-coding RNAs implicated in ribosomal biogenesis and protein translation were not correlated to dosage changes of the corresponding genes, revealing a gene dosage-independent, post-transcriptional mechanism of regulation. In contrast, abundance of gene products implicated in post-transcriptional regulation itself correlated to corresponding gene dosage changes. Thus, RNA abundance during parasite adaptation is controled by direct and indirect gene dosage changes. We correlated differential expression of small nucleolar RNAs (snoRNAs) with changes in rRNA modification, providing first evidence that Leishmania fitness gain in culture may be controlled by post-transcriptional and epitranscriptomic regulation. Our findings propose a novel model for Leishmania fitness gain in culture, where differential regulation of mRNA stability and the generation of modified ribosomes may potentially filter deleterious from beneficial gene dosage effects and provide proteomic robustness to genetically heterogenous, adapting parasite populations. This model challenges the current, genome-centric approach to Leishmania epidemiology and identifies the Leishmania transcriptome and non-coding small RNome as potential novel sources for the discovery of biomarkers that may be associated with parasite phenotypic adaptation in clinical settings.

Genome instability plays a central yet poorly understood role in human disease. Gene amplifications and deletions drive cancer development, microbial infection and therapeutic failure. The molecular mechanisms that harness the deleterious effects of genome instability to generate beneficial phenotypes in pathogenic systems are unknown. Here we study this important open question in the protozoan parasite Leishmania that causes devastating human diseases termed leishmaniases. Leishmania parasites lack transcriptional control and instead exploit genome instability to adapt to their host environment. Analyzing in vitro adaptation of hamster-derived parasites via gene copy number (genomic level) and gene expression changes (transcriptomic and proteomic levels), we show that these parasites likely exploit small nucleolar RNAs (snoRNAs) to mitigate toxic effects of genome instability by post-transcriptional regulation and the establishment of modified ribosomes. Our findings propose non-coding RNAs as potential novel biomarkers with diagnostic and prognostic value that may be linked to changes in parasite tissue tropism or drug susceptibility. This novel insight into Leishmania adaptation will be likely applicable to other fast evolving eukaryotic systems with unstable genomes, such as fungi or cancer cells.

Targeted Deletion of Centrin in Leishmania braziliensis Using CRISPR-Cas9-Based Editing

Leishmania braziliensis is the main causative agent of Tegumentary Leishmaniasis in the Americas. However, difficulties related to genome manipulation, experimental infection, and parasite growth have so far limited studies with this species. CRISPR-Cas9-based technology has made genome editing more accessible, and here we have successfully employed the LeishGEdit approach to attenuate L. braziliensis. We generated a transgenic cell line expressing Cas9 and T7 RNA polymerase, which was employed for the targeted deletion of centrin, a calcium-binding cytoskeletal protein involved in the centrosome duplication in eukaryotes. Centrin-deficient Leishmania exhibit growth arrest at the amastigote stage. Whole-genome sequencing of centrin-deficient L. braziliensis (LbCen^−/−) did not indicate the presence of off-target mutations. In vitro, the growth rates of LbCen^−/− and wild-type promastigotes were similar, but axenic and intracellular LbCen^−/− amastigotes showed a multinucleated phenotype with impaired survival following macrophage infection. Upon inoculation into BALB/c mice, LbCen^−/− were detected at an early time point but failed to induce lesion formation, contrary to control animals, infected with wild-type L. braziliensis. A significantly lower parasite burden was also observed in mice inoculated with LbCen^−/−, differently from control mice. Given that centrin-deficient Leishmania sp. have become candidates for vaccine development, we propose that LbCen^−/− can be further explored for the purposes of immunoprophylaxis against American Tegumentary Leishmaniasis.