Comparative genomic analysis of three Leishmania species that cause diverse human disease

Genomic and epidemiological evidence for the emergence of a L. infantum/L. donovani hybrid with unusual epidemiology in northern Italy

Leishmania (L.) infantum is one of the main causative agents of animal and human leishmaniasis across many endemic areas in South America, Europe, North Africa, and Asia. Despite its clinical significance, little is known about the genetic diversity of L. infantum circulating in a given endemic area. Here, we investigate this important open question by applying a comparative genomics approach to seven L. infantum isolates from different hosts and Italian regions, including the northern part of the country (Emilia-Romagna, RER), Sicily, and Sardinia, as an initial attempt to explore the breadth of parasite genetic heterogeneity in Italy. Additionally, microsatellite analysis was carried out to compare the isolates from RER with other 70 L. infantum strains from the same region as well as 65 strains belonging to the L. donovani complex from other countries. We revealed important karyotypic instability and identified strain-specific changes in gene dosage, which affected important virulence factors such as amastins and surface antigen-like proteins. Single nucleotide polymorphism-based clustering analysis of these genomes together with over 80 publicly available L. infantum and L. donovani genomes placed the Italian isolates into three geographically distinct clusters within the Mediterranean basin and uncovered three isolates clustering with putative L. infantum/L. donovani hybrids isolated in Cyprus. As judged by microsatellite profiling, these hybrid isolates are representative of a sub-population of parasites circulating in northern Italy that preferentially infect humans but not dogs. Our results place Italy at the crossroads of L. infantum infection in the Mediterranean and call attention to the public health risk represented by the introduction of non-European Leishmania species.IMPORTANCEThis study closes important knowledge gaps with respect to Leishmania (L.) infantum genetic heterogeneity in a given endemic country, as exemplified here for Italy, and reveals genetic hybridization as a main cause for re-emerging human leishmaniasis in northern Italy. The observed high diversity of Leishmania parasites on the Italian peninsula suggests different geographical origins, with genomic adaptation to various ecologies affecting both pathogenicity and transmission potential. This is documented by the discovery of a putative L. infantum/L. donovani hybrid strain, which has been shown to preferentially infect humans but not dogs. Our results provide important information to health authorities, which need to consider the public health risk represented by the introduction of new Leishmania species into EU countries due to population displacement or travel from countries where exotic/allochthonous parasite species are endemic.

Exploring the functionality and conservation of Alba proteins in Trypanosoma cruzi: A focus on biological diversity and RNA binding ability

Trypanosoma cruzi is the causative agent of Chagas disease, as well as a trypanosomatid parasite with a complex biological cycle that requires precise mechanisms for regulating gene expression. In Trypanosomatidae, gene regulation occurs mainly at the mRNA level through the recognition of cis elements by RNA-binding proteins (RBPs). Alba family members are ubiquitous DNA/RNA-binding proteins with representatives in trypanosomatid parasites functionally related to gene expression regulation. Although T. cruzi possesses two groups of Alba proteins (Alba1/2 and Alba30/40), their functional role remains poorly understood. Thus, herein, a characterization of T. cruzi Alba (TcAlba) proteins was undertaken. Physicochemical, structural, and phylogenetic analysis of TcAlba showed features compatible with RBPs, such as hydrophilicity, RBP domains/motifs, and evolutionary conservation of the Alba-domain, mainly regarding other trypanosomatid Alba. However, in silico RNA interaction analysis of T. cruzi Alba proteins showed that TcAlba30/40 proteins, but not TcAlba1/2, would directly interact with the assayed RNA molecules, suggesting that these two groups of TcAlba proteins have different targets. Given the marked differences existing between both T. cruzi Alba groups (TcAlba1/2 and TcAlba30/40), regarding sequence divergence, RNA binding potential, and life-cycle expression patterns, we suggest that they would be involved in different biological processes.

Epidemiological surveillance of Leishmania infantum in wild lagomorphs in Spanish Mediterranean ecosystems

Wild lagomorphs play a key epidemiological role as reservoirs of Leishmania infantum, causative agent of the largest outbreak of human leishmaniosis in Europe to date. A large-scale survey study was conducted on wild rabbit (Oryctolagus cuniculus) and Iberian hare (Lepus granatensis) populations in Spanish Mediterranean ecosystems to evaluate the exposure of L. infantum and investigate potential risk factors associated with exposure to this zoonotic parasite. Between 2018 and 2021, a total of 631 wild lagomorphs (471 wild rabbits and 160 Iberian hares) were collected in Andalusia (southern Spain) and tested for antibodies against L. infantum using the indirect fluorescent antibody test (IFAT). Spleen samples from 563 of the wild lagomorphs sampled (441 wild rabbits and 122 Iberian hares) were also evaluated by real-time quantitative PCR (qPCR) for detection of Leishmania kDNA. Exposure to L. infantum (positive by IFAT and/or qPCR) was detected in 56.4 % (356/631; 95 %CI: 52.3-60.3) of the lagomorphs analyzed. Anti-Leishmania antibodies were found in 12.8 % (81/631; 95 %CI: 10.2-15.5) of the animals, and L. infantum kDNA was detected in 59.0 % (332/563; 95 %CI: 54.9-63.0) of the spleen samples tested. Phylogenetic analysis revealed high homology (99.9-100 %) between L. infantum sequences obtained and strains previously isolated from humans in Spain. While apparent seroprevalence was significantly higher in Iberian hares (19.4 %; 95 %CI: 13.3-25.5) compared to wild rabbits (10.6 %; 95 %CI: 7.9-13.4), no significant differences in prevalence were found between wild rabbits (61.0 %; 95 %CI: 56.5-65.6) and Iberian hares (51.6 %; 95 %CI: 42.8-60.5). At least one positive animal was found on 64.8 % (70/108) of the hunting grounds sampled, and a high-risk spatial cluster (P < 0.001) was also identified in central Andalusia. The multivariable analysis identified bioclimatic level (meso-Mediterranean climate) and the presence of goats on hunting grounds as risk factors potentially associated with L. infantum exposure in wild lagomorphs. This study shows high, widespread exposure, but heterogeneous distribution of L. infantum in wild lagomorph populations in Mediterranean ecosystems in southern Spain. The results point to the need to promote integrated surveillance programs for the detection of Leishmania spp. in wild lagomorphs in order to establish effective control measures against human leishmaniosis under a One Health approach.

A panel of recombinant Leishmania donovani cell surface and secreted proteins identifies LdBPK_323600.1 as a serological marker of symptomatic infection

Visceral leishmaniasis is a deadly infectious disease and is one of the world's major neglected health problems. Because the symptoms of infection are similar to other endemic diseases, accurate diagnosis is crucial for appropriate treatment. Definitive diagnosis using splenic or bone marrow aspirates is highly invasive, and so, serological assays are preferred, including the direct agglutination test (DAT) or rK39 strip test. These tests, however, are either difficult to perform in the field (DAT) or lack specificity in some endemic regions (rK39), making the development of new tests a research priority. The availability of Leishmania spp. genomes presents an opportunity to identify new diagnostic targets. Here, we use genome data and a mammalian protein expression system to create a panel of 93 proteins consisting of the extracellular ectodomains of the Leishmania donovani cell surface and secreted proteins. We use these panel and sera from murine experimental infection models and natural human and canine infections to identify new candidates for serological diagnosis. We observed a concordance between the most immunoreactive antigens in different host species and transmission settings. The antigen encoded by the LdBPK_323600.1 gene can diagnose Leishmania infections with high sensitivity and specificity in patient cohorts from different endemic regions including Bangladesh and Ethiopia. In longitudinal sampling of treated patients, we observed reductions in immunoreactivity to LdBPK_323600.1 suggesting it could be used to diagnose treatment success. In summary, we have identified new antigens that could contribute to improved serological diagnostic tests to help control the impact of this deadly tropical infectious disease.

IMPORTANCE

Visceral leishmaniasis is fatal if left untreated with patients often displaying mild and non-specific symptoms during the early stages of infection making accurate diagnosis important. Current methods for diagnosis require highly trained medical staff to perform highly invasive biopsies of the liver or bone marrow which pose risks to the patient. Less invasive molecular tests are available but can suffer from regional variations in their ability to accurately diagnose an infection. To identify new diagnostic markers of visceral leishmaniasis, we produced and tested a panel of 93 proteins identified from the genome of the parasite responsible for this disease. We found that the pattern of host antibody reactivity to these proteins was broadly consistent across naturally acquired infections in both human patients and dogs, as well as experimental rodent infections. We identified a new protein called LdBPK_323600.1 that could accurately diagnose visceral leishmaniasis infections in humans.

Molecular tools to regulate gene expression in Trypanosoma cruzi

