Ecological divergence and hybridization of Neotropical Leishmania parasites

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.

The Andes as a semi-permeable geographical barrier: Genetic connectivity between structured populations in a widespread spider

Geographical barriers like mountain ranges impede genetic exchange among populations, promoting diversification. The effectiveness of these barriers in limiting gene flow varies between lineages due to each species' dispersal modes and capacities. Our understanding of how the Andes orogeny contributes to species diversification comes from well-studied vertebrates and a few arthropods and plants, neglecting organisms unable to fly or walk long distances. Some arachnids, such as Gasteracantha cancriformis, have been hypothesized to disperse long distances via ballooning (i.e. using their silk to interact with the wind). Yet, we do not know how the environment and geography shape its genetic diversity. Therefore, we tested whether the Andes contributed to the diversification of G. cancriformis acting as an absolute or semi-permeable barrier to genetic connectivity between populations of this spider at opposite sides of the mountain range. We sampled thousands of loci across the distribution of the species and implemented population genetics, phylogenetic, and landscape genetic analyses. We identified two genetically distinct groups structured by the Central Andes, and a third less structured group in the Northern Andes that shares ancestry with the previous two. This structure is largely explained by the altitude along the Andes, which decreases in some regions, possibly facilitating cross-Andean dispersal and gene flow. Our findings support that altitude in the Andes plays a major role in structuring populations in South America, but the strength of this barrier can be overcome by organisms with long-distance dispersal modes together with altitudinal depressions.

Rapid detection, quantification and speciation of Leishmania using real-time PCR and DNA sequencing at the rRNA Internal Transcribed Spacer 2

The ability to discriminate infection between closely related Leishmania species within the Viannia species complex, specifically L. braziliensis, L. guyanensis and L. panamensis is critical to inform the clinical diagnosis and determine the most efficacious treatment modality. We designed a nested primer set targeting the rRNA Internal Transcribed Spacer 2 (ITS2), located on Chromosome 27, to distinguish among all human infective Leishmania species. Separate nested and single primer pairs were developed for conventional and quantitative PCR approaches respectively. Species-specific single nucleotide polymorphisms and indels located within the PCR products were identified by Sanger sequencing. This single locus approach provides a sensitive and specific platform to identify the species of Leishmania causing infection.

High-throughput analysis of the Trypanosoma cruzi minicirculome (mcDNA) unveils structural variation and functional diversity

Trypanosoma cruzi causes Chagas disease and has a unique extranuclear genome enclosed in a structure called the kinetoplast, which contains circular genomes known as maxi- and minicircles. While the structure and function of maxicircles are well-understood, many aspects of minicircles remain to be discovered. Here, we performed a high-throughput analysis of the minicirculome (mcDNA) in 50 clones isolated from Colombia's diverse T. cruzi I populations. Results indicate that mcDNA comprises four diverse subpopulations with different structures, lengths, and numbers of interspersed semi-conserved (previously termed ultra-conserved regions mHCV) and hypervariable (mHVPs) regions. Analysis of mcDNA ancestry and inter-clone differentiation indicates the interbreeding of minicircle sequence classes is placed along diverse strains and hosts. These results support evidence of the multiclonal dynamics and random bi-parental segregation. Finally, we disclosed the guide RNA repertoire encoded by mcDNA at a clonal scale, and several attributes of its abundance and function are discussed.

Comparative genomics of Leishmania donovani progeny from genetic crosses in two sand fly species and impact on the diversity of diagnostic and vaccine candidates

Sand fly transmitted Leishmania species are responsible for severe, wide ranging, visceral and cutaneous leishmaniases. Genetic exchange can occur among natural Leishmania populations and hybrids can now be produced experimentally, with limitations. Feeding Phlebotomus orientalis or Phlebotomus argentipes on two strains of Leishmania donovani yielded hybrid progeny, selected using double drug resistance and fluorescence markers. Fluorescence activated cell sorting of cultured clones derived from these hybrids indicated diploid progeny. Multilocus sequence typing of the clones showed hybridisation and nuclear heterozygosity, although with inheritance of single haplotypes in a kinetoplastid target. Comparative genomics showed diversity of clonal progeny between single chromosomes, and extraordinary heterozygosity across all 36 chromosomes. Diversity between progeny was seen for the HASPB antigen, which has been noted previously as having implications for design of a therapeutic vaccine. Genomic diversity seen among Leishmania strains and hybrid progeny is of great importance in understanding the epidemiology and control of leishmaniasis. As an outcome of this study we strongly recommend that wider biological archives of different Leishmania species from endemic regions should be established and made available for comparative genomics. However, in parallel, performance of genetic crosses and genomic comparisons should give fundamental insight into the specificity, diversity and limitations of candidate diagnostics, vaccines and drugs, for targeted control of leishmaniasis.

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.

Leishmania genetic exchange is mediated by IgM natural antibodies

Host factors that mediate Leishmania genetic exchange are not well defined. Here we demonstrate that natural IgM (IgMn)1-4 antibodies mediate parasite genetic exchange by inducing the transient formation of a spherical parasite clump that promotes parasite fusion and hybrid formation. We establish that IgMn from Leishmania-free animals binds to the surface of Leishmania parasites to induce significant changes in the expression of parasite transcripts and proteins. Leishmania binding to IgMn is partially lost after glycosidase treatment, although parasite surface phosphoglycans, including lipophosphoglycan, are not required for IgMn-induced parasite clumping. Notably, the transient formation of parasite clumps is essential for Leishmania hybridization in vitro. In vivo, we observed a 12-fold increase in hybrid formation in sand flies provided a second blood meal containing IgMn compared with controls. Furthermore, the generation of recombinant progeny from mating hybrids and parental lines were only observed in sand flies provided with IgMn. Both in vitro and in vivo IgM-induced Leishmania crosses resulted in full genome hybrids that show equal patterns of biparental contribution. Leishmania co-option of a host natural antibody to facilitate mating in the insect vector establishes a new paradigm of parasite-host-vector interdependence that contributes to parasite diversity and fitness by promoting genetic exchange.

Unveiling drug-tolerant and persister-like cells in Leishmania braziliensis lines derived from patients with cutaneous leishmaniasis

