Phylogenomic reconstruction supports supercontinent origins for Leishmania

Mucosal leishmaniasis is associated with the Leishmania RNA virus and inappropriate cutaneous leishmaniasis treatment

BACKGROUND

Mucosal leishmaniasis (ML) is a severe clinical form of leishmaniasis that is characterized by the destruction of the nasal and/or the oral mucosae and appears as a late complication in 5% to 10% of cutaneous leishmaniasis (CL) cases produced by species belonging to Leishmania (Viannia) subgenus. Some strains of Leishmania spp. carry an RNA virus known as Leishmania RNA virus (LRV) that may contribute to the appearance of ML.

METHODS

To examine the role of LRV type 1 (LRV1) as a risk factor associated with ML, a retrospective case-control study involving 103 patients was conducted. Cases were defined as patients with ML (n = 33), and controls corresponded to patients with CL and without mucosal lesions (n = 70). Clinical data were recorded from the patient's medical records. Cryopreserved biopsies were used to detect LRV1 and identify Leishmania species.

RESULTS

The frequency of LRV1 in the 103 patients was 16.5% (95% CI,10.4-25.12) being higher in samples from cases [33.33% (95% CI,18.89-51.76) than from controls [8.57% (95% CI, 3.82-18.10)]. L. (V.) braziliensis was identified in 63.6% of cases and 55.7% of the controls. Multivariate logistic regression indicated that infection with Leishmania spp. carrying LRV1 (OR = 6.30; 95% CI,1.52-26.10, p = 0.011) acts as risk factors for ML occurrence, while the completed treatment for the cutaneous event decreases the risk of ML (OR = 0.039; 95% CI, 0.01-0.12, p < 0.0001).

CONCLUSIONS

Our data support the association between LRV1 and ML occurrence and emphasize the effect of completed treatment for CL in preventing ML.

Development of species-specific multiplex PCR for Leishmania identification

Leishmaniasis is a diverse group of clinical diseases caused by protozoan parasites of the Leishmania genus. Species-specific identification of Leishmania spp. is challenging due to the high number of different pathogenic species that sometimes co-circulate in the same foci, hampering efforts to effectively control the disease. Multiplex PCR is an attractive alternative for rapid differentiation of Leishmania species with high sensitivity and specificity. We aimed to generate a panel of primers optimized for a multiplex PCR assay capable of identifying different Leishmania species in a single reaction. Species-specific primers were designed based on genomic data using the TipMT tooL. Potential non-specific amplifications of other trypanosomatids as well as human, dog, and sandfly hosts were first evaluated in silico using the Primer-Blast tooL. Species-specific primers for Leishmania amazonensis, Leishmania braziliensis, Leishmania donovani, Leishmania infantum, Leishmania mexicana and for the Leishmania guyanensis complex were tested in vitro. The primers have a limit of detection ranging from 1 to 0.01 ng of parasite gDNA using the same annealing temperature of 66 °C. The primers were specific for their targets when tested against 13 species of Leishmania, six trypanosomatids, and Babesia sp., and to detect the target species in a prepared pool with gDNA of six pathogenic Leishmania species. The designed primers were optimized for multiplex PCR, enabling species-specific identification of all five Leishmania species and one species complex. This new primer set could allow for efficient, fast, and reliable identification of Leishmania parasites.

Antigenic reactivity of Leishmania (Viannia) lainsoni axenic amastigote proved to be a suitable alternative for optimizing Montenegro skin test

BACKGROUND

Laboratory diagnosis of American cutaneous leishmaniasis (ACL) requires a tool amenable to the epidemiological status of ACL in Brazil. Montenegro skin test (MST), an efficient immunological tool used for laboratory diagnosis of ACL, induces delayed-type hypersensitivity (DTH) response to the promastigote antigens of Leishmania; however, human immune responses against infection are modulated by the amastigote of the parasite. Leishmania (V.) lainsoni induces strong cellular immunity in humans; therefore, the antigenic reactivity of its axenic amastigote (AMA antigen) to MST was evaluated for the laboratory diagnosis of ACL.

METHODS

Among 70 individuals examined, 60 had a laboratory-confirmed diagnosis of ACL; 53 had localized cutaneous leishmaniasis (LCL), and 7 had mucosal leishmaniasis (ML). Patients were treated at the Evandro Chagas Institute's leishmaniasis clinic, Pará State, Brazil. Ten healthy individuals with no history of ACL (control group) were also examined. Leishmania (V.) braziliensis promastigote antigen (PRO) was used to compare the reactivity with that of AMA antigen. Paired Student's t-test, kappa agreement, and Spearman test were used to evaluate the reactivity of AMA and PRO.

RESULTS

The mean reactivity of AMA in ACL patients was 19.4 mm ± 13.3, which was higher (P < 0.001) than that of PRO: 12.1 mm ± 8.1. MST reactivity according to the clinical forms revealed that AMA reactivity in LCL and ML, 18.8 mm ± 13.3 and 24.3 mm ± 13.7, was higher (P < 0.001) than that of PRO, 11.8 mm ± 8.2 and 14.6 mm ± 8.4, respectively.

CONCLUSION

AMA reactivity was higher than that of PRO, indicating that AMA is a promising alternative for optimizing MST in the laboratory diagnosis of ACL.

Extensive transcriptome data providing great efficacy in genetic research and adaptive gene discovery: a case study of Elymus sibiricus L. (Poaceae, Triticeae)

Genetic markers play a central role in understanding genetic diversity, speciation, evolutionary processes, and how species respond to environmental stresses. However, conventional molecular markers are less effective when studying polyploid species with large genomes. In this study, we compared gene expression levels in 101 accessions of Elymus sibiricus, a widely distributed allotetraploid forage species across the Eurasian continent. A total of 20,273 high quality transcriptomic SNPs were identified. In addition, 72,344 evolutionary information loci of these accessions of E. sibiricus were identified using genome skimming data in conjunction with the assembled composite genome. The population structure results suggest that transcriptome SNPs were more effective than SNPs derived from genome skimming data in revealing the population structure of E. sibiricus from different locations, and also outperformed gene expression levels. Compared with transcriptome SNPs, the investigation of population-specifically-expressed genes (PSEGs) using expression levels revealed a larger number of locally adapted genes mainly involved in the ion response process in the Sichuan, Inner Mongolia, and Xizang geographical groups. Furthermore, we performed the weighted gene co-expression network analysis (WGCNA) and successfully identified potential regulators of PSEGs. Therefore, for species lacking genomic information, the use of transcriptome SNPs is an efficient approach to perform population structure analysis. In addition, analyzing genes under selection through nucleotide diversity and genetic differentiation index analysis based on transcriptome SNPs, and exploring PSEG through expression levels is an effective method for analyzing locally adaptive genes.

Seasonal dynamics, Leishmania diversity, and nanopore-based metabarcoding of blood meal origins in Culicoides spp. in the newly emerging focus of leishmaniasis in Northern Thailand

BACKGROUND

Clinical cases of leishmaniasis caused by Leishmania (Mundinia) parasites have been increasingly reported in Southeast Asia, particularly Thailand. Recent evidence has shown that Leishmania (Mundinia) parasites successfully developed into infective metacyclic promastigotes in Culicoides biting midges, strongly supporting their putative role in disease transmission. However, Culicoides diversity, host preference, and Leishmania prevalence in endemic areas remain largely unknown.

METHODS

We investigated the seasonal dynamics, infection prevalence, and blood meal identification of Culicoides collected from the emerging focus of visceral leishmaniasis in Lampang Province, Northern Thailand, during 2021-2023. Midge samples were molecularly screened for Leishmania using SSU rRNA-qPCR and ITS1-PCR, followed by Sanger plasmid sequencing, and parasite haplotype diversity was analyzed. Host blood meal origins were comparatively identified using host-specific Cytb-PCRs and a nanopore-based metabarcoding approach.

RESULTS

A total of 501 parous and gravid females and 46 blood-engorged ones belonging to at least 17 species of five subgenera (Remmia, Trithecoides, Avaritia, Hoffmania, and Meijerehelea) and two species groups (Shortti and Calvipalpis) were collected with temporal differences in abundance. Leishmania was detected by SSU rRNA-qPCR in 31 samples of at least 11 midge species, consisting of Culicoides oxystoma, C. guttifer, C. orientalis, C. mahasarakhamense, C (Trithecoides) spp., C. innoxius, C. shortti, C. arakawae, C. sumatrae, C. actoni, and C. fulvus, with the overall infection prevalence of 5.7%. The latter six species represent the new records as putative leishmaniasis vectors in Northern Thailand. The ITS1-PCR and plasmid sequencing revealed that Leishmania martiniquensis was predominantly identified in all qPCR-positive species, whereas L. orientalis was identified only in three C. oxystoma samples. The most dominant haplotype of L. martiniquensis in Thailand was genetically intermixed with those from other geographical regions, confirming its globalization. Neutrality test statistics were also significantly negative on regional and country-wide scales, suggesting rapid population expansion or selective sweeps. Nanopore-based blood meal analysis revealed that most Culicoides species are mammalophilic, with peridomestic and wild mammals (cow, pig, deer, and goat-like species) and humans as hosts, while C. guttifer and C. mahasarakhamense fed preferentially on chickens.

CONCLUSIONS

This study revealed seasonal dynamics and sympatric circulation of L. martiniquensis and L. orientalis in different species of Culicoides. Evidence of human blood feeding was also demonstrated, implicating Culicoides as putative vectors of human leishmaniasis in endemic areas. Further research is therefore urgently needed to develop vector control strategies and assess the infection status of their reservoir hosts to effectively minimize disease transmission.

Genetic diversity and epidemiological insights into cutaneous leishmaniasis in Pakistan: a comprehensive study on clinical manifestations and molecular characterization of Leishmania species

Cutaneous leishmaniasis (CL) stands out as a significant vector-borne endemic in Pakistan. Despite the rising incidence of CL, the genetic diversity of Leishmania species in the country's endemic regions remains insufficiently explored. This study aims to uncover the genetic diversity and molecular characteristics of Leishmania species in CL-endemic areas of Baluchistan, Khyber Pakhtunkhwa (KPK), and Punjab in Pakistan. Clinical samples from 300 CL patients were put to microscopic examination, real-time ITS-1 PCR, and sequencing. Predominantly affecting males between 16 to 30 years of age, with lesions primarily on hands and faces, the majority presented with nodular and plaque types. Microscopic analysis revealed a positivity rate of 67.8%, while real-time PCR identified 60.98% positive cases, mainly L. tropica, followed by L. infantum and L. major. Leishmania major (p = 0.009) showed substantially greater variation in nucleotide sequences than L. tropica (p = 0.07) and L. infantum (p = 0.03). Nucleotide diversity analysis indicated higher diversity in L. major and L. infantum compared to L. tropica. This study enhances our understanding of CL epidemiology in Pakistan, stressing the crucial role of molecular techniques in accurate species identification. The foundational data provided here emphasizes the necessity for future research to investigate deeper into genetic diversity and its implications for CL control at both individual and community levels.

Genomes of Endotrypanum monterogeii from Panama and Zelonia costaricensis from Brazil: Expansion of Multigene Families in Leishmaniinae Parasites That Are Close Relatives of Leishmania spp

The Leishmaniinae subfamily of the Trypanosomatidae contains both genus Zelonia (monoxenous) and Endotrypanum (dixenous). They are amongst the nearest known relatives of Leishmania, which comprises many human pathogens widespread in the developing world. These closely related lineages are models for the genomic biology of monoxenous and dixenous parasites. Herein, we used comparative genomics to identify the orthologous groups (OGs) shared among 26 Leishmaniinae species to investigate gene family expansion/contraction and applied two phylogenomic approaches to confirm relationships within the subfamily. The Endotrypanum monterogeii and Zelonia costaricensis genomes were assembled, with sizes of 29.9 Mb and 38.0 Mb and 9.711 and 12.201 predicted protein-coding genes, respectively. The genome of E. monterogeii displayed a higher number of multicopy cell surface protein families, including glycoprotein 63 and glycoprotein 46, compared to Leishmania spp. The genome of Z. costaricensis presents expansions of BT1 and amino acid transporters and proteins containing leucine-rich repeat domains, as well as a loss of ABC-type transporters. In total, 415 and 85 lineage-specific OGs were identified in Z. costaricensis and E. monterogeii. The evolutionary relationships within the subfamily were confirmed using the supermatrix (3384 protein-coding genes) and supertree methods. Overall, this study showed new expansions of multigene families in monoxenous and dixenous parasites of the subfamily Leishmaniinae.

First evidence of experimental genetic hybridization between cutaneous and visceral strains of Leishmania donovani within its natural vector Phlebotomus argentipes

Leishmaniasis is a neglected tropical disease caused by protozoan parasites of genus Leishmania, and transmitted by different species of Phlebotomine sand flies. More than 20 species of Leishmania are known to cause disease in humans and other animals. Leishmania donovani species complex is known to have a vast diversity of clinical manifestations in humans, but underlying mechanisms for such diversity are yet unknown. Long believed to be strictly asexual, Leishmania have been shown to undergo a cryptic sexual cycle inside its sandfly vector. Natural populations of hybrid parasites have been associated with the rise of atypical clinical outcomes in the Indian subcontinent (ISC). However, formal demonstration of genetic crossing in the major endemic sandfly species in the ISC remain unexplored. Here, we investigated the ability of two distinct variants of L. donovani associated with strikingly different forms of the disease to undergo genetic exchange inside its natural vector, Phlebotomus argentipes. Clinical isolates of L. donovani either from a Sri Lankan cutaneous leishmaniasis (CL) patient or an Indian visceral leishmaniasis (VL) patient were genetically engineered to express different fluorescent proteins and drug-resistance markers and subsequently used as parental strains in experimental sandfly co-infection. After 8 days of infection, sand flies were dissected and midgut promastigotes were transferred into double drug-selective media. Two double drug-resistant, dual fluorescent hybrid cell lines were recovered, which after cloning and whole genome sequencing, were shown to be full genomic hybrids. This study provides the first evidence of L. donovani hybridization within its natural vector Ph. argentipes.

DESCRIPTION, BIOLOGY, AND MEDICAL SIGNIFICANCE OF LEISHMANIA (MUNDINIA) CHANCEI N. SP. (KINETOPLASTEA: TRYPANOSOMATIDAE) FROM GHANA AND LEISHMANIA (MUNDINIA) PROCAVIENSIS N. SP. (KINETOPLASTEA: TRYPANOSOMATIDAE) FROM NAMIBIA

Genetic and phylogenetic analysis was performed on 2 isolates of Leishmania using DNA sequence data from the RNA polymerase II large subunit gene and the ribosomal protein L23a intergenic sequence. This showed the isolates to represent 2 new species within the subgenus Leishmania (Mundinia). The addition of Leishmania (Mundinia) chancei and Leishmania (Mundinia) procaviensis creates a total of 6 named species to date within this recently described subgenus of parasitic protozoa, containing both human pathogens and nonpathogens. Their widespread geographical distribution, basal phylogenetic position within the genus Leishmania, and probable non-sand fly vectors make these L. (Mundinia) species of significant medical and biological interest.

Comparative Genomic Analyses of New and Old World Viscerotropic Leishmanine Parasites: Further Insights into the Origins of Visceral Leishmaniasis Agents

Visceral leishmaniasis (VL), also known as kala-azar, is an anthropozoonotic disease affecting human populations on five continents. Aetiologic agents belong to the Leishmania (L.) donovani complex. Until the 1990s, three leishmanine parasites comprised this complex: L. (L.) donovani Laveran & Mesnil 1903, L. (L.) infantum Nicolle 1908, and L. (L.) chagasi Lainson & Shaw 1987 (=L. chagasi Cunha & Chagas 1937). The VL causal agent in the New World (NW) was previously identified as L. (L.) chagasi. After the development of molecular characterization, however, comparisons between L. (L.) chagasi and L. (L.) infantum showed high similarity, and L. (L.) chagasi was then regarded as synonymous with L. (L.) infantum. It was, therefore, suggested that L. (L.) chagasi was not native to the NW but had been introduced from the Old World by Iberian colonizers. However, in light of ecological evidence from the NW parasite’s enzootic cycle involving a wild phlebotomine vector (Lutzomyia longipalpis) and a wild mammal reservoir (the fox, Cerdocyon thous), we have recently analyzed by molecular clock comparisons of the DNA polymerase alpha subunit gene the whole-genome sequence of L. (L.) infantum chagasi of the most prevalent clinical form, atypical dermal leishmaniasis (ADL), from Honduras (Central America) with that of the same parasite from Brazil (South America), as well as those of L. (L.) donovani (India) and L. (L.) infantum (Europe), which revealed that the Honduran parasite is older ancestry (382,800 ya) than the parasite from Brazil (143,300 ya), L. (L.) donovani (33,776 ya), or L. (L.) infantum (13,000 ya). In the present work, we have now amplified the genomic comparisons among these leishmanine parasites, exploring mainly the variations in the genome for each chromosome, and the number of genomic SNPs for each chromosome. Although the results of this new analysis have confirmed a high genomic similarity (~99%) among these parasites [except L. (L.) donovani], the Honduran parasite revealed a single structural variation on chromosome 17, and the highest frequency of genomic SNPs (more than twice the number seen in the Brazilian one), which together to its extraordinary ancestry (382,800 ya) represent strong evidence that L. (L.) chagasi/L. (L.) infantum chagasi is, in fact, native to the NW, and therefore with valid taxonomic status. Furthermore, the Honduran parasite, the most ancestral viscerotropic leishmanine parasite, showed genomic and clinical taxonomic characteristics compatible with a new Leishmania species causing ADL in Central America.

Leishmania tarentolae: A new frontier in the epidemiology and control of the leishmaniases

Leishmaniasis (or the leishmaniases), classified as a neglected tropical parasitic disease, is found in parts of the tropics, subtropics and southern Europe. Leishmania parasites are transmitted by the bite of phlebotomine sand flies and million cases of human infection occur annually. Leishmania tarentolae has been historically considered a non-pathogenic protozoan of reptiles, which has been studied mainly for its potential biotechnological applications. However, some strains of L. tarentolae appear to be transiently infective to mammals. In areas where leishmaniasis is endemic, recent molecular diagnostics and serological positivity to L. tarentolae in humans and dogs have spurred interest in the interactions between these mammalian hosts, reptiles and Leishmania infantum, the main aetiologic agent of human and canine leishmaniasis. In this review, we discuss the systematics and biology of L. tarentolae in the insect vectors and the vertebrate hosts and address questions about evolution of reptilian leishmaniae. Furthermore, we discuss the possible usefulness of L. tarentolae for new vaccination strategies.

Leishmania guyanensis M4147 as a new LRV1-bearing model parasite: Phosphatidate phosphatase 2-like protein controls cell cycle progression and intracellular lipid content

Leishmaniasis is a parasitic vector-borne disease caused by the protistan flagellates of the genus Leishmania. Leishmania (Viannia) guyanensis is one of the most common causative agents of the American tegumentary leishmaniasis. It has previously been shown that L. guyanensis strains that carry the endosymbiotic Leishmania RNA virus 1 (LRV1) cause more severe form of the disease in a mouse model than those that do not. The presence of the virus was implicated into the parasite's replication and spreading. In this respect, studying the molecular mechanisms of cellular control of viral infection is of great medical importance. Here, we report ~30.5 Mb high-quality genome assembly of the LRV1-positive L. guyanensis M4147. This strain was turned into a model by establishing the CRISPR-Cas9 system and ablating the gene encoding phosphatidate phosphatase 2-like (PAP2L) protein. The orthologue of this gene is conspicuously absent from the genome of an unusual member of the family Trypanosomatidae, Vickermania ingenoplastis, a species with mostly bi-flagellated cells. Our analysis of the PAP2L-null L. guyanensis showed an increase in the number of cells strikingly resembling the bi-flagellated V. ingenoplastis, likely as a result of the disruption of the cell cycle, significant accumulation of phosphatidic acid, and increased virulence compared to the wild type cells.

Catalase impairs Leishmania mexicana development and virulence

Catalase is one of the most abundant enzymes on Earth. It decomposes hydrogen peroxide, thus protecting cells from dangerous reactive oxygen species. The catalase-encoding gene is conspicuously absent from the genome of most representatives of the family Trypanosomatidae. Here, we expressed this protein from the Leishmania mexicana Β-TUBULIN locus using a novel bicistronic expression system, which relies on the 2A peptide of Teschovirus A. We demonstrated that catalase-expressing parasites are severely compromised in their ability to develop in insects, to be transmitted and to infect mice, and to cause clinical manifestation in their mammalian host. Taken together, our data support the hypothesis that the presence of catalase is not compatible with the dixenous life cycle of Leishmania, resulting in loss of this gene from the genome during the evolution of these parasites.

Analyses of Leishmania-LRV Co-Phylogenetic Patterns and Evolutionary Variability of Viral Proteins

Leishmania spp. are important pathogens causing a vector-borne disease with a broad range of clinical manifestations from self-healing ulcers to the life-threatening visceral forms. Presence of Leishmania RNA virus (LRV) confers survival advantage to these parasites by suppressing anti-leishmanial immunity in the vertebrate host. The two viral species, LRV1 and LRV2 infect species of the subgenera Viannia and Leishmania, respectively. In this work we investigated co-phylogenetic patterns of leishmaniae and their viruses on a small scale (LRV2 in L. major) and demonstrated their predominant coevolution, occasionally broken by intraspecific host switches. Our analysis of the two viral genes, encoding the capsid and RNA-dependent RNA polymerase (RDRP), revealed them to be under the pressure of purifying selection, which was considerably stronger for the former gene across the whole tree. The selective pressure also differs between the LRV clades and correlates with the frequency of interspecific host switches. In addition, using experimental (capsid) and predicted (RDRP) models we demonstrated that the evolutionary variability across the structure is strikingly different in these two viral proteins.

Genomics of Trypanosomatidae: Where We Stand and What Needs to Be Done?

Trypanosomatids are easy to cultivate and they are (in many cases) amenable to genetic manipulation. Genome sequencing has become a standard tool routinely used in the study of these flagellates. In this review, we summarize the current state of the field and our vision of what needs to be done in order to achieve a more comprehensive picture of trypanosomatid evolution. This will also help to illuminate the lineage-specific proteins and pathways, which can be used as potential targets in treating diseases caused by these parasites.

The Maze Pathway of Coevolution: A Critical Review over the Leishmania and Its Endosymbiotic History

The description of the genus Leishmania as the causative agent of leishmaniasis occurred in the modern age. However, evolutionary studies suggest that the origin of Leishmania can be traced back to the Mesozoic era. Subsequently, during its evolutionary process, it achieved worldwide dispersion predating the breakup of the Gondwana supercontinent. It is assumed that this parasite evolved from monoxenic Trypanosomatidae. Phylogenetic studies locate dixenous Leishmania in a well-supported clade, in the recently named subfamily Leishmaniinae, which also includes monoxenous trypanosomatids. Virus-like particles have been reported in many species of this family. To date, several Leishmania species have been reported to be infected by Leishmania RNA virus (LRV) and Leishbunyavirus (LBV). Since the first descriptions of LRVs decades ago, differences in their genomic structures have been highlighted, leading to the designation of LRV1 in L. (Viannia) species and LRV2 in L. (Leishmania) species. There are strong indications that viruses that infect Leishmania spp. have the ability to enhance parasitic survival in humans as well as in experimental infections, through highly complex and specialized mechanisms. Phylogenetic analyses of these viruses have shown that their genomic differences correlate with the parasite species infected, suggesting a coevolutionary process. Herein, we will explore what has been described in the literature regarding the relationship between Leishmania and endosymbiotic Leishmania viruses and what is known about this association that could contribute to discussions about the worldwide dispersion of Leishmania.

Genome Analysis of Endotrypanum and Porcisia spp., Closest Phylogenetic Relatives of Leishmania, Highlights the Role of Amastins in Shaping Pathogenicity

While numerous genomes of Leishmania spp. have been sequenced and analyzed, an understanding of the evolutionary history of these organisms remains limited due to the unavailability of the sequence data for their closest known relatives, Endotrypanum and Porcisia spp., infecting sloths and porcupines. We have sequenced and analyzed genomes of three members of this clade in order to fill this gap. Their comparative analyses revealed only minute differences from Leishmaniamajor genome in terms of metabolic capacities. We also documented that the number of genes under positive selection on the Endotrypanum/Porcisia branch is rather small, with the flagellum-related group of genes being over-represented. Most significantly, the analysis of gene family evolution revealed a substantially reduced repertoire of surface proteins, such as amastins and biopterin transporters BT1 in the Endotrypanum/Porcisia species when compared to amastigote-dwelling Leishmania. This reduction was especially pronounced for δ-amastins, a subfamily of cell surface proteins crucial in the propagation of Leishmania amastigotes inside vertebrate macrophages and, apparently, dispensable for Endotrypanum/Porcisia, which do not infect such cells.

Molecular dating for phylogenies containing a mix of populations and species by using Bayesian and RelTime approaches

Simultaneous molecular dating of population and species divergences is essential in many biological investigations, including phylogeography, phylodynamics and species delimitation studies. In these investigations, multiple sequence alignments consist of both intra- and interspecies samples (mixed samples). As a result, the phylogenetic trees contain interspecies, interpopulation and within-population divergences. Bayesian relaxed clock methods are often employed in these analyses, but they assume the same tree prior for both inter- and intraspecies branching processes and require specification of a clock model for branch rates (independent vs. autocorrelated rates models). We evaluated the impact of a single tree prior on Bayesian divergence time estimates by analysing computer-simulated data sets. We also examined the effect of the assumption of independence of evolutionary rate variation among branches when the branch rates are autocorrelated. Bayesian approach with coalescent tree priors generally produced excellent molecular dates and highest posterior densities with high coverage probabilities. We also evaluated the performance of a non-Bayesian method, RelTime, which does not require the specification of a tree prior or a clock model. RelTime's performance was similar to that of the Bayesian approach, suggesting that it is also suitable to analyse data sets containing both populations and species variation when its computational efficiency is needed.

A review of the systematics, species identification and diagnostics of the Trypanosomatidae using the maxicircle kinetoplast DNA: from past to present

The Trypanosomatid family are a diverse and widespread group of protozoan parasites that belong to the higher order class Kinetoplastida. Containing predominantly monoxenous species (i.e. those having only a single host) that are confined to invertebrate hosts, this class is primarily known for its pathogenic dixenous species (i.e. those that have two hosts), serving as the aetiological agents of the important neglected tropical diseases including leishmaniasis, American trypanosomiasis (Chagas disease) and human African trypanosomiasis. Over the past few decades, a multitude of studies have investigated the diversity, classification and evolutionary history of the trypanosomatid family using different approaches and molecular targets. The mitochondrial-like DNA of the trypanosomatid parasites, also known as the kinetoplast, has emerged as a unique taxonomic and diagnostic target for exploring the evolution of this diverse group of parasitic eukaryotes. This review discusses recent advancements and important developments that have made a significant impact in the field of trypanosomatid systematics and diagnostics in recent years.

Proteins involved in the biosynthesis of lipophosphoglycan in Leishmania: a comparative genomic and evolutionary analysis

BACKGROUND

Leishmania spp. are digenetic parasites capable of infecting humans and causing a range of diseases collectively known as leishmaniasis. The main mechanisms involved in the development and permanence of this pathology are linked to evasion of the immune response. Crosstalk between the immune system and particularities of each pathogenic species is associated with diverse disease manifestations. Lipophosphoglycan (LPG), one of the most important molecules present on the surface of Leishmania parasites, is divided into four regions with high molecular variability. Although LPG plays an important role in host-pathogen and vector-parasite interactions, the distribution and phylogenetic relatedness of the genes responsible for its synthesis remain poorly explored. The recent availability of full genomes and transcriptomes of Leishmania parasites offers an opportunity to leverage insight on how LPG-related genes are distributed and expressed by these pathogens.

RESULTS

Using a phylogenomics-based framework, we identified a catalog of genes involved in LPG biosynthesis across 22 species of Leishmania from the subgenera Viannia and Leishmania, as well as 5 non-Leishmania trypanosomatids. The evolutionary relationships of these genes across species were also evaluated. Nine genes related to the production of the glycosylphosphatidylinositol (GPI)-anchor were highly conserved among compared species, whereas 22 genes related to the synthesis of the repeat unit presented variable conservation. Extensive gain/loss events were verified, particularly in genes SCG1-4 and SCA1-2. These genes act, respectively, on the synthesis of the side chain attached to phosphoglycans and in the transfer of arabinose residues. Phylogenetic analyses disclosed evolutionary patterns reflective of differences in host specialization, geographic origin and disease manifestation.

CONCLUSIONS

The multiple gene gain/loss events identified by genomic data mining help to explain some of the observed intra- and interspecies variation in LPG structure. Collectively, our results provide a comprehensive catalog that details how LPG-related genes evolved in the Leishmania parasite specialization process.

Kinetic and structural studies of Trypanosoma and Leishmania phosphofructokinases show evolutionary divergence and identify AMP as a switch regulating glycolysis versus gluconeogenesis

Trypanosomatids possess glycosome organelles that contain much of the glycolytic machinery, including phosphofructokinase (PFK). We present kinetic and structural data for PFK from three human pathogenic trypanosomatids, illustrating intriguing differences that may reflect evolutionary adaptations to differing ecological niches. The activity of Leishmania PFK – to a much larger extent than Trypanosoma PFK – is reliant on AMP for activity regulation, with 1 mm AMP increasing the L. infantum PFK (LiPFK) kcat/K_0.5^F6P value by 10‐fold, compared to only a 1.3‐ and 1.4‐fold increase for T. cruzi and T. brucei PFK, respectively. We also show that Leishmania PFK melts at a significantly lower (> 15 °C) temperature than Trypanosoma PFKs and that addition of either AMP or ATP results in a marked stabilization of the protein. Sequence comparisons of Trypanosoma spp. and Leishmania spp. show that divergence of the two genera involved amino acid substitutions that occur in the enzyme’s ‘reaching arms’ and ‘embracing arms’ that determine tetramer stability. The dramatic effects of AMP on Leishmania activity compared with the Trypanosoma PFKs may be explained by differences between the T‐to‐R equilibria for the two families, with the low‐melting Leishmania PFK favouring the flexible inactive T‐state in the absence of AMP. Sequence comparisons along with the enzymatic and structural data presented here also suggest there was a loss of AMP‐dependent regulation in Trypanosoma species rather than gain of this characteristic in Leishmania species and that AMP acts as a key regulator in Leishmania governing the balance between glycolysis and gluconeogenesis.

Comparative genomics of Leishmania (Mundinia)

Background

Trypanosomatids of the genus Leishmania are parasites of mammals or reptiles transmitted by bloodsucking dipterans. Many species of these flagellates cause important human diseases with clinical symptoms ranging from skin sores to life-threatening damage of visceral organs. The genus Leishmania contains four subgenera: Leishmania, Sauroleishmania, Viannia, and Mundinia. The last subgenus has been established recently and remains understudied, although Mundinia contains human-infecting species. In addition, it is interesting from the evolutionary viewpoint, representing the earliest branch within the genus and possibly with a different type of vector. Here we analyzed the genomes of L. (M.) martiniquensis, L. (M.) enriettii and L. (M.) macropodum to better understand the biology and evolution of these parasites.

Results

All three genomes analyzed were approximately of the same size (~ 30 Mb) and similar to that of L. (Sauroleishmania) tarentolae, but smaller than those of the members of subgenera Leishmania and Viannia, or the genus Endotrypanum (~ 32 Mb). This difference was explained by domination of gene losses over gains and contractions over expansions at the Mundinia node, although only a few of these genes could be identified. The analysis predicts significant changes in the Mundinia cell surface architecture, with the most important ones relating to losses of LPG-modifying side chain galactosyltransferases and arabinosyltransferases, as well as β-amastins. Among other important changes were gene family contractions for the oxygen-sensing adenylate cyclases and FYVE zinc finger-containing proteins.

Conclusions

We suggest that adaptation of Mundinia to different vectors and hosts has led to alternative host-parasite relationships and, thereby, made some proteins redundant. Thus, the evolution of genomes in the genus Leishmania and, in particular, in the subgenus Mundinia was mainly shaped by host (or vector) switches.

Evolutionary Insight into the Trypanosomatidae Using Alignment-Free Phylogenomics of the Kinetoplast

Advancements in next-generation sequencing techniques have led to a substantial increase in the genomic information available for analyses in evolutionary biology. As such, this data requires the exponential growth in bioinformatic methods and expertise required to understand such vast quantities of genomic data. Alignment-free phylogenomics offer an alternative approach for large-scale analyses that may have the potential to address these challenges. The evolutionary relationships between various species within the trypanosomatid family, specifically members belonging to the genera Leishmania and Trypanosoma have been extensively studies over the last 30 years. However, there is a need for a more exhaustive analysis of the Trypanosomatidae, summarising the evolutionary patterns amongst the entire family of these important protists. The mitochondrial DNA of the trypanosomatids, better known as the kinetoplast, represents a valuable taxonomic marker given its unique presence across all kinetoplastid protozoans. The aim of this study was to validate the reliability and robustness of alignment-free approaches for phylogenomic analyses and its applicability to reconstruct the evolutionary relationships between the trypanosomatid family. In the present study, alignment-free analyses demonstrated the strength of these methods, particularly when dealing with large datasets compared to the traditional phylogenetic approaches. We present a maxicircle genome phylogeny of 46 species spanning the trypanosomatid family, demonstrating the superiority of the maxicircle for the analysis and taxonomic resolution of the Trypanosomatidae.

Study of VIPER and TATE in kinetoplastids and the evolution of tyrosine recombinase retrotransposons

Background

Kinetoplastids are a flagellated group of protists, including some parasites, such as Trypanosoma and Leishmania species, that can cause diseases in humans and other animals. The genomes of these species enclose a fraction of retrotransposons including VIPER and TATE, two poorly studied transposable elements that encode a tyrosine recombinase (YR) and were previously classified as DIRS elements. This study investigated the distribution and evolution of VIPER and TATE in kinetoplastids to understand the relationships of these elements with other retrotransposons.

Results

We observed that VIPER and TATE have a discontinuous distribution among Trypanosomatidae, with several events of loss and degeneration occurring during a vertical transfer evolution. We were able to identify the terminal repeats of these elements for the first time, and we showed that these elements are potentially active in some species, including T. cruzi copies of VIPER. We found that VIPER and TATE are strictly related elements, which were named in this study as VIPER-like. The reverse transcriptase (RT) tree presented a low resolution, and the origin and relationships among YR groups remain uncertain. Conversely, for RH, VIPER-like grouped with Hepadnavirus, whereas for YR, VIPER-like sequences constituted two different clades that are closely allied to Crypton. Distinct topologies among RT, RH and YR trees suggest ancient rearrangements/exchanges in domains and a modular pattern of evolution with putative independent origins for each ORF.

Conclusions

Due to the presence of both elements in Bodo saltans, a nontrypanosomatid species, we suggested that VIPER and TATE have survived and remained active for more than 400 million years or were reactivated during the evolution of the host species. We did not find clear evidence of independent origins of VIPER-like from the other YR retroelements, supporting the maintenance of the DIRS group of retrotransposons. Nevertheless, according to phylogenetic findings and sequence structure obtained by this study and other works, we proposed separating DIRS elements into four subgroups: DIRS-like, PAT-like, Ngaro-like, and VIPER-like.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0175-2) contains supplementary material, which is available to authorized users.

