A battery of 12 microsatellite markers for genetic analysis of the Leishmania (Viannia) guyanensis complex

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.

Tree-based unrooted nonbinary phylogenetic networks

Phylogenetic networks are a generalisation of phylogenetic trees that allow for more complex evolutionary histories that include hybridisation-like processes. It is of considerable interest whether a network can be considered 'tree-like' or not, which leads to the introduction of tree-based networks in the rooted, binary context. Tree-based networks are those networks which can be constructed by adding additional edges into a phylogenetic tree, called the base tree. Previous extensions have considered extending to the binary, unrooted case and the nonbinary, rooted case. In this paper, we extend tree-based networks to the context of unrooted, nonbinary networks in three ways, depending on the types of additional edges that are permitted. Further, we study fully tree-based networks which are phylogenetic networks in which every embedded tree is a base tree. We also extend this concept to unrooted, nonbinary, phylogenetic networks and classify the resulting networks. Finally, we derive some results on the colourability of tree-based networks, which can be useful to determine whether a network is tree-based.

Distinct genetic profiles of Leishmania (Viannia) braziliensis associate with clinical variations in cutaneous-leishmaniasis patients from an endemic area in Brazil

American tegumentary leishmaniasis (ATL) samples obtained from the lesions of patients with typical (n = 25, 29%), atypical (n = 60, 69%) or both (n = 2%) clinical manifestations were analysed by multilocus enzyme electrophoresis, hsp70 restriction-fragment length polymorphism (PCR-RFLP), hsp70 sequencing and phylogenetics methods. The hsp70 PCR-RFLP analysis revealed two different profiles whose the most samples differed from those expected for Leishmania braziliensis and the other Leishmania species tested: of 39 samples evaluated, two (5%) had a restriction profile corresponding to L. braziliensis, and 37 (95%) had a restriction profile corresponding to a variant pattern. A 1300-bp hsp70 gene fragment was sequenced to aid in parasite identification and a phylogenetic analysis was performed including 26 consensus sequences from the ATL patient's samples and comparing to other Leishmania and trypanosomatids species. The dendrogram allowed to observe a potential population structure of L. braziliensis complex in the studied region, emphasizing that the majority of clinical samples presented a variant genetic profile. Of interest, the L. braziliensis diversity was associated with different clinical manifestations whose parasites with hsp70 variant profile were associated with atypical lesions. The results may be helpful to improve the diagnosis, treatment and control measures of the ATL in endemic areas.

Genetic Diversity and Population Structure of Leishmania infantum from Southeastern France: Evaluation Using Multi-Locus Microsatellite Typing

In the south of France, Leishmania infantum is responsible for numerous cases of canine leishmaniasis (CanL), sporadic cases of human visceral leishmaniasis (VL) and rare cases of cutaneous and muco-cutaneous leishmaniasis (CL and MCL, respectively). Several endemic areas have been clearly identified in the south of France including the Pyrénées-Orientales, Cévennes (CE), Provence (P), Alpes-Maritimes (AM) and Corsica (CO). Within these endemic areas, the two cities of Nice (AM) and Marseille (P), which are located 150 km apart, and their surroundings, concentrate the greatest number of French autochthonous leishmaniasis cases. In this study, 270 L. infantum isolates from an extended time period (1978–2011) from four endemic areas, AM, P, CE and CO, were assessed using Multi-Locus Microsatellite Typing (MLMT). MLMT revealed a total of 121 different genotypes with 91 unique genotypes and 30 repeated genotypes. Substantial genetic diversity was found with a strong genetic differentiation between the Leishmania populations from AM and P. However, exchanges were observed between these two endemic areas in which it seems that strains spread from AM to P. The genetic differentiations in these areas suggest strong epidemiological structuring. A model-based analysis using STRUCTURE revealed two main populations: population A (consisting of samples primarily from the P and AM endemic areas with MON-1 and non-MON-1 strains) and population B consisting of only MON-1 strains essentially from the AM endemic area. For four patients, we observed several isolates from different biological samples which provided insight into disease relapse and re-infection. These findings shed light on the transmission dynamics of parasites in humans. However, further data are required to confirm this hypothesis based on a limited sample set. This study represents the most extensive population analysis of L. infantum strains using MLMT conducted in France.