Introduction

Resistance against anti-Leishmania drugs (DR) has been studied for years, giving important insights into long-term adaptations of these parasites to drugs, through genetic modifications. However, microorganisms can also survive lethal drug exposure by entering into temporary quiescence, a phenomenon called drug tolerance (DT), which is rather unexplored in Leishmania.

Methods

We studied a panel of nine Leishmania braziliensis strains highly susceptible to potassium antimonyl tartrate (PAT), exposed promastigotes to lethal PAT pressure, and compared several cellular and molecular parameters distinguishing DT from DR.

Results and discussion

We demonstrated in vitro that a variable proportion of cells remained viable, showing all the criteria of DT and not of DR: i) signatures of quiescence, under drug pressure: reduced proliferation and significant decrease of rDNA transcription; ii) reversibility of the phenotype: return to low IC50 after removal of drug pressure; and iii) absence of significant genetic differences between exposed and unexposed lineages of each strain and absence of reported markers of DR. We found different levels of quiescence and DT among the different L. braziliensis strains. We provide here a new in-vitro model of drug-induced quiescence and DT in Leishmania. Research should be extended in vivo, but the current model could be further exploited to support R&D, for instance, to guide the screening of compounds to overcome the quiescence resilience of the parasite, thereby improving the therapy of leishmaniasis.

First evidence of experimental genetic hybridization between cutaneous and visceral strains of Leishmania donovani within its natural vector Phlebotomus argentipes

Leishmaniasis is a neglected tropical disease caused by protozoan parasites of genus Leishmania, and transmitted by different species of Phlebotomine sand flies. More than 20 species of Leishmania are known to cause disease in humans and other animals. Leishmania donovani species complex is known to have a vast diversity of clinical manifestations in humans, but underlying mechanisms for such diversity are yet unknown. Long believed to be strictly asexual, Leishmania have been shown to undergo a cryptic sexual cycle inside its sandfly vector. Natural populations of hybrid parasites have been associated with the rise of atypical clinical outcomes in the Indian subcontinent (ISC). However, formal demonstration of genetic crossing in the major endemic sandfly species in the ISC remain unexplored. Here, we investigated the ability of two distinct variants of L. donovani associated with strikingly different forms of the disease to undergo genetic exchange inside its natural vector, Phlebotomus argentipes. Clinical isolates of L. donovani either from a Sri Lankan cutaneous leishmaniasis (CL) patient or an Indian visceral leishmaniasis (VL) patient were genetically engineered to express different fluorescent proteins and drug-resistance markers and subsequently used as parental strains in experimental sandfly co-infection. After 8 days of infection, sand flies were dissected and midgut promastigotes were transferred into double drug-selective media. Two double drug-resistant, dual fluorescent hybrid cell lines were recovered, which after cloning and whole genome sequencing, were shown to be full genomic hybrids. This study provides the first evidence of L. donovani hybridization within its natural vector Ph. argentipes.

Deciphering anti-infectious compounds from Peruvian medicinal Cordoncillos extract library through multiplexed assays and chemical profiling

High prevalence of parasitic or bacterial infectious diseases in some world areas is due to multiple reasons, including a lack of an appropriate health policy, challenging logistics and poverty. The support to research and development of new medicines to fight infectious diseases is one of the sustainable development goals promoted by World Health Organization (WHO). In this sense, the traditional medicinal knowledge substantiated by ethnopharmacology is a valuable starting point for drug discovery. This work aims at the scientific validation of the traditional use of Piper species ("Cordoncillos") as firsthand anti-infectious medicines. For this purpose, we adapted a computational statistical model to correlate the LCMS chemical profiles of 54 extracts from 19 Piper species to their corresponding anti-infectious assay results based on 37 microbial or parasites strains. We mainly identified two groups of bioactive compounds (called features as they are considered at the analytical level and are not formally isolated). Group 1 is composed of 11 features being highly correlated to an inhibiting activity on 21 bacteria (principally Gram-positive strains), one fungus (C. albicans), and one parasite (Trypanosoma brucei gambiense). The group 2 is composed of 9 features having a clear selectivity on Leishmania (all strains, both axenic and intramacrophagic). Bioactive features in group 1 were identified principally in the extracts of Piper strigosum and P. xanthostachyum. In group 2, bioactive features were distributed in the extracts of 14 Piper species. This multiplexed approach provided a broad picture of the metabolome as well as a map of compounds putatively associated to bioactivity. To our knowledge, the implementation of this type of metabolomics tools aimed at identifying bioactive compounds has not been used so far.

Genome Analysis of Triploid Hybrid Leishmania Parasite from the Neotropics

We discovered a hybrid Leishmania parasite in Costa Rica that is genetically similar to hybrids from Panama. Genome analyses demonstrated the hybrid is triploid and identified L. braziliensis and L. guyanensis-related strains as parents. Our findings highlight the existence of poorly sampled Leishmania (Viannia) variants infectious to humans.

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?

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.

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.

Genome diversity of Leishmania aethiopica

Leishmania aethiopica is a zoonotic Old World parasite transmitted by Phlebotomine sand flies and causing cutaneous leishmaniasis in Ethiopia and Kenya. Despite a range of clinical manifestations and a high prevalence of treatment failure, L. aethiopica is one of the most neglected species of the Leishmania genus in terms of scientific attention. Here, we explored the genome diversity of L. aethiopica by analyzing the genomes of twenty isolates from Ethiopia. Phylogenomic analyses identified two strains as interspecific hybrids involving L. aethiopica as one parent and L. donovani and L. tropica respectively as the other parent. High levels of genome-wide heterozygosity suggest that these two hybrids are equivalent to F1 progeny that propagated mitotically since the initial hybridization event. Analyses of allelic read depths further revealed that the L. aethiopica - L. tropica hybrid was diploid and the L. aethiopica - L. donovani hybrid was triploid, as has been described for other interspecific Leishmania hybrids. When focusing on L. aethiopica, we show that this species is genetically highly diverse and consists of both asexually evolving strains and groups of recombining parasites. A remarkable observation is that some L. aethiopica strains showed an extensive loss of heterozygosity across large regions of the nuclear genome, which likely arose from gene conversion/mitotic recombination. Hence, our prospection of L. aethiopica genomics revealed new insights into the genomic consequences of both meiotic and mitotic recombination in Leishmania.

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.