Developing molecular strategies to manipulate gene expression in trypanosomatids is challenging, particularly with respect to the unique gene expression mechanisms adopted by these unicellular parasites, such as polycistronic mRNA transcription and multi-gene families. In the case of Trypanosoma cruzi (T. cruzi), the causative agent of Chagas Disease, the lack of RNA interference machinery further complicated functional genetic studies important for understanding parasitic biology and developing biomarkers and potential therapeutic targets. Therefore, alternative methods of performing knockout and/or endogenous labelling experiments were developed to identify and understand the function of proteins for survival and interaction with the host. In this review, we present the main tools for the genetic manipulation of T. cruzi, focusing on the Clustered Regularly Interspaced Short Palindromic Repeats Cas9-associated system technique widely used in this organism. Moreover, we highlight the importance of using these tools to elucidate the function of uncharacterized and glycosylated proteins. Further developments of these technologies will allow the identification of new biomarkers, therapeutic targets and potential vaccines against Chagas disease with greater efficiency and speed.

Ancestral aneuploidy and stable chromosomal duplication resulting in differential genome structure and gene expression control in trypanosomatid parasites

Aneuploidy is widely observed in both unicellular and multicellular eukaryotes, usually associated with adaptation to stress conditions. Chromosomal duplication stability is a tradeoff between the fitness cost of having unbalanced gene copies and the potential fitness gained from increased dosage of specific advantageous genes. Trypanosomatids, a family of protozoans that include species that cause neglected tropical diseases, are a relevant group to study aneuploidies. Their life cycle has several stressors that could select for different patterns of chromosomal duplications and/or losses, and their nearly universal use of polycistronic transcription increases their reliance on gene expansion/contraction, as well as post-transcriptional control as mechanisms for gene expression regulation. By evaluating the data from 866 isolates covering seven trypanosomatid genera, we have revealed that aneuploidy tolerance is an ancestral characteristic of trypanosomatids but has a reduced occurrence in a specific monophyletic clade that has undergone large genomic reorganization and chromosomal fusions. We have also identified an ancient chromosomal duplication that was maintained across these parasite's speciation, named collectively as the trypanosomatid ancestral supernumerary chromosome (TASC). TASC has most genes in the same coding strand, is expressed as a disomic chromosome (even having four copies), and has increased potential for functional variation, but it purges highly deleterious mutations more efficiently than other chromosomes. The evidence of stringent control over gene expression in this chromosome suggests that these parasites have adapted to mitigate the fitness cost associated with this ancient chromosomal duplication.

Effect of the aggregated protein dye YAT2150 on Leishmania parasite viability

The problems associated with the drugs currently used to treat leishmaniasis, including resistance, toxicity, and the high cost of some formulations, call for the urgent identification of new therapeutic agents with novel modes of action. The aggregated protein dye YAT2150 has been found to be a potent antileishmanial compound, with a half-maximal inhibitory concentration (IC50) of approximately 0.5 µM against promastigote and amastigote stages of Leishmania infantum. The encapsulation in liposomes of YAT2150 significantly improved its in vitro IC50 to 0.37 and 0.19 µM in promastigotes and amastigotes, respectively, and increased the half-maximal cytotoxic concentration in human umbilical vein endothelial cells to >50 µM. YAT2150 became strongly fluorescent when binding intracellular protein deposits in Leishmania cells. This fluorescence pattern aligns with the proposed mode of action of this drug in the malaria parasite Plasmodium falciparum, the inhibition of protein aggregation. In Leishmania major, YAT2150 rapidly reduced ATP levels, suggesting an alternative antileishmanial mechanism. To the best of our knowledge, this first-in-class compound is the only one described so far having significant activity against both Plasmodium and Leishmania, thus being a potential drug for the treatment of co-infections of both parasites.

Molecular Mechanisms of Persistence in Protozoan Parasites

Protozoan parasites are known for their remarkable capacity to persist within the bodies of vertebrate hosts, which frequently results in prolonged infections and the recurrence of diseases. Understanding the molecular mechanisms that underlie the event of persistence is of paramount significance to develop innovative therapeutic approaches, given that these pathways still need to be thoroughly elucidated. The present article provides a comprehensive overview of the latest developments in the investigation of protozoan persistence in vertebrate hosts. The focus is primarily on the function of persisters, their formation within the host, and the specific molecular interactions between host and parasite while they persist. Additionally, we examine the metabolomic, transcriptional, and translational changes that protozoan parasites undergo during persistence within vertebrate hosts, focusing on major parasites such as Plasmodium spp., Trypanosoma spp., Leishmania spp., and Toxoplasma spp. Key findings of our study suggest that protozoan parasites deploy several molecular and physiological strategies to evade the host immune surveillance and sustain their persistence. Furthermore, some parasites undergo stage differentiation, enabling them to acclimate to varying host environments and immune challenges. More often, stressors such as drug exposure were demonstrated to impact the formation of protozoan persisters significantly. Understanding the molecular mechanisms regulating the persistence of protozoan parasites in vertebrate hosts can reinvigorate our current insights into host-parasite interactions and facilitate the development of more efficacious disease therapeutics.

Complete assembly, annotation of virulence genes and CRISPR editing of the genome of Leishmania amazonensis PH8 strain

We report the sequencing and assembly of the PH8 strain of Leishmania amazonensis one of the etiological agents of leishmaniasis. After combining data from long Pacbio reads, short Illumina reads and synteny with the Leishmania mexicana genome, the sequence of 34 chromosomes with 8317 annotated genes was generated. Multigene families encoding three virulence factors, A2, amastins and the GP63 metalloproteases, were identified and compared to their annotation in other Leishmania species. As they have been recently recognized as virulence factors essential for disease establishment and progression of the infection, we also identified 14 genes encoding proteins involved in parasite iron and heme metabolism and compared to genes from other Trypanosomatids. To follow these studies with a genetic approach to address the role of virulence factors, we tested two CRISPR-Cas9 protocols to generate L. amazonensis knockout cell lines, using the Miltefosine transporter gene as a proof of concept.

Parasite Detection in Visceral Leishmaniasis Samples by Dye-Based qPCR Using New Gene Targets of Leishmania infantum and Crithidia

Visceral leishmaniasis (VL) is a neglected disease considered a serious public health problem, especially in endemic countries. Several studies have discovered monoxenous trypanosomatids (Leptomonas and Crithidia) in patients with VL. In different situations of leishmaniasis, investigations have examined cases of co-infection between Leishmania spp. and Crithidia spp. These coinfections have been observed in a wide range of vertebrate hosts, indicating that they are not rare. Diagnostic techniques require improvements and more robust tools to accurately detect the causative agent of VL. This study aimed to develop a real-time quantitative dye-based PCR (qPCR) assay capable of distinguishing Leishmania infantum from Crithidia-related species and to estimate the parasite load in samples of VL from humans and animals. The primer LinJ31_2420 targets an exclusive phosphatase of L. infantum; the primer Catalase_LVH60-12060_1F targets the catalase gene of Crithidia. Therefore, primers were designed to detect L. infantum and Crithidia sp. LVH60A (a novel trypanosomatid isolated from VL patients in Brazil), in samples related to VL. These primers were considered species-specific, based on sequence analysis using genome data retrieved from the TriTryp database and the genome assembling of Crithidia sp. LVH60A strain, in addition to experimental and clinical data presented herein. This novel qPCR assay was highly accurate in identifying and quantifying L. infantum and Crithidia sp. LVH60A in samples obtained experimentally (in vitro and in vivo) or collected from hosts (humans, dogs, cats, and vectors). Importantly, the screening of 62 cultured isolates from VL patients using these primers surprisingly revealed that 51 parasite cultures were PCR+ for Crithidia sp. In addition, qPCR assays identified the co-infection of L. infantum with Crithidia sp. LVH60A in two new VL cases in Brazil, confirming the suspicion of co-infection in a previously reported case of fatal VL. We believe that the species-specific genes targeted in this study can be helpful for the molecular diagnosis of VL, as well as for elucidating suspected co-infections with monoxenous-like trypanosomatids, which is a neglected fact of a neglected disease.

Leishmania species: A narrative review on surface proteins with structural aspects involved in host-pathogen interaction

In tropical and subtropical regions of the world, leishmaniasis is endemic and causes a range of clinical symptoms in people, from severe tegumentary forms (such as cutaneous, mucocutaneous, and diffuse leishmaniasis) to lethal visceral forms. The protozoan parasite of the genus Leishmania causes leishmaniasis, which is still a significant public health issue, according to the World Health Organization 2022. The public's worry about the neglected tropical disease is growing as new foci of the illness arise, which are exacerbated by alterations in behavior, changes in the environment, and an enlarged range of sand fly vectors. Leishmania research has advanced significantly during the past three decades in a few different avenues. Despite several studies on Leishmania, many issues, such as illness control, parasite resistance, parasite clearance, etc., remain unresolved. The key virulence variables that play a role in the pathogenicity-host-pathogen relationship of the parasite are comprehensively discussed in this paper. The important Leishmania virulence factors, such as Kinetoplastid Membrane Protein-11 (KMP-11), Leishmanolysin (GP63), Proteophosphoglycan (PPG), Lipophosphoglycan (LPG), Glycosylinositol Phospholipids (GIPL), and others, have an impact on the pathophysiology of the disease and enable the parasite to spread the infection. Leishmania infection may arise from virulence factors; they are treatable with medications or vaccinations more promptly and might greatly shorten the duration of treatment. Additionally, our research sought to present a modeled structure of a few putative virulence factors that might aid in the development of new chemotherapeutic approaches for the treatment of leishmaniasis. The predicted virulence protein's structure is utilized to design novel drugs, therapeutic targets, and immunizations for considerable advantage from a higher understanding of the host immune response.