In the south of France, the parasite Leishmania infantum is responsible for diseases that primarily affect dogs but can also impact humans. Several endemic areas have been clearly identified in the south of France including the Pyrénées-Orientales, Cévennes (CE), Provence (P), Alpes-Maritimes (AM) and Corsica (CO). In this study, 270 L. infantum isolates from four endemic areas, AM, P, CE and CO, were assessed using Multi-Locus Microsatellite Typing (MLMT), a tool applied for population genetic studies. MLMT revealed a strong genetic differentiation between the Leishmania populations from AM and P with exchanges observed between these two endemic areas. For four patients, the occurrence of disease relapses and re-infections was examined. These findings shed light on the transmission dynamics of parasites in humans. This study represents the most extensive population analysis of L. infantum isolates using MLMT conducted in France.

Genome-wide discovery and development of polymorphic microsatellites from Leishmania panamensis parasites circulating in central Panama

Background

The parasite Leishmania panamensis is the main cause of leishmaniasis in Panama. The disease is largely uncontrolled, with a rising incidence and no appropriate control measures. While microsatellites are considered some of the best genetic markers to study population genetics and molecular epidemiology in these and other parasites, none has been developed for L. panamensis.

Findings

Here we have developed and tested a new panel of microsatellites for this species, based on high-throughput genome-wide screening. The new set of microsatellites is composed of seventeen loci, mainly spanning trinucleotide or longer motifs. We have evaluated the sensitivity and specificity of the panel based on a sample of 27 isolates obtained from cutaneous leishmaniasis patients from central Panama and also several reference species from both L. (Leishmania) and L. (Viannia) subgenera. The genetic equilibrium was assessed both intra- and inter-loci, while the reproductive mode was evaluated using several tests. The new SSR panel shows high polymorphism and sensitivity, as well as good specificity. The preliminary data described here for L. panamensis suggest extensive departure from Hardy-Weinberg proportions, significant linkage disequilibrium and strong deficit of heterozygotes. Several recombination tests involving multilocus linkage disequilibrium and a phylogenetic approach allowed rejection of frequent recombination in our dataset.

Conclusions

The genome-wide strategy described here proved to be useful to identify and test new polymorphic SSR loci in Leishmania. The new panel of polymorphic microsatellites is a valuable contribution to the existing molecular markers for the study of genetic structure and other aspects of this important species.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-1153-2) contains supplementary material, which is available to authorized users.

A primer for Leishmania population genetic studies

Leishmaniases remain a major public health problem. Despite the development of elaborate experimental techniques and sophisticated statistical tools, how these parasites evolve, adapt themselves to new environmental compartments and hosts, and develop resistance to new drugs remains unclear. Leishmania parasites constitute a complex model from a biological, ecological, and epidemiological point of view but also with respect to their genetics and phylogenetics. With this in view, we seek to outline the criteria, caveats, and confounding factors to be considered for Leishmania population genetic studies. We examine how the taxonomic complexity, heterozygosity, intraspecific and interspecific recombination, aneuploidy, and ameiotic recombination of Leishmania intersect with population genetic studies of this parasite.

Genotype profile of Leishmania major strains isolated from tunisian rodent reservoir hosts revealed by multilocus microsatellite typing

Zoonotic cutaneous leishmaniasis (ZCL) caused by Leishmania (L.) major parasites represents a major health problem with a large spectrum of clinical manifestations. Psammomys (P.) obesus and Meriones (M.) shawi represent the most important host reservoirs of these parasites in Tunisia. We already reported that infection prevalence is different between these two rodent species. We aimed in this work to evaluate the importance of genetic diversity in L. major parasites isolated from different proven and suspected reservoirs for ZCL. Using the multilocus microsatellites typing (MLMT), we analyzed the genetic diversity among strains isolated from (i) P. obesus (n = 31), (ii) M. shawi (n = 8) and (iii) Mustela nivalis (n = 1), captured in Sidi Bouzid, an endemic region for ZCL located in the Center of Tunisia. Studied strains present a new homogeneous genotype profile so far as all tested markers and showed no polymorphism regardless of the parasite host-reservoir origin. This lack of genetic diversity among these L. major isolates is the first genetic information on strains isolated from Leishmania reservoirs hosts in Tunisia. This result indicates that rodent hosts are unlikely to exert a selective pressure on parasites and stresses on the similarity of geographic and ecological features in this study area. Overall, these results increase our knowledge among rodent reservoir hosts and L. major parasites interaction.

Polymorphisms and ambiguous sites present in DNA sequences of Leishmania clones: looking closer