An overview of the trypanosomatid (Kinetoplastida: Trypanosomatidae) parasites infecting several mammal species in Colombia

BACKGROUND

Trypanosomatids are among the most critical parasites for public health due to their impact on human, animal, and plant health. Diseases associated with these pathogens manifest mainly in poor and vulnerable populations, where social, environmental, and biological factors modulate the case incidence and geographical distribution.

METHODS

We used Sanger and amplicon-based next-generation sequencing (NGS) in samples from different mammals to identify trypanosomatid infections in several departments in Colombia. A total of 174 DNA samples (18 humans, 83 dogs, and 73 wild mammals) were analyzed by conventional PCR using a fragment of the heat shock protein 70 (Hsp70) gene and Sanger sequenced the positive samples. Twenty-seven samples were sent for amplicon-based NGS using the same gene fragment. Data obtained were used to perform diversity analyses.

RESULTS

One hundred and thirteen samples were positive for PCR by Hsp70 fragment; these corresponded to 22.1% Leishmania spp., 18.6% L. amazonensis, 9.7% L. braziliensis, 14.2% L. infantum, 8% L. panamensis, and 27.4% Trypanosoma cruzi. Comparison of the identified species by the two sequencing technologies used resulted in 97% concordance. Alpha and beta diversity indices were significant, mainly for dogs; there was an interesting index of coinfection events in the analyzed samples: different Leishmania species and the simultaneous presence of T. cruzi and even T. rangeli in one of the samples analyzed. Moreover, a low presence of L. braziliensis was observed in samples from wild mammals. Interestingly, to our knowledge, this is the first report of Leishmania detection in Hydrochaeris hydrochaeris (capybara) in Colombia.

CONCLUSIONS

The Hsp70 fragment used in this study is an optimal molecular marker for trypanosomatid identification in many hosts and allows the identification of different species in the same sample when amplicon-based sequencing is used. However, the use of this fragment for molecular diagnosis through conventional PCR should be carefully interpreted because of this same capacity to identify several parasites. This point is of pivotal importance in highly endemic countries across South America because of the co-circulation of different genera from the Trypanosomatidae family. The findings show an interesting starting point for One Health approaches in which coevolution and vector-host interactions can be studied.

Identification of a conserved maxicircle and unique minicircles as part of the mitochondrial genome of Leishmania martiniquensis strain PCM3 in Thailand

BACKGROUND

The mitochondrial DNA of trypanosomatids, including Leishmania, is known as kinetoplast DNAs (kDNAs). The kDNAs form networks of hundreds of DNA circles that are evidently interlocked and require complex RNA editing. Previous studies showed that kDNA played a role in drug resistance, adaptation, and survival of Leishmania. Leishmania martiniquensis is one of the most frequently observed species in Thailand, and its kDNAs have not been illustrated.

METHODS

This study aimed to extract the kDNA sequences from Illumina short-read and PacBio long-read whole-genome sequence data of L. martiniquensis strain PCM3 priorly isolated from the southern province of Thailand. A circular maxicircle DNA was reconstructed by de novo assembly using the SPAdes program, while the minicircle sequences were retrieved and assembled by the rKOMIC tool. The kDNA contigs were confirmed by blasting to the NCBI database, followed by comparative genomic and phylogenetic analysis.

RESULTS

We successfully constructed the complete circular sequence of the maxicircle (19,008 bp) and 214 classes of the minicircles from L. martiniquensis strain PCM3. The genome comparison and annotation showed that the maxicircle structure of L. martiniquensis strain PCM3 was similar to those of L. enriettii strain LEM3045 (84.29%), L. arabica strain LEM1108 (82.79%), and L. tarentolae (79.2%). Phylogenetic analysis also showed unique evolution of the minicircles of L. martiniquensis strain PCM3 from other examined Leishmania species.

CONCLUSIONS

This was the first report of the complete maxicircle and 214 minicircles of L. martiniquensis strain PCM3 using integrated whole-genome sequencing data. The information will be helpful for further improvement of diagnosis methods and monitoring genetic diversity changes of this parasite.

Deep kinetoplast genome analyses result in a novel molecular assay for detecting Trypanosoma brucei gambiense-specific minicircles

The World Health Organization targeted Trypanosoma brucei gambiense (Tbg) human African trypanosomiasis for elimination of transmission by 2030. Sensitive molecular markers that specifically detect Tbg type 1 (Tbg1) parasites will be important tools to assist in reaching this goal. We aim at improving molecular diagnosis of Tbg1 infections by targeting the abundant mitochondrial minicircles within the kinetoplast of these parasites. Using Next-Generation Sequencing of total cellular DNA extracts, we assembled and annotated the kinetoplast genome and investigated minicircle sequence diversity in 38 animal- and human-infective trypanosome strains. Computational analyses recognized a total of 241 Minicircle Sequence Classes as Tbg1-specific, of which three were shared by the 18 studied Tbg1 strains. We developed a minicircle-based assay that is applicable on animals and as specific as the TgsGP-based assay, the current golden standard for molecular detection of Tbg1. The median copy number of the targeted minicircle was equal to eight, suggesting our minicircle-based assay may be used for the sensitive detection of Tbg1 parasites. Annotation of the targeted minicircle sequence indicated that it encodes genes essential for the survival of the parasite and will thus likely be preserved in natural Tbg1 populations, the latter ensuring the reliability of our novel diagnostic assay.

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.

Organization of minicircle cassettes and guide RNA genes in Trypanosoma brucei

