Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2011 Dec;21(12):2143-56.
doi: 10.1101/gr.123430.111. Epub 2011 Oct 28.

Whole genome sequencing of multiple Leishmania donovani clinical isolates provides insights into population structure and mechanisms of drug resistance

Tim Downing  1 , Hideo Imamura, Saskia Decuypere, Taane G Clark, Graham H Coombs, James A Cotton, James D Hilley, Simonne de Doncker, Ilse Maes, Jeremy C Mottram, Mike A Quail, Suman Rijal, Mandy Sanders, Gabriele Schönian, Olivia Stark, Shyam Sundar, Manu Vanaerschot, Christiane Hertz-Fowler, Jean-Claude Dujardin, Matthew Berriman
Affiliations

Whole genome sequencing of multiple Leishmania donovani clinical isolates provides insights into population structure and mechanisms of drug resistance

Tim Downing et al. Genome Res. 2011 Dec.

Abstract

Visceral leishmaniasis is a potentially fatal disease endemic to large parts of Asia and Africa, primarily caused by the protozoan parasite Leishmania donovani. Here, we report a high-quality reference genome sequence for a strain of L. donovani from Nepal, and use this sequence to study variation in a set of 16 related clinical lines, isolated from visceral leishmaniasis patients from the same region, which also differ in their response to in vitro drug susceptibility. We show that whole-genome sequence data reveals genetic structure within these lines not shown by multilocus typing, and suggests that drug resistance has emerged multiple times in this closely related set of lines. Sequence comparisons with other Leishmania species and analysis of single-nucleotide diversity within our sample showed evidence of selection acting in a range of surface- and transport-related genes, including genes associated with drug resistance. Against a background of relative genetic homogeneity, we found extensive variation in chromosome copy number between our lines. Other forms of structural variation were significantly associated with drug resistance, notably including gene dosage and the copy number of an experimentally verified circular episome present in all lines and described here for the first time. This study provides a basis for more powerful molecular profiling of visceral leishmaniasis, providing additional power to track the drug resistance and epidemiology of an important human pathogen.

PubMed Disclaimer

Figures

Figure 1.
Figure 1.
The phylogenomic context of functional variation in clinical L. donovani lines. (A) A neighbor-joining phylogenomic tree of Leishmania species. The relative similarity is shown (1% was equivalent to 108,000 mutations). The first column indicates the fraction of the known genomes orthologously aligned to L. donovani and the second is the number of SNPs identified (Table 2). (B) Principle component analysis (PCA) of genomic SNP variation in the 17 L. donovani lines resistant (red) and susceptible (blue) to SGG. The two most significant PCs distinguished three main groups, accounting for 63.3% of the total variation. The PCA L. donovani line numbers correspond to those shown in C, a median-joining phylogenetic network of genome-wide nonsynonymous sites variation (396 SNPs) for samples resistant (red) and susceptible (blue) to SSG from India (beginning with BHU) and Nepal (BPK). The branch lengths displayed are proportional to the number of differences between lines. The position of the ancestral node for L. infantum (black) was shortened: For comparison, the L. donovani:L. infantum branch length was 138 times that between the most divergent L. donovani (1, BHU568/0cl1; and 10, BPK035/0cl1).
Figure 2.
Figure 2.
The neutrality of coding-site allele-frequency spectra. The allele-frequency spectrum for SNPs at nonsynonymous (black line) and synonymous (dashed line) sites among the 17 samples for the genome. Allele frequencies were computed using read-depth coverage values and were folded because the selective process may act on both the ancestral alleles inferred from L. infantum as well as the derived version with the L. donovani population. The ratio of nonsynonymous to synonymous SNPs approaches one as the allele frequencies increase, suggesting that higher frequency variants are subject to stronger purifying selection to remove deleterious alleles.
Figure 3.
Figure 3.
Aneuploidy in natural populations of L. donovani. The heatmap shows the copy-number status of the 36 chromosomes for the 17 lines as disomic (1, blue), trisomic (1.5, green), tetrasomic (2, yellow), and pentasomic (2.5, ivory). The color key shows the normalized chromosome read-depth and the distribution frequency.
Figure 4.
Figure 4.
Copy-number variation spanning MAPK-locus on chromosome 36. Copy number per cell in L. donovani (red, SSG-R; blue, SSG-S) was much higher for the MAPK locus than the chromosome 36 median for L. donovani (black dashed line) and L. infantum (black). The amplified region is marked in gray and their white edges indicate the location of two direct repeats. Genes within and adjacent to the amplified region are given in black and light gray, respectively. The MAPK (mitogen-activated protein kinase) gene is labeled as LMPK (LdBPK_366760). Recombination events may occur through a highly homologous 220-bp region starting at the start codon of the tartrate-sensitive acid phosphatase (LdBPK_366740) and histidine secretory acid phosphatase (LdBPK_366770). LdBPK_366750 encodes a protein of unknown function. Two RIME (L. donovani ribosomal mobile elements; light-blue) flank the duplicated region. Read-depth values were computed using 1-kb sliding windows, and error bars are shown for each kilobase.
See this image and copyright information in PMC

Similar articles

See all similar articles

Cited by

See all "Cited by" articles

References

    1. Adaui V, Castillo D, Zimic M, Gutiérrez A, Decuypere S, Vanaerschot M, De Doncker S, Llanos-Cuentas A, Arévalo J, Dujardin JC 2011. Comparison of gene expression patterns among Leishmania braziliensis clinical isolates showing a different in vitro susceptibility to pentavalent antimony. Parasitology 138: 183–193 - PubMed
    1. Alam MZ, Kuhls K, Schweynoch C, Sundar S, Rijal S, Shamsuzzaman AK, Raju BV, Salotra P, Dujardin JC, Schönian G 2009. Multilocus microsatellite typing (MLMT) reveals genetic homogeneity of Leishmania donovani strains in the Indian subcontinent. Infect Genet Evol 9: 24–31 - PubMed
    1. Alcolea PJ, Alonso A, Gómez MJ, Sánchez-Gorostiaga A, Moreno-Paz M, González-Pastor E, Toraño A, Parro V, Larraga V 2010. Temperature increase prevails over acidification in gene expression modulation of amastigote differentiation in Leishmania infantum. BMC Genomics 11: 31 doi: 10.1186/1471-2164-11-31 - PMC - PubMed
    1. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ 1990. Basic local alignment search tool. J Mol Biol 215: 403–410 - PubMed
    1. Assefa S, Keane TM, Otto TD, Newbold C, Berriman M 2009. ABACAS: algorithm-based automatic contiguation of assembled sequences. Bioinformatics 25: 1968–1969 - PMC - PubMed

Publication types

Associated data