A primer for Leishmania population genetic studies

GenAPoPop 1.0: A user-friendly software to analyse genetic diversity and structure from partially clonal and selfed autopolyploid organisms

Autopolyploidy is quite common in most clades of eukaryotes. The emergence of sequence-based genotyping methods with individual and marker tags now enables confident allele dosage, overcoming the main obstacle to the democratization of the population genetic approaches when studying ecology and evolution of autopolyploid populations and species. Reproductive modes, including clonality, selfing and allogamy, have deep consequences on the ecology and evolution of population and species. Analysing genetic diversity and its dynamics over generations is one efficient way to infer the relative importance of clonality, selfing and allogamy in populations. GenAPoPop is a user-friendly solution to compute the specific corpus of population genetic indices, including indices about genotypic diversity, needed to analyse partially clonal, selfed and allogamous polysomic populations genotyped with confident allele dosage. It also easily provides the posterior probabilities of quantitative reproductive modes in autopolyploid populations genotyped at two-time steps and a graphical representation of the minimum spanning trees of the genetic distances between polyploid individuals, facilitating the interpretation of the genetic coancestry between individuals in hierarchically structured populations. GenAPoPop complements the previously existing solutions, including SPAGEDI and POLYGENE, to use genotypings to study the ecology and evolution of autopolyploid populations. It was specially developed with a simple graphical interface and workflow, and comes with a simulator to facilitate practical courses and teaching of population genetics for autopolyploid populations.

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.

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.

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

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

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.

Reproduction in Trypanosomatids: Past and Present

Simple Summary

The reproduction of trypanosomatids is a fundamental issue for host–parasite interaction, and its biological importance lies in knowing how these species acquire new defense mechanisms against the countermeasures imposed by the host, which is consistent with the theory of the endless race or the Red Queen hypothesis for the existence of meiotic sex. Moreover, the way these species re-produce may also be at the origin of novel and more virulent clades and is relevant from a thera-peutic or vaccination point of view, as sex may contribute to increased tolerance and even to the rapid acquisition of drug resistance mechanisms. Kinetoplastids are single-celled organisms, many of them being responsible for important parasitic diseases, globally termed neglected diseases, which are endemic in low-income countries. Leishmaniasis, African (sleeping sickness) and American trypanosomiasis (Chagas disease) caused by trypanosomatids are among the most ne-glected tropical scourges related to poverty and poor health systems. The reproduction of these microorganisms has long been considered to be clonal due to population genetic observations. However, there is increasing evidence of true sex and genetic exchange events under laboratory conditions. We would like to highlight the importance of this topic in the field of host/parasite in-terplay, virulence, and drug resistance.

Abstract

Diseases caused by trypanosomatids (Sleeping sickness, Chagas disease, and leishmaniasis) are a serious public health concern in low-income endemic countries. These diseases are produced by single-celled parasites with a diploid genome (although aneuploidy is frequent) organized in pairs of non-condensable chromosomes. To explain the way they reproduce through the analysis of natural populations, the theory of strict clonal propagation of these microorganisms was taken as a rule at the beginning of the studies, since it partially justified their genomic stability. However, numerous experimental works provide evidence of sexual reproduction, thus explaining certain naturally occurring events that link the number of meiosis per mitosis and the frequency of mating. Recent techniques have demonstrated genetic exchange between individuals of the same species under laboratory conditions, as well as the expression of meiosis specific genes. The current debate focuses on the frequency of genomic recombination events and its impact on the natural parasite population structure. This paper reviews the results and techniques used to demonstrate the existence of sex in trypanosomatids, the inheritance of kinetoplast DNA (maxi- and minicircles), the impact of genetic exchange in these parasites, and how it can contribute to the phenotypic diversity of natural populations.