In genetic studies of Leishmania parasites, co-dominant markers are chosen for their ability to detect heterozygous polymorphisms, to infer the occurrence of inbreeding and to resolve genetic variability. The majority of DNA sequence based reports perform conventional dye terminator cycle sequencing where perfectly ambiguous sites or double peaks in the chromatogram are interpreted as heterozygous strains. However, molecular peculiarities of the parasite such as aneuploidy, mixed populations and homologous recombination advise that data from regular DNA sequence analysis should be carefully evaluated. We report here a closer look at ambiguous sites observed in 6pgd DNA sequences obtained for a multilocus sequence analysis project on Leishmania (Viannia) strains. After comparing 286 DNA sequences from biological and molecular clones of six L. (Viannia) strains we could distinguish events that contribute to genetic variation in Leishmania (recombination, mutation, chromosomal mosaics). Also, the results suggest how diversity might not be completely revealed through regular DNA sequence analysis and demonstrate the importance for molecular epidemiology research to be aware of such possibilities while choosing samples for studies.

A call to arms: on refining Plasmodium vivax microsatellite marker panels for comparing global diversity

BACKGROUND

Microsatellite (MS) markers have become an important tool for studying the population diversity, evolutionary history and multiplicity of infection (MOI) of malaria parasite infections. MS are typically selected on the basis of being highly polymorphic. However, it is known that the polymorphic potential (mutability) of each marker can vary as much as two orders of magnitude, which radically changes how diversity is represented in the genome from one marker to the next. Over the past decade, approximately 240 Plasmodium vivax MS have been published, comprising nine major panels of markers. Inconsistent usage of each panel has resulted in a surfeit of descriptive genetic diversity data that are largely incomparable between populations. The objective of this study was to statistically evaluate the quality of individual MS markers in order to validate a refined panel of markers that will provide a balanced picture of P. vivax population diversity.

METHODS

All previously published data, including genetic diversity indices, MS parameters, and population parameters, were assembled from 18 different global studies into a flat file to facilitate statistical analysis and modelling using JMP® Genomics 6.0 (SAS Institute Inc, Cary, NC, USA). Statistical modeling was employed to down-select markers with extreme variation among the mean number of alleles, expected heterozygosity, maximum repeat length and/or chromosomal location of the repeat. Individual MS were analysed by step-down whole model linear regression and standard least squares fit models, both stratified by annual parasite incidence to identify MS markers with values significantly different from the mean.

RESULTS

Of the 42 MS under evaluation in this study, 18 (nine high priority) were identified as ideal candidates for measuring population diversity between global regions, while five (two high priority) additional markers were identified as candidates for MOI studies.

CONCLUSIONS

MS diversity was found to be a function of endemicity and motif structure. Evaluation of individual MS permitted the assembly of a refined panel of markers that can be reliably utilized in the field to compare population structures between global regions.

Population structure and evidence for both clonality and recombination among Brazilian strains of the subgenus Leishmania (Viannia)

Background/Objectives: Parasites of the subgenus Leishmania (Viannia) cause varying clinical symptoms ranging from cutaneous leishmaniases (CL) with single or few lesions, disseminated CL (DL) with multiple lesions to disfiguring forms of mucocutaneous leishmaniasis (MCL). In this population genetics study, 37 strains of L. (V.) guyanensis, 63 of L. (V.) braziliensis, four of L. (V.) shawi, six of L. (V.) lainsoni, seven of L. (V.) naiffi, one each of L. (V.) utingensis and L. (V.) lindenbergi, and one L. (V.) lainsoni/L. naiffi hybrid from different endemic foci in Brazil were examined for variation at 15 hyper-variable microsatellite markers.

Methodology/Principal findings: The multilocus microsatellite profiles obtained for the 120 strains were analysed using both model- and distance-based methods. Significant genetic diversity was observed for all L. (Viannia) strains studied. The two cluster analysis approaches identified two principal genetic groups or populations, one consisting of strains of L. (V.) guyanensis from the Amazon region and the other of strains of L. (V.) braziliensis isolated along the Atlantic coast of Brazil. A third group comprised a heterogeneous assembly of species, including other strains of L. braziliensis isolated from the north of Brazil, which were extremely polymorphic. The latter strains seemed to be more closely related to those of L. (V.) shawi, L. (V.) naiffi, and L. (V.) lainsoni, also isolated in northern Brazilian foci. The MLMT approach identified an epidemic clone consisting of 13 strains of L. braziliensis from Minas Gerais, but evidence for recombination was obtained for the populations of L. (V.) braziliensis from the Atlantic coast and for L. (V.) guyanensis.

Conclusions/Significance: Different levels of recombination versus clonality seem to occur within the subgenus L. (Viannia). Though clearly departing from panmixia, sporadic, but long-term sustained recombination might explain the tremendous genetic diversity and limited population structure found for such L. (Viannia) strains.

