Comparative Life Cycle Transcriptomics Revises Leishmania mexicana Genome Annotation and Links a Chromosome Duplication with Parasitism of Vertebrates

Innovation through instability? Genome (dis)organisation in Trypanosoma cruzi

Chagas disease affects millions globally and is caused by Trypanosoma cruzi, a parasite with a remarkable ability to adapt and persist in complex environments. Emerging evidence suggests that this adaptability may be driven by its extensive genome instability yet the underlying mechanisms remain poorly understood. Drawing insights from related trypanosomatids, we explore how processes such as aneuploidy, copy number variations and genetic rearrangements can be strategically exploited for immune evasion, host adaptation and drug resistance, and consider whether similar processes shape the adaptive potential of T. cruzi. Unravelling the fundamental biology governing genome instability in this parasite will undoubtedly advance our understanding of its evolution, pathogenesis and resilience - crucial steps towards reducing the burden of Chagas disease.

Tubulin detyrosination shapes Leishmania cytoskeletal architecture and virulence

Tubulin detyrosination has been implicated in various human disorders and is important for regulating microtubule dynamics. While in most organisms this modification is restricted to α-tubulin, in trypanosomatid parasites, it occurs on both α- and β-tubulin. Here, we show that in Leishmania, a single vasohibin (LmVASH) enzyme is responsible for differential kinetics of α- and β-tubulin detyrosination. LmVASH knockout parasites, which are completely devoid of detyrosination, show decreased levels of glutamylation and exhibit a strongly diminished pathogenicity in mice, correlating with decreased proliferation in macrophages. Reduced virulence is associated with altered morphogenesis and flagellum remodeling in detyrosination-deficient amastigotes. Flagellum shortening in the absence of detyrosination is caused by hyperactivity of a microtubule-depolymerizing Kinesin-13 homolog, demonstrating its function as a key reader of the trypanosomatid-tubulin code. Taken together, our work establishes the importance of tubulin detyrosination in remodeling the microtubule-based cytoskeleton required for efficient proliferation in the mammalian host. This highlights tubulin detyrosination as a potential target for therapeutic action against leishmaniasis.

A novel strain of Leishmania braziliensis harbors not a toti- but a bunyavirus

Leishmania is a genus of the family Trypanosomatidae that unites obligatory parasitic flagellates causing a variety of vector-borne diseases collectively called leishmaniasis. The symptoms range from relatively innocuous skin lesions to complete failures of visceral organs. The disease is exacerbated if a parasite harbors Leishmania RNA viruses (LRVs) of the family Pseudototiviridae. Screening a novel isolate of L. braziliensis, we revealed that it possesses not a toti-, but a bunyavirus of the family Leishbuviridae. To the best of our knowledge, this is a very first discovery of a bunyavirus infecting a representative of the Leishmania subgenus Viannia. We suggest that these viruses may serve as potential factors of virulence in American leishmaniasis and encourage researchers to test leishmanial strains for the presence of not only LRVs, but also other RNA viruses.

LeishGEM Team. LeishGEM: genome-wide deletion mutant fitness and protein localisations in Leishmania

LeishGEM is a genome-wide functional annotation community resource for Leishmania mexicana, where deletion mutant growth in vitro and in vivo is measured and protein localisation is determined by endogenous tagging and LOPIT-DC (localisation of organelle proteins by isotope tagging with differential centrifugation) spatial proteomics. Data are being made available pre-publication via http://leishgem.org which allows data-driven identification of the mechanisms for Leishmania parasitism.

Formation of functional E3 ligase complexes with UBC2 and UEV1 of Leishmania mexicana

In eukaryotic cells, molecular fate and cellular responses are shaped by multicomponent enzyme systems which reversibly attach ubiquitin and ubiquitin-like modifiers to target proteins. The extent of the ubiquitin proteasome system in Leishmania mexicana and its importance for parasite survival has recently been established through deletion mutagenesis and life-cycle phenotyping studies. The ubiquitin conjugating E2 enzyme UBC2, and the E2 enzyme variant UEV1, with which it forms a stable complex in vitro, were shown to be essential for the differentiation of promastigote parasites to the infectious amastigote form. To investigate further, we used immunoprecipitation of Myc-UBC2 or Myc-UEV1 to identify interacting proteins in L. mexicana promastigotes. The interactome of UBC2 comprises multiple ubiquitin-proteasome components including UEV1 and four RING E3 ligases, as well as potential substrates predicted to have roles in carbohydrate metabolism and intracellular trafficking. The smaller UEV1 interactome comprises six proteins, including UBC2 and shared components of the UBC2 interactome consistent with the presence of intracellular UBC2-UEV1 complexes. Recombinant RING1, RING2 and RING4 E3 ligases were shown to support ubiquitin transfer reactions involving the E1, UBA1a, and UBC2 to available substrate proteins or to unanchored ubiquitin chains. These studies define additional components of a UBC2-dependent ubiquitination pathway shown previously to be essential for promastigote to amastigote differentiation.

Refinement of Leishmania donovani Genome Annotations in the Light of Ribosome-Protected mRNAs Fragments (Ribo-Seq Data)

Advances in next-generation sequencing methodologies have facilitated the assembly of an ever-increasing number of genomes. Gene annotations are typically conducted via specialized software, but the most accurate results require additional manual curation that incorporates insights derived from functional and bioinformatic analyses (e.g., transcriptomics, proteomics, and phylogenetics). In this study, we improved the annotation of the Leishmania donovani (strain HU3) genome using publicly available data from the deep sequencing of ribosome-protected mRNA fragments (Ribo-Seq). As a result of this analysis, we uncovered 70 previously non-annotated protein-coding genes and improved the annotation of around 600 genes. Additionally, we present evidence for small upstream open reading frames (uORFs) in a significant number of transcripts, indicating their potential role in the translational regulation of gene expression. The bioinformatics pipelines developed for these analyses can be used to improve the genome annotations of other organisms for which Ribo-Seq data are available. The improvements provided by these studies will bring us closer to the ultimate goal of a complete and accurately annotated L. donovani genome and will enhance future transcriptomics, proteomics, and genetics studies.

Kinetoplast Genome of Leishmania spp. Is under Strong Purifying Selection

Instability is an intriguing characteristic of many protist genomes, and trypanosomatids are not an exception in this respect. Some regions of trypanosomatid genomes evolve fast. For instance, the trypanosomatid mitochondrial (kinetoplast) genome consists of fairly conserved maxicircle and minicircle molecules that can, nevertheless, possess high nucleotide substitution rates between closely related strains. Recent experiments have demonstrated that rapid laboratory evolution can result in the non-functionality of multiple genes of kinetoplast genomes due to the accumulation of mutations or loss of critical genomic components. An example of a loss of critical components is the reported loss of entire minicircle classes in Leishmania tarentolae during laboratory cultivation, which results in an inability to generate some correctly encoded genes. In the current work, we estimated the evolutionary rates of mitochondrial and nuclear genome regions of multiple natural Leishmania spp. We analyzed synonymous and non-synonymous substitutions and, rather unexpectedly, found that the coding regions of kinetoplast maxicircles are among the most variable regions of both genomes. In addition, we demonstrate that synonymous substitutions greatly predominate among maxicircle coding regions and that most maxicircle genes show signs of purifying selection. These results imply that maxicircles in natural Leishmania populations remain functional despite their high mutation rate.