Leishmania tarentolae: Taxonomic classification and its application as a promising biotechnological expression host

In this review, we summarize the current knowledge concerning the eukaryotic protozoan parasite Leishmania tarentolae, with a main focus on its potential for biotechnological applications. We will also discuss the genus, subgenus, and species-level classification of this parasite, its life cycle and geographical distribution, and similarities and differences to human-pathogenic species, as these aspects are relevant for the evaluation of biosafety aspects of L. tarentolae as host for recombinant DNA/protein applications. Studies indicate that strain LEM-125 but not strain TARII/UC of L. tarentolae might also be capable of infecting mammals, at least transiently. This could raise the question of whether the current biosafety level of this strain should be reevaluated. In addition, we will summarize the current state of biotechnological research involving L. tarentolae and explain why this eukaryotic parasite is an advantageous and promising human recombinant protein expression host. This summary includes overall biotechnological applications, insights into its protein expression machinery (especially on glycoprotein and antibody fragment expression), available expression vectors, cell culture conditions, and its potential as an immunotherapy agent for human leishmaniasis treatment. Furthermore, we will highlight useful online tools and, finally, discuss possible future applications such as the humanization of the glycosylation profile of L. tarentolae or the expression of mammalian recombinant proteins in amastigote-like cells of this species or in amastigotes of avirulent human-pathogenic Leishmania species.