Mitochondrial DNA of protists of order Kinetoplastida comprises thousands of interlinked circular molecules arranged in a network. There are two types of molecules called minicircles and maxicircles. Minicircles encode guide RNA (gRNA) genes whose transcripts mediate post-transcriptional editing of maxicircle encoded genes. Minicircles are diverse. The human sleeping sickness parasite Trypanosoma brucei has one of the most diverse sets of minicircle classes of all studied trypanosomatids with hundreds of different classes, each encoding one to four genes mainly within cassettes framed by 18 bp inverted repeats. A third of cassettes have no identifiable gRNA genes even though their sequence structures are similar to cassettes with identifiable genes. Only recently have almost all minicircle classes for some subspecies and isolates of T. brucei been sequenced and annotated with corresponding verification of gRNA expression by small-RNA transcriptome data. These data sets provide a rich resource for understanding the structure of minicircle classes, cassettes and gRNA genes and their transcription. Here, we provide a statistical description of the functionality, expression status, structure and sequence of gRNA genes in a differentiation-competent, laboratory-adapted strain of T. brucei We obtain a clearer definition of what is a gRNA gene. Our analysis supports the idea that many, if not all, cassettes without an identifiable gRNA gene contain decaying remnants of once functional gRNA genes. Finally, we report several new, unexplained discoveries such as the association between cassette position on the minicircle and gene expression and functionality, and the association between gene initiation sequence and anchor position.

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.

Genetic diversity and population structure of Leishmania (Viannia) braziliensis in the Peruvian jungle

Background

Human cutaneous leishmaniasis caused by Leishmania (Viannia) braziliensis is highly prevalent in the Peruvian jungle, where it affects military forces deployed to fight against drug trafficking and civilian people that migrate from the highland to the lowland jungle for economic activities such as mining, agriculture, construction, and chestnut harvest.

We explored the genetic diversity and population structure of 124 L. (V.) braziliensis isolates collected from the highland (Junín, Cusco, and Ayacucho) and lowland Peruvian jungle (Loreto, Ucayali, and Madre de Dios). All samples were genotyped using Multilocus Microsatellite Typing (MLMT) of ten highly polymorphic markers.

Principal findings

High polymorphism and genetic diversity were found in Peruvian isolates of L. (V.) braziliensis. Most markers are not in Hardy-Weinberg equilibrium; this deviation is most likely caused by local inbreeding, as shown by the positive F_IS values. Linkage Disequilibrium in subpopulations was not strong, suggesting the reproduction was not strictly clonal. Likewise, for the first time, two genetic clusters of this parasite were determined, distributed in both areas of the Peruvian jungle, which suggested a possible recent colonization event of the highland jungle from the lowland jungle.

Conclusions

L. (V.) braziliensis exhibits considerable genetic diversity with two different clusters in the Peruvian jungle. Migration analysis suggested a colonization event between geographical areas of distribution. Although no human migration was observed at the time of sampling, earlier displacement of humans, reservoirs, or vectors could have been responsible for the parasite spread in both regions.

L. (V.) braziliensis is widespread in the Peruvian jungle region. In this region, the departments of Cusco and Madre de Dios account for a large number of patients that get infected while working in the virgin forest. For the first time, we described ample genetic diversity among Peruvian L. (V.) braziliensis isolates with new alleles that were not previously reported in South America. In addition, two different genetic clusters or subpopulations of L. (V.) braziliensis in the Amazon Rainforest of Peru were described. This finding reveals the important distribution of parasite populations and suggests a possible colonization event between ecoregions of the highland and lowland Peruvian jungle independently of recent human migration.

Experimental Hybridization in Leishmania: Tools for the Study of Genetic Exchange

Despite major advances over the last decade in our understanding of Leishmania reproductive strategies, the sexual cycle in Leishmania has defied direct observation and remains poorly investigated due to experimental constraints. Here, we summarize the findings and conclusions drawn from genetic analysis of experimental hybrids generated in sand flies and highlight the recent advances in generating hybrids in vitro. The ability to hybridize between culture forms of different species and strains of Leishmania should invite more intensive investigation of the mechanisms underlying genetic exchange and provide a rich source of recombinant parasites for future genetic analyses.

Microevolution of Trypanosoma cruzi reveals hybridization and clonal mechanisms driving rapid genome diversification

Protozoa and fungi are known to have extraordinarily diverse mechanisms of genetic exchange. However, the presence and epidemiological relevance of genetic exchange in Trypanosoma cruzi, the agent of Chagas disease, has been controversial and debated for many years. Field studies have identified both predominantly clonal and sexually recombining natural populations. Two of six natural T. cruzi lineages (TcV and TcVI) show hybrid mosaicism, using analysis of single-gene locus markers. The formation of hybrid strains in vitro has been achieved and this provides a framework to study the mechanisms and adaptive significance of genetic exchange. Using whole genome sequencing of a set of experimental hybrids strains, we have confirmed that hybrid formation initially results in tetraploid parasites. The hybrid progeny showed novel mutations that were not attributable to either (diploid) parent showing an increase in amino acid changes. In long-term culture, up to 800 generations, there was a variable but gradual erosion of progeny genomes towards triploidy, yet retention of elevated copy number was observed at several core housekeeping loci. Our findings indicate hybrid formation by fusion of diploid T. cruzi, followed by sporadic genome erosion, but with substantial potential for adaptive evolution, as has been described as a genetic feature of other organisms, such as some fungi.

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.

Signatures of hybridization in Trypanosoma brucei

Genetic exchange among disease-causing micro-organisms can generate progeny that combine different pathogenic traits. Though sexual reproduction has been described in trypanosomes, its impact on the epidemiology of Human African Trypanosomiasis (HAT) remains controversial. However, human infective and non-human infective strains of Trypanosoma brucei circulate in the same transmission cycles in HAT endemic areas in subsaharan Africa, providing the opportunity for mating during the developmental cycle in the tsetse fly vector. Here we investigated inheritance among progeny from a laboratory cross of T. brucei and then applied these insights to genomic analysis of field-collected isolates to identify signatures of past genetic exchange. Genomes of two parental and four hybrid progeny clones with a range of DNA contents were assembled and analysed by k-mer and single nucleotide polymorphism (SNP) frequencies to determine heterozygosity and chromosomal inheritance. Variant surface glycoprotein (VSG) genes and kinetoplast (mitochondrial) DNA maxi- and minicircles were extracted from each genome to examine how each of these components was inherited in the hybrid progeny. The same bioinformatic approaches were applied to an additional 37 genomes representing the diversity of T. brucei in subsaharan Africa and T. evansi. SNP analysis provided evidence of crossover events affecting all 11 pairs of megabase chromosomes and demonstrated that polyploid hybrids were formed post-meiotically and not by fusion of the parental diploid cells. VSGs and kinetoplast DNA minicircles were inherited biparentally, with approximately equal numbers from each parent, whereas maxicircles were inherited uniparentally. Extrapolation of these findings to field isolates allowed us to distinguish clonal descent from hybridization by comparing maxicircle genotype to VSG and minicircle repertoires. Discordance between maxicircle genotype and VSG and minicircle repertoires indicated inter-lineage hybridization. Significantly, some of the hybridization events we identified involved human infective and non-human infective trypanosomes circulating in the same geographic areas.