ULK4 and Fused/STK36 interact to mediate assembly of a motile flagellum

Unc-51-like kinase (ULK) family serine-threonine protein kinase homologues have been linked to the function of motile cilia in diverse species. Mutations in Fused/STK36 and ULK4 in mice resulted in hydrocephalus and other phenotypes consistent with ciliary defects. How either protein contributes to the assembly and function of motile cilia is not well understood. Here we studied the phenotypes of ULK4 and Fused gene knockout (KO) mutants in the flagellated protist Leishmania mexicana. Both KO mutants exhibited a variety of structural defects of the flagellum cytoskeleton. Biochemical approaches indicate spatial proximity of these proteins and indicate a direct interaction between the N-terminus of LmxULK4 and LmxFused. Both proteins display a dispersed localization throughout the cell body and flagellum, with enrichment near the flagellar base and tip. The stable expression of LmxULK4 was dependent on the presence of LmxFused. Fused/STK36 was previously shown to localize to mammalian motile cilia, and we demonstrate here that ULK4 also localizes to the motile cilia in mouse ependymal cells. Taken together these data suggest a model where the pseudokinase ULK4 is a positive regulator of the kinase Fused/ STK36 in a pathway required for stable assembly of motile cilia.

Ribosome Specialization in Protozoa Parasites

Ribosomes, in general, are viewed as constitutive macromolecular machines where protein synthesis takes place; however, this view has been recently challenged, supporting the hypothesis of ribosome specialization and opening a completely new field of research. Recent studies have demonstrated that ribosomes are heterogenous in their nature and can provide another layer of gene expression control by regulating translation. Heterogeneities in ribosomal RNA and ribosomal proteins that compose them favor the selective translation of different sub-pools of mRNAs and functional specialization. In recent years, the heterogeneity and specialization of ribosomes have been widely reported in different eukaryotic study models; however, few reports on this topic have been made on protozoa and even less on protozoa parasites of medical importance. This review analyzes heterogeneities of ribosomes in protozoa parasites highlighting the specialization in their functions and their importance in parasitism, in the transition between stages in their life cycle, in the change of host and in response to environmental conditions.

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.

Cell-to-flagellum attachment and surface architecture in kinetoplastids

A key morphological feature of kinetoplastid parasites is the position and length of flagellum attachment to the cell body. This lateral attachment is mediated by the flagellum attachment zone (FAZ), a large complex cytoskeletal structure, which is essential for parasite morphogenesis and pathogenicity. Despite the complexity of the FAZ only two transmembrane proteins, FLA1 and FLA1BP, are known to interact and connect the flagellum to the cell body. Across the different kinetoplastid species, each only has a single FLA/FLABP pair, except in Trypanosoma brucei and Trypanosoma congolense where there has been an expansion of these genes. Here, we focus on the selection pressure behind the evolution of the FLA/FLABP proteins and the likely impact this will have on host-parasite interactions.

Genome sequence of Leishmania mexicana MNYC/BZ/62/M379 expressing Cas9 and T7 RNA polymerase

We present the genome sequence of Leishmania mexicana MNYC/BZ/62/M379 modified to express Cas9 and T7 RNA-polymerase, revealing high similarity to the reference genome (MHOM/GT2001/U1103). Through RNAseq-based annotation of coding sequences and untranslated regions, we provide primer sequences for construct and sgRNA template generation for CRISPR-assisted gene deletion and endogenous tagging.

Genome sequence of Leishmania mexicana MNYC/BZ/62/M379 expressing Cas9 and T7 RNA polymerase

We present the genome sequence of Leishmania mexicana MNYC/BZ/62/M379 modified to express Cas9 and T7 RNA-polymerase, revealing high similarity to the reference genome (MHOM/GT2001/U1103). Through RNAseq-based annotation of coding sequences and untranslated regions, we provide primer sequences for construct and sgRNA template generation for CRISPR-assisted gene deletion and endogenous tagging.

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.

Discovery of Leishmania Druggable Serine Proteases by Activity-Based Protein Profiling

Leishmaniasis are a group of diseases caused by parasitic protozoa of the genus Leishmania. Current treatments are limited by difficult administration, high cost, poor efficacy, toxicity, and growing resistance. New agents, with new mechanisms of action, are urgently needed to treat the disease. Although extensively studied in other organisms, serine proteases (SPs) have not been widely explored as antileishmanial drug targets. Herein, we report for the first time an activity-based protein profiling (ABPP) strategy to investigate new therapeutic targets within the SPs of the Leishmania parasites. Active-site directed fluorophosphonate probes (rhodamine and biotin-conjugated) were used for the detection and identification of active Leishmania serine hydrolases (SHs). Significant differences were observed in the SHs expression levels throughout the Leishmania life cycle and between different Leishmania species. Using iTRAQ-labelling-based quantitative proteomic mass spectrometry, we identified two targetable SPs in Leishmania mexicana: carboxypeptidase LmxM.18.0450 and prolyl oligopeptidase LmxM.36.6750. Druggability was ascertained by selective inhibition using the commercial serine protease inhibitors chymostatin, lactacystin and ZPP, which represent templates for future anti-leishmanial drug discovery programs. Collectively, the use of ABPP method complements existing genetic methods for target identification and validation in Leishmania.

Bromodomain factor 5 is an essential regulator of transcription in Leishmania

Leishmania are unicellular parasites that cause human and animal diseases. Like other kinetoplastids, they possess large transcriptional start regions (TSRs) which are defined by histone variants and histone lysine acetylation. Cellular interpretation of these chromatin marks is not well understood. Eight bromodomain factors, the reader modules for acetyl-lysine, are found across Leishmania genomes. Using L. mexicana, Cas9-driven gene deletions indicate that BDF1–5 are essential for promastigotes. Dimerisable, split Cre recombinase (DiCre)-inducible gene deletion of BDF5 show it is essential for both promastigotes and murine infection. ChIP-seq identifies BDF5 as enriched at TSRs. XL-BioID proximity proteomics shows the BDF5 landscape is enriched for BDFs, HAT2, proteins involved in transcriptional activity, and RNA processing; revealing a Conserved Regulators of Kinetoplastid Transcription (CRKT) Complex. Inducible deletion of BDF5 causes global reduction in RNA polymerase II transcription. Our results indicate the requirement of Leishmania to interpret histone acetylation marks through the bromodomain-enriched CRKT complex for normal gene expression and cellular viability.

Leishmania use large (5–10 kb) transcriptional start regions, where the chromatin is highly enriched for acetylated histones, to drive the expression of polycistronic gene arrays. Here the authors show bromodomain-containing protein BDF5 is enriched at transcriptional start sites and its depletion leads to cell death in vitro and in murine infections, and they identify its interactors.

In vivo transcriptional analysis of mice infected with Leishmania major unveils cellular heterogeneity and altered transcriptomic profiling at single-cell resolution