Preliminary Information of Iranian Lizard Leishmania Promastigote Transcriptome Sequencing by Next Generation Sequencing (NGS) Method

Background

A lizard Leishmania has been isolated from a lizard (Agama agilis) in Iran. Its genome sequence has not been determined, so far.

Methods

The study was done at Shahid Beheshti University of Medical Sciences, Tehran, Iran in 2017-2023. Leishmania promastigotes were cultured in RPMI1640 culture medium and collected at logarithmic phase by centrigugation. Parasite RNA was extracted by the Qiagene standard kit and its quantity and quality was determined and sequenced by NGS method with Illumina PE machine at BGI Company (China).

Results

The number of 8316 mRNA, 83 tRNA, 63 rRNA, 83 ncRNA, 5 snRNA, 1039 snoRNA, 36 region, and 3 repeat regions, 8343 CDS, 9597 Exon and 9292 Genes were identified in promastigote of Iranian lizard Leishmania.

Conclusion

Genomic elements of Iranian lizards Leishmania (with unique characteristics) were determined and identified by NGS system.

Identification of Immunodominant Proteins of the Leishmania (Viannia) naiffi SubProteome as Pan-Specific Vaccine Targets against Leishmaniasis

Leishmaniasis is a wide-spectrum disease caused by parasites from Leishmania genus. A well-modulated immune response that is established after the long-lasting clinical cure of leishmaniasis can represent a standard requirement for a vaccine. Previous studies demonstrated that Leishmania (Viannia) naiffi causes benign disease and its antigens induce well-modulated immune responses in vitro. In this work we aimed to identify the immunodominant proteins present in the soluble extract of L. naiffi (sLnAg) as candidates for composing a pan-specific anti-leishmaniasis vaccine. After immunoblotting using cured patients of cutaneous leishmaniasis sera and proteomics approaches, we identified a group of antigenic proteins from the sLnAg. In silico analyses allowed us to select mildly similar proteins to the host; in addition, we evaluated the binding potential and degree of promiscuity of the protein epitopes to HLA molecules and to B-cell receptors. We selected 24 immunodominant proteins from a sub-proteome with 328 proteins. Homology analysis allowed the identification of 13 proteins with the most orthologues among seven Leishmania species. This work demonstrated the potential of these proteins as promising vaccine targets capable of inducing humoral and cellular pan-specific immune responses in humans, which may in the future contribute to the control of leishmaniasis.

Hit-to-lead optimization of a pyrazinylpiperazine series against Leishmania infantum and Leishmania braziliensis

An early hit-to-lead optimization of a novel pyrazinylpiperazine series against L. infantum and L. braziliensis has been performed after an extensive SAR focusing on the benzoyl fragment of hit (4). Deletion of the meta-Cl of (4) led to the obtention of the para-hydroxyl derivative (12), on which the design of most monosubstituted derivatives of the SAR was based. Further optimization of the series, involving disubstituted benzoyl fragments and the hydroxyl substituent of (12), allowed the obtention of a total of 15 compounds with increased antileishmanial potency (IC₅₀ < 10 μM), nine of which displayed activity in the low micromolar range (IC₅₀ < 5 μM). This optimization ultimately identified the ortho, meta-dihydroxyl derivative (46) as an early lead for this series (IC₅₀ (L. infantum) = 2.8 μM, IC₅₀ (L. braziliensis) = 0.2 μM). Additional assessment of some selected compounds against other trypanosomatid parasites revealed that this series is selective towards Leishmania parasites, and in silico ADMET predictions revealed satisfactory profiles for these compounds, allowing further lead optimization of the pyrazinylpiperazine class against Leishmania.

Leishmania infantum (JPCM5) Transcriptome, Gene Models and Resources for an Active Curation of Gene Annotations

Leishmania infantum is one of the causative agents of visceral leishmaniases, the most severe form of leishmaniasis. An improved assembly for the L. infantum genome was published five years ago, yet delineation of its transcriptome remained to be accomplished. In this work, the transcriptome annotation was attained by a combination of both short and long RNA-seq reads. The good agreement between the results derived from both methodologies confirmed that transcript assembly based on Illumina RNA-seq and further delimitation according to the positions of spliced leader (SAS) and poly-A (PAS) addition sites is an adequate strategy to annotate the transcriptomes of Leishmania, a procedure previously used for transcriptome annotation in other Leishmania species and related trypanosomatids. These analyses also confirmed that the Leishmania transcripts boundaries are relatively slippery, showing extensive heterogeneity at the 5′- and 3′-ends. However, the use of RNA-seq reads derived from the PacBio technology (referred to as Iso-Seq) allowed the authors to uncover some complex transcription patterns occurring at particular loci that would be unnoticed by the use of short RNA-seq reads alone. Thus, Iso-Seq analysis provided evidence that transcript processing at particular loci would be more dynamic than expected. Another noticeable finding was the observation of a case of allelic heterozygosity based on the existence of chimeric Iso-Seq reads that might be generated by an event of intrachromosomal recombination. In addition, we are providing the L. infantum gene models, including both UTRs and CDS regions, that would be helpful for undertaking whole-genome expression studies. Moreover, we have built the foundations of a communal database for the active curation of both gene/transcript models and functional annotations for genes and proteins.

The glycoprotein gp63- a potential pan drug target for developing new antileishmanial agents

Leishmaniasis is a tropical parasitic disease caused by Leishmania spp. They cause several presentations of illness ranging from cutaneous leishmaniasis to visceral leishmaniasis. The current arsenal of drugs to treat leishmaniasis is limited, and drug resistance further impedes the problem. Therefore, it is necessary to revisit the available information to identify an alternative or new target for treatment. The glycoprotein 63 (gp63), is a potential anti-leishmanial target that plays a significant role in host-pathogen interaction and virulence. Many studies are ongoing to develop gp63 inhibitors or use it as a vaccine target. In this review, we will discuss the potential of gp63 as a drug target. This review summarises the studies focusing on gp63 as a drug target and its inhibitors identified using in silico approaches.

Immunoproteomics approach for the discovery of antigens applied to the diagnosis of canine visceral leishmaniasis

In the present study, an immunoproteomic approach using Leishmania infantum parasites isolated from naturally infected dogs from an endemic region of the disease, was carried out to identify new antigens to be used in the diagnosis of canine visceral leishmaniasis (CVL). Protein extracts, obtained from parasites isolated from asymptomatic (CanLA) and symptomatic (CanLS) dogs, were used to perform the two-dimensional gels. Western Blotting assays were carried out by employing a pool of sera from dogs with visceral leishmaniasis (CanLA or CanLS), healthy dogs from an endemic area, or dogs with similar diseases associated with cross-reactions (babesiosis and ehrlichiosis). With these results, it was possible to exclude the spots that showed a cross-reactivity of the sera from groups of healthy dogs, and those with babesiosis or ehrlichiosis. Taken together, 20 proteins were identified, 15 of which have already been described in the literature and 5 of which are hypothetical. An immunogenomic screen strategy was applied to identify conserved linear B-cell epitopes in the identified hypothetical proteins. Two peptides were synthesized and tested in ELISA experiments as a proof of concept for the validation of our immunoproteomics findings. The results demonstrated that the antigens presented sensitivity and specificity values ranging from 81.93% to 97.59% and 78.14 to 85.12%, respectively. As a comparative antigen, a preparation of a Leishmania extract showed sensitivity and specificity values of 75.90% and 74.88%, respectively. The present study was able to identify proteins capable of being used for the serodiagnosis of canine visceral leishmaniasis.

The paradigm of intracellular parasite survival and drug resistance in leishmanial parasite through genome plasticity and epigenetics: Perception and future perspective

Leishmania is an intracellular, zoonotic, kinetoplastid eukaryote with more than 1.2 million cases all over the world. The leishmanial chromosomes are divided into polymorphic chromosomal ends, conserved central domains, and antigen-encoding genes found in telomere-proximal regions. The genome flexibility of chromosomal ends of the leishmanial parasite is known to cause drug resistance and intracellular survival through the evasion of host defense mechanisms. Therefore, in this review, we discuss the plasticity of Leishmania genome organization which is the primary cause of drug resistance and parasite survival. Moreover, we have not only elucidated the causes of such genome plasticity which includes aneuploidy, epigenetic factors, copy number variation (CNV), and post-translation modification (PTM) but also highlighted their impact on drug resistance and parasite survival.

Revisiting epidemiology of leishmaniasis in central Asia: lessons learnt

In this work we reviewed historical and recent data on Leishmania spp. infection combining data collected in Turkmenistan, Uzbekistan, Kazakhstan, Kyrgyzstan, Iran, China and Mongolia. We specifically focused on a complex of co-existing species (Leishmania major, Leishmania turanica and Leishmania gerbilli) sharing the same animal reservoirs and vectors. In addition, we analysed the presence of dsRNA viruses in these species and discussed future research directions to identify species-specific traits, which may determine susceptibility of different Leishmania spp. to viral infection.