Nuclear and kinetoplast DNA analyses reveal genetically complex Leishmania strains with hybrid and mito-nuclear discordance in Peru

Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) analysis of the mannose phosphate isomerase (mpi) gene was applied to 134 skin samples collected from patients with cutaneous leishmaniasis (CL) in Peru for identification of the infecting parasite at the species level, and the results were compared with those of cytochrome b (cyt b) gene sequencing obtained in previous studies. Although most results (121/134) including 4 hybrids of Leishmania (Viannia) braziliensis and L. (V.) peruviana corresponded to those obtained in the previous study, PCR-RFLP analyses revealed the distribution of putative hybrid strains between L. (V.) peruviana and L. (V.) lainsoni in two samples, which has never been reported. Moreover, parasite strains showing discordance between kinetoplast and nuclear genes (kDNA and nDNA), so-called mito-nuclear discordance, were identified in 11 samples. Of these, six strains had the kDNAs of L. (V.) braziliensis or L. (V.) peruviana and nDNAs of L. (V.) guyanensis, and three strains had the kDNAs of L. (V.) shawi and nDNAs of L. (V.) braziliensis. The rest were identified as mito-nuclear discordance strains having kDNAs of L. (V.) braziliensis or L. (V.) peruviana and nDNAs of L. (V.) lainsoni, and kDNAs of L. (V.) lainsoni and nDNAs of L. (V.) braziliensis. The results demonstrate that Leishmania strains in Peru are genetically more complex than previously considered.

Protozoan parasites of the genus Leishmania are able to undergo genetic exchange during their growth. The previous description of hybrids in Peru and the recent discovery of unexpected genetically complex strains having characteristics of both hybrid and mito-nuclear discordance in its neighbouring country (Ecuador) with a similar eco-epidemiological situation led us to consider that the genetic structure of Leishmania strains in Peru is more complicated than previously thought. In an effort to revise the data on Leishmania strain dispersion in Peru and to search for evidence of genetic recombination, the present study was conducted. A polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) analysis targeting the mannose phosphate isomerase (mpi) gene sequence was performed to identify the infecting parasite at the species level in 134 skin samples collected from patients with cutaneous leishmaniasis (CL) in Peru, and the results were compared with those of cytochrome b (cyt b) gene sequencing obtained in previous studies. Most results (121/134) including 4 hybrids between L. (V.) braziliensis and L. (V.) peruviana showed agreement between PCR-RFLP of the mpi gene and cyt b gene sequence analysis; however, 13 of 134 samples revealed the distribution of strains with hybrids and mito-nuclear discordance. The results demonstrate that genetically complex Leishmania strains are present in Peru. These findings indicate that Leishmania strain dispersion in Peru is genetically more complex than previously considered. Further prospective studies including larger samples and the isolation of parasite strains are required to update the available data.

Ecological divergence and hybridization of Neotropical Leishmania parasites

Parasites are interesting models for studying speciation processes because they have a high potential for specialization, thanks to the intimate ecological association with their hosts and vectors. Yet little is known about the circumstances under which new parasite lineages emerge. Here we studied the genome diversity of parasites of the Leishmania braziliensis species complex that inhabit both Amazonian and Andean biotas in Peru. We identify three major parasite lineages that occupy particular ecological niches and show that these emerged during forestation changes over the past 150,000 y. We furthermore discovered that meiotic recombination between Amazonian and Andean lineages resulted in full-genome hybrids presenting mixed mitochondrial genomes, providing insights into the genetic consequences of hybridization in parasitic protozoa.

The tropical Andes are an important natural laboratory to understand speciation in many taxa. Here we examined the evolutionary history of parasites of the Leishmania braziliensis species complex based on whole-genome sequencing of 67 isolates from 47 localities in Peru. We first show the origin of Andean Leishmania as a clade of near-clonal lineages that diverged from admixed Amazonian ancestors, accompanied by a significant reduction in genome diversity and large structural variations implicated in host–parasite interactions. Within the Andean species, patterns of population structure were strongly associated with biogeographical origin. Molecular clock and ecological niche modeling suggested that the history of diversification of the Andean lineages is limited to the Late Pleistocene and intimately associated with habitat contractions driven by climate change. These results suggest that changes in forestation over the past 150,000 y have influenced speciation and diversity of these Neotropical parasites. Second, genome-scale analyses provided evidence of meiotic-like recombination between Andean and Amazonian Leishmania species, resulting in full-genome hybrids. The mitochondrial genome of these hybrids consisted of homogeneous uniparental maxicircles, but minicircles originated from both parental species. We further show that mitochondrial minicircles—but not maxicircles—show a similar evolutionary pattern to the nuclear genome, suggesting that compatibility between nuclear-encoded mitochondrial genes and minicircle-encoded guide RNA genes is essential to maintain efficient respiration. By comparing full nuclear and mitochondrial genome ancestries, our data expand our appreciation on the genetic consequences of diversification and hybridization in parasitic protozoa.

Development of a Multilocus sequence typing (MLST) scheme for Pan-Leishmania