Whole Genome Shotgun Phylogenomics Resolves the Pattern and Timing of Swallowtail Butterfly Evolution

Evolutionary relationships have remained unresolved in many well-studied groups, even though advances in next-generation sequencing and analysis, using approaches such as transcriptomics, anchored hybrid enrichment, or ultraconserved elements, have brought systematics to the brink of whole genome phylogenomics. Recently, it has become possible to sequence the entire genomes of numerous nonbiological models in parallel at reasonable cost, particularly with shotgun sequencing. Here, we identify orthologous coding sequences from whole-genome shotgun sequences, which we then use to investigate the relevance and power of phylogenomic relationship inference and time-calibrated tree estimation. We study an iconic group of butterflies—swallowtails of the family Papilionidae—that has remained phylogenetically unresolved, with continued debate about the timing of their diversification. Low-coverage whole genomes were obtained using Illumina shotgun sequencing for all genera. Genome assembly coupled to BLAST-based orthology searches allowed extraction of 6621 orthologous protein-coding genes for 45 Papilionidae species and 16 outgroup species (with 32% missing data after cleaning phases). Supermatrix phylogenomic analyses were performed with both maximum-likelihood (IQ-TREE) and Bayesian mixture models (PhyloBayes) for amino acid sequences, which produced a fully resolved phylogeny providing new insights into controversial relationships. Species tree reconstruction from gene trees was performed with ASTRAL and SuperTriplets and recovered the same phylogeny. We estimated gene site concordant factors to complement traditional node-support measures, which strengthens the robustness of inferred phylogenies. Bayesian estimates of divergence times based on a reduced data set (760 orthologs and 12% missing data) indicate a mid-Cretaceous origin of Papilionoidea around 99.2 Ma (95% credibility interval: 68.6–142.7 Ma) and Papilionidae around 71.4 Ma (49.8–103.6 Ma), with subsequent diversification of modern lineages well after the Cretaceous-Paleogene event. These results show that shotgun sequencing of whole genomes, even when highly fragmented, represents a powerful approach to phylogenomics and molecular dating in a group that has previously been refractory to resolution.