Leishmania parasites cause cutaneous leishmaniasis (CL), a disease characterized by disfiguring, ulcerative skin lesions. Both parasite and host gene expression following infection with various Leishmania species has been investigated in vitro, but global transcriptional analysis following L. major infection in vivo is lacking. Thus, we conducted a comprehensive transcriptomic profiling study combining bulk RNA sequencing (RNA-Seq) and single-cell RNA sequencing (scRNA-Seq) to identify global changes in gene expression in vivo following L. major infection. Bulk RNA-Seq analysis revealed that host immune response pathways like the antigen processing and presentation pathway were significantly enriched amongst differentially expressed genes (DEGs) upon infection, while ribosomal pathways were significantly downregulated in infected mice compared to naive controls. scRNA-Seq analyses revealed cellular heterogeneity including distinct resident and recruited cell types in the skin following murine L. major infection. Within the individual immune cell types, several DEGs indicative of many interferon induced GTPases and antigen presentation molecules were significantly enhanced in the infected ears including macrophages, resident macrophages, and inflammatory monocytes. Ingenuity Pathway Analysis of scRNA-Seq data indicated the antigen presentation pathway was increased with infection, while EIF2 signaling is the top downregulated pathway followed by eIF4/p70S6k and mTOR signaling in multiple cell types including macrophages, blood and lymphatic endothelial cells. Altogether, this transcriptomic profile highlights known recruitment of myeloid cells to lesions and recognizes a potential role for EIF2 signaling in murine L. major infection in vivo.

Transcriptome Analysis of Intracellular Amastigotes of Clinical Leishmania infantum Lines from Therapeutic Failure Patients after Infection of Human Macrophages

Leishmaniasis is considered to be one of the most neglected tropical diseases affecting humans and animals around the world. Due to the absence of an effective vaccine, current treatment is based on chemotherapy. However, the continuous appearance of drug resistance and therapeutic failure (TF) lead to an early obsolescence of treatments. Identification of the factors that contribute to TF and drug resistance in leishmaniasis will constitute a useful tool for establishing future strategies to control this disease. In this manuscript, we evaluated the transcriptomic changes in the intracellular amastigotes of the Leishmania infantum parasites isolated from patients with leishmaniasis and TF at 96 h post-infection of THP-1 cells. The adaptation of the parasites to their new environment leads to expression alterations in the genes involved mainly in the transport through cell membranes, energy and redox metabolism, and detoxification. Specifically, the gene that codes for the prostaglandin f2α synthase seems to be relevant in the pathogenicity and TF since it appears substantially upregulated in all the L. infantum lines. Overall, our results show that at the late infection timepoint, the transcriptome of the parasites undergoes significant changes that probably improve the survival of the Leishmania lines in the host cells, contributing to the TF phenotype as well as drug therapy evasion.

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.

Novel Insights Into Leishmania (Viannia) braziliensis In Vitro Fitness Guided by Temperature Changes Along With Its Subtilisins and Oligopeptidase B

Proteases are virulence factors with a recognized impact on the Leishmania spp. life cycle. This study considers a set of analyses measuring phenotypic factors of L. (V.) braziliensis clinical isolates as promastigotes growth curves, murine peritoneal macrophages infection, inflammatory mediators production, and serine proteases gene expression (subtilisin 13: S13, subtilisin 28: S28, oligopeptidase B: OPB) assessing these isolates’ fitness on in vitro conditions. Parasites had different behavior during the early growth phase from day zero to day three, and all isolates reached the stationary growth phase between days four and seven. Macrophages infection showed two tendencies, one of decreased infection rate and number of parasites per macrophage (Infection Index <1000) and another with a constant infection index (≥1400). TNF-α (≥10 pg/mL) detected in infections by 75% of isolates, IL-6 (≥80 pg/mL) by 30% of isolates and low levels of NO (≥0.01µM) in almost all infections. Gene expression showed higher values of S13 (≥2RQ) in the intracellular amastigotes of all the isolates evaluated. On the contrary, S28 expression was low (≤1RQ) in all isolates. OPB expression was different between promastigotes and intracellular amastigotes, being significantly higher (≥2RQ) in the latter form of 58% of the isolates. Predictive structural assays of S13 and OPB were performed to explore temperature influence on gene expression and the encoded proteases. Gene expression data is discussed based on in silico predictions of regulatory regions that show plasticity in the linearity index of secondary structures of S13 and OPB 3’-untranslated regions of mRNA, dependent on temperature changes. While hairpin structures suggest an active region of mRNA for both genes above 26°C, pseudoknot structure found in S13 is an indication of a particular profile of this gene at mammalian host temperatures (37°C). Furthermore, the predicted 3D structures are in accordance with the influence of these temperatures on the catalytic site stability of both enzymes, favoring their action over peptide substrates. Data gathered here suggest that L. (V.) braziliensis serine proteases can be influenced by the temperature conditions affecting parasite fitness throughout its life cycle.

Absence of DEATH kinesin is fatal for Leishmania mexicana amastigotes

Kinesins are motor proteins present in organisms from protists to mammals playing important roles in cell division, intracellular organisation and flagellum formation and maintenance. Leishmania mexicana is a protozoan parasite of the order Kinetoplastida causing human cutaneous leishmaniasis. Kinetoplastida genome sequence analyses revealed a large number of kinesins showing sequence and structure homology to eukaryotic kinesins. Here, we investigate the L. mexicana kinesin LmxKIN29 (LmxM.29.0350), also called DEATH kinesin. The activated MAP kinase LmxMPK3, a kinase affecting flagellum length in Leishmania, is able to phosphorylate recombinant full length LmxKIN29 at serine 554. Insect promastigote LmxKIN29 Leishmania null mutants showed no obvious phenotype. However, in mouse infection experiments, the null mutants were unable to cause the disease, whereas LmxKIN29 add-backs and single allele knockouts caused footpad lesions. Localisation using promastigotes expressing GFP-tagged LmxKIN29 revealed that the kinesin is predominantly found in between the nucleus and the flagellar pocket, while in dividing cells the GFP-fusion protein was found at the anterior and posterior ends of the cells indicating a role in cytokinesis. The inability to cause lesions in infected animals and the amino acid sequence divergence from mammalian kinesins suggests that LmxKIN29 is a potential drug target against leishmaniasis.

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.

Application of single-cell transcriptomics to kinetoplastid research

Kinetoplastid parasites are responsible for both human and animal diseases across the globe where they have a great impact on health and economic well-being. Many species and life cycle stages are difficult to study due to limitations in isolation and culture, as well as to their existence as heterogeneous populations in hosts and vectors. Single-cell transcriptomics (scRNA-seq) has the capacity to overcome many of these difficulties, and can be leveraged to disentangle heterogeneous populations, highlight genes crucial for propagation through the life cycle, and enable detailed analysis of host–parasite interactions. Here, we provide a review of studies that have applied scRNA-seq to protozoan parasites so far. In addition, we provide an overview of sample preparation and technology choice considerations when planning scRNA-seq experiments, as well as challenges faced when analysing the large amounts of data generated. Finally, we highlight areas of kinetoplastid research that could benefit from scRNA-seq technologies.

Advances in Understanding Leishmania Pathobiology: What Does RNA-Seq Tell Us?

Leishmaniasis is a vector-borne disease caused by a protozoa parasite from over 20 Leishmania species. The clinical manifestations and the outcome of the disease vary greatly. Global RNA sequencing (RNA-Seq) analyses emerged as a powerful technique to profile the changes in the transcriptome that occur in the Leishmania parasites and their infected host cells as the parasites progresses through their life cycle. Following the bite of a sandfly vector, Leishmania are transmitted to a mammalian host where neutrophils and macrophages are key cells mediating the interactions with the parasites and result in either the elimination the infection or contributing to its proliferation. This review focuses on RNA-Seq based transcriptomics analyses and summarizes the main findings derived from this technology. In doing so, we will highlight caveats in our understanding of the parasite's pathobiology and suggest novel directions for research, including integrating more recent data highlighting the role of the bacterial members of the sandfly gut microbiota and the mammalian host skin microbiota in their potential role in influencing the quantitative and qualitative aspects of leishmaniasis pathology.