Cutaneous leishmaniasis (CL) constitutes a significant public health problem in all federal states of Brazil. Most cases are caused by parasites of the subgenus Leishmania (Viannia) which can cause a variety of clinical symptoms ranging from single or few lesions, disseminated CL with multiple lesions, to disfiguring forms of mucocutaneous leishmaniasis. This study has used a multilocus microsatellite typing approach for exploring the genetic diversity and population structure among 120 strains representing different subgenus L. (Viannia) species and different Brazilian CL foci. Genetic diversity within the subgenus was much higher than expected, especially within L. (V.) braziliensis, L. (V.) shawi, L. (V.) naiffi, and L. (V.) lainsoni which were all from the north of Brazil. These strains could not be assigned to well-defined populations, but presented a rather loosely associated group. Strains of L. (V.) braziliensis isolated along the Atlantic coast of Brazil and strains of L. (V.) guyanensis formed, however, two clearly separated populations exhibiting remarkable levels of sexual exchange. The latter finding is in contrast to previous studies suggesting clonal modes of propagation or inbreeding for natural populations of Leishmania parasites and might explain the genetic heterogeneity and limited population structure for Brazilian strains of subgenus L. (Viannia) observed in this study.

New insights on taxonomy, phylogeny and population genetics of Leishmania (Viannia) parasites based on multilocus sequence analysis

The Leishmania genus comprises up to 35 species, some with status still under discussion. The multilocus sequence typing (MLST)--extensively used for bacteria--has been proposed for pathogenic trypanosomatids. For Leishmania, however, a detailed analysis and revision on the taxonomy is still required. We have partially sequenced four housekeeping genes--glucose-6-phosphate dehydrogenase (G6PD), 6-phosphogluconate dehydrogenase (6PGD), mannose phosphate isomerase (MPI) and isocitrate dehydrogenase (ICD)--from 96 Leishmania (Viannia) strains and assessed their discriminatory typing capacity. The fragments had different degrees of diversity, and are thus suitable to be used in combination for intra- and inter-specific inferences. Species-specific single nucleotide polymorphisms were detected, but not for all species; ambiguous sites indicating heterozygosis were observed, as well as the putative homozygous donor. A large number of haplotypes were detected for each marker; for 6PGD a possible ancestral allele for L. (Viannia) was found. Maximum parsimony-based haplotype networks were built. Strains of different species, as identified by multilocus enzyme electrophoresis (MLEE), formed separated clusters in each network, with exceptions. NeighborNet of concatenated sequences confirmed species-specific clusters, suggesting recombination occurring in L. braziliensis and L. guyanensis. Phylogenetic analysis indicates L. lainsoni and L. naiffi as the most divergent species and does not support L. shawi as a distinct species, placing it in the L. guyanensis cluster. BURST analysis resulted in six clonal complexes (CC), corresponding to distinct species. The L. braziliensis strains evaluated correspond to one widely geographically distributed CC and another restricted to one endemic area. This study demonstrates the value of systematic multilocus sequence analysis (MLSA) for determining intra- and inter-species relationships and presents an approach to validate the species status of some entities. Furthermore, it contributes to the phylogeny of L. (Viannia) and might be helpful for epidemiological and population genetics analysis based on haplotype/diplotype determinations and inferences.

The genetic structure of Leishmania infantum populations in Brazil and its possible association with the transmission cycle of visceral leishmaniasis

Leishmania infantum is the etiologic agent of visceral leishmaniasis (VL) in the Americas, Mediterranean basin and West and Central Asia. Although the geographic structure of L. infantum populations from the Old World have been described, few studies have addressed the population structure of this parasite in the Neotropical region. We employed 14 microsatellites to analyze the population structure of the L. infantum strains isolated from humans and dogs from most of the Brazilian states endemic for VL and from Paraguay. The results indicate a low genetic diversity, high inbreeding estimates and a depletion of heterozygotes, which together indicate a predominantly clonal breeding system, but signs of sexual events are also present. Three populations were identified from the clustering analysis, and they were well supported by F statistics inferences and partially corroborated by distance-based. POP1 (111 strains) was observed in all but one endemic area. POP2 (31 strains) is also well-dispersed, but it was the predominant population in Mato Grosso (MT). POP3 (31 strains) was less dispersed, and it was observed primarily in Mato Grosso do Sul (MS). Strains originated from an outbreak of canine VL in Southern Brazil were grouped in POP1 with those from Paraguay, which corroborates the hypothesis of dispersal from Northeastern Argentina and Paraguay. The distribution of VL in MS seems to follow the west-east construction of the Bolivia-Brazil pipeline from Corumbá municipality. This may have resulted in a strong association of POP3 and Lutzomyia cruzi, which is the main VL vector in Corumbá, and a dispersion of this population in this region that was shaped by human interference. This vector also occurs in MT and may influence the structure of POP2. This paper presents significant advances in the understanding of the population structure of L. infantum in Brazil and its association with eco-epidemiological aspects of VL.