Review on the Drug Intolerance and Vaccine Development for the Leishmaniasis

Leishmaniasis is one of the Neglected Tropical Diseases (NTDs), a zoonotic disease of vector-borne nature that is caused by a protozoan parasite Leishmania. This parasite is transmitted by the vector sandfly into the human via a bite. Visceral leishmaniasis (VL), also called kala-azar, is the most fatal among the types of leishmaniasis, with high mortality mostly spread in the East Africa and South Asia regions. WHO report stated that approximately 3.3 million disabilities occur every year due to the disease along with approximately 50,000 annual deaths. The real matter of concern is that there is no particular effective medicine/vaccine available against leishmaniasis to date except a few approved drugs and chemotherapy for the infected patient. The current selection of small compounds was constrained, and their growing drug resistance had been a major worry. Additionally, the serious side effects on humans of the available therapy or drugs have made it essential to discover efficient and low-cost methods to speed up the development of new drugs against leishmaniasis. Ideally, the vaccine could be a low risk and effective alternative for both CL and VL and elicit long-lasting immunity against the disease. There are a number of vaccine candidates at various stages of clinical development and preclinical stage. However, none has successfully passed all clinical trials. But, the successful development and approval of commercially available vaccines for dogs against canine leishmaniasis (CanL) provides evidence that it can be possible for humans in distant future. In the present article, the approaches used for the development of vaccines for leishmaniasis are discussed and the progress being made is briefly reviewed.

Polyunsaturated fatty acids alter the formation of lipid droplets and eicosanoid production in Leishmania promastigotes

BACKGROUND

The knowledge about eicosanoid metabolism and lipid droplet (LD) formation in the Leishmania is very limited and new approaches are needed to identify which bioactive molecules are produced of them.

OBJECTIVES

Herein, we compared LDs and eicosanoids biogenesis in distinct Leishmania species which are etiologic agents of different clinical forms of leishmaniasis.

METHODS

For this, promastigotes of Leishmania amazonensis, L. braziliensis and L. infantum were stimulated with polyunsaturated fatty acids (PUFA) and LD and eicosanoid production was evaluated. We also compared mutations in structural models of human-like cyclooxygenase-2 (GP63) and prostaglandin F synthase (PGFS) proteins, as well as the levels of these enzymes in parasite cell extracts.

FINDINGS

PUFAs modulate the LD formation in L. braziliensis and L. infantum. Leishmania spp with equivalent tissue tropism had same protein mutations in GP63 and PGFS. No differences in GP63 production were observed among Leishmania spp, however PGFS production increased during the parasite differentiation. Stimulation with arachidonic acid resulted in elevated production of hydroxyeicosatetraenoic acids compared to prostaglandins.

MAIN CONCLUSIONS

Our data suggest LD formation and eicosanoid production are distinctly modulated by PUFAS dependent of Leishmania species. In addition, eicosanoid-enzyme mutations are more similar between Leishmania species with same host tropism.

An RNA Interference (RNAi) Toolkit and Its Utility for Functional Genetic Analysis of Leishmania (Viannia)

RNA interference (RNAi) is a powerful tool whose efficacy against a broad range of targets enables functional genetic tests individually or systematically. However, the RNAi pathway has been lost in evolution by a variety of eukaryotes including most Leishmania sp. RNAi was retained in species of the Leishmania subgenus Viannia, and here we describe the development, optimization, and application of RNAi tools to the study of L. (Viannia) braziliensis (Lbr). We developed vectors facilitating generation of long-hairpin or "stem-loop" (StL) RNAi knockdown constructs, using Gateway^TM site-specific recombinase technology. A survey of applications of RNAi in L. braziliensis included genes interspersed within multigene tandem arrays such as quinonoid dihydropteridine reductase (QDPR), a potential target or modulator of antifolate sensitivity. Other tests include genes involved in cell differentiation and amastigote proliferation (A600), and essential genes of the intraflagellar transport (IFT) pathway. We tested a range of stem lengths targeting the L. braziliensis hypoxanthine-guanine phosphoribosyl transferase (HGPRT) and reporter firefly luciferase (LUC) genes and found that the efficacy of RNAi increased with stem length, and fell off greatly below about 128 nt. We used the StL length dependency to establish a useful 'hypomorphic' approach not possible with other gene ablation strategies, with shorter IFT140 stems yielding viable cells with compromised flagellar morphology. We showed that co-selection for RNAi against adenine phosphoryl transferase (APRT1) using 4-aminopyrazolpyrimidine (APP) could increase the efficacy of RNAi against reporter constructs, a finding that may facilitate improvements in future work. Thus, for many genes, RNAi provides a useful tool for studying Leishmania gene function with some unique advantages.

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.

Two Leishmania species separation targeting the ITS-rDNA and Cyt b genes by developing and evaluating HRM- qPCR

BACKGROUND

Incidence of Cutaneous Leishmaniasis as an infectious and neglected disease is increasing, for the diagnosis of which several traditional methods and conventional PCR techniques have been developed, employing different genes for species identification.

METHODS

Leishmania parasites were sampled, DNA was extracted, and new specific and sensitive primers were designed. Two ITS-rDNA and Cyt b genes were targeted by qPCR using the High- Resolution Melting method to identify Leishmania parasites. The standard curves were drawn, compared, and identified by high-resolution melting curve analysis.

RESULTS

Melting temperature and Cycle of Threshold of ITS-rDNA was higher than Cyt b but Cyt b was more sensitive than ITS-rDNA when Leishmania major and Leishmania tropica were analyzed and evaluated. By aligning melt curves, normalizing fluorescence curves, and difference plotting melt curves, each Leishmania species was distinguished easily. L. major and L. tropica were separated at 83.6 °C and 84.7 °C, respectively, with less than 0.9 °C of temperature difference. Developing sensitivity and specificity of real-time PCR based on EvaGreen could detect DNA concentration to less than one pmol.

CONCLUSIONS

Precise identification of Leishmania parasites is crucial for strategies of disease control. Real-time PCR using EvaGreen provides rapid, highly sensitive, and specific detection of parasite's DNA. The modified High-Resolution Melting could determine unique curves and was able to detect single nucleotide polymorphisms according to small differences in the nucleotide content of Leishmania parasites.

Thiophene derivatives activity against the protozoan parasite Leishmania infantum

Treatments against leishmaniasis are limited and the development of new molecules is crucial. One class of developmental drug that has shown activity against the parasite Leishmania are thiophene derivatives. Here we synthetized thirty-eight novel thiophene compounds and characterized their activity and potential for resistance against L. infantum. Half of the molecules had an EC₅₀ in the low micromolar range, the piperidine derivatives being more potent than the tetramethylpyran derivatives. Resistance was challenging to select for, and resistant cells could only be raised against one (GC1-19) of the four most active compounds. Using chemogenomic screens we show that a gene conversion event at the ABCG2 locus as well as the overexpression of a tryparedoxin peroxidase are responsible for a weak but significant resistance to the GC1-19 drug candidate. Together, our results suggest that thiophene is a scaffold of interest for further drug development against leishmaniasis.

Three types of Leishmania mexicana amastigotes: Proteome comparison by quantitative proteomic analysis

Leishmania is the unicellular parasite transmitted by phlebotomine sand fly bite. It exists in two different forms; extracellular promastigotes, occurring in the gut of sand flies, and intracellular, round-shaped amastigotes residing mainly in vertebrate macrophages. As amastigotes originating from infected animals are often present in insufficient quality and quantity, two alternative types of amastigotes were introduced for laboratory experiments: axenic amastigotes and amastigotes from macrophages infected in vitro. Nevertheless, there is very little information about the degree of similarity/difference among these three types of amastigotes on proteomic level, whose comparison is crucial for assessing the suitability of using alternative types of amastigotes in experiments. In this study, L. mexicana amastigotes obtained from lesion of infected BALB/c mice were proteomically compared with alternatively cultivated amastigotes (axenic and macrophage-derived ones). Amastigotes of all three types were isolated, individually treated and analysed by LC-MS/MS proteomic analysis with quantification using TMT10-plex isobaric labeling. Significant differences were observed in the abundance of metabolic enzymes, virulence factors and proteins involved in translation and condensation of DNA. The most pronounced differences were observed between axenic amastigotes and lesion-derived amastigotes, macrophage-derived amastigotes were mostly intermediate between axenic and lesion-derived ones.

The lipidome of Crithidia fasiculataand its plasticity