Gene Annotation and Transcriptome Delineation on a De Novo Genome Assembly for the Reference Leishmania major Friedlin Strain

Leishmania major is the main causative agent of cutaneous leishmaniasis in humans. The Friedlin strain of this species (LmjF) was chosen when a multi-laboratory consortium undertook the objective of deciphering the first genome sequence for a parasite of the genus Leishmania. The objective was successfully attained in 2005, and this represented a milestone for Leishmania molecular biology studies around the world. Although the LmjF genome sequence was done following a shotgun strategy and using classical Sanger sequencing, the results were excellent, and this genome assembly served as the reference for subsequent genome assemblies in other Leishmania species. Here, we present a new assembly for the genome of this strain (named LMJFC for clarity), generated by the combination of two high throughput sequencing platforms, Illumina short-read sequencing and PacBio Single Molecular Real-Time (SMRT) sequencing, which provides long-read sequences. Apart from resolving uncertain nucleotide positions, several genomic regions were reorganized and a more precise composition of tandemly repeated gene loci was attained. Additionally, the genome annotation was improved by adding 542 genes and more accurate coding-sequences defined for around two hundred genes, based on the transcriptome delimitation also carried out in this work. As a result, we are providing gene models (including untranslated regions and introns) for 11,238 genes. Genomic information ultimately determines the biology of every organism; therefore, our understanding of molecular mechanisms will depend on the availability of precise genome sequences and accurate gene annotations. In this regard, this work is providing an improved genome sequence and updated transcriptome annotations for the reference L. major Friedlin strain.

Axenic amastigotes of Leishmania species as a suitable model for in vitro studies

Leishmania spp. are etiological agents of infection diseases, which in some cases can be fatal. The main forms of their biological cycle, promastigotes and amastigotes, can be maintained in vitro. While promastigotes are easier to maintain, amastigotes are more complex and can be obtained through different ways, including infection assays of tissues or in vitro cells, and differentiation from promastigotes to axenic amastigotes. Several protocols have been proposed for in vitro differentiation for at least 12 Leishmania spp. of both subgenera, Leishmania and Viannia. In this review we propose a critical summary of axenic amastigotes induction, as well as the impact of these strategies on metabolic pathways and regulatory networks analyzed by omics approaches. The parameters used by different research groups show considerable variations in temperature, pH and induction stages, as highlighted here for Leishmania (Viannia) braziliensis. Therefore, a consensus on strategies for inducing amastigogenesis is necessary to improve accuracy and even define stage-specific biomarkers. In fact, the axenic amastigote model has contributed to elucidate several aspects of the parasite cycle, however, since it does not reproduce the intracellular environment, its use requires several precautions. In addition, we present a discussion about using axenic amastigotes for drug screening, suggesting the need of a more sensitive methodology to verify cell viability in these tests. Collectively, this review explores the advantages and limitations found in studies with axenic amastigotes, done for more than 30 years, and discuss the gaps that impair their use as a suitable model for in vitro studies.

Protein Serine/Threonine Phosphatase Type 2C of Leishmania mexicana

Protein phosphorylation and dephosphorylation are increasingly recognized as important processes for regulating multiple physiological mechanisms. Phosphorylation is carried out by protein kinases and dephosphorylation by protein phosphatases. Phosphoprotein phosphatases (PPPs), one of three families of protein serine/threonine phosphatases, have great structural diversity and are involved in regulating many cell functions. PP2C, a type of PPP, is found in Leishmania, a dimorphic protozoan parasite and the causal agent of leishmaniasis. The aim of this study was to clone, purify, biochemically characterize and quantify the expression of PP2C in Leishmania mexicana (LmxPP2C). Recombinant LmxPP2C dephosphorylated a specific threonine (with optimal activity at pH 8) in the presence of the manganese divalent cation (Mn⁺²). LmxPP2C activity was inhibited by sanguinarine (a specific inhibitor) but was unaffected by protein tyrosine phosphatase inhibitors. Western blot analysis indicated that anti-LmxPP2C antibodies recognized a molecule of 45.2 kDa. Transmission electron microscopy with immunodetection localized LmxPP2C in the flagellar pocket and flagellum of promastigotes but showed poor staining in amastigotes. Interestingly, LmxPP2C belongs to the ortholog group OG6_142542, which contains only protozoa of the family Trypanosomatidae. This suggests a specific function of the enzyme in the flagellar pocket of these microorganisms.

Genomics and functional genomics in Leishmania and Trypanosoma cruzi: statuses, challenges and perspectives

The availability of Trypanosomatid genomic data in public databases has opened myriad experimental possibilities that have contributed to a more comprehensive understanding of the biology of these parasites and their interactions with hosts. In this review, after brief remarks on the history of the Trypanosoma cruzi and Leishmania genome initiatives, we present an overview of the relevant contributions of genomics, transcriptomics and functional genomics, discussing the primary obstacles, challenges, relevant achievements and future perspectives of these technologies.

VSG mRNA levels are regulated by the production of functional VSG protein

The bloodstream form of Trypanosoma brucei persists in mammalian hosts through a population survival strategy depending on antigenic variation of a cell surface coat composed of the variant surface glycoprotein (VSG). The integrity of the VSG coat is essential and blocking its synthesis results in a cell division cycle arrest just prior to cytokinesis. This observation indicates that VSG levels are monitored and that the cell has mechanisms to respond to a disruption of synthesis. Here, the regulation of VSG mRNA levels has been investigated by first measuring VSG mRNA copy number, and second using ectopic expression of VSG transgenes containing premature termination codons. The findings are that (i) VSG mRNA copy number varies with the identity of the VSG and (ii) a pathway detects synthesis of non-functional VSG protein and results in an increase in VSG mRNA levels.

Leishmania differentiation requires ubiquitin conjugation mediated by a UBC2-UEV1 E2 complex