Sexual reproduction allows genes from different individuals to be mixed up in the offspring. This is particularly important for disease-causing microbes, because new combinations of harmful traits can arise, potentially leading to more severe outbreaks of disease. Tsetse-transmitted trypanosomes are single-celled parasites that cause severe human and livestock diseases in tropical Africa. During their developmental cycle in the tsetse fly, trypanosomes can mate and produce hybrid trypanosomes, which have one set of chromosomes from each parent. But polyploid hybrids, with more than one set of chromosomes from one or both parents, are often observed too. Here we have investigated how these polyploid hybrids are formed by comparing the genomes of hybrid progeny with those of their parents. Analysis of the large, paired chromosomes of both diploid and polyploid hybrids showed frequent crossovers, which are the hallmark of meiosis, the special form of division that produces haploid gametes. This indicates that the polyploids were formed after meiosis rather than by fusion of the parental diploid cells. We also investigated the inheritance of two other features of trypanosomes: the large family of variant surface glycoprotein (VSG) genes, and the mitochondrial (kinetoplast) DNA. Hybrid clones had inherited about half the VSG genes from each parent, and also showed biparental inheritance of one component of the kinetoplast DNA, the minicircles. We assessed the relatedness of field-collected trypanosomes by comparing their VSG and minicircle repertoires, together with maxicircle genotype. While most isolates shared few VSGs or minicircles, a group of mostly human-infective strains from Uganda had a large proportion of their repertoires in common. Most of these trypanosomes were probably related by clonal descent, but we also identified that some were hybrids by the mismatch between their maxicircle genotype and their VSG and minicircle repertoires. These signals of hybridization were also detected in some of the other field-collected isolates, suggesting that genetic exchange is widespread in nature. Significantly, the hybridization events involved human infective and non-human infective trypanosomes circulating in the same geographic areas, providing a mechanism for the generation of new, potentially more pathogenic, trypanosome strains causing human disease.

Culture of Cutaneous Leishmania from Skin Biopsy Specimens

This protocol describes the culture of Leishmania parasites from skin biopsy samples of patients with cutaneous lesions. The use of antibiotics to prevent bacterial contamination of these cultures increases the ability of researchers to collect isolates for various research purposes, including genetic analysis and in vitro and in vivo experiments. © Published 2022. This article is a U.S. Government work and is in the public domain in the USA. Basic Protocol: Culture of Leishmania from skin biopsy specimens.

Impact of Genetic Diversity and Genome Plasticity of Leishmania spp. in Treatment and the Search for Novel Chemotherapeutic Targets

Leishmaniasis is one of the major public health concerns in Latin America, Africa, Asia, and Europe. The absence of vaccines for human use and the lack of effective vector control programs make chemotherapy the main strategy to control all forms of the disease. However, the high toxicity of available drugs, limited choice of therapeutic agents, and occurrence of drug-resistant parasite strains are the main challenges related to chemotherapy. Currently, only a small number of drugs are available for leishmaniasis treatment, including pentavalent antimonials (Sb^V), amphotericin B and its formulations, miltefosine, paromomycin sulphate, and pentamidine isethionate. In addition to drug toxicity, therapeutic failure of leishmaniasis is a serious concern. The occurrence of drug-resistant parasites is one of the causes of therapeutic failure and is closely related to the diversity of parasites in this genus. Owing to the enormous plasticity of the genome, resistance can occur by altering different metabolic pathways, demonstrating that resistance mechanisms are multifactorial and extremely complex. Genetic variability and genome plasticity cause not only the available drugs to have limitations, but also make the search for new drugs challenging. Here, we examined the biological characteristics of parasites that hinder drug discovery.

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.

Morphological and genomic characterisation of the Schistosoma hybrid infecting humans in Europe reveals admixture between Schistosoma haematobium and Schistosoma bovis

Schistosomes cause schistosomiasis, the world's second most important parasitic disease after malaria in terms of public health and social-economic impacts. A peculiar feature of these dioecious parasites is their ability to produce viable and fertile hybrid offspring. Originally only present in the tropics, schistosomiasis is now also endemic in southern Europe. Based on the analysis of two genetic markers the European schistosomes had previously been identified as hybrids between the livestock- and the human-infective species Schistosoma bovis and Schistosoma haematobium, respectively. Here, using PacBio long-read sequencing technology we performed genome assembly improvement and annotation of S. bovis, one of the parental species for which no satisfactory genome assembly was available. We then describe the whole genome introgression levels of the hybrid schistosomes, their morphometric parameters (eggs and adult worms) and their compatibility with two European snail strains used as vectors (Bulinus truncatus and Planorbarius metidjensis). Schistosome-snail compatibility is a key parameter for the parasites life cycle progression, and thus the capability of the parasite to establish in a given area. Our results show that this Schistosoma hybrid is strongly introgressed genetically, composed of 77% S. haematobium and 23% S. bovis origin. This genomic admixture suggests an ancient hybridization event and subsequent backcrosses with the human-specific species, S. haematobium, before its introduction in Corsica. We also show that egg morphology (commonly used as a species diagnostic) does not allow for accurate hybrid identification while genetic tests do.

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.

Overcoming roadblocks in the development of vaccines for leishmaniasis

INTRODUCTION

The leishmaniases represent a group of parasitic diseases caused by infection with one of several species of Leishmania parasites. Disease presentation varies because of differences in parasite and host genetics and may be influenced by additional factors such as host nutritional status or co-infection. Studies in experimental models of Leishmania infection, vaccination of companion animals and human epidemiological data suggest that many forms of leishmaniasis could be prevented by vaccination, but no vaccines are currently available for human use.

