Plasmodium vivax genomic surveillance in the Peruvian Amazon with Pv AmpliSeq assay.

BACKGROUND: Plasmodium vivax is the most predominant malaria species in Latin America, constituting 71.5% of malaria cases in 2021. With several countries aiming for malaria elimination, it is crucial to prioritize effectiveness of national control programs by optimizing the utilization of available resources and strategically implementing necessary changes. To support this, there is a need for innovative approaches such as genomic surveillance tools that can investigate changes in transmission intensity, imported cases and sources of reintroduction, and can detect molecular markers associated with drug resistance. METHODOLOGY/PRINCIPAL FINDINGS: Here, we apply a modified highly-multiplexed deep sequencing assay: Pv AmpliSeq v2 Peru. The tool targets a newly developed 41-SNP Peru barcode for parasite population analysis within Peru, the 33-SNP vivaxGEN-geo panel for country-level classification, and 11 putative drug resistance genes. It was applied to 230 samples from the Peruvian Amazon (2007-2020), generating baseline surveillance data. We observed a heterogenous P. vivax population with high diversity and gene flow in peri-urban areas of Maynas province (Loreto region) with a temporal drift using all SNPs detected by the assay (nSNP = 2909). In comparison, in an indigenous isolated area, the parasite population was genetically differentiated (FST = 0.07-0.09) with moderate diversity and high relatedness between isolates in the community. In a remote border community, a clonal P. vivax cluster was identified, with distinct haplotypes in drug resistant genes and ama1, more similar to Brazilian isolates, likely representing an introduction of P. vivax from Brazil at that time. To test its applicability for Latin America, we evaluated the SNP Peru barcode in P. vivax genomes from the region and demonstrated the capacity to capture local population clustering at within-country level. CONCLUSIONS/SIGNIFICANCE: Together this data shows that P. vivax transmission is heterogeneous in different settings within the Peruvian Amazon. Genetic analysis is a key component for regional malaria control, offering valuable insights that should be incorporated into routine surveillance.

Population genomic evidence of structured and connected Plasmodium vivax populations under host selection in Latin America.

Pathogen genomic epidemiology has the potential to provide a deep understanding of population dynamics, facilitating strategic planning of interventions, monitoring their impact, and enabling timely responses, and thereby supporting control and elimination efforts of parasitic tropical diseases. Plasmodium vivax, responsible for most malaria cases outside Africa, shows high genetic diversity at the population level, driven by factors like sub-patent infections, a hidden reservoir of hypnozoites, and early transmission to mosquitoes. While Latin America has made significant progress in controlling Plasmodium falciparum, it faces challenges with residual P. vivax. To characterize genetic diversity and population structure and dynamics, we have analyzed the largest collection of P. vivax genomes to date, including 1474 high-quality genomes from 31 countries across Asia, Africa, Oceania, and America. While P. vivax shows high genetic diversity globally, Latin American isolates form a distinctive population, which is further divided into sub-populations and occasional clonal pockets. Genetic diversity within the continent was associated with the intensity of transmission. Population differentiation exists between Central America and the North Coast of South America, vs. the Amazon Basin, with significant gene flow within the Amazon Basin, but limited connectivity between the Northwest Coast and the Amazon Basin. Shared genomic regions in these parasite populations indicate adaptive evolution, particularly in genes related to DNA replication, RNA processing, invasion, and motility - crucial for the parasite's survival in diverse environments. Understanding these population-level adaptations is crucial for effective control efforts, offering insights into potential mechanisms behind drug resistance, immune evasion, and transmission dynamics.

Investigation of an airport-associated cluster of falciparum malaria in Frankfurt, Germany, 2022.