Since the description of the Leishmania genus, its identification and organization have been a challenge. A high number of molecular markers have been developed to resolve phylogenetic differences at the species level and for addressing key epidemiological and population genetics questions. Based on Multilocus enzyme electrophoresis (MLEE), Multilocus sequence typing (MLST) schemes have been developed using different gene candidates. From 38 original gene targets proposed by other authors, 27 of them were chosen. In silico selection was made by analyzing free access genomic sequence data of 33 Leishmania species, one Paraleishmania representative, and one outgroup, in order to select the best 15 loci. De novo amplifications and primers redesign of these 15 genes were analyzed over a panel of 20 reference strains and isolates. Phylogenetic analysis was made at every step. Two MLST schemes were selected. The first one was based on the analysis of three-gene fragments, and it is suitable for species assignment as well as basic phylogenetic studies. By the addition of seven-genes, an approach based on the analysis of ten-gene fragments was also proposed. This is the first work that two optimized MLST schemes have been suggested, validated against a phylogenetically diverse panel of Leishmania isolates. MLST is potentially a powerful phylogenetic approach, and most probably the new gold standard for Leishmania spp. characterization.

Revisiting FIS, FST, Wahlund Effects, and Null Alleles

Null alleles and Wahlund effects are well known causes of heterozygote deficits in empirical population genetics studies as compared to Hardy-Weinberg genotypic expectations. Some authors have theoretically studied the relationship of Wright's FIS computed from subsamples displaying a Wahlund effect and FST before the Wahlund effect, as can occasionally be obtained from populations of long-lived organisms. In the 2 subsample case, a positive relationship between these 2 parameters across loci would represent a signature of Wahlund effects. Nevertheless, for most organisms, getting 2 independent subsamples of the same cohort and population, one with a Wahlund effect and the other without, is almost never achieved and most of the time, empirical population geneticists only collect a single sample, with or without a Wahlund effect, or with or without null alleles. Another issue is that null allele increase FIS and FST altogether and thus may also create such correlation. In this article, I show that, for organisms collected in a single sample, which corresponds to the most common situation, Wahlund effects and null alleles affect the values of both FIS and FST though in the opposite direction. I also show that Wahlund effect produces no or weak positive correlation between the 2 F-statistics, while null alleles generate a strong positive correlation between them. Variation of these F-statistics is small and even minimized for FST under Wahlund effects as compared to null alleles. I finally propose a determination key to interpret data with heterozygote deficits.

Dispersion of Leishmania (Leishmania) infantum in central-southern Brazil: Evidence from an integrative approach

Leishmania (Leishmania) infantum is the zoonotic agent of visceral leishmaniasis (VL), a disease with a global distribution. The transmission scenario of VL has been undergoing changes worldwide, with the biologic cycle invading urbanized areas and dispersing the parasites into other previously free areas. The epidemiological cycle in Brazil has dispersed from the Northeast to other regions of the country. In this study, an integrative approach, including genotyping Brazilian strains of L. (L.) infantum for 14 microsatellite markers and reviewing historical records of the disease, was used to assess dispersion routes throughout central-southern Brazil. Our results support three L. (L.) infantum dispersion routes: A) dispersion from Bolivia to the states of Mato Grosso, Mato Grosso do Sul and São Paulo via the Bolivia-Brazil gas pipeline from 1998 to 2005; B) VL dispersion from Paraguay to the Brazilian side of the triple border (Foz do Iguaçu and Santa Terezinha de Itaipu) during after 2012; and C) emergence of a new L. (L.) infantum cluster in western Santa Catarina State and its dispersion to southern Paraná State (municipality of Pato Branco), after 2013. Hypotheses regarding possible entries of Leishmania (L.) infantum into the area of the triple border are presented and discussed. Understanding how VL has dispersed is vital to the development of control measures for this disease and to avoid future dispersion events.

The dispersion of visceral leishmaniasis is an enigma. The State of Paraná, in southern Brazil, borders the states of São Paulo and Mato Grosso, which have experienced LV epidemics over the past 20 years. Therefore, we expected that the disease would enter this state through the contiguity of epidemics from other regions following by "ghost shadows". However, in 2012, the vectors of the parasite were reported in the western region (Foz do Iguaçu) of Paraná state, far from the epidemic regions. In the cross-sectional study, 23.8% of the dogs were infected, which is more than the eyes can see, showing an unexpected scenario where the disease was already widespread in the city. Now the question was: where does the life cycle element came from? In this study, we used genetic markers to understand the dispersion of Leishmania infantum throughout central-southern Brazil. Our results showed two possible agent inputs in the Paraná state, one coming from Paraguay and, another from Santa Catarina state. When we verify our results we perceived the monitoring importance of the distribution of these agents by diverse hypotheses, not only those that the scientific literature presents. Another relevant factor is always to be attentive to the environmental and socioeconomic events that can provide this dispersion.

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

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

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

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

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

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