Fully automated sequence alignment methods are comparable to, and much faster than, traditional methods in large data sets: an example with hepatitis B virus

Aligning sequences for phylogenetic analysis (multiple sequence alignment; MSA) is an important, but increasingly computationally expensive step with the recent surge in DNA sequence data. Much of this sequence data is publicly available, but can be extremely fragmentary (i.e., a combination of full genomes and genomic fragments), which can compound the computational issues related to MSA. Traditionally, alignments are produced with automated algorithms and then checked and/or corrected "by eye" prior to phylogenetic inference. However, this manual curation is inefficient at the data scales required of modern phylogenetics and results in alignments that are not reproducible. Recently, methods have been developed for fully automating alignments of large data sets, but it is unclear if these methods produce alignments that result in compatible phylogenies when compared to more traditional alignment approaches that combined automated and manual methods. Here we use approximately 33,000 publicly available sequences from the hepatitis B virus (HBV), a globally distributed and rapidly evolving virus, to compare different alignment approaches. Using one data set comprised exclusively of whole genomes and a second that also included sequence fragments, we compared three MSA methods: (1) a purely automated approach using traditional software, (2) an automated approach including by eye manual editing, and (3) more recent fully automated approaches. To understand how these methods affect phylogenetic results, we compared resulting tree topologies based on these different alignment methods using multiple metrics. We further determined if the monophyly of existing HBV genotypes was supported in phylogenies estimated from each alignment type and under different statistical support thresholds. Traditional and fully automated alignments produced similar HBV phylogenies. Although there was variability between branch support thresholds, allowing lower support thresholds tended to result in more differences among trees. Therefore, differences between the trees could be best explained by phylogenetic uncertainty unrelated to the MSA method used. Nevertheless, automated alignment approaches did not require human intervention and were therefore considerably less time-intensive than traditional approaches. Because of this, we conclude that fully automated algorithms for MSA are fully compatible with older methods even in extremely difficult to align data sets. Additionally, we found that most HBV diagnostic genotypes did not correspond to evolutionarily-sound groups, regardless of alignment type and support threshold. This suggests there may be errors in genotype classification in the database or that HBV genotypes may need a revision.