Airport malaria is uncommon but increasing in Europe and often difficult to diagnose. We describe the clinical, epidemiological and environmental investigations of a cluster of airport malaria cases and measures taken in response. Three Frankfurt International Airport employees without travel histories to malaria-endemic areas were diagnosed with Plasmodium falciparum malaria in Germany in 2022. Two cases were diagnosed within 1 week, and the third one after 10 weeks. Two cases had severe disease, all three recovered fully. The cases worked in separate areas and no specific location for the transmissions could be identified. No additional cases were detected among airport employees. In June and July, direct flights from Equatorial Guinea, Nigeria and Angola and one parcel originating in Ghana arrived at Frankfurt airport. No vector-competent mosquitoes could be trapped to identify the source of the outbreak. Whole genome sequencing of P. falciparum genomes showed a high genetic relatedness between samples of the three cases and suggested the geographical origin closest to Ghana. A diagnosis of airport malaria should prompt appropriate and comprehensive outbreak investigations to identify the source and to prevent severe forms of falciparum malaria.

Source Tracing of Leishmania donovani in Emerging Foci of Visceral Leishmaniasis, Western Nepal.

We sequenced Leishmania donovani genomes in blood samples collected in emerging foci of visceral leishmaniasis in western Nepal. We detected lineages very different from the preelimination main parasite population, including a new lineage and a rare one previously reported in eastern Nepal. Our findings underscore the need for genomic surveillance.

Molecular surveillance of Plasmodium falciparum drug-resistance markers in Vietnam using multiplex amplicon sequencing (2000-2016).

Emergence and spread of Plasmodium falciparum resistance to artemisinin-based combination therapies (ACT) is a major challenge for Greater Mekong Subregion countries in their goal to eliminate malaria by 2030. Tools to efficiently monitor drug resistance beyond resource-demanding therapeutic efficacy studies are necessary. A custom multiplex amplicon sequencing assay based on Illumina technology was designed to target the marker of partial resistance to artemisinin (K13), five candidate modulators of artemisinin resistance, the marker of resistance to chloroquine (crt), and four neutral microsatellite loci. The assay was used to genotype 635 P. falciparum-positive blood samples collected across seven provinces of Vietnam and one of Cambodia between 2000 and 2016. Markers of resistance to artemisinin partner-drugs piperaquine (copy number of plasmepsin-2) and mefloquine (copy number of multidrug-resistance 1) were determined by qPCR. Parasite population structure was further assessed using a 101-SNP barcode. Validated mutations of artemisinin partial resistance in K13 were found in 48.1% of samples, first detection was in 2000, and by 2015 prevalence overcame > 50% in Central Highlands and Binh Phuoc province. K13-C580Y variant became predominant country-wide, quickly replacing an outbreak of K13-I543T in Central Highlands. Mutations in candidate artemisinin resistance modulator genes paralleled the trends of K13 mutants, whereas resistance to piperaquine and mefloquine remained low (≈ 10%) by 2015-2016. Genomic tools applied to malaria surveillance generate comprehensive information on dynamics of drug resistance and population structure and reflect drug efficacy profiles from in vivo studies.

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

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

Malaria Molecular Surveillance in the Peruvian Amazon with a Novel Highly Multiplexed Plasmodium falciparum AmpliSeq Assay.

Molecular surveillance for malaria has great potential to support national malaria control programs (NMCPs). To bridge the gap between research and implementation, several applications (use cases) have been identified to align research, technology development, and public health efforts. For implementation at NMCPs, there is an urgent need for feasible and cost-effective tools. We designed a new highly multiplexed deep sequencing assay (Pf AmpliSeq), which is compatible with benchtop sequencers, that allows high-accuracy sequencing with higher coverage and lower cost than whole-genome sequencing (WGS), targeting genomic regions of interest. The novelty of the assay is its high number of targets multiplexed into one easy workflow, combining population genetic markers with 13 nearly full-length resistance genes, which is applicable for many different use cases. We provide the first proof of principle for hrp2 and hrp3 deletion detection using amplicon sequencing. Initial sequence data processing can be performed automatically, and subsequent variant analysis requires minimal bioinformatic skills using any tabulated data analysis program. The assay was validated using a retrospective sample collection (n = 254) from the Peruvian Amazon between 2003 and 2018. By combining phenotypic markers and a within-country 28-single-nucleotide-polymorphism (SNP) barcode, we were able to distinguish different lineages with multiple resistance haplotypes (in dhfr, dhps, crt and mdr1) and hrp2 and hrp3 deletions, which have been increasing in recent years. We found no evidence to suggest the emergence of artemisinin (ART) resistance in Peru. These findings indicate a parasite population that is under drug pressure but is susceptible to current antimalarials and demonstrate the added value of a highly multiplexed molecular tool to inform malaria strategies and surveillance systems. IMPORTANCE While the power of next-generation sequencing technologies to inform and guide malaria control programs has become broadly recognized, the integration of genomic data for operational incorporation into malaria surveillance remains a challenge in most countries where malaria is endemic. The main obstacles include limited infrastructure, limited access to high-throughput sequencing facilities, and the need for local capacity to run an in-country analysis of genomes at a large-enough scale to be informative for surveillance. In addition, there is a lack of standardized laboratory protocols and automated analysis pipelines to generate reproducible and timely results useful for relevant stakeholders. With our standardized laboratory and bioinformatic workflow, malaria genetic surveillance data can be readily generated by surveillance researchers and malaria control programs in countries of endemicity, increasing ownership and ensuring timely results for informed decision- and policy-making.