Post-translational modifications such as ubiquitination are important for orchestrating the cellular transformations that occur as the Leishmania parasite differentiates between its main morphological forms, the promastigote and amastigote. 2 E1 ubiquitin-activating (E1), 13 E2 ubiquitin-conjugating (E2), 79 E3 ubiquitin ligase (E3) and 20 deubiquitinating cysteine peptidase (DUB) genes can be identified in the Leishmania mexicana genome but, currently, little is known about the role of E1, E2 and E3 enzymes in this parasite. Bar-seq analysis of 23 E1, E2 and HECT/RBR E3 null mutants generated in promastigotes using CRISPR-Cas9 revealed numerous loss-of-fitness phenotypes in promastigote to amastigote differentiation and mammalian infection. The E2s UBC1/CDC34, UBC2 and UEV1 and the HECT E3 ligase HECT2 are required for the successful transformation from promastigote to amastigote and UBA1b, UBC9, UBC14, HECT7 and HECT11 are required for normal proliferation during mouse infection. Of all ubiquitination enzyme null mutants examined in the screen, Δubc2 and Δuev1 exhibited the most extreme loss-of-fitness during differentiation. Null mutants could not be generated for the E1 UBA1a or the E2s UBC3, UBC7, UBC12 and UBC13, suggesting these genes are essential in promastigotes. X-ray crystal structure analysis of UBC2 and UEV1, orthologues of human UBE2N and UBE2V1/UBE2V2 respectively, reveal a heterodimer with a highly conserved structure and interface. Furthermore, recombinant L. mexicana UBA1a can load ubiquitin onto UBC2, allowing UBC2-UEV1 to form K63-linked di-ubiquitin chains in vitro. Notably, UBC2 can cooperate in vitro with human E3s RNF8 and BIRC2 to form non-K63-linked polyubiquitin chains, showing that UBC2 can facilitate ubiquitination independent of UEV1, but association of UBC2 with UEV1 inhibits this ability. Our study demonstrates the dual essentiality of UBC2 and UEV1 in the differentiation and intracellular survival of L. mexicana and shows that the interaction between these two proteins is crucial for regulation of their ubiquitination activity and function.

The post-translational modification of proteins is key for allowing Leishmania parasites to transition between the different life cycle stages that exist in its insect vector and mammalian host. In particular, components of the ubiquitin system are important for the transformation of Leishmania from its insect (promastigote) to mammalian (amastigote) stage and normal infection in mice. However, little is known about the role of the enzymes that generate ubiquitin modifications in Leishmania. Here we characterise 28 enzymes of the ubiquitination pathway and show that many are required for life cycle progression or mouse infection by this parasite. Two proteins, UBC2 and UEV1, were selected for further study based on their importance in the promastigote to amastigote transition. We demonstrate that UBC2 and UEV1 form a heterodimer capable of carrying out ubiquitination and that the structural basis for this activity is conserved between Leishmania, Saccharomyces cerevisiae and humans. We also show that the interaction of UBC2 with UEV1 alters the nature of the ubiquitination activity performed by UBC2. Overall, we demonstrate the important role that ubiquitination enzymes play in the life cycle and infection process of Leishmania and explore the biochemistry underlying UBC2 and UEV1 function.

Dual transcriptome analysis reveals differential gene expression modulation influenced by Leishmania arginase and host genetic background

The outcome of Leishmania infection is strongly influenced by the host's genetic background. BALB/c mice are susceptible to Leishmania infection, while C57BL/6 mice show discrete resistance. Central to the fate of the infection is the availability of l-arginine and the related metabolic processes in the host and parasite. Depending on l-arginine availability, nitric oxide synthase 2 (NOS2) of the host cell produces nitric oxide (NO) controlling the parasite growth. On the other hand, Leishmania can also use host l-arginine for the production of polyamines through its own arginase activity, thus favouring parasite replication. Considering RNA-seq data, we analysed the dual modulation of host and parasite gene expression of BALB/c or C57BL/6 mouse bone marrow-derived macrophages (BMDMs) after 4 h of infection with Leishmania amazonensis wild-type (La-WT) or L. amazonensis arginase knockout (La-arg-). We identified 12 641 host transcripts and 8282 parasite transcripts by alignment analysis with the respective Mus musculus and L. mexicana genomes. The comparison of BALB/c_La-arg-versus BALB/c_La-WT revealed 233 modulated transcripts, with most related to the immune response and some related to the amino acid transporters and l-arginine metabolism. In contrast, the comparison of C57BL/6_La-arg-vs. C57BL/6_La-WT revealed only 30 modulated transcripts, including some related to the immune response but none related to amino acid transport or l-arginine metabolism. The transcriptome profiles of the intracellular amastigote revealed 94 modulated transcripts in the comparison of La-arg-_BALB/c vs. La-WT_BALB/c and 45 modulated transcripts in the comparison of La-arg-_C57BL/6 vs. La-WT_C57BL/6. Taken together, our data present new insights into the impact of parasite arginase activity on the orchestration of the host gene expression modulation, including in the immune response and amino acid transport and metabolism, mainly in susceptible BALB/c-infected macrophages. Moreover, we show how parasite arginase activity affects parasite gene expression modulation, including amino acid uptake and amastin expression.

Identification of gene expression profiles in Leishmania major infection by integrated bioinformatics analyses

Gene expression profiling in mouse models of leishmaniasis has given useful information to understand the molecular pathways active in lesions and to discover new diagnostic/therapeutic targets. Although the host response plays a critical role in protection from leishmaniasis and promoting disease severity, there are still unexplained aspects in the mechanism of non-healing cutaneous lesions, which need biomarkers for both targeted- therapy and diagnosis. To address this, transcriptional profiling of the skin lesions obtained from BALB/c mice infected with Leishmania major and healthy skin from naïve mice were evaluated by bioinformatics analysis, and then the results were validated by Revers Transcriptase-PCR. Five genes among the up-regulated differentially expressed genes named FCGR4, CCL4, CXCL9, Arg1 and IL-1β were found to have relatively high diagnostic value for CL due to L. major. Pathway analysis revealed that Triggering Receptor Expressed On Myeloid Cells 1 (TREM1) signaling pathways are active in cutaneous lesions, providing new insights for the understanding and treatment of leishmaniasis.

Frequent Recombination Events in Leishmania donovani: Mining Population Data

The Leishmania donovani species complex consists of all L. donovani and L. infantum strains mainly responsible for visceral leishmaniasis (VL). It was suggested that genome rearrangements in Leishmania spp. occur very often, thus enabling parasites to adapt to the different environmental conditions. Some of these rearrangements may be directly linked to the virulence or explain the reduced efficacy of antimonial drugs in some isolates. In the current study, we focused on a large-scale analysis of putative gene conversion events using publicly available datasets. Previous population study of L. donovani suggested that population variability of L. donovani is relatively low, however the authors used masking procedures and strict read selection criteria. We decided to re-analyze DNA-seq data without masking sequences, because we were interested in the most dynamic fraction of the genome. The majority of samples have an excess of putative gene conversion/recombination events in the noncoding regions, however we found an overall excess of putative intrachromosomal gene conversion/recombination in the protein coding genes, compared to putative interchromosomal gene conversion/recombination events.

Bar-seq strategies for the LeishGEdit toolbox

The number of fully sequenced genomes increases steadily but the function of many genes remains unstudied. To accelerate dissection of gene function in Leishmania spp. and other kinetoplastids we previously developed a streamlined pipeline for CRISPR-Cas9 gene editing, which we termed LeishGEdit. To facilitate high-throughput mutant screens we have adapted this pipeline by barcoding mutants with unique 17-nucleotide barcodes, allowing loss-of-function screens in mixed populations. Here we present primer design and analysis tools that facilitate these bar-seq strategies. We have developed a standalone easy-to-use pipeline to design CRISPR primers suitable for the LeishGEdit toolbox for any given genome and have generated a list of 14,995 barcodes. Barcodes and oligo sequences are now accessible through our website www.leishgedit.net allowing researchers to pursue bar-seq experiments in all currently available TriTrypDB genomes (release 41). This will streamline CRISPR bar-seq assays in kinetoplastids, enabling pooled mutant screens across the community.

Distinct gene expression patterns in vector-residing Leishmania infantum identify parasite stage-enriched markers

BACKGROUND