Multifaceted population structure and reproductive strategy in Leishmania donovani complex in one Sudanese village

Leishmania species of the subgenus Leishmania and especially L. donovani are responsible for a large proportion of visceral leishmaniasis cases. The debate on the mode of reproduction and population structure of Leishmania parasites remains opened. It has been suggested that Leishmania parasites could alternate different modes of reproduction, more particularly clonality and frequent recombinations either between related individuals (endogamy) or between unrelated individuals (outcrossing) within strongly isolated subpopulations. To determine whether this assumption is generalized to other species, a population genetics analysis within Leishmania donovani complex strains was conducted within a single village. The results suggest that a mixed-mating reproduction system exists, an important heterogeneity of subsamples and the coexistence of several genetic entities in Sudanese L. donovani. Indeed, results showed significant genetic differentiation between the three taxa (L. donovani, L. infantum and L. archibaldi) and between the human or canine strains of such taxa, suggesting that there may be different imbricated transmission cycles involving either dogs or humans. Results also are in agreement with an almost strict specificity of L. donovani stricto sensu to human hosts. This empirical study demonstrates the complexity of population structure in the genus Leishmania and the need to pursue such kind of analyses at the smallest possible spatio-temporal and ecological scales.

Leishmaniases are a serious public health problem, especially in developing countries, caused by Leishmania parasites and transmitted by sandfly bites. More information is needed on the population biology of these pathogens for diagnostic and epidemiological inquiries and for drug and vaccine elaboration. For studies dealing with the population genetics, exploring the genetic patterns of such organisms at microgeographic scales is fundamental. In this context, we made a population genetic study, based on 20 microsatellite loci, on 61 strains of Leishmania donovani complex collected in a Sudanese village, Babar El Fugara, during the epidemic of 1996–2000. Results showed that considering the whole sample as a single population was not adequate because of the coexistence of several genetic entities and a genetic differentiation between the human or canine strains. In addition, our findings suggested that clonality may have a strong impact on the L. donovani complex, unlike other Leishmania species. This study demonstrates the need to pursue population genetics studies in Leishmania species from sampling designs that control maximum possible confounding factors and to elaborate such kinds of analyses at the smallest possible spatio-temporal and ecological scales.

Reproductive strategies and population structure in Leishmania: substantial amount of sex in Leishmania Viannia guyanensis

Leishmania species of the subgenus Viannia and especially Leishmania Viannia guyanensis are responsible for a large proportion of New World leishmaniasis cases. Since a recent publication on Leishmania Viannia braziliensis, the debate on the mode of reproduction of Leishmania parasites has been reopened. A predominant endogamic reproductive mode (mating with relatives), together with strong Wahlund effects (sampling of strains from heterogeneous subpopulations), was indeed evidenced. To determine whether this hypothesis can be generalized to other Leishmania Viannia species, we performed a population genetic study on 153 human strains of L. (V.) guyanensis from French Guiana based on 12 microsatellite loci. The results revealed important homozygosity and very modest linkage disequilibrium, which is in agreement with a high level of sexual recombination and substantial endogamy. These results also revealed a significant isolation by distance with relatively small neighbourhoods and hence substantial viscosity of Leishmania populations in French Guiana. These results are of epidemiological relevance and suggest a major role for natural hosts and/or vectors in parasite strain diffusion across the country as compared to human hosts.

Microsatellite marker development in the protozoan parasite Perkinsus olseni

The analysis of an enriched partial genomic library and of public expressed sequence tag (EST) resources allowed the characterization of the first microsatellite loci in the protozoan parasite Perkinsus olseni. Clonal cultures from laboratory isolates derived from infected clams Ruditapes decussatus (from Spain), R. philippinarum (from Spain and Japan), and Austrovenus stutchburyi (from New Zealand) were used for the characterization of 12 microsatellites. Low variation was detected at most loci, with the number of alleles at polymorphic loci ranging from 2 to 7 (average 3.20 +/- 0.51) and gene diversity from 0.11 to 0.79 (average 0.40 +/- 0.07). Preliminary results show that (1) isolates of P. olseni are diploid cells, and (2) multiple infections can occur within a single host. Eight of the loci analyzed successfully cross-amplified in the congeneric species P. mediterraneus. These microsatellite markers will be useful to analyze in detail the population genetic structure of P. olseni, crucial for the efficient management of this parasitic disease.