Elucidating in vitro and in vivo phenotypic behaviour of L. infantum/L. major natural hybrids

The clinical manifestation and course of Leishmania infections depend on factors such as species, virulence and host-immunity. Although trypanosomatids are considered to have clonal propagation, genetic hybridization has produced successful natural hybrid lineages. Hybrids displaying strong selective advantages may have an impact on pathogenesis and the eco-epidemiology of leishmaniasis. Thus, characterization of phenotypic properties of Leishmania hybrids could bring significant insight into the biology, infectivity, pathogenicity and transmission dynamics of these atypical strains. The present study focuses on phenotypic features and survival capacity of Leishmania infantum/Leishmania major hybrid isolates as compared with representative putative parental species, L. infantum and L. major. In vitro assays (growth kinetics, susceptibility to different conditions) and in vivo infection (parasite detection and histopathological alterations) showed that hybrids present higher growth capacity and decreased susceptibility to reactive oxygen species. Furthermore, evaluation of infected spleen tissue suggests that hybrids induce a stronger immune reaction than their putative parents, leading to the development of white pulp hyperplasia in B-lymphocyte compartments. Overall, these hybrids have shown high plasticity in terms of their general behaviour within the different phenotypic parameters, suggesting that they might have acquired genetic features conferring different mechanisms to evade host cells.

Evidence of Subpopulations with Distinct Biological Features Within a Leishmania (Viannia) braziliensis Strain

The present study demonstrates that the Leishmania (Viannia) braziliensis strain MCAN/BR/1998/R619 is composed of multiple subpopulations with measurable distinctions. Single parasites were separated from a culture of promastigotes in stationary phase by cell sorting and then cultivated as subpopulations. Subsequently, these subpopulations were evaluated for features of in vitro growth, infectivity to murine macrophages and proteinase gene expression. The first evidence of distinct characteristics was observed during the in vitro cultivation of isolated subpopulations, as distinct clusters of patterns were formed among the cultures, indicating the existence of quantifiable fluctuations in metrics. Further, when infecting murine macrophages, the subpopulations induced distinct patterns of production of immune response mediators. While some subpopulations mainly induced the production of IL-1β, IL-6 and TNF-α, others induced the production of IL-12p70 and nitric oxide. Finally, amastigotes of these subpopulations had higher expression of proteinase genes than promastigotes. Additionally, cysteine proteinase, serine proteinase, metalloproteinase and aspartic proteinases were differentially expressed in promastigote and amastigote forms. These data suggest the existence of distinct profiles for the L. (V.) braziliensis MCAN/BR/1998/R619 strain and subpopulations that could drive the success of parasite adaptation to the environments that they inhabit.

Induction of allopurinol resistance in Leishmania infantum isolated from dogs

Resistance to allopurinol in zoonotic canine leishmaniasis has been recently shown to be associated with disease relapse in naturally-infected dogs. However, information regarding the formation of resistance and its dynamics is lacking. This study describes the successful in-vitro induction of allopurinol resistance in Leishmania infantum cultured under increasing drug pressure. Allopurinol susceptibility and growth rate of induced parasites were monitored over 23 weeks and parasite clones were tested at selected time points and compared to their parental lines, both as promastigotes and as amastigotes. Allopurinol resistance was formed in strains from two parasite stocks producing a 20-fold rise in IC₅₀ along three distinct growth phases. In addition, characteristic differential clustering of single nucleotide polymorphisms (SNP) was found in drug sensitive and resistant parasite clones. Results confirm that genetic polymorphism, as well as clonal heterogeneity, contribute to in-vitro resistance to allopurinol, which is likely to occur in natural infection.