Leishmaniasis is a vector-borne neglected disease. Inside the natural sand fly vector, the promastigote forms of Leishmania undergo a series of extracellular developmental stages to reach the infectious stage, the metacyclic promastigote. There is limited information regarding the expression profile of L. infantum developmental stages inside the sand fly vector, and molecular markers that can distinguish the different parasite stages are lacking.

METHODOLOGY/PRINCIPAL FINDINGS

We performed RNAseq on unaltered midguts of the sand fly Lutzomyia longipalpis after infection with L. infantum parasites. RNAseq was carried out at various time points throughout parasite development. Principal component analysis separated the transcripts corresponding to the different Leishmania promastigote stages, the procyclic, nectomonad, leptomonad and metacyclics. Importantly, there were a significant number of differentially expressed genes when comparing the sequential development of the various Leishmania stages in the sand fly. There were 836 differentially expressed (DE) genes between procyclic and long nectomonad promastigotes; 113 DE genes between nectomonad and leptomonad promastigotes; and 302 DE genes between leptomonad and metacyclic promastigotes. Most of the DE genes do not overlap across stages, highlighting the uniqueness of each Leishmania stage. Furthermore, the different stages of Leishmania parasites exhibited specific transcriptional enrichment across chromosomes. Using the transcriptional signatures exhibited by distinct Leishmania stages during their development in the sand fly midgut, we determined the genes predominantly enriched in each stage, identifying multiple potential stage-specific markers for L. infantum.

CONCLUSIONS

Overall, these findings demonstrate the transcriptional plasticity of the Leishmania parasite inside the sand fly vector and provide a repertoire of potential stage-specific markers for further development as molecular tools for epidemiological studies.

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.

LAX28 is required for the stable assembly of the inner dynein arm f complex, and the tether and tether head complex in Leishmania flagella

Motile eukaryotic flagella beat through coordinated activity of dynein motor proteins; however, the mechanisms of dynein coordination and regulation are incompletely understood. The inner dynein arm (IDA) f complex (also known as the I1 complex), and the tether and tether head (T/TH) complex are thought to be key regulators of dynein action but, unlike the IDA f complex, T/TH proteins remain poorly characterised. Here, we characterised T/TH-associated proteins in the protist Leishmania mexicana Proteome analysis of axonemes from null mutants for the CFAP44 T/TH protein showed that they lacked the IDA f protein IC140 and a novel 28-kDa axonemal protein, LAX28. Sequence analysis identified similarities between LAX28 and the uncharacterised human sperm tail protein TEX47, both sharing features with sensory BLUF-domain-containing proteins. Leishmania lacking LAX28, CFAP44 or IC140 retained some motility, albeit with reduced swimming speed and directionality and a propensity for flagellar curling. Expression of tagged proteins in different null mutant backgrounds showed that the axonemal localisation of LAX28 requires CFAP44 and IC140, and the axonemal localisations of CFAP44 and IC140 both depend on LAX28. These data demonstrate a role for LAX28 in motility and show mutual dependencies of IDA f and T/TH-associated proteins for axonemal assembly in Leishmania.

Analysis of the Physiological and Metabolic State of Leishmania Using Heavy Water Labeling

This protocol describes the use of heavy water (²H₂O) labeling to determine the growth rate and metabolic state of Leishmania parasites in culture and in infected animals. In vitro labeling studies are undertaken by cultivating defined parasite developmental stages in standard medium supplemented with 5% ²H₂O, resulting in the incorporation of deuterium (²H) into a range of metabolic precursors used in macromolecule (DNA, RNA, protein, lipid, and glycan) synthesis. The rate of turnover of different parasite macromolecules can subsequently be determined by analysis of deuterium enrichment in the different constituents of these macromolecules by gas chromatography-mass spectrometry (GC-MS). To measure the growth rate and physiological state of parasite stages in lesion tissue, infected mice were provided with 9% ²H₂O in their drinking water for various periods of time and ²H-enrichment in the macromolecular constituents of isolated lesion-derived parasite stages determined by GC-MS. This protocol provides quantitative information on key cellular processes, such as replication (DNA turnover), transcription (RNA turnover), translation (protein turnover), membrane biogenesis (lipid turnover), and central carbon metabolism (glycan turnover) that define the growth state and phenome of different parasite stages in vitro and in vivo. This approach can be used to assess the impact of host immune responses on parasite growth and physiology (using different Leishmania strains/species, mouse lines), characterize different parasite populations during chronic and acute infections, and assess parasite responses to drug treatments. It is also broadly applicable to other microbial pathogens.

Differential immune response modulation in early Leishmania amazonensis infection of BALB/c and C57BL/6 macrophages based on transcriptome profiles

The fate of Leishmania infection can be strongly influenced by the host genetic background. In this work, we describe gene expression modulation of the immune system based on dual global transcriptome profiles of bone marrow-derived macrophages (BMDMs) from BALB/c and C57BL/6 mice infected with Leishmania amazonensis. A total of 12,641 host transcripts were identified according to the alignment to the Mus musculus genome. Differentially expressed genes (DEGs) profiling revealed a differential modulation of the basal genetic background between the two hosts independent of L. amazonensis infection. In addition, in response to early L. amazonensis infection, 10 genes were modulated in infected BALB/c vs. non-infected BALB/c macrophages; and 127 genes were modulated in infected C57BL/6 vs. non-infected C57BL/6 macrophages. These modulated genes appeared to be related to the main immune response processes, such as recognition, antigen presentation, costimulation and proliferation. The distinct gene expression was correlated with the susceptibility and resistance to infection of each host. Furthermore, upon comparing the DEGs in BMDMs vs. peritoneal macrophages, we observed no differences in the gene expression patterns of Jun, Fcgr1 and Il1b, suggesting a similar activation trends of transcription factor binding, recognition and phagocytosis, as well as the proinflammatory cytokine production in response to early L. amazonensis infection. Analysis of the DEG profile of the parasite revealed only one DEG among the 8,282 transcripts, indicating that parasite gene expression in early infection does not depend on the host genetic background.

Mapping benznidazole resistance in trypanosomatids and exploring evolutionary histories of nitroreductases and ABCG transporter protein sequences

The nitro-heterocyclic compound benznidazole (BZ) is the first-line drug for the treatment of Chagas disease, caused by the protozoan Trypanosoma cruzi. However, therapeutic failures are common for reasons that include the influences of parasite and host genetics, the effects of toxicity on adherence to treatment, and difficulties in demonstrating parasitological cure. To obtain information on the origin of the resistance to BZ and eliminate from the scenery the participation of the host, initially we mapped the susceptibility to the drug in thirteen species of seven genera of the family Trypanosomatidae. We verified that all Trypanosoma species are sensitive to low concentrations of the drug (IC₅₀ 2.7 to 25 µM) while Non-Trypanosoma species are highly resistant to these concentrations. The two groups of parasites correspond to the major phylogenetic lineages of trypanosomatids. Next, we searched in the trypanosomatid genome databases homologs of two type-I nitroreductases (NTR-1 and OYE) and an ABC transporter (ABCG1) that have been associated with BZ resistance in T. cruzi. The predicted proteins were characterized regarding domains and used for phylogenetic analyses. Homologous NTR-1 genes were found in all trypanosomatids investigated and the structural characteristics of the enzyme suggest that it may be functional. OYE genes were absent in BZ-sensitive African trypanosomes, which excludes the participation of this enzyme in BZ bio-activation. Two copies of ABCG1 genes were observed in most BZ resistant species, while Trypanosoma species exhibit only one copy per haploid genome. Functional studies are required to verify the involvement of these genes in BZ resistance. In addition, since multiple mechanisms can contribute to BZ susceptibility, our study poses a range of organisms highly resistant to BZ in which these aspects can be investigated. Preliminary studies on BZ uptake indicate marked differences between BZ-sensitive and BZ-resistant species.