Phylogenetic Studies

Phylogenetics is an important component of the systems biology approach. Knowledge about evolution of the genus Leishmania is essential to understand various aspects of basic biology of these parasites, such as parasite-host or parasite-vector relationships, biogeography, or epidemiology. Here, we present a comprehensive guideline for performing phylogenetic studies based on DNA sequence data, but with principles that can be adapted to protein sequences or other molecular markers. It is presented as a compilation of the most commonly used genetic targets for phylogenetic studies of Leishmania, including their respective primers for amplification and references, as well as details of PCR assays. Guidelines are, then, presented to choose the best targets in relation to the types of samples under study. Finally, and importantly, instructions are given to obtain optimal sequences, alignments, and datasets for the subsequent data analysis and phylogenetic inference. Different bioinformatics methods and software for phylogenetic inference are presented and explained. This chapter aims to provide a compilation of methods and generic guidelines to conduct phylogenetics of Leishmania for nonspecialists.

Leishmania (Mundinia) orientalis n. sp. (Trypanosomatidae), a parasite from Thailand responsible for localised cutaneous leishmaniasis

BACKGROUND

Leishmaniasis is an emerging disease in Thailand with an unknown incidence or prevalence. Although the number of properly characterized and clinically confirmed cases is about 20, it is suspected that this low number masks a potentially high prevalence, with clinical disease typically manifesting itself against an immunocompromised background, but with a substantial number of subclinical or cured cases of infection. To date leishmaniasis in Thailand has been mainly ascribed to two taxa within the recently erected subgenus Mundinia Shaw, Camargo & Teixeira, 2016, Leishmania (Mundinia) martiniquensis Desbois, Pratlong & Dedet, 2014 and a species that has not been formally described prior to this study.