Visceral leishmaniasis caused by the parasite Leishmania infantum is a neglected tropical disease transmitted from animal hosts to humans by sand fly bites. This potentially fatal disease affects thousands of people annually and threatens millions who live in disease risk areas. Domestic dogs are considered as the main reservoir of this parasite which can also cause a severe chronic canine disease. Allopurinol is the main drug used for long term treatment of this disease but it often does not eliminate infection in dogs. We have recently demonstrated that allopurinol resistant parasites can be isolated from naturally infected dogs that have developed clinical recurrence of disease during allopurinol treatment. In this study we aimed to see if resistance can be induced in susceptible parasite strains isolated from sick dogs by growing them in increasing drug concentrations under laboratory conditions. The changes in allopurinol susceptibility were measured and the impact of drug on parasite growth was monitored over 23 weeks. Induction of resistance was successful producing parasites 20-folds less susceptible to the drug. The pattern of change in drug susceptibility suggests that a genetic change is responsible for the increased resistance which is likely to mimic the formation of resistance in dogs.

Comparative analyses of whole genome sequences of Leishmania infantum isolates from humans and dogs in northeastern Brazil

The genomic sequences of 20 Leishmania infantum isolates collected in northeastern Brazil were compared with each other and with the available genomic sequences of 29 L. infantum/donovani isolates from Nepal and Turkey. The Brazilian isolates were obtained in the early 1990s or since 2009 from patients with visceral or non-ulcerating cutaneous leishmaniasis, asymptomatic humans, or dogs with visceral leishmaniasis. Two isolates were from the blood and bone marrow of the same visceral leishmaniasis patient. All 20 genomic sequences display 99.95% identity with each other and slightly less identity with a reference L. infantum genome from a Spanish isolate. Despite the high identity, analysis of individual differences among the 32 million base pair genomes showed sufficient variation to allow the isolates to be clustered based on the primary sequence. A major source of variation detected was in chromosome somy, with only four of the 36 chromosomes being predominantly disomic in all 49 isolates examined. In contrast, chromosome 31 was predominantly tetrasomic/pentasomic, consistent with its regions of synteny on two different disomic chromosomes of Trypanosoma brucei. In the Brazilian isolates, evidence for recombination was detected in 27 of the 36 chromosomes. Clustering analyses suggested two populations, in which two of the five older isolates from the 1990s clustered with a majority of recent isolates. Overall the analyses do not suggest individual sequence variants account for differences in clinical outcome or adaptation to different hosts. For the first known time, DNA of isolates from asymptomatic subjects were sequenced. Of interest, these displayed lower diversity than isolates from symptomatic subjects, an observation that deserves further investigation with additional isolates from asymptomatic subjects.

Reproduction in Leishmania: A focus on genetic exchange

One key process of the life cycle of pathogens is their mode of reproduction. Indeed, this fundamental biological process conditions the multiplication and the transmission of genes and thus the propagation of diseases in the environment. Reproductive strategies of protozoan parasites have been a subject of debate for many years, principally due to the difficulty in making direct observations of sexual reproduction (i.e. genetic recombination). Traditionally, these parasites were considered as characterized by a preeminent clonal structure. Nevertheless, with the development of elaborate culture experiments, population genetics and evolutionary and population genomics, several studies suggested that most of these pathogens were also characterized by constitutive genetic recombination events. In this opinion, we focused on Leishmania parasites, pathogens responsible of leishmaniases, a major public health issue. We first discuss the evolutionary advantages of a mixed mating reproductive strategy, then we review the evidence of genetic exchange, and finally we detail available tools to detect naturally occurring genetic recombination in Leishmania parasites and more generally in protozoan parasites.

Leishmania infections: Molecular targets and diagnosis

Progress in the diagnosis of leishmaniases depends on the development of effective methods and the discovery of suitable biomarkers. We propose firstly an update classification of Leishmania species and their synonymies. We demonstrate a global map highlighting the geography of known endemic Leishmania species pathogenic to humans. We summarize a complete list of techniques currently in use and discuss their advantages and limitations. The available data highlights the benefits of molecular markers in terms of their sensitivity and specificity to quantify variation from the subgeneric level to species complexes, (sub) species within complexes, and individual populations and infection foci. Each DNA-based detection method is supplied with a comprehensive description of markers and primers and proposal for a classification based on the role of each target and primer in the detection, identification and quantification of leishmaniasis infection. We outline a genome-wide map of genes informative for diagnosis that have been used for Leishmania genotyping. Furthermore, we propose a classification method based on the suitability of well-studied molecular markers for typing the 21 known Leishmania species pathogenic to humans. This can be applied to newly discovered species and to hybrid strains originating from inter-species crosses. Developing more effective and sensitive diagnostic methods and biomarkers is vital for enhancing Leishmania infection control programs.

Is Predominant Clonal Evolution a Common Evolutionary Adaptation to Parasitism in Pathogenic Parasitic Protozoa, Fungi, Bacteria, and Viruses?