P-type transport ATPases in Leishmania and Trypanosoma

P-type ATPases are critical to the maintenance and regulation of cellular ion homeostasis and membrane lipid asymmetry due to their ability to move ions and phospholipids against a concentration gradient by utilizing the energy of ATP hydrolysis. P-type ATPases are particularly relevant in human pathogenic trypanosomatids which are exposed to abrupt and dramatic changes in their external environment during their life cycles. This review describes the complete inventory of ion-motive, P-type ATPase genes in the human pathogenic Trypanosomatidae; eight Leishmania species (L. aethiopica, L. braziliensis, L. donovani, L. infantum, L. major, L. mexicana, L. panamensis, L. tropica), Trypanosoma cruzi and three Trypanosoma brucei subspecies (Trypanosoma brucei brucei TREU927, Trypanosoma brucei Lister strain 427, Trypanosoma brucei gambiense DAL972). The P-type ATPase complement in these trypanosomatids includes the P_1B (metal pumps), P_2A (SERCA, sarcoplasmic-endoplasmic reticulum calcium ATPases), P_2B (PMCA, plasma membrane calcium ATPases), P_2D (Na⁺ pumps), P_3A (H⁺ pumps), P₄ (aminophospholipid translocators), and P_5B (no assigned specificity) subfamilies. These subfamilies represent the P-type ATPase transport functions necessary for survival in the Trypanosomatidae as P-type ATPases for each of these seven subfamilies are found in all Leishmania and Trypanosoma species included in this analysis. These P-type ATPase subfamilies are correlated with current molecular and biochemical knowledge of their function in trypanosomatid growth, adaptation, infectivity, and survival.

The confounding effects of high genetic diversity on the determination and interpretation of differential gene expression analysis in the parasitic nematode Haemonchus contortus

Differential expression analysis between parasitic nematode strains is commonly used to implicate candidate genes in anthelmintic resistance or other biological functions. We have tested the hypothesis that the high genetic diversity of an organism such as Haemonchus contortus could complicate such analyses. First, we investigated the extent to which sequence polymorphism affects the reliability of differential expression analysis between the genetically divergent H. contortus strains MHco3(ISE), MHco4(WRS) and MHco10(CAVR). Using triplicates of 20 adult female worms from each population isolated under parallel experimental conditions, we found that high rates of sequence polymorphism in RNAseq reads were associated with lower efficiency read mapping to gene models under default TopHat2 parameters, leading to biased estimates of inter-strain differential expression. We then showed it is possible to largely compensate for this bias by optimising the read mapping single nucleotide polymorphism (SNP) allowance and filtering out genes with particularly high single nucleotide polymorphism rates. Once the sequence polymorphism biases were removed, we then assessed the genuine transcriptional diversity between the strains, finding ≥824 differentially expressed genes across all three pairwise strain comparisons. This high level of inter-strain transcriptional diversity not only suggests substantive inter-strain phenotypic variation but also highlights the difficulty in reliably associating differential expression of specific genes with phenotypic differences. To provide a practical example, we analysed two gene families of potential relevance to ivermectin drug resistance; the ABC transporters and the ligand-gated ion channels (LGICs). Over half of genes identified as differentially expressed using default TopHat2 parameters were shown to be an artifact of sequence polymorphism differences. This work illustrates the need to account for sequence polymorphism in differential expression analysis. It also demonstrates that a large number of genuine transcriptional differences can occur between H. contortus strains and these must be considered before associating the differential expression of specific genes with phenotypic differences between strains.

Transcriptional responses of Leishmania (Leishmania) amazonensis in the presence of trivalent sodium stibogluconate

Background

In the last decade, resistance to antimonials has become a serious problem due to the emergence of drug-resistant strains. Therefore, understanding the mechanisms used by Leishmania parasites to survive under drug pressure is essential, particularly for species of medical-veterinary importance such as L. amazonensis.

Methods

Here, we used RNA-seq technology to analyse transcriptome profiles and identify global changes in gene expression between antimony-resistant and -sensitive L. amazonensis promastigotes.

Results

A total of 723 differentially expressed genes were identified between resistant and sensitive lines. Comparative transcriptomic analysis revealed that genes encoding proteins involved in metabolism (fatty acids) and stress response, as well as those associated with antimony resistance in other Leishmania species, were upregulated in the antimony-resistant line. Most importantly, we observed upregulation of genes encoding autophagy proteins, suggesting that in the presence of trivalent stibogluconate (Sb^III) L. amazonensis can activate these genes either as a survival strategy or to induce cell death, as has been observed in other parasites.

Conclusions

This work identified global transcriptomic changes in an in vitro-adapted strain in response to Sb^III. Our results provide relevant information to continue understanding the mechanism used by parasites of the subgenus Leishmania (L. amazonensis) to generate an antimony-resistant phenotype.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3603-8) contains supplementary material, which is available to authorized users.

The mRNA-bound Proteome of Leishmania mexicana: Novel Genetic Insight into an Ancient Parasite

Leishmania parasite infections, termed the leishmaniases, cause significant global infectious disease burden. The lifecycle of the parasite embodies three main stages that require precise coordination of gene regulation to survive environmental shifts between sandfly and mammalian hosts. Constitutive transcription in kinetoplastid parasites means that gene regulation is overwhelmingly reliant on post-transcriptional mechanisms, yet strikingly few Leishmania trans-regulators are known. Using optimized crosslinking and deep, quantified mass spectrometry, we present a comprehensive analysis of 1400 mRNA binding proteins (mRBPs) and whole cell proteomes from the three main Leishmania lifecycle stages. Supporting the validity, although the crosslinked RBPome is magnitudes more enriched, the protein identities of the crosslinked and non-crosslinked RBPomes were nearly identical. Moreover, multiple candidate RBPs were endogenously tagged and found to associate with discrete mRNA target pools in a stage-specific manner. Results indicate that in L. mexicana parasites, mRNA levels are not a strong predictor of the whole cell expression or RNA binding potential of encoded proteins. Evidence includes a low correlation between transcript and corresponding protein expression and stage-specific variation in protein expression versus RNA binding potential. Unsurprisingly, RNA binding protein enrichment correlates strongly with relative replication efficiency of the specific lifecycle stage. Our study is the first to quantitatively define and compare the mRBPome of multiple stages in kinetoplastid parasites. It provides novel, in-depth insight into the trans-regulatory mRNA:Protein (mRNP) complexes that drive Leishmania parasite lifecycle progression.

Genetic dissection of a Leishmania flagellar proteome demonstrates requirement for directional motility in sand fly infections