RESULTS

A case of simple cutaneous leishmaniasis was diagnosed in a patient from Nan Province, Thailand. Molecular analysis of parasites derived from a biopsy sample revealed this to be a new species of Leishmania Ross, 1908, which has been named as Leishmania (Mundinia) orientalis Bates & Jariyapan n. sp. A formal description is provided, and this new taxon supercedes some isolates from the invalid taxon "Leishmania siamensis". A summary of all known cases of leishmaniasis with a corrected species identification is provided.

CONCLUSIONS

Three species of parasites are now known to cause leishmaniasis is Thailand, L. martiniquensis and L. orientalis n. sp. in the subgenus Mundinia, which contains the type-species Leishmania enriettii Muniz & Medina, 1948, and a single case of Leishmania infantum Nicolle, 1908. This study now enables epidemiological and other investigations into the biology of these unusual parasites to be conducted. It is recommended that the use of the taxonomically invalid name "L. siamensis" should be discontinued.

Leishmania naiffi and Leishmania guyanensis reference genomes highlight genome structure and gene evolution in the Viannia subgenus

The unicellular protozoan parasite Leishmania causes the neglected tropical disease leishmaniasis, affecting 12 million people in 98 countries. In South America, where the Viannia subgenus predominates, so far only L. (Viannia) braziliensis and L. (V.) panamensis have been sequenced, assembled and annotated as reference genomes. Addressing this deficit in molecular information can inform species typing, epidemiological monitoring and clinical treatment. Here, L. (V.) naiffi and L. (V.) guyanensis genomic DNA was sequenced to assemble these two genomes as draft references from short sequence reads. The methods used were tested using short sequence reads for L. braziliensis M2904 against its published reference as a comparison. This assembly and annotation pipeline identified 70 additional genes not annotated on the original M2904 reference. Phylogenetic and evolutionary comparisons of L. guyanensis and L. naiffi with 10 other Viannia genomes revealed four traits common to all Viannia: aneuploidy, 22 orthologous groups of genes absent in other Leishmania subgenera, elevated TATE transposon copies and a high NADH-dependent fumarate reductase gene copy number. Within the Viannia, there were limited structural changes in genome architecture specific to individual species: a 45 Kb amplification on chromosome 34 was present in all bar L. lainsoni, L. naiffi had a higher copy number of the virulence factor leishmanolysin, and laboratory isolate L. shawi M8408 had a possible minichromosome derived from the 3' end of chromosome 34. This combination of genome assembly, phylogenetics and comparative analysis across an extended panel of diverse Viannia has uncovered new insights into the origin and evolution of this subgenus and can help improve diagnostics for leishmaniasis surveillance.

Codon usage bias controls mRNA and protein abundance in trypanosomatids

Protein abundance differs from a few to millions of copies per cell. Trypanosoma brucei presents an excellent model for studies on codon bias and differential gene expression because transcription is broadly unregulated and uniform across the genome. T. brucei is also a major human and animal protozoal pathogen. Here, an experimental assessment, using synthetic reporter genes, revealed that GC3 codons have a major positive impact on both mRNA and protein abundance. Our estimates of relative expression, based on coding sequences alone (codon usage and sequence length), are within 2-fold of the observed values for the majority of measured cellular mRNAs (n > 7000) and proteins (n > 2000). Our estimates also correspond with expression measures from published transcriptome and proteome datasets from other trypanosomatids. We conclude that codon usage is a key factor affecting global relative mRNA and protein expression in trypanosomatids and that relative abundance can be effectively estimated using only protein coding sequences.

Evolutionary analyses of myosin genes in trypanosomatids show a history of expansion, secondary losses and neofunctionalization

Myosins are motor proteins that comprise a large and diversified family important for a broad range of functions. Two myosin classes, I and XIII, were previously assigned in Trypanosomatids, based mainly on the studies of Trypanosoma cruzi, T. brucei and Leishmania major, and important human pathogenic species; seven orphan myosins were identified in T. cruzi. Our results show that the great variety of T. cruzi myosins is also present in some closely related species and in Bodo saltans, a member of an early divergent branch of Kinetoplastida. Therefore, these myosins should no longer be considered "orphans". We proposed the classification of a kinetoplastid-specific myosin group into a new class, XXXVI. Moreover, our phylogenetic data suggest that a great repertoire of myosin genes was present in the last common ancestor of trypanosomatids and B. saltans, mainly resulting from several gene duplications. These genes have since been predominantly maintained in synteny in some species, and secondary losses explain the current distribution. We also found two interesting genes that were clearly derived from myosin genes, demonstrating that possible redundant or useless genes, instead of simply being lost, can serve as raw material for the evolution of new genes and functions.

High-throughput sequencing of kDNA amplicons for the analysis of Leishmania minicircles and identification of Neotropical species

Leishmania kinetoplast DNA contains thousands of small circular molecules referred to as kinetoplast DNA (kDNA) minicercles. kDNA minicircles are the preferred targets for sensitive Leishmania detection, because they are present in high copy number and contain conserved sequence blocks in which polymerase chain reaction (PCR) primers can be designed. On the other hand, the heterogenic nature of minicircle networks has hampered the use of this peculiar genomic region for strain typing. The characterization of Leishmania minicirculomes used to require isolation and cloning steps prior to sequencing. Here, we show that high-throughput sequencing of single minicircle PCR products allows bypassing these laborious laboratory tasks. The 120 bp long minicircle conserved region was amplified by PCR from 18 Leishmania strains representative of the major species complexes found in the Neotropics. High-throughput sequencing of PCR products enabled recovering significant numbers of distinct minicircle sequences from each strain, reflecting minicircle class diversity. Minicircle sequence analysis revealed patterns that are congruent with current hypothesis of Leishmania relationships. Then, we show that a barcoding-like approach based on minicircle sequence comparisons may allow reliable identifications of Leishmania spp. This work opens up promising perspectives for the study of kDNA minicercles and a variety of applications in Leishmania research.

A pilot study on fingerprinting Leishmania species from the Old World using Fourier transform infrared spectroscopy

Leishmania species are protozoan parasites and the causative agents of leishmaniasis, a vector borne disease that imposes a large health burden on individuals living mainly in tropical and subtropical regions. Different Leishmania species are responsible for the distinct clinical patterns, such as cutaneous, mucocutaneous, and visceral leishmaniasis, with the latter being potentially fatal if left untreated. For this reason, it is important to perform correct species identification and differentiation. Fourier transform infrared spectroscopy (FTIR) is an analytical spectroscopic technique increasingly being used as a potential tool for identification of microorganisms for diagnostic purposes. By employing mid-infrared (MIR) spectral data, it is not only possible to assess the chemical structures but also to achieve differentiation supported by multivariate statistic analysis. This work comprises a pilot study on differentiation of Leishmania species of the Old World (L. major, L. tropica, L. infantum, and L. donovani) as well as hybrids of distinct species by using vibrational spectroscopic fingerprints. Films of intact Leishmania parasites and their deoxyribonucleic acid (DNA) were characterized comparatively with respect to their biochemical nature and MIR spectral patterns. The strains’ hyperspectral datasets were multivariately examined by means of variance-based principal components analysis (PCA) and distance-based hierarchical cluster analysis (HCA). With the implementation of MIR spectral datasets we show that a phenotypic differentiation of Leishmania at species and intra-species level is feasible. Thus, FTIR spectroscopy can be further exploited for building up spectral databases of Leishmania parasites in view of high-throughput analysis of clinical specimens.