Novel highly-multiplexed AmpliSeq targeted assay for Plasmodium vivax genetic surveillance use cases at multiple geographical scales.

Although the power of genetic surveillance tools has been acknowledged widely, there is an urgent need in malaria endemic countries for feasible and cost-effective tools to implement in national malaria control programs (NMCPs) that can generate evidence to guide malaria control and elimination strategies, especially in the case of Plasmodium vivax. Several genetic surveillance applications ('use cases') have been identified to align research, technology development, and public health efforts, requiring different types of molecular markers. Here we present a new highly-multiplexed deep sequencing assay (Pv AmpliSeq). The assay targets the 33-SNP vivaxGEN-geo panel for country-level classification, and a newly designed 42-SNP within-country barcode for analysis of parasite dynamics in Vietnam and 11 putative drug resistance genes in a highly multiplexed NGS protocol with easy workflow, applicable for many different genetic surveillance use cases. The Pv AmpliSeq assay was validated using: 1) isolates from travelers and migrants in Belgium, and 2) routine collections of the national malaria control program at sentinel sites in Vietnam. The assay targets 229 amplicons and achieved a high depth of coverage (mean 595.7 ± 481) and high accuracy (mean error-rate of 0.013 ± 0.007). P. vivax parasites could be characterized from dried blood spots with a minimum of 5 parasites/µL and 10% of minority-clones. The assay achieved good spatial specificity for between-country prediction of origin using the 33-SNP vivaxGEN-geo panel that targets rare alleles specific for certain countries and regions. A high resolution for within-country diversity in Vietnam was achieved using the designed 42-SNP within-country barcode that targets common alleles (median MAF 0.34, range 0.01-0.49. Many variants were detected in (putative) drug resistance genes, with different predominant haplotypes in the pvmdr1 and pvcrt genes in different provinces in Vietnam. The capacity of the assay for high resolution identity-by-descent (IBD) analysis was demonstrated and identified a high rate of shared ancestry within Gia Lai Province in the Central Highlands of Vietnam, as well as between the coastal province of Binh Thuan and Lam Dong. Our approach performed well in geographically differentiating isolates at multiple spatial scales, detecting variants in putative resistance genes, and can be easily adjusted to suit the needs in other settings in a country or region. We prioritize making this tool available to researchers and NMCPs in endemic countries to increase ownership and ensure data usage for decision-making and malaria policy.

Two fatal autochthonous cases of airport malaria, Belgium, 2020.

We report an outbreak investigation of two fatal cases of autochthonous Plasmodium falciparum malaria that occurred in Belgium in September 2020. Various hypotheses of the potential source of infection were investigated. The most likely route of transmission was through an infectious exotic Anopheles mosquito that was imported via the international airport of Brussels or the military airport Melsbroek and infected the cases who lived at 5 km from the airports. Based on genomic analysis of the parasites collected from the two cases, the most likely origin of the Plasmodium was Gabon or Cameroon. Further, the parasites collected from the two Belgian patients were identical by descent, which supports the assumption that the two infections originated from the bite of the same mosquito, during interrupted feeding. Although airport malaria remains a rare event, it has significant implications, particularly for the patient, as delayed or missed diagnosis of the cause of illness often results in complications and mortality. Therefore, to prevent such severe or fatal outcomes, we suggest a number of public health actions including increased awareness among health practitioners, especially those working in the vicinity of airports, and increased surveillance of exotic mosquito species at airports.

Genomic and Phenotypic Characterization of Experimentally Selected Resistant Leishmania donovani Reveals a Role for Dynamin-1-Like Protein in the Mechanism of Resistance to a Novel Antileishmanial Compound.

The implementation of prospective drug resistance (DR) studies in the research-and-development (R&D) pipeline is a common practice for many infectious diseases but not for neglected tropical diseases (NTDs). Here, we explored and demonstrated the importance of this approach using as paradigms Leishmania donovani, the etiological agent of visceral leishmaniasis (VL), and TCMDC-143345, a promising compound of the GlaxoSmithKline (GSK) "Leishbox" to treat VL. We experimentally selected resistance to TCMDC-143345 in vitro and characterized resistant parasites at the genomic and phenotypic levels. We found that it took more time to develop resistance to TCMDC-143345 than to other drugs in clinical use and that there was no cross-resistance to these drugs, suggesting a new and unique mechanism. By whole-genome sequencing, we found two mutations in the gene encoding the L. donovani dynamin-1-like protein (LdoDLP1) that were fixed at the highest drug pressure. Through phylogenetic analysis, we identified LdoDLP1 as a family member of the dynamin-related proteins, a group of proteins that impacts the shapes of biological membranes by mediating fusion and fission events, with a putative role in mitochondrial fission. We found that L. donovani lines genetically engineered to harbor the two identified LdoDLP1 mutations were resistant to TCMDC-143345 and displayed altered mitochondrial properties. By homology modeling, we showed how the two LdoDLP1 mutations may influence protein structure and function. Taken together, our data reveal a clear involvement of LdoDLP1 in the adaptation/reduced susceptibility of L. donovani to TCMDC-143345. IMPORTANCE Humans and their pathogens are continuously locked in a molecular arms race during which the eventual emergence of pathogen drug resistance (DR) seems inevitable. For neglected tropical diseases (NTDs), DR is generally studied retrospectively once it has already been established in clinical settings. We previously recommended to keep one step ahead in the host-pathogen arms race and implement prospective DR studies in the R&D pipeline, a common practice for many infectious diseases but not for NTDs. Here, using Leishmania donovani, the etiological agent of visceral leishmaniasis (VL), and TCMDC-143345, a promising compound of the GSK Leishbox to treat VL, as paradigms, we experimentally selected resistance to the compound and proceeded to genomic and phenotypic characterization of DR parasites. The results gathered in the present study suggest a new DR mechanism involving the L. donovani dynamin-1-like protein (LdoDLP1) and demonstrate the practical relevance of prospective DR studies.

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.

Cyclical Patterns Affect Microbial Dynamics in the Water Basin of a Nuclear Research Reactor.

The BR2 nuclear research reactor in Mol, Belgium, runs in successive phases of operation (cycles) and shutdown, whereby a water basin surrounding the reactor vessel undergoes periodic changes in physico-chemical parameters such as flow rate, temperature, and radiation. The aim of this study was to explore the microbial community in this unique environment and to investigate its long-term dynamics using a 16S rRNA amplicon sequencing approach. Results from two sampling campaigns spanning several months showed a clear shift in community profiles: cycles were mostly dominated by two Operational Taxonomic Units (OTUs) assigned to unclassified Gammaproteobacterium and Pelomonas, whereas shutdowns were dominated by an OTU assigned to Methylobacterium. Although 1 year apart, both campaigns showed similar results, indicating that the system remained stable over this 2-year period. The community shifts were linked with changes in physico-chemical parameters by Non-metric Multidimensional Scaling (NMDS) and correlation analyses. In addition, radiation was hypothesized to cause a decrease in cell number, whereas temperature had the opposite effect. Chemoautotrophic use of H2 and dead cell recycling are proposed to be used as a strategies for nutrient retrieval in this extremely oligotrophic environment.

Phenotypic and Genetic Characterization of Temperature-Induced Mutagenesis and Mortality in Cupriavidus metallidurans.

Cupriavidus metallidurans strains display a decreased viability when incubated in rich medium at a temperature of 37°C compared to their normal growth temperature of 30°C, a phenomenon coined "temperature-induced mortality and mutagenesis" (TIMM). To scrutinize this aberrant phenotype further, the contributions of specific inducers and protective agents were determined. Different growth media, including lysogeny broth (LB) and Schatz, and components, including casamino acids, in particular amino acids (proline, cysteine, glycine, glutamine, leucine, histidine and phenylalanine) and ammonium, were found to induce TIMM at 37°C. Sorbitol was found to counteract TIMM. Furthermore, although TIMM is well conserved within the C. metallidurans species, multiple and strain-specific TIMM inducers exist. Twenty-nine percent of the TIMM survivors inherited resistance to TIMM. Whole-genome sequencing of two resistant derivatives revealed an important role of an uncharacterized oxidoreductase, indicating putative metabolic poisoning when grown in high-concentration nitrogen-containing media at 37°C.

Microbiota-based markers predictive of development of Clostridioides difficile infection.

Antibiotic-induced modulation of the intestinal microbiota can lead to Clostridioides difficile infection (CDI), which is associated with considerable morbidity, mortality, and healthcare-costs globally. Therefore, identification of markers predictive of CDI could substantially contribute to guiding therapy and decreasing the infection burden. Here, we analyze the intestinal microbiota of hospitalized patients at increased CDI risk in a prospective, 90-day cohort-study before and after antibiotic treatment and at diarrhea onset. We show that patients developing CDI already exhibit significantly lower diversity before antibiotic treatment and a distinct microbiota enriched in Enterococcus and depleted of Ruminococcus, Blautia, Prevotella and Bifidobacterium compared to non-CDI patients. We find that antibiotic treatment-induced dysbiosis is class-specific with beta-lactams further increasing enterococcal abundance. Our findings, validated in an independent prospective patient cohort developing CDI, can be exploited to enrich for high-risk patients in prospective clinical trials, and to develop predictive microbiota-based diagnostics for management of patients at risk for CDI.

Adaptation of Cupriavidus metallidurans CH34 to Toxic Zinc Concentrations Involves an Uncharacterized ABC-Type Transporter.

Cupriavidus metallidurans CH34 is a well-studied metal-resistant ß-proteobacterium and contains a battery of genes participating in metal metabolism and resistance. Here, we generated a mutant (CH34ZnR) adapted to high zinc concentrations in order to study how CH34 could adaptively further increase its resistance against this metal. Characterization of CH34ZnR revealed that it was also more resistant to cadmium, and that it incurred seven insertion sequence-mediated mutations. Among these, an IS1088 disruption of the glpR gene (encoding a DeoR-type transcriptional repressor) resulted in the constitutive expression of the neighboring ATP-binding cassette (ABC)-type transporter. GlpR and the adjacent ABC transporter are highly similar to the glycerol operon regulator and ATP-driven glycerol importer of Rhizobium leguminosarum bv. viciae VF39, respectively. Deletion of glpR or the ABC transporter and complementation of CH34ZnR with the parental glpR gene further demonstrated that loss of GlpR function and concomitant derepression of the adjacent ABC transporter is pivotal for the observed resistance phenotype. Importantly, addition of glycerol, presumably by glycerol-mediated attenuation of GlpR activity, also promoted increased zinc and cadmium resistance in the parental CH34 strain. Upregulation of this ABC-type transporter is therefore proposed as a new adaptation route towards metal resistance.

Residential exposure to air pollution and access to neighborhood greenspace in relation to hair cortisol concentrations during the second and third trimester of pregnancy.

BACKGROUND: Exposure to air pollution during pregnancy has been associated with adverse pregnancy outcomes in studies worldwide, other studies have described beneficial effects of residential greenspace on pregnancy outcomes. The biological mechanisms that underlie these associations are incompletely understood. A biological stress response, which implies release of cortisol, may underlie associations of air pollution exposure and access to neighborhood greenspaces with health. METHODS: We explored residential exposure to air pollution and residential access to neighborhood greenspaces in relation to hair cortisol concentrations of participants in a prospective pregnancy cohort study in Flanders, Belgium. Hair samples were collected at the end of the second pregnancy trimester (n = 133) and shortly after delivery (n = 81). Cortisol concentrations were measured in 3-cm scalp-near hair sections, to reflect second and third pregnancy trimester cortisol secretion. We estimated long-term (3 months before sampling) residential exposure to fine particulate matter (PM2.5), nitrogen dioxide (NO2) and black carbon (BC), assessed residential distance to major roads and residential access to neighborhood greenspaces (NHGS). Associations between residential exposures and hair cortisol concentrations were studied using linear regression models while adjusting for season of sampling. RESULTS: Three-month mean residential NO2 and BC concentrations were positively associated with third pregnancy trimester hair cortisol concentrations (p = 0.008 and p = 0.017). Access to a large NHGS (10 ha or more within 800 m from residence) was negatively associated with third trimester hair cortisol concentrations (p = 0.019). Access to a large NHGS significantly moderated the association between residential proximity to major roads and second trimester hair cortisol concentrations (p = 0.021). Residential distance to major roads was negatively associated with second trimester hair cortisol concentrations of participants without access to a large NHGS (p = 0.003). The association was not significant for participants with access to a large NHGS. The moderation tended towards significance in the third pregnancy trimester (p < 0.10). CONCLUSIONS: Our findings suggest a positive association between long-term residential exposure to air pollution and biological stress during pregnancy, residential access to neighborhood greenspaces may moderate the association. Further research is needed to confirm our results. TRIAL REGISTRATION: The IPANEMA study is registered under number  NCT02592005 at clinicaltrials.gov .

Direct Meta-Analyses Reveal Unexpected Microbial Life in the Highly Radioactive Water of an Operating Nuclear Reactor Core.

The pools of nuclear reactor facilities constitute harsh environments for life, bathed with ionizing radiation, filled with demineralized water and containing toxic radioactive elements. The very few studies published to date have explored water pools used to store spent nuclear fuels. Due to access restrictions and strong handling constraints related to the high radioactivity level, nothing is presently known about life in water pools that directly cool nuclear cores. In this work, we investigated the microbial communities in the cooling pool of the French Osiris nuclear reactor using direct meta-omics approaches, namely, DNA metabarcoding and proteotyping based on 16S ribosomal RNA gene sequencing and on peptide analysis, respectively. We identified 25 genera in the highly radioactive core water supply during operation with radionuclide activity higher than 3 × 109 Bq/m3. The prevailing genera Variovorax and Sphingomonas at operation were supplanted by Methylobacterium, Asanoa, and Streptomyces during shutdown. Variovorax might use dihydrogen produced by water radiolysis as an energy source.

Evaluation of whole genome amplification and bioinformatic methods for the characterization of Leishmania genomes at a single cell level.

Here, we report a pilot study paving the way for further single cell genomics studies in Leishmania. First, the performances of two commercially available kits for Whole Genome Amplification (WGA), PicoPLEX and RepliG were compared on small amounts of Leishmania donovani DNA, testing their ability to preserve specific genetic variations, including aneuploidy levels and SNPs. We show here that the choice of WGA method should be determined by the planned downstream genetic analysis, PicoPLEX and RepliG performing better for aneuploidy and SNP calling, respectively. This comparison allowed us to evaluate and optimize corresponding bio-informatic methods. As PicoPLEX was shown to be the preferred method for studying single cell aneuploidy, this method was applied in a second step, on single cells of L. braziliensis, which were sorted by fluorescence activated cell sorting (FACS). Even sequencing depth was achieved in 28 single cells, allowing accurate somy estimation. A dominant karyotype with three aneuploid chromosomes was observed in 25 cells, while two different minor karyotypes were observed in the other cells. Our method thus allowed the detection of aneuploidy mosaicism, and provides a solid basis which can be further refined to concur with higher-throughput single cell genomic methods.