Crithidia fasiculata belongs to the trypanosomatidae order of protozoan parasites, bearing close relation to other kinetoplastid parasites such as Trypanosoma brucei and Leishmania spp. As an early diverging lineage of eukaryotes, the study of kinetoplastid parasites has provided unique insights into alternative mechanisms to traditional eukaryotic metabolic pathways. Crithidia are a monogenetic parasite for mosquito species and have two distinct lifecycle stages both taking place in the mosquito gut. These consist of a motile choanomastigote form and an immotile amastigote form morphologically similar to amastigotes in Leishmania. Owing to their close relation to Leishmania, Crithidia are a growing research tool, with continuing interest in its use as a model organism for kinetoplastid research with the added benefit that they are non-pathogenic to humans and can be grown with no special equipment or requirements for biological containment. Although comparatively little research has taken place on Crithidia, similarities to other kinetoplast species has been shown in terms of energy metabolism and genetics. Crithidia also show similarities to kinetoplastids in their production of the monosaccharide D-arabinopyranose similar to Leishmania, which is incorporated into a lipoarabinogalactan a major cell surface GPI-anchored molecule. Additionally, Crithidia have been used as a eukaryotic expression system to express proteins from other kinetoplastids and potentially other eukaryotes including human proteins allowing various co- and post-translational protein modifications to the recombinant proteins. Despite the obvious usefulness and potential of this organism very little is known about its lipid metabolism. Here we describe a detailed lipidomic analyses and demonstrate the possible placidity of Crithidia’s lipid metabolis. This could have important implications for biotechnology approaches and how other kinetoplastids interact with, and scavenge nutrients from their hosts.

Uncovering Pseudogenes and Intergenic Protein-coding Sequences in TriTryps' Genomes

Trypanosomatids belong to a remarkable group of unicellular, parasitic organisms of the order Kinetoplastida, an early diverging branch of the phylogenetic tree of eukaryotes, exhibiting intriguing biological characteristics affecting gene expression (intronless polycistronic transcription, trans-splicing, and RNA editing), metabolism, surface molecules, and organelles (compartmentalization of glycolysis, variation of the surface molecules, and unique mitochondrial DNA), cell biology and life cycle (phagocytic vacuoles evasion and intricate patterns of cell morphogenesis). With numerous genomic-scale data of several trypanosomatids becoming available since 2005 (genomes, transcriptomes, and proteomes), the scientific community can further investigate the mechanisms underlying these unusual features and address other unexplored phenomena possibly revealing biological aspects of the early evolution of eukaryotes. One fundamental aspect comprises the processes and mechanisms involved in the acquisition and loss of genes throughout the evolutionary history of these primitive microorganisms. Here, we present a comprehensive in silico analysis of pseudogenes in three major representatives of this group: Leishmania major, Trypanosoma brucei, and Trypanosoma cruzi. Pseudogenes, DNA segments originating from altered genes that lost their original function, are genomic relics that can offer an essential record of the evolutionary history of functional genes, as well as clues about the dynamics and evolution of hosting genomes. Scanning these genomes with functional proteins as proxies to reveal intergenic regions with protein-coding features, relying on a customized threshold to distinguish statistically and biologically significant sequence similarities, and reassembling remnant sequences from their debris, we found thousands of pseudogenes and hundreds of open reading frames, with particular characteristics in each trypanosomatid: mutation profile, number, content, density, codon bias, average size, single- or multi-copy gene origin, number and type of mutations, putative primitive function, and transcriptional activity. These features suggest a common process of pseudogene formation, different patterns of pseudogene evolution and extant biological functions, and/or distinct genome organization undertaken by those parasites during evolution, as well as different evolutionary and/or selective pressures acting on distinct lineages.

Genomic diversity and genetic variation of Leishmania panamensis within its endemic range

Species belonging to the Leishmania (Viannia) subgenus are important causative agents of cutaneous and mucocutaneous leishmaniasis in Central and South America. These parasites possess several distinctive biological features that are influenced by their genetics, population structure, and genome instability. To date, several studies have revealed varying degrees of genetic diversity within Leishmania species. Particularly, in species of the L. (Viannia) subgenus, a generalized high intraspecific genetic diversity has been reported, although, conflicting conclusions have been drawn using different molecular techniques. Despite being the most common Leishmania species circulating in Panama and Colombia, few studies have analyzed clinical samples of Leishmania panamensis using whole-genome sequencing, and their restricted number of samples has limited the information they can provide to understand the population structure of L. panamensis. Here, we used next generation sequencing (NGS) to explore the genetic diversity of L. panamensis within its endemic range, analyzing data from 43 isolates of Colombian and Panamanian origin. Our results show the occurrence of three well-defined geographically correlated groups, and suggests the possible occurrence of additional phylogeographic groups. Furthermore, these results support the existence of a mixed mode of reproduction in L. panamensis, with varying frequencies of events of genetic recombination occurring primarily within subpopulations of closely related strains. This study offers important insights into the population genetics and reproduction mode of L. panamensis, paving the way to better understand their population structure and the emergence and maintenance of key eco-epidemiological traits.

DNA methylation alterations caused by Leishmania infection may generate a microenvironment prone to tumour development

DNA methylation is an epigenetic signature consisting of a methyl group at the 5’ cytosine of CpG dinucleotides. Modifications in DNA methylation pattern have been detected in cancer and infectious diseases and may be associated with gene expression changes. In cancer development DNA methylation aberrations are early events whereas in infectious diseases these epigenetic changes may be due to host/pathogen interaction. In particular, in leishmaniasis, a parasitic disease caused by the protozoan Leishmania, DNA methylation alterations have been detected in macrophages upon infection with Leishmania donovani and in skin lesions from patients with cutaneous leishmaniasis. Interestingly, different types of cancers, such as cutaneous malignant lesions, lymphoma and hepatocellular carcinoma, have been diagnosed in patients with a history of leishmaniasis. In fact, it is known that there exists an association between cancer and infectious diseases. Leishmania infection may increase susceptibility to develop cancer, but the mechanisms involved are not entirely clear. Considering these aspects, in this review we discuss the hypothesis that DNA methylation alterations induced by Leishmania may trigger tumorigenesis in long term infection since these epigenetic modifications may enhance and accumulate during chronic leishmaniasis.

Impact of Visceral Leishmaniasis on Local Organ Metabolism in Hamsters

Leishmania is an intracellular parasite with different species pathogenic to humans and causing the disease leishmaniasis. Leishmania donovani causes visceral leishmaniasis (VL) that manifests as hepatosplenomegaly, fever, pancytopenia and hypergammaglobulinemia. If left without treatment, VL can cause death, especially in immunocompromised people. Current treatments have often significant adverse effects, and resistance has been reported in some countries. Determining the metabolites perturbed during VL can lead us to find new treatments targeting disease pathogenesis. We therefore compared metabolic perturbation between L. donovani-infected and uninfected hamsters across organs (spleen, liver, and gut). Metabolites were extracted, analyzed by liquid chromatography-mass spectrometry, and processed with MZmine and molecular networking to annotate metabolites. We found few metabolites commonly impacted by infection across all three sites, including glycerophospholipids, ceramides, acylcarnitines, peptides, purines and amino acids. In accordance with VL symptoms and parasite tropism, we found a greater overlap of perturbed metabolites between spleen and liver compared to spleen and gut, or liver and gut. Targeting pathways related to these metabolite families would be the next focus that can lead us to find more effective treatments for VL.

Making the Most of Its Short Reads: A Bioinformatics Workflow for Analysing the Short-Read-Only Data of Leishmania orientalis (Formerly Named Leishmania siamensis) Isolate PCM2 in Thailand

Simple Summary

Leishmaniasis is a parasitic disease caused by flagellated protozoa of the genus Leishmania. Multiple genome sequencing platforms have been employed to complete Leishmania genomes at the expense of high cost. This study proposes an integrative bioinformatic workflow for assembling only the short-read data of Leishmania orientalis isolate PCM2 from Thailand and produce an acceptable-quality genome for further genomic analysis. This workflow gives extensive information required for identifying strain-specific markers and virulence-associated genes useful for drug and vaccine development before a more exhaustive and expensive investigation.

Abstract

Background: Leishmania orientalis (formerly named Leishmania siamensis) has been neglected for years in Thailand. The genomic study of L. orientalis has gained much attention recently after the release of the first high-quality reference genome of the isolate LSCM4. The integrative approach of multiple sequencing platforms for whole-genome sequencing has proven effective at the expense of considerably expensive costs. This study presents a preliminary bioinformatic workflow including the use of multi-step de novo assembly coupled with the reference-based assembly method to produce high-quality genomic drafts from the short-read Illumina sequence data of L. orientalis isolate PCM2. Results: The integrating multi-step de novo assembly by MEGAHIT and SPAdes with the reference-based method using the L. enriettii genome and salvaging the unmapped reads resulted in the 30.27 Mb genomic draft of L. orientalis isolate PCM2 with 3367 contigs and 8887 predicted genes. The results from the integrated approach showed the best integrity, coverage, and contig alignment when compared to the genome of L. orientalis isolate LSCM4 collected from the northern province of Thailand. Similar patterns of gene ratios and frequency were observed from the GO biological process annotation. Fifty GO terms were assigned to the assembled genomes, and 23 of these (accounting for 61.6% of the annotated genes) showed higher gene counts and ratios when results from our workflow were compared to those of the LSCM4 isolate. Conclusions: These results indicated that our proposed bioinformatic workflow produced an acceptable-quality genome of L. orientalis strain PCM2 for functional genomic analysis, maximising the usage of the short-read data. This workflow would give extensive information required for identifying strain-specific markers and virulence-associated genes useful for drug and vaccine development before a more exhaustive and expensive investigation.

Tissue tropism: Is it an intrinsic characteristic of Leishmania species?