Graphical abstract

Revealing the mystery of metabolic adaptations using a genome scale model of Leishmania infantum

Human macrophage phagolysosome and sandfly midgut provide antagonistic ecological niches for Leishmania parasites to survive and proliferate. Parasites optimize their metabolism to utilize the available inadequate resources by adapting to those environments. Lately, a number of metabolomics studies have revived the interest to understand metabolic strategies utilized by the Leishmania parasite for optimal survival within its hosts. For the first time, we propose a reconstructed genome-scale metabolic model for Leishmania infantum JPCM5, the analyses of which not only captures observations reported by metabolomics studies in other Leishmania species but also divulges novel features of the L. infantum metabolome. Our results indicate that Leishmania metabolism is organized in such a way that the parasite can select appropriate alternatives to compensate for limited external substrates. A dynamic non-essential amino acid motif exists within the network that promotes a restricted redistribution of resources to yield required essential metabolites. Further, subcellular compartments regulate this metabolic re-routing by reinforcing the physiological coupling of specific reactions. This unique metabolic organization is robust against accidental errors and provides a wide array of choices for the parasite to achieve optimal survival.

The evolution of trypanosomatid taxonomy

Trypanosomatids are protozoan parasites of the class Kinetoplastida predominately restricted to invertebrate hosts (i.e. possess a monoxenous life-cycle). However, several genera are pathogenic to humans, animals and plants, and have an invertebrate vector that facilitates their transmission (i.e. possess a dixenous life-cycle). Phytomonas is one dixenous genus that includes several plant pathogens transmitted by phytophagous insects. Trypanosoma and Leishmania are dixenous genera that infect vertebrates, including humans, and are transmitted by hematophagous invertebrates. Traditionally, monoxenous trypanosomatids such as Leptomonas were distinguished from morphologically similar dixenous species based on their restriction to an invertebrate host. Nonetheless, this criterion is somewhat flawed as exemplified by Leptomonas seymouri which reportedly infects vertebrates opportunistically. Similarly, Novymonas and Zelonia are presumably monoxenous genera yet sit comfortably in the dixenous clade occupied by Leishmania. The isolation of Leishmania macropodum from a biting midge (Forcipomyia spp.) rather than a phlebotomine sand fly calls into question the exclusivity of the Leishmania-sand fly relationship, and its suitability for defining the Leishmania genus. It is now accepted that classic genus-defining characteristics based on parasite morphology and host range are insufficient to form the sole basis of trypanosomatid taxonomy as this has led to several instances of paraphyly. While improvements have been made, resolution of evolutionary relationships within the Trypanosomatidae is confounded by our incomplete knowledge of its true diversity. The known trypanosomatids probably represent a fraction of those that exist and isolation of new species will help resolve relationships in this group with greater accuracy. This review incites a dialogue on how our understanding of the relationships between certain trypanosomatids has shifted, and discusses new knowledge that informs the present taxonomy of these important parasites.

What pre-Columbian mummies could teach us about South American leishmaniases?

A recent report on the taxonomic profile of the human gut microbiome in pre-Columbian mummies (Santiago-Rodriguez et al. 2016) gives for the first time evidence of the presence of Leishmania DNA (sequences similar to Leishmania donovani according to the authors) that can be reminiscent of visceral leishmaniasis during the pre-Columbian era. It is commonly assumed that Leishmania infantum, the etiological agent of American visceral leishmaniasis (AVL) was introduced into the New World by the Iberian conquest. This finding is really surprising and must be put into perspective with what is known from an AVL epidemiological and historical point of view. Beside L. infantum, there are other species that are occasionally reported to cause AVL in the New World. Among these, L. colombiensis is present in the region of pre-Columbian mummies studied. Other explanations for these findings include a more ancient introduction of a visceral species of Leishmania from the Old World or the existence of a yet unidentified endemic species causing visceral leishmaniasis in South America. Unfortunately, very few molecular data are known about this very long pre-Columbian period concerning the circulating species of Leishmania and their diversity in America.

The genome of Leishmania adleri from a mammalian host highlights chromosome fission in Sauroleishmania

Control of pathogens arising from humans, livestock and wild animals can be enhanced by genome-based investigation. Phylogenetically classifying and optimal construction of these genomes using short sequence reads are key to this process. We examined the mammal-infecting unicellular parasite Leishmania adleri belonging to the lizard-infecting Sauroleishmania subgenus. L. adleri has been associated with cutaneous disease in humans, but can be asymptomatic in wild animals. We sequenced, assembled and investigated the L. adleri genome isolated from an asymptomatic Ethiopian rodent (MARV/ET/75/HO174) and verified it as L. adleri by comparison with other Sauroleishmania species. Chromosome-level scaffolding was achieved by combining reference-guided with de novo assembly followed by extensive improvement steps to produce a final draft genome with contiguity comparable with other references. L. tarentolae and L. major genome annotation was transferred and these gene models were manually verified and improved. This first high-quality draft Leishmania adleri reference genome is also the first Sauroleishmania genome from a non-reptilian host. Comparison of the L. adleri HO174 genome with those of L. tarentolae Parrot-TarII and lizard-infecting L. adleri RLAT/KE/1957/SKINK-7 showed extensive gene amplifications, pervasive aneuploidy, and fission of chromosomes 30 and 36. There was little genetic differentiation between L. adleri extracted from mammals and reptiles, highlighting challenges for leishmaniasis surveillance.

The history of leishmaniasis

In this review article the history of leishmaniasis is discussed regarding the origin of the genus Leishmania in the Mesozoic era and its subsequent geographical distribution, initial evidence of the disease in ancient times, first accounts of the infection in the Middle Ages, and the discovery of Leishmania parasites as causative agents of leishmaniasis in modern times. With respect to the origin and dispersal of Leishmania parasites, the three currently debated hypotheses (Palaearctic, Neotropical and supercontinental origin, respectively) are presented. Ancient documents and paleoparasitological data indicate that leishmaniasis was already widespread in antiquity. Identification of Leishmania parasites as etiological agents and sand flies as the transmission vectors of leishmaniasis started at the beginning of the 20^th century and the discovery of new Leishmania and sand fly species continued well into the 21^st century. Lately, the Syrian civil war and refugee crises have shown that leishmaniasis epidemics can happen any time in conflict areas and neighbouring regions where the disease was previously endemic.

Leishmania infections: Molecular targets and diagnosis

Progress in the diagnosis of leishmaniases depends on the development of effective methods and the discovery of suitable biomarkers. We propose firstly an update classification of Leishmania species and their synonymies. We demonstrate a global map highlighting the geography of known endemic Leishmania species pathogenic to humans. We summarize a complete list of techniques currently in use and discuss their advantages and limitations. The available data highlights the benefits of molecular markers in terms of their sensitivity and specificity to quantify variation from the subgeneric level to species complexes, (sub) species within complexes, and individual populations and infection foci. Each DNA-based detection method is supplied with a comprehensive description of markers and primers and proposal for a classification based on the role of each target and primer in the detection, identification and quantification of leishmaniasis infection. We outline a genome-wide map of genes informative for diagnosis that have been used for Leishmania genotyping. Furthermore, we propose a classification method based on the suitability of well-studied molecular markers for typing the 21 known Leishmania species pathogenic to humans. This can be applied to newly discovered species and to hybrid strains originating from inter-species crosses. Developing more effective and sensitive diagnostic methods and biomarkers is vital for enhancing Leishmania infection control programs.

Comparative genomics of canine-isolated Leishmania (Leishmania) amazonensis from an endemic focus of visceral leishmaniasis in Governador Valadares, southeastern Brazil

Leishmaniasis is a highly diverse group of diseases caused by kinetoplastid of the genus Leishmania. These parasites are taxonomically diverse, with human pathogenic species separated into two subgenera according to their development site inside the alimentary tract of the sand fly insect vector. The disease encompasses a variable spectrum of clinical manifestations with tegumentary or visceral symptoms. Among the causative species in Brazil, Leishmania (Leishmania) amazonensis is an important etiological agent of human cutaneous leishmaniasis that accounts for more than 8% of all cases in endemic regions. L. (L.) amazonensis is generally found in the north and northeast regions of Brazil. Here, we report the first isolation of L. (L.) amazonensis from dogs with clinical manifestations of visceral leishmaniasis in Governador Valadares, an endemic focus in the southeastern Brazilian State of Minas Gerais where L. (L.) infantum is also endemic. These isolates were characterized in terms of SNPs, chromosome and gene copy number variations, confirming that they are closely related to a previously sequenced isolate obtained in 1973 from the typical Northern range of this species. The results presented in this article will increase our knowledge of L. (L.) amazonensis-specific adaptations to infection, parasite survival and the transmission of this Amazonian species in a new endemic area of Brazil.