The protozoan parasite Leishmania possesses a single flagellum, which is remodelled during the parasite’s life cycle from a long motile flagellum in promastigote forms in the sand fly to a short immotile flagellum in amastigotes residing in mammalian phagocytes. This study examined the protein composition and in vivo function of the promastigote flagellum. Protein mass spectrometry and label free protein enrichment testing of isolated flagella and deflagellated cell bodies defined a flagellar proteome for L. mexicana promastigote forms (available via ProteomeXchange with identifier PXD011057). This information was used to generate a CRISPR-Cas9 knockout library of 100 mutants to screen for flagellar defects. This first large-scale knockout screen in a Leishmania sp. identified 56 mutants with altered swimming speed (52 reduced and 4 increased) and defined distinct mutant categories (faster swimmers, slower swimmers, slow uncoordinated swimmers and paralysed cells, including aflagellate promastigotes and cells with curled flagella and disruptions of the paraflagellar rod). Each mutant was tagged with a unique 17-nt barcode, providing a simple barcode sequencing (bar-seq) method for measuring the relative fitness of L. mexicana mutants in vivo. In mixed infections of the permissive sand fly vector Lutzomyia longipalpis, paralysed promastigotes and uncoordinated swimmers were severely diminished in the fly after defecation of the bloodmeal. Subsequent examination of flies infected with a single paralysed mutant lacking the central pair protein PF16 or an uncoordinated swimmer lacking the axonemal protein MBO2 showed that these promastigotes did not reach anterior regions of the fly alimentary tract. These data show that L. mexicana need directional motility for successful colonisation of sand flies.

Leishmania are protozoan parasites, transmitted between mammals by the bite of phlebotomine sand flies. Promastigote forms in the sand fly have a long flagellum, which is motile and used for anchoring the parasites to prevent clearance with the digested blood meal remnants. To dissect flagellar functions and their importance in life cycle progression, we generated here a comprehensive list of >300 flagellar proteins and produced a CRISPR-Cas9 gene knockout library of 100 mutant Leishmania. We studied their behaviour in vitro before examining their fate in the sand fly Lutzomyia longipalpis. Measuring mutant swimming speeds showed that about half behaved differently compared to the wild type: a few swam faster, many slower and some were completely paralysed. We also found a group of uncoordinated swimmers. To test whether flagellar motility is required for parasite migration from the fly midgut to the foregut from where they reach the next host, we infected sand flies with a mixed mutant population. Each mutant carried a unique tag and tracking these tags up to nine days after infection showed that paralysed and uncoordinated Leishmania were rapidly lost from flies. These data indicate that directional swimming is important for successful colonisation of sand flies.

Genome and transcriptome analyses of Leishmania spp.: opening Pandora's box

In the last 30 years, significant advances in genetic manipulation tools along with complete genome and transcriptome sequencing have advanced our understanding of the biology of Leishmania parasites and their interplay with the sand fly and mammalian hosts. High-throughput sequencing in association with CRISPR/Cas9 have prepared the ground for significant advances. Given the richness of the progress made over the last decade, in this article, we focused on the most recent contributions of genome-wide and transcriptome analyses of Leishmania spp., which permit the comparison of life cycle stages, the evaluation of different strains and species in their natural niches and in the field and the simultaneously comparison of the gene expression profiles of parasites and hosts.

Complete assembly of the Leishmania donovani (HU3 strain) genome and transcriptome annotation

Leishmania donovani is a unicellular parasite that causes visceral leishmaniasis, a fatal disease in humans. In this study, a complete assembly of the genome of L. donovani is provided. Apart from being the first published genome of this strain (HU3), this constitutes the best assembly for an L. donovani genome attained to date. The use of a combination of sequencing platforms enabled to assemble, without any sequence gap, the 36 chromosomes for this species. Additionally, based on this assembly and using RNA-seq reads derived from poly-A + RNA, the transcriptome for this species, not yet available, was delineated. Alternative SL addition sites and heterogeneity in the poly-A addition sites were commonly observed for most of the genes. After a complete annotation of the transcriptome, 2,410 novel transcripts were defined. Additionally, the relative expression for all transcripts present in the promastigote stage was determined. Events of cis-splicing have been documented to occur during the maturation of the transcripts derived from genes LDHU3_07.0430 and LDHU3_29.3990. The complete genome assembly and the availability of the gene models (including annotation of untranslated regions) are important pieces to understand how differential gene expression occurs in this pathogen, and to decipher phenotypic peculiarities like tissue tropism, clinical disease, and drug susceptibility.

Comparative transcriptomics in Leishmania braziliensis: disclosing differential gene expression of coding and putative noncoding RNAs across developmental stages

Leishmaniasis is a worldwide public health problem caused by protozoan parasites of the genus Leishmania. Leishmania braziliensis is the most important species responsible for tegumentary leishmaniases in Brazil. An understanding of the molecular mechanisms underlying the success of this parasite is urgently needed. An in-depth study on the modulation of gene expression across the life cycle stages of L. braziliensis covering coding and noncoding RNAs (ncRNAs) was missing and is presented herein. Analyses of differentially expressed (DE) genes revealed that most prominent differences were observed between the transcriptomes of insect and mammalian proliferative forms (6,576 genes). Gene ontology (GO) analysis indicated stage-specific enriched biological processes. A computational pipeline and 5 ncRNA predictors allowed the identification of 11,372 putative ncRNAs. Most of the DE ncRNAs were found between the transcriptomes of insect and mammalian proliferative stages (38%). Of the DE ncRNAs, 295 were DE in all three stages and displayed a wide range of lengths, chromosomal distributions and locations; many of them had a distinct expression profile compared to that of their protein-coding neighbors. Thirty-five putative ncRNAs were submitted to northern blotting analysis, and one or more hybridization-positive signals were observed in 22 of these ncRNAs. This work presents an overview of the L. braziliensis transcriptome and its adjustments throughout development. In addition to determining the general features of the transcriptome at each life stage and the profile of protein-coding transcripts, we identified and characterized a variety of noncoding transcripts. The novel putative ncRNAs uncovered in L. braziliensis might be regulatory elements to be further investigated.

Genomic instability at the locus of sterol C24-methyltransferase promotes amphotericin B resistance in Leishmania parasites

Amphotericin B is an increasingly important tool in efforts to reduce the global disease burden posed by Leishmania parasites. With few other chemotherapeutic options available for the treatment of leishmaniasis, the potential for emergent resistance to this drug is a considerable threat. Here we characterised four novel amphotericin B-resistant Leishmania mexicana lines. All lines exhibited altered sterol biosynthesis, and hypersensitivity to pentamidine. Whole genome sequencing demonstrated resistance-associated mutation of the sterol biosynthesis gene sterol C5-desaturase in one line. However, in three out of four lines, RNA-seq revealed loss of expression of sterol C24-methyltransferase (SMT) responsible for drug resistance and altered sterol biosynthesis. Additional loss of the miltefosine transporter was associated with one of those lines. SMT is encoded by two tandem gene copies, which we found to have very different expression levels. In all cases, reduced overall expression was associated with loss of the 3' untranslated region of the dominant gene copy, resulting from structural variations at this locus. Local regions of sequence homology, between the gene copies themselves, and also due to the presence of SIDER1 retrotransposon elements that promote multi-gene amplification, correlate to these structural variations. Moreover, in at least one case loss of SMT expression was not associated with loss of virulence in primary macrophages or in vivo. Whilst such repeat sequence-mediated instability is known in Leishmania genomes, its presence associated with resistance to a major antileishmanial drug, with no evidence of associated fitness costs, is a significant concern.