The genus Leishmania comprises a wide range of species, some of which are pathogenic to humans and each of which has a different tissue preference, resulting in one of the three clinical forms of human leishmaniasis: visceral, cutaneous, or mucocutaneous. Although, all pathogenic species are deposited intradermally in the mammalian host upon an infectious sand fly bite, only the viscerotropic strains can leave the skin and reach the internal organs. We assume that Leishmania tissue tropism is not only the result of Leishmania genetic determinism but is also governed by the interaction of the parasite with different vectorial and human host elements. To shed light on these elements and key steps determining the course of the infection, we describe throughout this review the disease's progression from the early stages of infection taking place in the skin to the late stages succeeding in the parasite's visceral dissemination. Hence, we address the question of Leishmania tropism, through providing relevant hypotheses and answers gathered from the literature.

Identification of Novel Inhibitor of Enoyl-Acyl Carrier Protein Reductase (InhA) Enzyme in Mycobacterium tuberculosis from Plant-Derived Metabolites: An In Silico Study

Mycobacterium tuberculosis (M.tb.) enoyl-acyl carrier protein (ACP) reductase (InhA) is validated as a useful target for tuberculosis therapy and is considered an attractive enzyme to drug discovery. This study aimed to identify the novel inhibitor of the InhA enzyme, a potential target of M.tb. involved in the type II fatty acid biosynthesis pathway that controls mycobacterial cell envelope synthesis. We compiled 80 active compounds from Ruta graveolens and citrus plants belonging to the Rutaceae family for pharmacokinetics and molecular docking analyses. The chemical structures of the 80 phytochemicals and the 3D structure of the target protein were retrieved from the PubChem database and RCSB Protein Data Bank, respectively. The evaluation of druglikeness was performed based on Lipinski’s Rule of Five, while the computed phytochemical properties and molecular descriptors were used to predict the ADMET of the compounds. Amongst these, 11 pharmacokinetically-screened compounds were further examined by performing molecular docking analysis with an InhA target using AutoDock 4.2. The docking results showed that gravacridonediol, a major glycosylated natural alkaloid from Ruta graveolens, might possess a promising inhibitory potential against InhA, with a binding energy (B.E.) of −10.80 kcal/mole and inhibition constant (Ki) of 600.24 nM. These contrast those of the known inhibitor triclosan, which has a B.E. of −6.69 kcal/mole and Ki of 12.43 µM. The binding efficiency of gravacridonediol was higher than that of the well-known inhibitor triclosan against the InhA target. The present study shows that the identified natural compound gravacridonediol possesses drug-like properties and also holds promise in inhibiting InhA, a key target enzyme of M.tb.

Reconstitution of Mycobacterium marinum Nonhomologous DNA End Joining Pathway in Leishmania

In mammalian cells, DNA double-strand breaks (DSBs) are mainly repaired by nonhomologous end joining (NHEJ) pathway. Ku (a heterodimer formed by Ku70 and Ku80 proteins) and DNA ligase IV are the core NHEJ factors. Ku could also be involved in other cellular processes, including telomere length regulation, DNA replication, transcription, and translation control. Leishmania, an early branching eukaryote and the causative agent of leishmaniasis, has no functional NHEJ pathway due to its lack of DNA ligase IV and other NHEJ factors but retains Ku70 and Ku80 proteins. In this study, we generated Leishmania donovani Ku70 disruption mutants and Ku70 and Ku80 double gene (Ku70/80) disruption mutants. We found that Leishmania Ku is still involved in DSB repair, possibly through its binding to DNA ends to block and slowdown 5′ end resections and Ku-Ku or other protein interactions. Depending on location of a DSB between the direct repeat genomic sequences, Leishmania Ku could have an inhibiting effect, no effect or a promoting effect on the DSB repair mediated by single strand annealing (SSA), the most frequently used DSB repair pathway in Leishmania. Ku70/80 proteins are also required for the healthy proliferation of Leishmania cells. Interestingly, unlike in Trypanosoma brucei and L. mexicana, Ku70/80 proteins are dispensable for maintaining the normal lengths of telomeres in L. donovani. We also show it is possible to reconstitute the two components (Ku and Ligase D) NHEJ pathway derived from Mycobacterium marinum in Leishmania. This improved DSB repair fidelity and efficiency in Leishmania and sets up an example that the bacterial NHEJ pathway can be successfully reconstructed in an NHEJ-deficient eukaryotic parasite.

IMPORTANCE Nonhomologous end joining (NHEJ) is the most efficient double-stranded DNA break (DSB) repair pathway in mammalian cells. In contrast, the protozoan parasite Leishmania has no functional NHEJ pathway but retains the core NHEJ factors of Ku70 and Ku80 proteins. In this study, we found that Leishmania Ku heterodimers are still participating in DSB repair possibly through blocking 5′ end resections and Ku-Ku protein interactions. Depending on the DSB location, Ku could have an inhibiting or promoting effect on DSB repair mediated by the single-strand annealing repair pathway. Ku is also required for the normal growth of the parasite but surprisingly dispensable for maintaining the telomere lengths. Further, we show it is possible to introduce Mycobacterium marinum NHEJ pathway into Leishmania. Understanding DSB repair mechanisms of Leishmania may improve the CRISPR gene targeting specificity and efficiency and help identify new drug targets for this important human parasite.

Development of an immunogen containing CD4+/CD8+ T-cell epitopes for the prophylaxis of tegumentary leishmaniasis

Abstract

Tegumentary leishmaniasis (TL) is a disease of high severity and incidence in Brazil, and Leishmania braziliensis is its main etiological agent. The inefficiency of control measures, such as high toxicity and costs of current treatments and the lack of effective immunoprophylactic strategies, makes the development of vaccines indispensable and imminent. In this light, the present work developed a gene encoding multiple T-cell (CD4⁺/CD8⁺) epitope, derived from conserved proteins found in Leishmania species and associated with TL, to generate a chimeric protein (rMEP/TL) and compose a vaccine formulation. For this, six T-cell epitopes were selected by immunoinformatics approaches from proteins present in the amastigote stage and associated with host-parasite interactions. The following formulations were then tested in an L. braziliensis murine infection model: rMEP/TL in saline or associated with MPLA-PHAD^®. Our data revealed that, after immunization (three doses; 14-day intervals) and subsequent challenging, rMEP/TL and rMEP/TL + MPLA-vaccinated mice showed an increased production of key immunological biomarkers of protection, such as IgG_2a, IgG_2a/IgG₁, NO, CD4⁺, and CD8⁺ T-cells with IFN-γ and TNF-α production, associated with a reduction in CD4⁺IL-10⁺ and CD8⁺IL-10⁺ T-cells. Vaccines also induced the development of central (CD44^highCD62L^high) and effector (CD44^highCD62L^low) memory of CD4⁺ and CD8⁺ T-cells. These findings, associated with the observation of lower rates of parasite burdens in the vaccinated groups, when compared to the control groups, suggest that immunization with rMEP/TL and, preferably, associated with an adjuvant, may be considered an effective tool to prevent TL.

Key points

• Rational design approaches for vaccine development.

• Central and effector memory of CD4+ and CD8+ T-cells.

• Vaccine comprised of rMEP/TL plus MPLA as an effective tool to prevent TL.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00253-022-12033-7.

Replacement of Leishmania (Leishmania) infantum Populations in an Endemic Focus of Visceral Leishmaniasis in Brazil

Visceral leishmaniasis is an important global health problem with an estimated of 50,000 to 90,000 new cases per year. VL is the most serious form of leishmaniasis as it can be fatal in 95% of the cases if it remains untreated. VL is a particularly acute problem in Brazil which contributed with 97% of all cases reported in 2020 in the Americas. In this country, VL affects mainly the poorest people in both urban and rural areas and continues to have a high mortality rate estimated around 8.15%. Here, we performed a temporal parasite population study using whole genome sequence data from a set of 34 canine isolates sampled in 2008, 2012 and 2015 from a re-emergent focus in Southeastern Brazil. Our study found the presence of two distinct sexual subpopulations that corresponded to two isolation periods. These subpopulations diverged hundreds of years ago with no apparent gene flow between them suggesting a process of rapid replacement during a two-year period. Sequence comparisons and analysis of nucleotide diversity also showed evidence of balancing selection acting on transport-related genes and antigenic families. To our knowledge this is the first population genomic study showing a turn-over of parasite populations in an endemic region for leishmaniasis. The complexity and rapid adaptability of these parasites pose new challenges to control activities and demand more integrated approaches to understand this disease in New World foci.

Paving the Way: Contributions of Big Data to Apicomplexan and Kinetoplastid Research

In the age of big data an important question is how to ensure we make the most out of the resources we generate. In this review, we discuss the major methods used in Apicomplexan and Kinetoplastid research to produce big datasets and advance our understanding of Plasmodium, Toxoplasma, Cryptosporidium, Trypanosoma and Leishmania biology. We debate the benefits and limitations of the current technologies, and propose future advancements that may be key to improving our use of these techniques. Finally, we consider the difficulties the field faces when trying to make the most of the abundance of data that has already been, and will continue to be, generated.

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.