Isolation of Novel Trypanosomatid, Zelonia australiensis sp. nov. (Kinetoplastida: Trypanosomatidae) Provides Support for a Gondwanan Origin of Dixenous Parasitism in the Leishmaniinae

The genus Leishmania includes approximately 53 species, 20 of which cause human leishmaniais; a significant albeit neglected tropical disease. Leishmaniasis has afflicted humans for millennia, but how ancient is Leishmania and where did it arise? These questions have been hotly debated for decades and several theories have been proposed. One theory suggests Leishmania originated in the Palearctic, and dispersed to the New World via the Bering land bridge. Others propose that Leishmania evolved in the Neotropics. The Multiple Origins theory suggests that separation of certain Old World and New World species occurred due to the opening of the Atlantic Ocean. Some suggest that the ancestor of the dixenous genera Leishmania, Endotrypanum and Porcisia evolved on Gondwana between 90 and 140 million years ago. In the present study a detailed molecular and morphological characterisation was performed on a novel Australian trypanosomatid following its isolation in Australia’s tropics from the native black fly, Simulium (Morops) dycei Colbo, 1976. Phylogenetic analyses were conducted and confirmed this parasite as a sibling to Zelonia costaricensis, a close relative of Leishmania previously isolated from a reduviid bug in Costa Rica. Consequently, this parasite was assigned the name Zelonia australiensis sp. nov. Assuming Z. costaricensis and Z. australiensis diverged when Australia and South America became completely separated, their divergence occurred between 36 and 41 million years ago at least. Using this vicariance event as a calibration point for a phylogenetic time tree, the common ancestor of the dixenous genera Leishmania, Endotrypanum and Porcisia appeared in Gondwana approximately 91 million years ago. Ultimately, this study contributes to our understanding of trypanosomatid diversity, and of Leishmania origins by providing support for a Gondwanan origin of dixenous parasitism in the Leishmaniinae.

The genus Leishmania includes approximately 53 species, 20 of which cause human leishmaniais, a significant disease that has afflicted humans for millennia. But how ancient is Leishmania and where did it arise? Some suggest Leishmania originated in the Palearctic. Others suggest it appeared in the Neotropics. The Multiple Origins theory proposes that separation of certain Old World and Neotropical species occurred following the opening of the Atlantic. Others suggest that an ancestor to the Euleishmania and Paraleishmania appeared on Gondwana 90 to 140 million years ago (MYA). We performed a detailed molecular and morphological characterisation of a novel Australian trypanosomatid. This parasite is a sibling to the Neotropical Zelonia costaricensis, a close relative of Leishmania, and designated as Zelonia australiensis sp. nov. Assuming Z. costaricensis and Z. australiensis split when Australia and South America separated, their divergence occurred between 36 and 41 MYA. Using this event as a calibration point for a phylogenetic time tree, an ancestor of the dixenous Leishmaniinae appeared in Gondwana ~ 91 MYA. This study contributes to our understanding of trypanosomatid diversity by describing a unique Australian trypanosomatid and to our understanding of Leishmania evolution by inferring a Gondwanan origin for dixenous parasitism in the Leishmaniinae.

An appraisal of the taxonomy and nomenclature of trypanosomatids presently classified as Leishmania and Endotrypanum

We propose a taxonomic revision of the dixenous trypanosomatids currently classified as Endotrypanum and Leishmania, including parasites that do not fall within the subgenera L. (Leishmania) and L. (Viannia) related to human leishmaniasis or L. (Sauroleishmania) formed by leishmanias of lizards: L. colombiensis, L. equatorensis, L. herreri, L. hertigi, L. deanei, L. enriettii and L. martiniquensis. The comparison of these species with newly characterized isolates from sloths, porcupines and phlebotomines from central and South America unveiled new genera and subgenera supported by past (RNA PolII gene) and present (V7V8 SSU rRNA, Hsp70 and gGAPDH) phylogenetic analyses of the organisms. The genus Endotrypanum is restricted to Central and South America, comprising isolates from sloths and transmitted by phlebotomines that sporadically infect humans. This genus is the closest to the new genus Porcisia proposed to accommodate the Neotropical porcupine parasites originally described as L. hertigi and L. deanei. A new subgenus Leishmania (Mundinia) is created for the L. enriettii complex that includes L. martiniquensis. The new genus Zelonia harbours trypanosomatids from Neotropical hemipterans placed at the edge of the Leishmania-Endotrypanum-Porcisia clade. Finally, attention is drawn to the status of L. siamensis and L. australiensis as nomem nudums.

Genomic and phylogenetic evidence of VIPER retrotransposon domestication in trypanosomatids

Transposable elements are important residents of eukaryotic genomes and eventually the host can domesticate them to serve cellular functions. We reported here a possible domestication event of the vestigial interposed retroelement (VIPER) in trypanosomatids. We found a large gene in a syntenic location in Leishmania braziliensis, L. panamensis, Leptomanas pyrrhocoris, and Crithidia fasciculata whose products share similarity in the C-terminal portion with the third protein of VIPER. No remnants of other VIPER regions surrounding the gene sequence were found. We hypothesise that the domestication event occurred more than 50 mya and the conservation of this gene suggests it might perform some function in the host species.

William Boog Leishman: parasitologist and politician

Famous for the discovery of the parasite, Leishmania, named after him, and the invention of Leishman's stain, William Boog Leishman should perhaps be better known for his work in military and public health, particularly the prevention of typhoid. Leishman was a Medical Officer in the British Army from 1887 until his death in 1926. His early research was on diseases affecting troops posted to stations within the British Empire. He saw cases of Leishmaniasis while stationed in India, and was able to identify the causative organism from his detailed records of his observations. Leishman's most important contribution to public health, however, was his work with typhoid, a major cause of morbidity and mortality in the army. Leishman planned experiments and the collection of data to demonstrate the efficacy of anti-typhoid inoculation and, using his considerable political skills, advocated the adoption of the vaccine. He planned for the inoculation of troops in an emergency so, when war broke out in 1914, the vaccine was available to save thousands of lives. Leishman's colleagues and mentors included Ronald Ross and Almroth Wright. Leishman was less outspoken than either Ross or Wright; this paper shows how the different contributions of the three men overlapped.

Molecular detection, identification and phylogenetic inference of Leishmania spp. in some desert lizards from Northwest China by using internal transcribed spacer 1 (ITS1) sequences

Leishmaniasis caused by Leishmania is still endemic in Northwest China. It has been thought that reptiles could be a reservoir for mammalian leishmaniasis. However, data are still scarce on natural infection of lizards with Leishmania spp. in China. The present study deals with detection, identification and phylogenetic inference of Leishmania parasites at species and intraspecies levels isolated from six desert lizard species from 10 geographical locations in Northwest China using amplification and sequencing of ITS-rDNA. In total, 83 haplotypes were found among 137 ITS1 sequences obtained from up to 64.6% of all captured lizards. Representative sequences of Leishmania available in GenBank were compiled for comparison with the obtained haplotypes. Tree-based species delimitation was achieved by using Bayesian phylogenitc analyses and maximum parsimony approach. Phylogenetic trees congruently supported that the haplotypes were found to belong to three Leishmania species including L. (sauroleishmania) sp., Leishmania tropica and Leishmania donovani complex. A network approach revealed paraphyletic populations of L. (sauroleishmania) sp. and L. tropica at intraspecies level regarding geographical origin and low host specificity. Chinese L. tropica from lizards showed significant heterogeneity as the obtained haplotypes were distributed in different clusters from other countries. Common ancestry was observed between some sequences of L. tropica from lizards and other sequence types from clinical samples from other countries. This may lend support to the potential reservoir role of lizards for human leishmaniasis. Our results appear to be the first molecular evidence for natural infection of lizards in Northwest China with reptilian Leishmania and mammalian Leishmania species. Desert lizards may be considered as putative reservoir hosts for Leishmania in China. Further studies on persistence of the Leishmania parasites in lizards and sandflies are recommended for the better understanding of their epidemiological involvement.