AREAS COVERED

We describe some of the existing roadblocks to the development and implementation of an effective leishmaniasis vaccine, based on a review of recent literature found on PubMed, BioRxiv and MedRxiv. In addition to discussing scientific unknowns that hinder vaccine candidate identification and selection, we explore gaps in knowledge regarding the commercial and public health value propositions underpinning vaccine development and provide a route map for future research and advocacy.

EXPERT OPINION

Despite significant progress, leishmaniasis vaccine development remains hindered by significant gaps in understanding that span the vaccine development pipeline. Increased coordination and adoption of a more holistic view to vaccine development will be required to ensure more rapid progress in the years ahead.

Development of an Amplicon-Based Next-Generation Sequencing Protocol to Identify Leishmania Species and Other Trypanosomatids in Leishmaniasis Endemic Areas

Trypanosomatid infections are an important public health threat affecting many low-income countries across the tropics, particularly in the Americas. Trypanosomatids can infect many vertebrate, invertebrate, and plant species and play an important role as human pathogens. Among these clinically relevant pathogens are species from the genera Leishmania and Trypanosoma. Mixed trypanosomatid infections remain a largely unexplored phenomenon. Herein, we describe the application of an amplicon-based next-generation sequencing (NGS) assay to detect and identify trypanosomatid species in mammalian reservoirs, human patients, and sand fly vectors throughout regions of Leishmania endemicity. Sixty-five samples from different departments of Colombia, including two samples from Venezuela, were analyzed: 49 samples from cutaneous leishmaniasis (CL) patients, 8 from sand flies, 2 from domestic reservoirs (Canis familiaris), and 6 from wild reservoirs (Phyllostomus hastatus). DNA from each sample served to identify the presence of trypanosomatids through conventional PCR using heat shock protein 70 (HSP70) gene as the target. PCR products underwent sequencing by Sanger sequencing and NGS, and trypanosomatid species were identified by using BLASTn against a reference database built from trypanosomatid-derived HSP70 sequences. The alpha and beta diversity indexes of amplicon sequence variants were calculated for each group. The results revealed the presence of mixed infections with more than two Leishmania species in 34% of CL samples analyzed. Trypanosoma cruzi was identified in samples from wild reservoirs, as well as in sand fly vectors. Coinfection events with three different Leishmania species were identified in domestic reservoirs. These findings depose the traditional paradigm of leishmaniasis as being a single-species-driven infection and redraw the choreography of host-pathogen interaction in the context of multiparasitism. Further research is needed to decipher how coinfections may influence disease progression. This knowledge is key to developing an integrated approach for diagnosis and treatment.

IMPORTANCE Traditionally, there has been a frequent, yet incorrect assumption that phlebotomine vectors, animal reservoirs, and human hosts are susceptible to Leishmania infection by a single parasite species. However, current evidence supports that these new vector-parasite-reservoir associations lend vectors and reservoirs greater permissiveness to certain Leishmania species, thus promoting the appearance of coinfection events, particularly in disease-endemic regions. The present study describes the application of an amplicon-based next-generation sequencing (NGS) assay to detect and identify trypanosomatid species in mammalian reservoirs, human patients, and sand fly vectors from regions of endemicity for leishmaniasis. This changes our understanding of the clinical course of leishmaniasis in areas of endemicity.

rKOMICS: an R package for processing mitochondrial minicircle assemblies in population-scale genome projects

Background

The advent of population-scale genome projects has revolutionized our biological understanding of parasitic protozoa. However, while hundreds to thousands of nuclear genomes of parasitic protozoa have been generated and analyzed, information about the diversity, structure and evolution of their mitochondrial genomes remains fragmentary, mainly because of their extraordinary complexity. Indeed, unicellular flagellates of the order Kinetoplastida contain structurally the most complex mitochondrial genome of all eukaryotes, organized as a giant network of homogeneous maxicircles and heterogeneous minicircles. We recently developed KOMICS, an analysis toolkit that automates the assembly and circularization of the mitochondrial genomes of Kinetoplastid parasites. While this tool overcomes the limitation of extracting mitochondrial assemblies from Next-Generation Sequencing datasets, interpreting and visualizing the genetic (dis)similarity within and between samples remains a time-consuming process.

Results

Here, we present a new analysis toolkit—rKOMICS—to streamline the analyses of minicircle sequence diversity in population-scale genome projects. rKOMICS is a user-friendly R package that has simple installation requirements and that is applicable to all 27 trypanosomatid genera. Once minicircle sequence alignments are generated, rKOMICS allows to examine, summarize and visualize minicircle sequence diversity within and between samples through the analyses of minicircle sequence clusters. We showcase the functionalities of the (r)KOMICS tool suite using a whole-genome sequencing dataset from a recently published study on the history of diversification of the Leishmania braziliensis species complex in Peru. Analyses of population diversity and structure highlighted differences in minicircle sequence richness and composition between Leishmania subspecies, and between subpopulations within subspecies.

Conclusion

The rKOMICS package establishes a critical framework to manipulate, explore and extract biologically relevant information from mitochondrial minicircle assemblies in tens to hundreds of samples simultaneously and efficiently. This should facilitate research that aims to develop new molecular markers for identifying species-specific minicircles, or to study the ancestry of parasites for complementary insights into their evolutionary history.

Hybridized Zoonotic Schistosoma Infections Result in Hybridized Morbidity Profiles: A Clinical Morbidity Study amongst Co-Infected Human Populations of Senegal

Hybridization of infectious agents is a major emerging public and veterinary health concern at the interface of evolution, epidemiology, and control. Whilst evidence of the extent of hybridization amongst parasites is increasing, their impact on morbidity remains largely unknown. This may be predicted to be particularly pertinent where parasites of animals with contrasting pathogenicity viably hybridize with human parasites. Recent research has revealed that viable zoonotic hybrids between human urogenital Schistosoma haematobium with intestinal Schistosoma species of livestock, notably Schistosoma bovis, can be highly prevalent across Africa and beyond. Examining human populations in Senegal, we found increased hepatic but decreased urogenital morbidity, and reduced improvement following treatment with praziquantel, in those infected with zoonotic hybrids compared to non-hybrids. Our results have implications for effective monitoring and evaluation of control programmes, and demonstrate for the first time the potential impact of parasite hybridizations on host morbidity.