Worldwide, over one million people are getting infected by the parasitic flagellates of the genus Leishmania annually leading to ~30,000 deaths. Notably, there is still no approved vaccine against human leishmaniases. A range of methods of forward and reverse genetics has recently been developed for several model Leishmania species. Unfortunately, these methods are often not transferrable to non-model species, which may be of even greater medical importance. Leishmania guyanensis is one of such cases. It frequently carries a symbiotic RNA virus that contributes to the development of a more aggressive form of leishmaniasis in an experimental murine model. In order to establish and optimize the system for genetic manipulations in L. guyanensis, we sequenced and annotated its genome. Next, we applied the CRISPR-Cas9 technology to target a gene of interest. This approach was validated by ablating a gene encoding a protein involved in lipid metabolism. In this work, we document that deletion of this gene leads to the disturbance of cell cycle and affects the ratio of critical intracellular lipids. We believe that our study will facilitate research into more effective treatment of leishmaniases.

Geographic Origin and Vertical Transmission of Leishmania infantum Parasites in Hunting Hounds, United States

Vertical transmission of leishmaniasis is common but is difficult to study against the background of pervasive vector transmission. We present genomic data from dogs in the United States infected with Leishmania infantum parasites; these infections have persisted in the apparent absence of vector transmission. We demonstrate that these parasites were introduced from the Old World separately and more recently than L. infantum from South America. The parasite population shows unusual genetics consistent with a lack of meiosis: a high level of heterozygous sites shared across all isolates and no decrease in linkage with genomic distance between variants. Our data confirm that this parasite population has been evolving with little or no sexual reproduction. This demonstration of vertical transmission has profound implications for the population genetics of Leishmania parasites. When investigating transmission in complex natural settings, considering vertical transmission alongside vector transmission is vital.

In silico analysis of Leishmania proteomes and protein-protein interaction network: Prioritizing therapeutic targets and drugs for repurposing to treat leishmaniasis

