Genomic and Metabolomic Polymorphism among Experimentally Selected Paromomycin-Resistant Leishmania donovani Strains.

Understanding the mechanism(s) underpinning drug resistance could lead to novel treatments to reverse the increased tolerance of a pathogen. In this study, paromomycin (PMM) resistance (PMMr) was induced in three Nepalese clinical strains of Leishmania donovani with different inherent susceptibilities to antimony (Sb) drugs by stepwise exposure of promastigotes to PMM. Exposure to PMM resulted in the production of mixed populations of parasites, even though a single cloned population was used at the start of selection. PMM 50% inhibitory concentration (IC50) values for PMMr parasites varied between 104 and 481 µM at the promastigote stage and 32 and 195 µM at the intracellular amastigote stage. PMM resistance was associated with increased resistance to nitric oxide at the amastigote stage but not the promastigote stage (P < 0.05). This effect was most marked in the Sb-resistant (Sbr) PMMr clone, in which PMM resistance was associated with a significant upregulation of glutathione compared to that in its wild type (P < 0.05), although there was no change in the regulation of trypanothione (detected in its oxidized form). Interestingly, PMMr strains showed an increase in either the keto acid derivative of isoleucine (Sb intermediate PMMr) or the 2-hydroxy acids derived from arginine and tyrosine (Sb susceptible PMMr and Sbr PMMr). These results are consistent with the recent finding that the upregulation of the branched-chain amino acid aminotransferase and d-lactate dehydrogenase is linked to PMMr In addition, we found that PMMr is associated with a significant increase in aneuploidy during PMM selection in all the strains, which could allow the rapid selection of genetic changes that confer a survival advantage.

The global diversity of Haemonchus contortus is shaped by human intervention and climate.

Haemonchus contortus is a haematophagous parasitic nematode of veterinary interest. We have performed a survey of its genome-wide diversity using single-worm whole genome sequencing of 223 individuals sampled from 19 isolates spanning five continents. We find an African origin for the species, together with evidence for parasites spreading during the transatlantic slave trade and colonisation of Australia. Strong selective sweeps surrounding the ß-tubulin locus, a target of benzimidazole anthelmintic drug, are identified in independent populations. These sweeps are further supported by signals of diversifying selection enriched in genes involved in response to drugs and other anthelmintic-associated biological functions. We also identify some candidate genes that may play a role in ivermectin resistance. Finally, genetic signatures of climate-driven adaptation are described, revealing a gene acting as an epigenetic regulator and components of the dauer pathway. These results begin to define genetic adaptation to climate in a parasitic nematode.

Modulation of Aneuploidy in Leishmania donovani during Adaptation to Different In Vitro and In Vivo Environments and Its Impact on Gene Expression.

Aneuploidy is usually deleterious in multicellular organisms but appears to be tolerated and potentially beneficial in unicellular organisms, including pathogens. Leishmania, a major protozoan parasite, is emerging as a new model for aneuploidy, since in vitro-cultivated strains are highly aneuploid, with interstrain diversity and intrastrain mosaicism. The alternation of two life stages in different environments (extracellular promastigotes and intracellular amastigotes) offers a unique opportunity to study the impact of environment on aneuploidy and gene expression. We sequenced the whole genomes and transcriptomes of Leishmania donovani strains throughout their adaptation to in vivo conditions mimicking natural vertebrate and invertebrate host environments. The nucleotide sequences were almost unchanged within a strain, in contrast to highly variable aneuploidy. Although high in promastigotes in vitro, aneuploidy dropped significantly in hamster amastigotes, in a progressive and strain-specific manner, accompanied by the emergence of new polysomies. After a passage through a sand fly, smaller yet consistent karyotype changes were detected. Changes in chromosome copy numbers were correlated with the corresponding transcript levels, but additional aneuploidy-independent regulation of gene expression was observed. This affected stage-specific gene expression, downregulation of the entire chromosome 31, and upregulation of gene arrays on chromosomes 5 and 8. Aneuploidy changes in Leishmania are probably adaptive and exploited to modulate the dosage and expression of specific genes; they are well tolerated, but additional mechanisms may exist to regulate the transcript levels of other genes located on aneuploid chromosomes. Our model should allow studies of the impact of aneuploidy on molecular adaptations and cellular fitness.IMPORTANCE Aneuploidy is usually detrimental in multicellular organisms, but in several microorganisms, it can be tolerated and even beneficial. Leishmania-a protozoan parasite that kills more than 30,000 people each year-is emerging as a new model for aneuploidy studies, as unexpectedly high levels of aneuploidy are found in clinical isolates. Leishmania lacks classical regulation of transcription at initiation through promoters, so aneuploidy could represent a major adaptive strategy of this parasite to modulate gene dosage in response to stressful environments. For the first time, we document the dynamics of aneuploidy throughout the life cycle of the parasite, in vitro and in vivo We show its adaptive impact on transcription and its interaction with regulation. Besides offering a new model for aneuploidy studies, we show that further genomic studies should be done directly in clinical samples without parasite isolation and that adequate methods should be developed for this.

