Genome wide comparison of Ethiopian Leishmania donovani strains reveals differences potentially related to parasite survival.

PLoS Genetics. 2018;14(1):e1007133 DOI 10.1371/journal.pgen.1007133


Journal Homepage

Journal Title: PLoS Genetics

ISSN: 1553-7390 (Print); 1553-7404 (Online)

Publisher: Public Library of Science (PLoS)

LCC Subject Category: Science: Biology (General): Genetics

Country of publisher: United States

Language of fulltext: English

Full-text formats available: PDF, HTML, XML



Arie Zackay

James A Cotton

Mandy Sanders

Asrat Hailu

Abedelmajeed Nasereddin

Alon Warburg

Charles L Jaffe


Peer review

Editorial Board

Instructions for authors

Time From Submission to Publication: 26 weeks


Abstract | Full Text

Leishmania donovani is the main cause of visceral leishmaniasis (VL) in East Africa. Differences between northern Ethiopia/Sudan (NE) and southern Ethiopia (SE) in ecology, vectors, and patient sensitivity to drug treatment have been described, however the relationship between differences in parasite genotype between these two foci and phenotype is unknown. Whole genomic sequencing (WGS) was carried out for 41 L. donovani strains and clones from VL and VL/HIV co-infected patients in NE (n = 28) and SE (n = 13). Chromosome aneuploidy was observed in all parasites examined with each isolate exhibiting a unique karyotype. Differences in chromosome ploidy or karyotype were not correlated with the geographic origin of the parasites. However, correlation between single nucleotide polymorphism (SNP) and geographic origin was seen for 38/41 isolates, separating the NE and SE parasites into two large groups. SNP restricted to NE and SE groups were associated with genes involved in viability and parasite resistance to drugs. Unique copy number variation (CNV) were also associated with NE and SE parasites, respectively. One striking example is the folate transporter (FT) family genes (LdBPK_100390, LdBPK_100400 and LdBPK_100410) on chromosome 10 that are single copy in all 13 SE isolates, but either double copy or higher in 39/41 NE isolates (copy number 2-4). High copy number (= 4) was also found for one Sudanese strain examined. This was confirmed by quantitative polymerase chain reaction for LdBPK_100400, the L. donovani FT1 transporter homolog. Good correlation (p = 0.005) between FT copy number and resistance to methotrexate (0.5 mg/ml MTX) was also observed with the haploid SE strains examined showing higher viability than the NE strains at this concentration. Our results emphasize the advantages of whole genome analysis to shed light on vital parasite processes in Leishmania.