Leishmaniasis is a serious world health problem and its current therapies have several limitations demanding to develop novel therapeutics for this disease. The present study aims to prioritize novel broad-spectrum targets using proteomics and protein-protein interaction network (PPIN) data for 11 Leishmania species. Proteome comparison and host non-homology analysis resulted in 3605 pathogen-specific conserved core proteins. Gene ontology analysis indicated their involvement in major molecular functions like DNA binding, transportation, dioxygenase, and catalytic activity. PPIN analysis of these core proteins identified eight hub proteins (viz., vesicle-trafficking protein (LBRM2903_190011800), ribosomal proteins S17 (LBRM2903_34004790) and L2 (LBRM2903_080008100), eukaryotic translation initiation factor 3 (LBRM2903_350086700), replication factor A (LBRM2903_150008000), U3 small nucleolar RNA-associated protein (LBRM2903_340025600), exonuclease (LBRM2903_200021800), and mitochondrial RNA ligase (LBRM2903_200074100)). Among the hub proteins, six were classified as drug targets and two as vaccine candidates. Further, druggability analysis indicated three hub proteins, namely eukaryotic translation initiation factor 3, ribosomal proteins S17 and L2 as druggable. Their three-dimensional structures were modelled and docked with the identified ligands (2-methylthio-N6-isopentenyl-adenosine-5'-monophosphate, artenimol and omacetaxine mepesuccinate). These ligands could be experimentally validated (in vitro and in vivo) and repurposed for the development of novel antileishmanial agents.

The History of Live Attenuated Centrin Gene-Deleted Leishmania Vaccine Candidates

Leishmaniasis, caused by an infection of the Leishmania protozoa, is a neglected tropical disease and a major health problem in tropical and subtropical regions of the world, with approximately 350 million people worldwide at risk and 2 million new cases occurring annually. Current treatments for leishmaniasis are not highly efficacious and are associated with high costs, especially in low- and middle-income endemic countries, and high toxicity. Due to a surge in the incidence of leishmaniases worldwide, the development of new strategies such as a prophylactic vaccine has become a high priority. However, the ability of Leishmania to undermine immune recognition has limited our efforts to design safe and efficacious vaccines against leishmaniasis. Numerous antileishmanial vaccine preparations based on DNA, subunit, and heat-killed parasites with or without adjuvants have been tried in several animal models but very few have progressed beyond the experimental stage. However, it is known that people who recover from Leishmania infection can be protected lifelong against future infection, suggesting that a successful vaccine requires a controlled infection to develop immunologic memory and subsequent long-term immunity. Live attenuated Leishmania parasites that are non-pathogenic and provide a complete range of antigens similarly to their wild-type counterparts could evoke such memory and, thus, would be effective vaccine candidates. Our laboratory has developed several live attenuated Leishmania vaccines by targeted centrin gene disruptions either by homologous recombination or, more recently, by using genome editing technologies involving CRISPR-Cas9. In this review, we focused on the sequential history of centrin gene-deleted Leishmania vaccine development, along with the characterization of its safety and efficacy. Further, we discussed other major considerations regarding the transition of dermotropic live attenuated centrin gene-deleted parasites from the laboratory to human clinical trials.

In-Depth Quantitative Proteomics Characterization of In Vitro Selected Miltefosine Resistance in Leishmania infantum

Visceral leishmaniasis (VL) is a neglected disease caused by Leishmania parasites. Although significant morbidity and mortality in tropical and subtropical regions of the world are associated with VL, the low investment for developing new treatment measures is chronic. Moreover, resistance and treatment failure are increasing for the main medications, but the emergence of resistance phenotypes is poorly understood at the protein level. Here, we analyzed the development of resistance to miltefosine upon experimental selection in a L. infantum strain. Time to miltefosine resistance emergence was ~six months and label-free quantitative mass-spectrometry-based proteomics analyses revealed that this process involves a remodeling of components of the membrane and mitochondrion, with significant increase in oxidative phosphorylation complexes, particularly on complex IV and ATP synthase, accompanied by increased energy metabolism mainly dependent on β-oxidation of fatty acids. Proteins canonically involved in ROS detoxification did not contribute to the resistant process whereas sterol biosynthesis enzymes could have a role in this development. Furthermore, changes in the abundance of proteins known to be involved in miltefosine resistance such as ABC transporters and phospholipid transport ATPase were detected. Together, our data show a more complete picture of the elements that make up the miltefosine resistance phenotype in L. infantum.

Comparative Draft Genomes of Leishmania orientalis Isolate PCM2 (Formerly Named Leishmania siamensis) and Leishmania martiniquensis Isolate PCM3 from the Southern Province of Thailand

Simple Summary

This study successfully sequenced the draft genomes of the southern isolates of Leishmania orientalis and Leishmania martiniquensis in Thailand. The comparison with the genomes of the northern isolates revealed species-level similarity with a level of genome and proteome variation, suggesting the different emerging strains. Comparing the proteins of these southern strains to those of the northern ones and 14 other Leishmania species showed six protein groups with numerous unique proteins: 53 for the southern strain PCM2 of L. orientalis and 97 for the strain PCM3 of L. martiniquensis. Some of these proteins were related to virulence, drug resistance, drug target, and stress response, which could be targeted for further experimental characterization. Therefore, the findings could initiate further genetic and population genomic investigation, and the close monitoring of L. orientalis and L. martiniquensis in Thailand and neighboring regions.

Abstract

(1) Background: Autochthonous leishmaniasis, a sandfly-borne disease caused by the protozoan parasites Leishmania orientalis (formerly named Leishmania siamensis) and Leishmania martiniquensis, has been reported for immunocompromised and immunocompetent patients in the southern province of Thailand. Apart from the recent genomes of the northern isolates, limited information is known on the emergence and genetics of these parasites. (2) Methods: This study sequenced and compared the genomes of L. orientalis isolate PCM2 and L. martiniquensis isolate PCM3 with those of the northern isolates and other 14 Leishmania species using short-read whole-genome sequencing methods and comparative bioinformatic analyses. (3) Results: The genomes of the southern isolates of L. orientalis and L. martiniquensis were 30.01 Mbp and 32.39 Mbp, and the comparison with the genomes of the northern isolates revealed species-level similarity with a level of genome and proteome variation, suggesting the different strains. Comparative proteome analysis showed six protein groups with 53 unique proteins for the strain PCM2 and 97 for the strain PCM3. Certain proteins were related to virulence, drug resistance, and stress response. (4) Conclusion: Therefore, the findings could indicate the need for more genetic and population genomic investigation, and the close monitoring of L. orientalis and L. martiniquensis in Thailand and neighboring regions.

Systematic identification of genes encoding cell surface and secreted proteins that are essential for in vitro growth and infection in Leishmania donovani

Leishmaniasis is an infectious disease caused by protozoan parasites belonging to the genus Leishmania for which there are no approved human vaccines. Infections localise to different tissues in a species-specific manner with the visceral form of the disease caused by Leishmania donovani and L. infantum being the most deadly in humans. Although Leishmania spp. parasites are predominantly intracellular, the visceral disease can be prevented in dogs by vaccinating with a complex mixture of secreted products from cultures of L. infantum promastigotes. With the logic that extracellular parasite proteins make good subunit vaccine candidates because they are directly accessible to vaccine-elicited host antibodies, here we attempt to discover proteins that are essential for in vitro growth and host infection with the goal of identifying subunit vaccine candidates. Using an in silico analysis of the Leishmania donovani genome, we identified 92 genes encoding proteins that are predicted to be secreted or externally anchored to the parasite membrane by a single transmembrane region or a GPI anchor. By selecting a transgenic L. donovani parasite that expresses both luciferase and the Cas9 nuclease, we systematically attempted to target all 92 genes by CRISPR genome editing and identified four that were required for in vitro growth. For fifty-five genes, we infected cohorts of mice with each mutant parasite and by longitudinally quantifying parasitaemia with bioluminescent imaging, showed that nine genes had evidence of an attenuated infection although all ultimately established an infection. Finally, we expressed two genes as full-length soluble recombinant proteins and tested them as subunit vaccine candidates in a murine preclinical infection model. Both proteins elicited significant levels of protection against the uncontrolled development of a splenic infection warranting further investigation as subunit vaccine candidates against this deadly infectious tropical disease.

Leishmaniasis is a parasitic infectious disease that is responsible for many tens of thousands of human deaths per year, primarily in impoverished parts of the world. Although there are drugs to treat this parasite infection, resistance is emerging and there are no approved human vaccines. Extracellular parasite proteins can make good vaccine targets because they are directly accessible to host antibodies; however, not all parasite surface proteins can elicit protective immune responses. With the goal of identifying new vaccine targets, we selected over 90 genes that encode parasite cell surface and secreted proteins and used the latest CRISPR gene editing technology to individually target them. Using these mutant parasites, we identified four genes required for parasite growth in the laboratory. We expressed two of the proteins as subunit vaccines and a preclinical infection model was used to determine if they could elicit protective immune responses. We found that two of our candidates were able to confer significant levels of protection in a murine model of visceral leishmaniasis. Our study will contribute to the search for a highly effective vaccine against visceral leishmaniasis to improve the lives of people living in some of the poorest regions on the planet.

Targeted Deletion of Centrin in Leishmania braziliensis Using CRISPR-Cas9-Based Editing

Leishmania braziliensis is the main causative agent of Tegumentary Leishmaniasis in the Americas. However, difficulties related to genome manipulation, experimental infection, and parasite growth have so far limited studies with this species. CRISPR-Cas9-based technology has made genome editing more accessible, and here we have successfully employed the LeishGEdit approach to attenuate L. braziliensis. We generated a transgenic cell line expressing Cas9 and T7 RNA polymerase, which was employed for the targeted deletion of centrin, a calcium-binding cytoskeletal protein involved in the centrosome duplication in eukaryotes. Centrin-deficient Leishmania exhibit growth arrest at the amastigote stage. Whole-genome sequencing of centrin-deficient L. braziliensis (LbCen^−/−) did not indicate the presence of off-target mutations. In vitro, the growth rates of LbCen^−/− and wild-type promastigotes were similar, but axenic and intracellular LbCen^−/− amastigotes showed a multinucleated phenotype with impaired survival following macrophage infection. Upon inoculation into BALB/c mice, LbCen^−/− were detected at an early time point but failed to induce lesion formation, contrary to control animals, infected with wild-type L. braziliensis. A significantly lower parasite burden was also observed in mice inoculated with LbCen^−/−, differently from control mice. Given that centrin-deficient Leishmania sp. have become candidates for vaccine development, we propose that LbCen^−/− can be further explored for the purposes of immunoprophylaxis against American Tegumentary Leishmaniasis.

Mutation Characteristics and Phylogenetic Analysis of Five Leishmania Clinical Isolates

Simple Summary

Leishmaniasis, a neglected tropical disease, is caused by infection with the Leishmania species, threatening millions of people in approximately 100 endemic countries. The emergence of antimony-resistant Leishmania strains have brought difficulties to the treatment and elimination of leishmaniasis. This study performed genome-wide resequencing and phylogenetic analysis of five isolates from the Leishmania donovani complex, focusing on finding mutations related to antimony resistance and virulence of the newly isolated Leishmania strain L_HCZ in 2016. By combining whole-genome sequencing and whole-genome phylogenetic analysis, Leishmania isolates L_801, L_9044 and L_Liu were identified as Leishmania donovani, and L_HCZ as Leishmania infantum. By discovering genome-wide single-nucleotide polymorphisms and structural variations, we identified mutations of drug resistance-related genes in the antimony-resistant Leishmania isolate L_HCZ. The new Leishmania isolate L_HCZ has strong virulence and strong drug resistance, which should be taken seriously by the relevant health departments and scientific researchers.

Abstract

Leishmaniasis is a neglected tropical disease threatening millions of people worldwide. The emergence of antimony-resistant Leishmania strains have brought difficulties to the treatment and elimination of leishmaniasis. This study performed genome sequencing, phylogenetic analysis and mutation analysis of five Leishmania clinical isolates, especially the Leishmania strain L_HCZ isolated in 2016, which shows strong virulence and antimony resistance. By phylogenetic analysis, four isolates (L_DD8, L_801, L_Liu and L_9044) were identified as Leishmania donovani, the isolate L_HCZ was identified as Leishmania infantum and the isolate L_DD8 as a standard strain of L.donovani. Genome-wide mutation analysis was applied to identify mutations related to the drug resistance and virulence of the newly isolated L_HCZ. Compared with the other four Leishmania isolates, L_HCZ had the most mutations in genes associated with antimony resistance, including the ABC transporter, ascorbate-dependent peroxidase, gamma–glutamylcysteine synthetase, glucose-6-phosphate 1-dehydrogenase, ATP-binding cassette protein subfamily A and multi-drug resistance protein-like genes. Among the genes associated with virulence, L_HCZ had the most mutations in cysteine peptidase A, cysteine peptidase B, cysteine peptidase C, heat-shock protein 70, gp63, acid phosphatase, kinesin k39, kinesin, phosphoglycan beta 1, amastin-like surface protein and amastin-like proteins. The mutations in L_HCZ might possibly contribute to its antimony resistance and strong virulence in clinical patients. Whole-genome resequencing has exhibited broad application prospects and may be put into clinical use in the future for parasite identifying and epidemiological investigations.

Impact of Genetic Diversity and Genome Plasticity of Leishmania spp. in Treatment and the Search for Novel Chemotherapeutic Targets

Leishmaniasis is one of the major public health concerns in Latin America, Africa, Asia, and Europe. The absence of vaccines for human use and the lack of effective vector control programs make chemotherapy the main strategy to control all forms of the disease. However, the high toxicity of available drugs, limited choice of therapeutic agents, and occurrence of drug-resistant parasite strains are the main challenges related to chemotherapy. Currently, only a small number of drugs are available for leishmaniasis treatment, including pentavalent antimonials (Sb^V), amphotericin B and its formulations, miltefosine, paromomycin sulphate, and pentamidine isethionate. In addition to drug toxicity, therapeutic failure of leishmaniasis is a serious concern. The occurrence of drug-resistant parasites is one of the causes of therapeutic failure and is closely related to the diversity of parasites in this genus. Owing to the enormous plasticity of the genome, resistance can occur by altering different metabolic pathways, demonstrating that resistance mechanisms are multifactorial and extremely complex. Genetic variability and genome plasticity cause not only the available drugs to have limitations, but also make the search for new drugs challenging. Here, we examined the biological characteristics of parasites that hinder drug discovery.

A profile of research on the parasitic trypanosomatids and the diseases they cause

The parasitic trypanosomatids cause lethal and debilitating diseases, the leishmaniases, Chagas disease, and the African trypanosomiases, with major impacts on human and animal health. Sustained research has borne fruit by assisting efforts to reduce the burden of disease and by improving our understanding of fundamental molecular and cell biology. But where has the research primarily been conducted, and which research areas have received the most attention? These questions are addressed below using publication and citation data from the past few decades.

Metabolic characterization and biomarkers screening for visceral leishmaniasis in golden hamsters

A better understanding of the changes in metabolic molecules during visceral leishmaniasis (VL) is essential to develop new strategies for diagnosis and treatment. Previous metabolomics studies on Leishmania have increased our knowledge of leishmaniasis and its causative pathogen. As these studies were mainly carried out in vitro, to go further, we conducted this global metabolomics analysis on the serum of golden hamsters. Serum samples were detected over a time course of 2, 4, 8 and 12 weeks post infection. Our results revealed that under extensively disturbed metabolomes between the infection group and controls, glycerophospholipid (GPL) metabolism was most affected over the infection time, followed by α-linoleic acid metabolism and arachidonic acid metabolism. Within GPLs, phosphatidylcholine (PC) and phosphatidylethanolamine (PE) were found to be significantly increased, while their enzyme-catalysed metabolites lysophosphatidylcholine (LPC) and lysophosphatidylethanolamine (LPE) showed no significant changes. Moreover, eight differential metabolites were selected. The ability of these metabolites to be used as a diagnostic biomarker panel was supported by receiver operating characteristic (ROC) analysis. Our findings revealed that GPL metabolism might play an important role in the response of the host to Leishmania infection. The metabolism of PC and PE might be crucial in the in vivo progression of VL. The panel of eight potential biomarkers might contribute to the diagnosis of VL.