Comparative genomics of Leishmania isolates from Brazil confirms the presence of Leishmania major in the Americas

Leishmania (Leishmania) major is an important agent of cutaneous leishmaniasis, having as a vector sandflies belonging to the genus Phlebotomus. Although this species has been described as restricted to the Old World, parasites similar to L. major have been isolated from South American patients who have never travelled abroad. These parasites were named "L. major-like", and several studies have been carried out to characterise them biochemically, molecularly, and biologically. However, the phylogenetic origin of these isolates is still unknown. In the present study we characterised three L. major-like isolates, named BH49, BH121 and BH129, using comparative genomics approaches. We evaluated the presence of gene and segmental duplications/deletions and the presence of aneuploidies that could explain the differences in infectivity observed in the BH49 and BH121 isolates. All isolates presented a pattern of mosaic aneuploidy and gene copy number variation, which are common in the genus Leishmania. Virulence factors such as phosphatases and peptidases were found to have increased gene copy numbers in the infective isolate, which could explain the difference in infectivity previously observed between BH121 and BH49. Phylogenetic analyses revealed that BH49, BH121 and BH129 L. major-like grouped with L. major isolates, and suggest they were imported from the Old World in at least two independent events. We suggest that new epidemiological inquiries should also evaluate L. major infections in South America, to assess the epidemiological importance of this species in the New World.

Tegumentary leishmaniasis by Leishmania braziliensis complex in Cochabamba, Bolivia including the presence of L. braziliensis outlier: Tegumentary leishmaniasis in Cochabamba, Bolivia

Leishmaniasis is caused by protozoans of the Leishmania genus, which includes more than 20 species capable of infecting humans worldwide. In the Americas, the most widespread specie is L. braziliensis, present in 18 countries including Bolivia. The taxonomic position of the L. braziliensis complex has been a subject of controversy, complicated further by the recent identification of a particular subpopulation named L. braziliensis atypical or outlier. The aim of this study was to carry out a systematic analysis of the L. braziliensis complex in Bolivia and to describe the associated clinical characteristics. Forty-one strains were analyzed by sequencing an amplified 1245 bp fragment of the hsp70 gene, which allowed its identification as: 24 (59%) L. braziliensis, 16 (39%) L. braziliensis outlier, and one (2%) L. peruviana. In a dendrogram constructed, L. braziliensis and L. peruviana are grouped in the same cluster, whilst L. braziliensis outlier appears in a separate branch. Sequence alignment allowed the identification of five non-polymorphic nucleotide positions (288, 297, 642, 993, and 1213) that discriminate L. braziliensis and L. peruviana from L. braziliensis outlier. Moreover, nucleotide positions 51 and 561 enable L. peruviana to be discriminated from the other two taxa. A greater diversity was observed in L. braziliensis outlier than in L. braziliensis-L. peruviana. The 41 strains came from 32 patients with tegumentary leishmaniasis, among which 22 patients (69%) presented cutaneous lesions (11 caused by L. braziliensis and 11 by L. braziliensis outlier) and 10 patients (31%) mucocutaneous lesions (eight caused by L. braziliensis, one by L. braziliensis outlier, and one by L. peruviana). Nine patients (28%) simultaneously provided two isolates, each from a separate lesion, and in each case the same genotype was identified in both. Treatment failure was observed in six patients infected with L. braziliensis and one patient with L. peruviana.

The sexual side of parasitic protists

Much of the vast evolutionary landscape occupied by Eukaryotes is dominated by protists. Though parasitism has arisen in many lineages, there are three main groups of parasitic protists of relevance to human and livestock health: the Apicomplexa, including the malaria parasite Plasmodium and coccidian pathogens of livestock such as Eimeria; the excavate flagellates, encompassing a diverse range of protist pathogens including trypanosomes, Leishmania, Giardia and Trichomonas; and the Amoebozoa, including pathogenic amoebae such as Entamoeba. These three groups represent separate, deep branches of the eukaryote tree, underlining their divergent evolutionary histories. Here, I explore what is known about sex in these three main groups of parasitic protists.

Microsatellite based molecular epidemiology of Leishmania infantum from re-emerging foci of visceral leishmaniasis in Armenia and pilot risk assessment by ecological niche modeling

BACKGROUND

Visceral leishmaniasis (VL) is re-emerging in Armenia since 1999 with 167 cases recorded until 2019. The objectives of this study were (i) to determine for the first time the genetic diversity and population structure of the causative agent of VL in Armenia; (ii) to compare these genotypes with those from most endemic regions worldwide; (iii) to monitor the diversity of vectors in Armenia; (iv) to predict the distribution of the vectors and VL in time and space by ecological niche modeling.

METHODOLOGY/PRINCIPAL FINDINGS

Human samples from different parts of Armenia previously identified by ITS-1-RFLP as L. infantum were studied by Multilocus Microsatellite Typing (MLMT). These data were combined with previously typed L. infantum strains from the main global endemic regions for population structure analysis. Within the 23 Armenian L. infantum strains 22 different genotypes were identified. The combined analysis revealed that all strains belong to the worldwide predominating MON1-population, however most closely related to a subpopulation from Southeastern Europe, Maghreb, Middle East and Central Asia. The three observed Armenian clusters grouped within this subpopulation with strains from Greece/Turkey, and from Central Asia, respectively. Ecological niche modeling based on VL cases and collected proven vectors (P. balcanicus, P. kandelakii) identified Yerevan and districts Lori, Tavush, Syunik, Armavir, Ararat bordering Georgia, Turkey, Iran and Azerbaijan as most suitable for the vectors and with the highest risk for VL transmission. Due to climate change the suitable habitat for VL transmission will expand in future all over Armenia.

CONCLUSIONS

Genetic diversity and population structure of the causative agent of VL in Armenia were addressed for the first time. Further genotyping studies should be performed with samples from infected humans, animals and sand flies from all active foci including the neighboring countries to understand transmission cycles, re-emergence, spread, and epidemiology of VL in Armenia and the entire Transcaucasus enabling epidemiological monitoring.