We propose that predominant clonal evolution (PCE) in microbial pathogens be defined as restrained recombination on an evolutionary scale, with genetic exchange scarce enough to not break the prevalent pattern of clonal population structure. The main features of PCE are (1) strong linkage disequilibrium, (2) the widespread occurrence of stable genetic clusters blurred by occasional bouts of genetic exchange ('near-clades'), (3) the existence of a "clonality threshold", beyond which recombination is efficiently countered by PCE, and near-clades irreversibly diverge. We hypothesize that the PCE features are not mainly due to natural selection but also chiefly originate from in-built genetic properties of pathogens. We show that the PCE model obtains even in microbes that have been considered as 'highly recombining', such as Neisseria meningitidis, and that some clonality features are observed even in Plasmodium, which has been long described as panmictic. Lastly, we provide evidence that PCE features are also observed in viruses, taking into account their extremely fast genetic turnover. The PCE model provides a convenient population genetic framework for any kind of micropathogen. It makes it possible to describe convenient units of analysis (clones and near-clades) for all applied studies. Due to PCE features, these units of analysis are stable in space and time, and clearly delimited. The PCE model opens up the possibility of revisiting the problem of species definition in these organisms. We hypothesize that PCE constitutes a major evolutionary strategy for protozoa, fungi, bacteria, and viruses to adapt to parasitism.

Taxonomy, diversity, temporal and geographical distribution of Cutaneous Leishmaniasis in Colombia: A retrospective study

Leishmaniases are tropical zoonotic diseases, caused by kinetoplastid parasites from the genus Leishmania. New World (NW) species are related to sylvatic cycles although urbanization processes have been reported in some South American Countries such as Colombia. Currently, few studies show the relative distribution of Leishmania species related to cutaneous Leishmaniasis (CL) in South America due to the lack of accurate surveillance and public health systems. Herein, we conducted a systematic estimation of the Leishmania species causing CL in Colombia from 1980 to 2001 via molecular typing and isoenzymes. A total of 327 Leishmania isolates from humans, sandflies and reservoirs were typed as L. panamensis 61.3% (201), L. braziliensis 27.1% (88), L. lainsoni 0.6% (2), L. guyanensis 0.9% (3), L. infantum chagasi 4% (12), L. equatoriensis 0.6% (2), L. mexicana 2.1% (8), L. amazonensis 2.8% (9) and L. colombiensis 0.6% (2). This is the first report of two new Leishmania species circulating in Colombia and suggests the need to convince the Colombian government about the need to deploy and standardize tools for the species identification to provide adequate management to individuals suffering this pathology.

The population genetics of Trypanosoma cruzi revisited in the light of the predominant clonal evolution model

Comparing the population structure of Trypanosoma cruzi with that of other pathogens, including parasitic protozoa, fungi, bacteria and viruses, shows that the agent of Chagas disease shares typical traits with many other species, related to a predominant clonal evolution (PCE) pattern: statistically significant linkage disequilibrium, overrepresented multilocus genotypes, near-clades (genetic subdivisions somewhat blurred by occasional genetic exchange/hybridization) and "Russian doll" patterns (PCE is observed, not only at the level of the whole species, but also, within the near-clades). Moreover, T. cruzi population structure exhibits linkage with the diversity of several strongly selected genes, with gene expression profiles, and with some major phenotypic traits. We discuss the evolutionary significance of these results, and their implications in terms of applied research (molecular epidemiology/strain typing, analysis of genes of interest, vaccine and drug design, immunological diagnosis) and of experimental evolution. Lastly, we revisit the long-term debate of describing new species within the T. cruzi taxon.

Multilocus microsatellite typing of Leishmania and clinical applications: a review

Microsatellite markers have been used for Leishmania genetic studies worldwide, giving useful insight into leishmaniasis epidemiology. Understanding the geographic distribution, dynamics of Leishmania populations, and disease epidemiology improved markedly with this tool. In endemic foci, the origins of antimony-resistant strains and multidrug treatment failures were explored with multilocus microsatellite typing (MLMT). High genetic variability was detected but no association between parasite genotypes and drug resistance was established. An association between MLMT profiles and clinical disease manifestations was highlighted in only three studies and this data needs further confirmation. At the individual level, MLMT provided information on relapse and reinfection when multiple leishmaniasis episodes occurred. This information could improve knowledge of epidemiology and guide therapeutic choices for active chronic visceral leishmaniasis, the disease form in some HIV-positive patients.