Haemonchus contortus: Genome Structure, Organization and Comparative Genomics.

One of the first genome sequencing projects for a parasitic nematode was that for Haemonchus contortus. The open access data from the Wellcome Trust Sanger Institute provided a valuable early resource for the research community, particularly for the identification of specific genes and genetic markers. Later, a second sequencing project was initiated by the University of Melbourne, and the two draft genome sequences for H. contortus were published back-to-back in 2013. There is a pressing need for long-range genomic information for genetic mapping, population genetics and functional genomic studies, so we are continuing to improve the Wellcome Trust Sanger Institute assembly to provide a finished reference genome for H. contortus. This review describes this process, compares the H. contortus genome assemblies with draft genomes from other members of the strongylid group and discusses future directions for parasite genomics using the H. contortus model.

In vitro selection of miltefosine resistance in promastigotes of Leishmania donovani from Nepal: genomic and metabolomic characterization.

In this study, we followed the genomic, lipidomic and metabolomic changes associated with the selection of miltefosine (MIL) resistance in two clinically derived Leishmania donovani strains with different inherent resistance to antimonial drugs (antimony sensitive strain Sb-S; and antimony resistant Sb-R). MIL-R was easily induced in both strains using the promastigote-stage, but a significant increase in MIL-R in the intracellular amastigote compared to the corresponding wild-type did not occur until promastigotes had adapted to 12.2 µM MIL. A variety of common and strain-specific genetic changes were discovered in MIL-adapted parasites, including deletions at the LdMT transporter gene, single-base mutations and changes in somy. The most obvious lipid changes in MIL-R promastigotes occurred to phosphatidylcholines and lysophosphatidylcholines and results indicate that the Kennedy pathway is involved in MIL resistance. The inherent Sb resistance of the parasite had an impact on the changes that occurred in MIL-R parasites, with more genetic changes occurring in Sb-R compared with Sb-S parasites. Initial interpretation of the changes identified in this study does not support synergies with Sb-R in the mechanisms of MIL resistance, though this requires an enhanced understanding of the parasite's biochemical pathways and how they are genetically regulated to be verified fully.

Parallel signatures of sequence evolution among hearing genes in echolocating mammals: an emerging model of genetic convergence.

Recent findings of sequence convergence in the Prestin gene among some bats and cetaceans suggest that parallel adaptations for high-frequency hearing have taken place during the evolution of echolocation. To determine if this gene is an exception, or instead similar processes have occurred in other hearing genes, we have examined Tmc1 and Pjvk, both of which are associated with non-syndromic hearing loss in mammals. These genes were amplified and sequenced from a number of mammalian species, including echolocating and non-echolocating bats and whales, and were analysed together with published sequences. Sections of both genes showed phylogenetic signals that conflicted with accepted species relationships, with coding regions uniting laryngeal echolocating bats in a monophyletic clade. Bayesian estimates of posterior probabilities of convergent and divergent substitutions provided more direct evidence of sequence convergence between the two groups of laryngeal echolocating bats as well as between echolocating bats and dolphins. We found strong evidence of positive selection acting on some echolocating bat species and echolocating cetaceans, contrasting with purifying selection on non-echolocating bats. Signatures of sequence convergence and molecular adaptation in two additional hearing genes suggest that the acquisition of high-frequency hearing has involved multiple loci.

Vertebrate phylogenomics: reconciled trees and gene duplications.

Ancient gene duplication events have left many traces in vertebrate genomes. Reconciled trees represent the differences between gene family trees and the species phylogeny those genes are sampled from, allowing us to both infer gene duplication events and estimate a species phylogeny from a sample of gene families. We show that analysis of 118 gene families yields a phylogeny of vertebrates largely in agreement with other data. We formulate the problem of locating episodes of gene duplication as a set cover problem: given a species tree in which each node has a set of gene duplications associated with it, the smallest set of species nodes whose union includes all gene duplications specifies the locations of gene duplication episodes. By generating a unique mapping from this cover set we can determine the minimal number of such episodes at each location. When applied to our data, this method reveals a complex history of gene duplications in vertebrate evolution that does not conform to the "2R" hypothesis.