Prevalence of Genetically Complex Leishmania Strains With Hybrid and Mito-Nuclear Discordance

Approximately 20 Leishmania species are known to cause cutaneous, mucocutaneous, and visceral disorders in humans. Identification of the causative species in infected individuals is important for appropriate treatment and a favorable prognosis because infecting species are known to be the major determinant of clinical manifestations and may affect treatments for leishmaniasis. Although Leishmania species have been conventionally identified by multilocus enzyme electrophoresis, genetic analysis targeting kinetoplast and nuclear DNA (kDNA and nDNA, respectively) is now widely used for this purpose. Recently, we conducted countrywide epidemiological studies of leishmaniasis in Ecuador and Peru to reveal prevalent species using PCR-RFLP targeting nDNA, and identified unknown hybrid parasites in these countries together with species reported previously. Furthermore, comparative analyses of kDNA and nDNA revealed the distribution of parasites with mismatches between these genes, representing the first report of mito-nuclear discordance in protozoa. The prevalence of an unexpectedly high rate (~10%) of genetically complex strains including hybrid strains, in conjunction with the observation of mito-nuclear discordance, suggests that genetic exchange may occur more frequently than previously thought in natural Leishmania populations. Hybrid Leishmania strains resulting from genetic exchanges are suggested to cause more severe clinical symptoms when compared with parental strains, and to have increased transmissibility by vectors of the parental parasite species. Therefore, it is important to clarify how such genetic exchange influences disease progression and transmissibility by sand flies in nature. In addition, our aim was to identify where and how the genetic exchange resulting in the formation of hybrid and mito-nuclear discordance occurs.

Models in parasite and pathogen evolution: Genomic analysis reveals predominant clonality and progressive evolution at all evolutionary scales in parasitic protozoa, yeasts and bacteria

The predominant clonal evolution (PCE) model of pathogenic microorganisms postulates that the impact of genetic recombination in those pathogens' natural populations is not enough to erase a persistent phylogenetic signal at all evolutionary scales from microevolution till geological times in the whole ecogeographical range of the species considered. We have tested this model with a set of representative parasitic protozoa, yeasts and bacteria in the light of the most recent genomic data. All surveyed species, including those that were considered as highly recombining, exhibit similar PCE patterns above and under the species level, from macro- to micro-evolutionary scales (Russian doll pattern), suggesting gradual evolution. To our knowledge, it is the first time that such a strong common evolutionary feature among very diverse pathogens has been evidenced. The implications of this model for basic biology and applied research are exposed. These implications include our knowledge on the pathogens' reproductive mode, their population structure, the possibility to type strain and to follow up epidemics (molecular epidemiology) and to revisit pathogens' taxonomy through a flexible use of the phylogenetic species concept (Cracraft, 1983).

Genomic Diversification, Structural Plasticity, and Hybridization in Leishmania (Viannia) braziliensis

Leishmania (Viannia) braziliensis is an important Leishmania species circulating in several Central and South American countries. Among Leishmania species circulating in Brazil, Argentina and Colombia, L. braziliensis has the highest genomic variability. However, genomic variability at the whole genome level has been only studied in Brazilian and Peruvian isolates; to date, no Colombian isolates have been studied. Considering that in Colombia, L. braziliensis is a species with great clinical and therapeutic relevance, as well as the role of genetic variability in the epidemiology of leishmaniasis, we analyzed and evaluated intraspecific genomic variability of L. braziliensis from Colombian and Bolivian isolates and compared them with Brazilian isolates. Twenty-one genomes were analyzed, six from Colombian patients, one from a Bolivian patient, and 14 Brazilian isolates downloaded from public databases. The results obtained of Phylogenomic analysis showed the existence of four well-supported clades, which evidenced intraspecific variability. The whole-genome analysis revealed structural variations in the somy, mainly in the Brazilian genomes (clade 1 and clade 3), low copy number variations, and a moderate number of single-nucleotide polymorphisms (SNPs) in all genomes analyzed. Interestingly, the genomes belonging to clades 2 and 3 from Colombia and Brazil, respectively, were characterized by low heterozygosity (~90% of SNP loci were homozygous) and regions suggestive of loss of heterozygosity (LOH). Additionally, we observed the drastic whole genome loss of heterozygosity and possible hybridization events in one genome belonging to clade 4. Unique/shared SNPs between and within the four clades were identified, revealing the importance of some of them in biological processes of L. braziliensis. Our analyses demonstrate high genomic variability of L. braziliensis in different regions of South America, mainly in Colombia and suggest that this species exhibits striking genomic diversity and a capacity of genomic hybridization; additionally, this is the first study to report whole-genome sequences of Colombian L. braziliensis isolates.

Application of CRISPR/Cas9-Based Reverse Genetics in Leishmania braziliensis: Conserved Roles for HSP100 and HSP23

The protozoan parasite Leishmania (Viannia) braziliensis (L. braziliensis) is the main cause of human tegumentary leishmaniasis in the New World, a disease affecting the skin and/or mucosal tissues. Despite its importance, the study of the unique biology of L. braziliensis through reverse genetics analyses has so far lagged behind in comparison with Old World Leishmania spp. In this study, we successfully applied a cloning-free, PCR-based CRISPR–Cas9 technology in L. braziliensis that was previously developed for Old World Leishmania major and New World L. mexicana species. As proof of principle, we demonstrate the targeted replacement of a transgene (eGFP) and two L. braziliensis single-copy genes (HSP23 and HSP100). We obtained homozygous Cas9-free HSP23- and HSP100-null mutants in L. braziliensis that matched the phenotypes reported previously for the respective L. donovani null mutants. The function of HSP23 is indeed conserved throughout the Trypanosomatida as L. majorHSP23 null mutants could be complemented phenotypically with transgenes from a range of trypanosomatids. In summary, the feasibility of genetic manipulation of L. braziliensis by CRISPR–Cas9-mediated gene editing sets the stage for testing the role of specific genes in that parasite’s biology, including functional studies of virulence factors in relevant animal models to reveal novel therapeutic targets to combat American tegumentary leishmaniasis.