KKT23 is a kinetoplastid-specific kinetochore protein that has a C-terminal GCN5-related histone acetyltransferase domain that acetylates the C-terminal tail of histone H2A. Here, we present the 1H, 13C and 15N resonance assignments for the C-terminal region of KKT23 (KKT23125-348) from Trypanosoma brucei in complex with known cofactors for acetyltransferases, acetyl coenzyme A and coenzyme A. These assignments provide the starting point for detailed investigation of the structure, dynamics and interactions of KKT23 in solution.

Doublet microtubules (DMTs) are flagellar components required for the protist Trichomonas vaginalis (Tv) to swim through the human genitourinary tract to cause trichomoniasis, the most common non-viral sexually transmitted disease. Lack of structures of Tv's DMT (Tv-DMT) has prevented structure-guided drug design to manage Tv infection. Here, we determine the 16 nm, 32 nm, 48 nm and 96 nm-repeat structures of native Tv-DMT at resolution ranging from 3.4 to 4.4 Å by cryogenic electron microscopy (cryoEM) and built an atomic model for the entire Tv-DMT. These structures show that Tv-DMT is composed of 30 different proteins, including the α- and β-tubulin, 19 microtubule inner proteins (MIPs) and 9 microtubule outer proteins. While the A-tubule of Tv-DMT is simplistic compared to DMTs of other organisms, the B-tubule of Tv-DMT features parasite-specific proteins, such as TvFAP40 and TvFAP35. Notably, TvFAP40 and TvFAP35 form filaments near the inner and outer junctions, respectively, and interface with stabilizing MIPs. This atomic model of the Tv-DMT highlights diversity of eukaryotic motility machineries and provides a structural framework to inform rational design of therapeutics against trichomoniasis.

Leishmania donovani is the causative agent of the fatal visceral leishmaniasis (VL) disease in humans in the tropical regions, mainly the Indian Subcontinent and Africa. We have previously described centrin1, a basal body associated cell division specific protein in this parasite important for the parasite's host intracellular stage. In this study, we identified a novel centrin1-binding protein called LdDRP through pull-down and MS/MS analysis, which is a homolog of the XPC protein of humans involved in DNA damage. The protein interaction with LdCen1 was also confirmed through peptide spectrum analysis against the UniProt database. Immunofluorescence analysis confirms that LdDRP is localized within the nucleus, suggesting the protein's possible role in DNA interaction. The overexpression of three LdDRP forms in the parasite, each fused with HA-tag (LdDRPF [full length] LdDRPN [only N-terminal], and LdDRPC [only C-terminal]), revealed that only LdDRPF and LdDRPC were able to support the retention of the parasite's shape and promote rapid division following the UV-damage recovery period. This was also correlated to the elevated expression level of both LdDRPC and LdCen1, by Western blot analysis soon after UV-C exposure in the parasites compared to control. The study emphasizes the role of the LdDRP, and its crucial domains involved in the DNA binding process, DNA damage response, and interaction with centrin, particularly in response to UV-C light-induced DNA damage.

Diplonemids are among the most abundant and species-rich protists in the oceans. Marine heterotrophic flagellates, including diplonemids, have been suggested to play important roles in global biogeochemical cycles. Diplonemids are also the sister taxon of kinetoplastids, home to trypanosomatid parasites of global health importance, and thus are informative about the evolution of kinetoplastid biology. However, the genomic and cellular complement that underpins diplonemids' highly successful lifestyle is underexplored. At the same time, our framework describing cellular processes may not be as broadly applicable as presumed, as it is largely derived from animal and fungal model organisms, a small subset of extant eukaryotic diversity. In addition to uniquely evolved machinery in animals and fungi, there exist components with sporadic (i.e., "patchy") distributions across other eukaryotes. A most intriguing subset are components ("jötnarlogs") stochastically present in a wide range of eukaryotes but lost in animal and/or fungal models. Such components are considered exotic curiosities but may be relevant to inferences about the complexity of the last eukaryotic common ancestor (LECA) and frameworks of modern cell biology. Here, we use comparative genomics and phylogenetics to comprehensively assess the membrane-trafficking system of diplonemids. They possess several proteins thought of as kinetoplastid specific, as well as an extensive set of patchy proteins, including jötnarlogs. Diplonemids apparently function with endomembrane machinery distinct from existing cell biological models but comparable with other free-living heterotrophic protists, highlighting the importance of including such exotic components when considering different models of ancient eukaryotic genomic complexity and the cell biology of non-opisthokont organisms.

Neglected tropical diseases such as Chagas disease, human African trypanosomiasis, leishmaniasis, and schistosomiasis have a significant global health impact in predominantly developing countries, although these diseases are spreading due to increased international travel and population migration. Drug repurposing with a focus on increasing antiparasitic potency and drug-like properties is a cost-effective and efficient route to the development of new therapies. Here we identify compounds that have potent activity against Trypanosoma cruzi and Leishmania donovani, and the latter were progressed into the murine model of infection. Despite the potent in vitro activity, there was no effect on parasitemia, necessitating further work to improve the pharmacokinetic properties of this series. Nonetheless, valuable insights have been obtained into the structure-activity and structure-property relationships of this compound series.

The exon junction complex (EJC) is a key player in metazoan mRNA quality control and is placed upstream of the exon-exon junction after splicing. Its inner core is composed of Magoh, Y14, eIF4AIII and BTZ and the outer core of proteins involved in mRNA splicing (CWC22), export (Yra1), translation (PYM) and nonsense mediated decay (NMD, UPF1/2/3). Trypanosoma brucei encodes only two genes with introns, but all mRNAs are processed by trans-splicing. The presence of three core EJC proteins and a potential BTZ homologue (Rbp25) in trypanosomes has been suggested to adapt of the EJC function to mark trans-spliced mRNAs. We analysed trypanosome EJC components and noticed major differences between eIF4AIII and Magoh/Y14: (i) whilst eIF4AIII is essential, knocking out both Magoh and Y14 elicits only a mild growth phenotype (ii) eIF4AIII localization is mostly nucleolar, while Magoh and Y14 are nucleolar and nucleoplasmic but excluded from the cytoplasm (iii) eIF4AIII associates with nucleolar proteins and the splicing factor CWC22, but not with Y14 or Magoh, while Magoh and Y14 associate with each other, but not with eIF4AIII, CWC22 or nucleolar proteins. Our data argue against the presence of a functional EJC in trypanosomes, but indicate that eIF4AIII adopted non-EJC related, essential functions, while Magoh and Y14 became redundant. Trypanosomes also possess homologues to the NMD proteins UPF1 and UPF2. Depletion of UPF1 causes only a minor reduction in growth and phylogenetic analyses show several independent losses of UPF1 and UPF2, as well as complete loss of UPF3 in the Kinetoplastida group, indicating that UPF1-dependent NMD is not essential. Regardless, we demonstrate that UPF1 depletion restores the mRNA levels of a PTC reporter. Altogether, we show that the almost intron-less trypanosomes are in the process of losing the canonical EJC/NMD pathways: Y14 and Magoh have become redundant and the still-functional UPF1-dependent NMD pathway is not essential.

Trypanosoma cruzi, the protozoan causing Chagas disease, is responsible for a neglected tropical disease affecting millions in Latin America. Its genome contains rapidly evolving multigene families, such as mucins (TcMUC), trans-sialidases (TS), and mucin-associated surface proteins (MASP), which are essential for parasite transmission and disease mechanisms. However, methodological challenges in genome assembly and annotation have limited the characterization of these gene families, particularly MASPs.

We developed a bioinformatic pipeline for the automatic identification, characterization, and annotation of MASPs directly from T. cruzi genome assemblies. This algorithm, based on a manually curated MASP database and HMM-based identification of MASP diagnostic motifs, enables the robust classification of these molecules into canonical MASPs, MASP-related molecules (mostly pseudogenes), and chimeric sequences combining MASPs and TcMUC/TS genes. Validation against a rigorously annotated dataset demonstrated high accuracy, and allowed us to reclassify misanotated sequences and, more crucially, to accurately identify previously unrecognized canonical MASPs and MASP chimeras. This algorithm was then used to explore the MASP repertoire in the genomes of 13 parasite strains from different evolutionary lineages, revealing patterns of diversity. For instance, TcI and TcII strains exhibited higher ratios of canonical MASP/MASP-related molecules and a greater abundance of MASP chimeras, suggesting that their genomes are under strong selective pressures towards maintaining a broader panel of full-length MASP genes at the expense of pseudogenes. On the contrary, structural features of canonical MASPs, MASP-related sequences, and MASP-chimeras were largely conserved across parasite genomes.

This novel pipeline automates the annotation of MASPs, a key surface protein family unique to T. cruzi, improving genome annotation and enabling robust comparative analyses. It provides an essential tool for exploring the evolutionary dynamics of multigene families in T. cruzi and could be extended to other gene families.

Kinetoplastid parasites cause diseases that threaten human and animal health. To survive transitions between vertebrate hosts and insect vectors, these parasites rely on precise regulation of gene expression to adapt to environmental changes. Since gene regulation in kinetoplastids is primarily post-transcriptional, developing efficient genetic tools for modifying genes at their endogenous loci while preserving regulatory mRNA elements is crucial for studying their complex biology. We present a CRISPR/Cas9-based tagging system that preserves untranslated regulatory elements and uses a viral 2A peptide from Thosea asigna to generate two separate proteins from a single transcript: a drug-selectable marker and a tagged protein of interest. This dual-function design maintains native control elements, allowing discrimination between regulation of transcript abundance, translational efficiency, and post-translational events. We validate the system by tagging six Trypanosoma brucei proteins and demonstrate (i) high-efficiency positive selection and separation of drug-selectable marker and target protein, (ii) preservation of regulatory responses to environmental cues like heat shock and iron availability, and (iii) maintenance of stage-specific regulation during developmental transitions. This versatile toolkit is applicable to all kinetoplastids amenable to CRISPR/Cas9 editing, providing a powerful reverse genetic tool for studying post-transcriptional regulation and protein function in organisms where post-transcriptional control is dominant.

Gene expression is a multi-step process that converts DNA-encoded information into proteins, involving RNA transcription, maturation, degradation, and translation. While transcriptional control is a major regulator of protein levels, the role of post-transcriptional processes such as RNA processing and degradation is less well understood due to the challenge of measuring their contributions individually. To address this challenge, we investigated the control of gene expression in Trypanosoma brucei, a unicellular parasite assumed to lack transcriptional control. Instead, mRNA levels in T. brucei are controlled by post-transcriptional processes, which enabled us to disentangle the contribution of both processes to total mRNA levels. In this study, we developed an efficient metabolic RNA labeling approach and combined ultra-short metabolic labeling with transient transcriptome sequencing (TT-seq) to confirm the long-standing assumption that RNA polymerase II transcription is unregulated in T. brucei. In addition, we established thiol (SH)-linked alkylation for metabolic sequencing of RNA (SLAM-seq) to globally quantify RNA processing rates and half-lives. Our data, combined with scRNA-seq data, indicate that RNA processing and stability independently affect total mRNA levels and contribute to the variability seen between individual cells in African trypanosomes.

The Trypanosoma brucei genome is structurally complex. Eleven megabase-sized chromosomes each comprise a transcribed core flanked by silent subtelomeres, housing thousands of Variant Surface Glycoprotein (VSG) genes. Additionally, hundreds of sub-megabase chromosomes contain 177 bp repeats of unknown function, and VSG transcription sites localise to many telomeres. DNA replication dynamics have only been described in the megabase chromosome cores, and in the single active VSG transcription site. Using a Nanopore genome assembly, we show that megabase chromosome subtelomeres display a paucity of replication initiation events relative to the core, correlating with increased instability. In addition, replication of the active VSG transcription site is shown to originate from the telomere, likely causing targeted VSG recombination. Lastly, we provide evidence that the 177 bp repeats act as conserved DNA replication origins, explaining submegabase chromosome stability. Compartmentalized DNA replication therefore explains how T. brucei balances stable genome transmission with localised instability driving immune evasion.

Biological information is encoded in a sequence of biochemical molecules such as nucleic acids and amino acids, and nanopore sequencing is a long-read sequencing technology capable of directly decoding these molecules into human-readable signals. The long reads from nanopore sequencing offer the advantage of obtaining contiguous information, which is particularly beneficial for decoding complex or repetitive regions in a genome. In this study, we investigated the efficacy of nanopore sequencing in decoding biological information from distinctive genomes in metagenomic samples, which pose significant challenges for traditional short-read sequencing technologies. Specifically, we sequenced blood and fecal samples from mice infected with Trypanosoma brucei, a unicellular protozoan known for its hypervariable and dynamic regions that help it evade host immunity. Such characteristics are also prevalent in other host-dependent parasites, such as bacteriophages. The taxonomic classification results showed a high proportion of nanopore reads identified as T. brucei in the infected blood samples, with no significant identification in the control blood samples and fecal samples. Furthermore, metagenomic de novo assembly of these nanopore reads yielded contigs that mapped to the reference genome of T. brucei in the infected blood samples with over 96 % accuracy. This exploratory work demonstrates the potential of nanopore sequencing for the challenging task of classifying and assembling hypervariable and dynamic genomes from metagenomic samples.

Microtubule-severing enzymes play essential roles in regulating diverse cellular processes, including mitosis and cytokinesis, by modulating microtubule dynamics. In the early branching protozoan parasite Trypanosoma brucei, microtubule-severing enzymes are involved in cytokinesis and flagellum length control during different life cycle stages, but none of them have been found to regulate mitosis in any life cycle form. Here, we report the biochemical and functional characterization of the microtubule-severing enzyme spastin in the procyclic form of T. brucei. We demonstrate that spastin catalyzes microtubule severing in vitro and ectopic overexpression of spastin disrupts spindle microtubules in vivo in trypanosome cells, leading to defective chromosome segregation. Knockdown of spastin impairs spindle integrity and disrupts chromosome alignment in metaphase and chromosome segregation in anaphase. We further show that the function of spastin requires the catalytic AAA-ATPase domain, the microtubule-binding domain, and the microtubule interacting and trafficking domain, and that the association of spastin with spindle depends on the microtubule-binding domain. Together, these results uncover an essential role for spastin in chromosome segregation by regulating spindle dynamics in this unicellular eukaryote.

Animal trypanosomiasis, such as nagana, surra, and dourine, represent a significant challenge to animal health and economic development, especially in tropical and subtropical regions where livestock production is an essential component of a country's economy. Despite advances in the control of human trypanosomiasis, animal diseases caused by several species of trypanosomes remain neglected. The lack of funding for the development of new treatments and vaccines contributes to sustaining the severe economic impacts these diseases have on the farming industry, especially in low-income rural areas. Recent advances in the understanding of the immune processes involved during infection have been essential for the development of new approaches towards disease control including vaccines. These new approaches must be part of integrated control programs, which must also include vector management and the awareness of good veterinary practices. Addressing the challenges posed by the control of animal trypanosomiasis requires collaborative and continuous efforts shared among scientists, governments, and the farming industry, if significant progress is to be made to mitigate the impact of these diseases. In this literature review, we discuss the main challenges for the development of vaccines for animal trypanosomiasis and the research underway, including the prospects for employing new vaccine platforms, such as an mRNA vaccine, vector-based vaccine, and CRISPR-attenuated parasite vaccine.

Transmission of genetic material from one generation to the next is a fundamental feature of all living cells. In eukaryotes, a macromolecular complex called the kinetochore plays crucial roles during chromosome segregation by linking chromosomes to spindle microtubules. Little is known about this process in evolutionarily diverse protists. Within the supergroup Discoba, Euglenozoa forms a speciose group of unicellular flagellates-kinetoplastids, euglenids, and diplonemids. Kinetoplastids have an unconventional kinetochore system, while euglenids have subunits that are conserved among most eukaryotes. For diplonemids, a group of extremely diverse and abundant marine flagellates, it remains unclear what kind of kinetochores are present. Here, we employed deep homology detection protocols using profile-versus-profile Hidden Markov Model searches and AlphaFold-based structural comparisons to detect homologies that might have been previously missed. Interestingly, we still could not detect orthologs for most of the kinetoplastid or canonical kinetochore subunits with few exceptions including a putative centromere-specific histone H3 variant (cenH3/CENP-A), the spindle checkpoint protein Mad2, the chromosomal passenger complex members Aurora and INCENP, and broadly conserved proteins like CLK kinase and the meiotic synaptonemal complex proteins SYCP2/3 that also function at kinetoplastid kinetochores. We examined the localization of five candidate kinetochore-associated proteins in the model diplonemid, Paradiplonema papillatum. PpCENP-A shows discrete dots in the nucleus, implying that it is likely a kinetochore component. PpMad2, PpCLKKKT10/19, PpSYCP2L1KKT17/18, and PpINCENP reside in the nucleus, but no clear kinetochore localization was observed. Altogether, these results point to the possibility that diplonemids evolved a hitherto unknown type of kinetochore system.

A macromolecular assembly called the kinetochore is essential for the segregation of genetic material during eukaryotic cell division. Therefore, characterization of kinetochores across species is essential for understanding the mechanisms involved in this key process across the eukaryotic tree of life. In particular, little is known about kinetochores in divergent protists such as Euglenozoa, a group of unicellular flagellates that includes kinetoplastids, euglenids, and diplonemids, the latter being a highly diverse and abundant component of marine plankton. While kinetoplastids have an unconventional kinetochore system and euglenids have a canonical one similar to traditional model eukaryotes, preliminary searches detected neither unconventional nor canonical kinetochore components in diplonemids. Here, we employed state-of-the-art deep homology detection protocols but still could not detect orthologs for the bulk of kinetoplastid-specific nor canonical kinetochore proteins in diplonemids except for a putative centromere-specific histone H3 variant. Our results suggest that diplonemids evolved kinetochores that do not resemble previously known ones.

In trypanosomatids, a group of unicellular eukaryotes that includes numerous important human parasites, cis-splicing has been previously reported for only two genes: a poly(A) polymerase and an RNA helicase. Conversely, trans-splicing, which involves the attachment of a spliced leader sequence, is observed for nearly every protein-coding transcript. So far, our understanding of splicing in this protistan group has stemmed from the analysis of only a few medically relevant species. In this study, we used an extensive dataset encompassing all described trypanosomatid genera to investigate the distribution of intron-containing genes and the evolution of splice sites.

We identified a new conserved intron-containing gene encoding an RNA-binding protein that is universally present in Kinetoplastea. We show that Perkinsela sp., a kinetoplastid endosymbiont of Amoebozoa, represents the first eukaryote completely devoid of cis-splicing, yet still preserving trans-splicing. We also provided evidence for reverse transcriptase-mediated intron loss in Kinetoplastea, extensive conservation of 5' splice sites, and the presence of non-coding RNAs within a subset of retained trypanosomatid introns.

All three intron-containing genes identified in Kinetoplastea encode RNA-interacting proteins, with a potential to fine-tune the expression of multiple genes, thus challenging the perception of cis-splicing in these protists as a mere evolutionary relic. We suggest that there is a selective pressure to retain cis-splicing in trypanosomatids and that this is likely associated with overall control of mRNA processing. Our study provides new insights into the evolution of introns and, consequently, the regulation of gene expression in eukaryotes.

Primary cilia are antenna-like sensory organelles that are evolutionarily conserved in nearly all modern eukaryotes, from the single-celled green alga, Chlamydomonas reinhardtii, to vertebrates and mammals. Cilia are microtubule-based cellular projections that have adapted to perform a broad range of species-specific functions, from cell motility to detection of light and the transduction of extracellular mechanical and chemical signals. These functions render cilia essential for organismal development and survival. The high conservation of cilia has allowed for discoveries in C. reinhardtii to inform our understanding of the basic biology of mammalian primary cilia, and to provide insight into the genetic etiology of ciliopathies. Over the last two decades, a growing number of studies has revealed that multiple aspects of ciliary homeostasis are regulated by the actin cytoskeleton, including centrosome migration and positioning, vesicle transport to the basal body, ectocytosis, and ciliary-mediated signaling. Here, we review actin regulation of ciliary homeostasis, and highlight conserved and divergent mechanisms in C. reinhardtii and mammalian cells. Further, we compare the disease manifestations of patients with ciliopathies to those with mutations in actin and actin-associated genes, and propose that primary cilia defects caused by genetic alteration of the actin cytoskeleton may underlie certain birth defects.

Faithful chromosome segregation in eukaryotes requires the assembly of a bipolar spindle and the faithful attachment of kinetochores to spindle microtubules, which are regulated by various spindle-associated proteins (SAPs) that play distinct functions in regulating spindle dynamics and microtubule-kinetochore attachment. The protozoan parasite Trypanosoma brucei employs evolutionarily conserved and kinetoplastid-specific proteins, including some kinetoplastid-specific nucleus- and spindle-associated proteins (NuSAPs), to regulate chromosome segregation. Here, we characterized NuSAP4 and its functional interplay with diverse SAPs in promoting chromosome segregation in T. brucei. NuSAP4 associates with the spindle during mitosis and concentrates at spindle poles where it interacts with SPB1 and MAP103. Knockdown of NuSAP4 impairs chromosome segregation by disrupting bipolar spindle assembly and spindle pole protein localization. These results uncover the mechanistic role of NuSAP4 in regulating chromosome segregation by promoting bipolar spindle assembly, and highlight the unusual features of mitotic regulation by spindle-associated proteins in this early divergent microbial eukaryote.

Trypanosomes are the extracellular protozoan parasites that cause human African trypanosomiasis disease in humans and nagana disease in animals. Tsetse flies act as a vector for the transmission of the disease in African countries. Animals infected with these parasites become useless or workless, and if not treated, disease can be fatal. There are many side effects associated with old treatments and some of them result in death in 5% of cases. There is a major surface glycoprotein in the parasite known as variant surface glycoprotein. The immune system of the host develops antibodies against this antigen but due to antigenic variation, parasites evade the immune response. Currently, no vaccine is available that provides complete protection. In murine models, only partial protection was observed using certain antigens. In order to develop vaccines against trypanosomes, molecular biology and immunology tools have been used. Immunization is the sole method for the control of disease because the eradication of the vector from endemic areas is an impossible task. Genetic vaccines can carry multiple genes encoding different antigens of the same parasite or different parasites. DNA immunization induces the activation of both cellular immune response and humoral immune response along with the generation of memory. This review highlights the importance of DNA vaccines and advances in the development of DNA vaccines against T. brucei.

Sets of electrophilic probes are generally prepared using a narrow toolkit of robust reactions, which tends to limit both their structural and functional diversity. A unified synthesis of skeletally-diverse sulfonyl fluorides was developed that relied upon photoredox-catalysed dehydrogenative couplings between hetaryl sulfonyl fluorides and hydrogen donor building blocks. A set of 32 diverse probes was prepared, and then screened against Trypanosoma brucei. Four of the probes were found to have sub-micromolar anti-trypanosomal activity. A chemical proteomic approach, harnessing an alkynylated analogue and broad-spectrum fluorophosphonate tools, provided insights into the observed anti-trypanosomal activity, which likely stems from covalent modification of multiple protein targets. It is envisaged that the unified diversity-oriented approach may enable the discovery of electrophilic probes that have value in the elucidation of biological and biomedical mechanisms.

The role of malate dehydrogenase (MDH) in the metabolism of various medically significant protozoan parasites is reviewed. MDH is an NADH-dependent oxidoreductase that catalyzes interconversion between oxaloacetate and malate, provides metabolic intermediates for both catabolic and anabolic pathways, and can contribute to NAD+/NADH balance in multiple cellular compartments. MDH is present in nearly all organisms; isoforms of MDH from apicomplexans (Plasmodium falciparum, Toxoplasma gondii, Cryptosporidium spp.), trypanosomatids (Trypanosoma brucei, T. cruzi) and anaerobic protozoans (Trichomonas vaginalis, Giardia duodenalis) are presented here. Many parasitic species have complex life cycles and depend on the environment of their hosts for carbon sources and other nutrients. Metabolic plasticity is crucial to parasite transition between host environments; thus, the regulation of metabolic processes is an important area to explore for therapeutic intervention. Common themes in protozoan parasite metabolism include emphasis on glycolytic catabolism, substrate-level phosphorylation, non-traditional uses of common pathways like tricarboxylic acid cycle and adapted or reduced mitochondria-like organelles. We describe the roles of MDH isoforms in these pathways, discuss unusual structural or functional features of these isoforms relevant to activity or drug targeting, and review current studies exploring the therapeutic potential of MDH and related genes. These studies show that MDH activity has important roles in many metabolic pathways, and thus in the metabolic transitions of protozoan parasites needed for success as pathogens.

Damage to the genetic material of the cell poses a universal threat to all forms of life. The DNA damage response is a coordinated cellular response to a DNA break, key to which is the phosphorylation signaling cascade. Identifying which proteins are phosphorylated is therefore crucial to understanding the mechanisms that underlie it. We have used stable isotopic labeling of amino acids in cell culture-based quantitative phosphoproteomics to profile changes in phosphorylation site abundance following double stranded DNA breaks, at two distinct loci in the genome of the single cell eukaryote Trypanosoma brucei. Here, we report on the T. brucei phosphoproteome following a single double-strand break at either a chromosome internal or subtelomeric locus, specifically the bloodstream form expression site. We detected >6500 phosphorylation sites, of which 211 form a core set of double-strand break responsive phosphorylation sites. Along with phosphorylation of canonical DNA damage factors, we have identified two novel phosphorylation events on histone H2A and found that in response to a chromosome internal break, proteins are predominantly phosphorylated, while a greater proportion of proteins dephosphorylated following a DNA break at a subtelomeric bloodstream form expression site. Our data represent the first DNA damage phosphoproteome and provides novel insights into repair at distinct chromosomal contexts in T. brucei.

Antigenic variation, using large genomic repertoires of antigen-encoding genes, allows pathogens to evade host antibody. Many pathogens, including the African trypanosome Trypanosoma brucei, extend their antigenic repertoire through genomic diversification. While evidence suggests that T. brucei depends on the generation of new variant surface glycoprotein (VSG) genes to maintain a chronic infection, a lack of experimentally tractable tools for studying this process has obscured its underlying mechanisms. Here, we present a highly sensitive targeted sequencing approach for measuring VSG diversification. Using this method, we demonstrate that a Cas9-induced DNA double-strand break within the VSG coding sequence can induce VSG recombination with patterns identical to those observed during infection. These newly generated VSGs are antigenically distinct from parental clones and thus capable of facilitating immune evasion. Together, these results provide insight into the mechanisms of VSG diversification and an experimental framework for studying the evolution of antigen repertoires in pathogenic microbes.

This systematic review evaluated the efficacy and safety of photodynamic therapy (PDT) in the management of cutaneous leishmaniasis (CL). The electronic search for identification of relevant studies, adhered to the PICOS (Population, Intervention, Comparator, Outcomes and Study type) framework, was conducted through PubMed, Google scholar, Dimensions, X-mol, and Semantic Scholar till December 2023. All types of studies reporting PDT in the management of CL with no language restriction were included. Methodological quality appraised of the selected studies was performed using Jadad index. Of the 317 identified studies, 21 reported PDT for the treatment of CL lesions, consisting of two randomized controlled trials (RCTs), four single-center open study, one case series and 14 case reports. Collectively, these studies presented a total of 304 patients with ages ranging from 1 to 82 years, undergoing varying number of PDT sessions (3-28) and follow-up durations spanning 4 weeks to 24 months. The CL lesions predominantly manifested on the exposed body areas, such as face, limbs, neck, ear and nose, and characterized with the use of clinical variables, such as plaques, papules, erythema and ulceration. PDT protocols differed in the photosensitizer type, incubation time, light source characteristics (e.g., wavelength, output power, and energy density), duration of light illumination, number of PDT sessions and their respective frequencies. Treatment response was assessed through the clinical presentation (i.e., at the baseline and after PDT completion) or by the absence of Leishmania parasites. Adverse effects comprised of pain, burning and tingling sensation experienced during PDT, followed by erythema, pigmentation changes and edema post-treatment. This systematic review revealed that PDT is an efficacious and safe modality for the treatment of CL, with mild and transient side effects.

Leishmaniasis, caused by the protozoan parasite poses a significant health burden globally. With a very few specific drugs, increased drug resistance it is important to look for drug repurposing along with the identification of pre-clinical candidates against visceral leishmaniasis. This study aims to identify potential drug candidates against visceral leishmaniasis by targeting leishmanial MAP kinases and screening FDA approved protein kinase inhibitors.

MAP kinases were identified from the Leishmania genome. 12 FDA approved protein kinase inhibitors were screened against Leishmania MAP kinases. Binding affinity, ADME and toxicity of identified drug candidates were profiled. The anti-proliferative effects and mechanism of action were assessed in Leishmania, including changes in cell morphology, flagellar length, cell cycle progression, reactive oxygen species (ROS) generation, and intra-macrophage parasitic burden.

23 MAP kinases were identified from the Leishmania genome. Sorafenib and imatinib emerged as repurposable drug candidates and demonstrated excellent anti-proliferative effects in Leishmania. Treatment with these inhibitors resulted in significant changes in cell morphology, flagellar length, and cell cycle arrest. Furthermore, sorafenib and imatinib promoted ROS generation and reduced intra-macrophage parasitic burden, and elicited anti-leishmanial activity in in vivo experimental VL models.

Collectively, these results imply involvement of MAP kinases in infectivity and survival of the parasite and can pave the avenue for repurposing sorafenib and imatinib as anti-leishmanial agents. These findings contribute to the exploration of new treatment options for visceral leishmaniasis, particularly in the context of emerging drug resistance.

Outer dynein arms (ODAs) are responsible for ciliary beating in eukaryotes. They are assembled in the cytoplasm and shipped by intraflagellar transport (IFT) before attachment to microtubule doublets via the docking complex. The LRRC56 protein has been proposed to contribute to ODAs maturation. Mutations or deletion of the LRRC56 gene lead to reduced ciliary motility in all species investigated so far, but with variable impact on dynein arm presence. Here, we investigated the role of LRRC56 in the protist Trypanosoma brucei, where its absence results in distal loss of ODAs, mostly in growing flagella. We show that LRRC56 is a transient cargo of IFT trains during flagellum construction and surprisingly, is required for efficient attachment of a subset of docking complex proteins present in the distal portion of the organelle. This relation is interdependent since the knockdown of the distal docking complex prevents LRRC56's association with the flagellum. Intriguingly, lrrc56-/- cells display shorter flagella whose maturation is delayed. Inhibition of cell division compensates for the distal ODAs absence thanks to the redistribution of the proximal docking complex, restoring ODAs attachment but not the flagellum length phenotype. This work reveals an unexpected connection between LRRC56 and the docking complex.

The hyper-modified DNA base J helps control termination of Pol II transcription at polycistronic transcription units (PTUs) in T. brucei and L. major , allowing epigenetic control of gene expression. The Telomere Repeat-containing RNA (TERRA) is synthesized in T. brucei by Pol I readthrough transcription of a telomeric PTU. While little is understood regarding TERRA synthesis and function, the hyper-modified DNA base J is highly enriched at telomeres in L. major promastigotes. We now show that TERRA is synthesized by Pol II in L. major and loss of base J leads to increased TERRA. For at least one site, the increased TERRA is by Pol II readthrough transcription from an adjacent PTU. Furthermore, Pol II readthrough defects and increased TERRA correlate with increased differentiation of promastigotes to the infectious metacyclic life stage and decreased cell viability. These results help explain the essential nature of base J in Leishmania and provide insight regarding epigenetic control of coding and non-coding RNA expression and parasite development during the life cycle of L. major .

Whole genome sequencing (WGS) plays an important role in the advanced characterization of pathogen transmission and is widely used in studies of major bacterial and viral diseases. Although protozoan parasites cause serious diseases in humans and animals, WGS data on them are relatively scarce due to the large genomes and lack of cultivation techniques for some. In this review, we have illustrated bioinformatic analyses of WGS data and their applications in studies of the genomic epidemiology of apicomplexan parasites. WGS has been used in outbreak detection and investigation, studies of pathogen transmission and evolution, and drug resistance surveillance and tracking. However, comparative analysis of parasite WGS data is still in its infancy, and available WGS data are mainly from a few genera of major public health importance, such as Plasmodium, Toxoplasma, and Cryptosporidium. In addition, the utility of third-generation sequencing technology for complete genome assembly at the chromosome level, studies of the biological significance of structural genomic variation, and molecular surveillance of pathogens has not been fully exploited. These issues require large-scale WGS of various protozoan parasites of public health and veterinary importance using both second- and third-generation sequencing technologies.

Faithful transmission of genetic material is crucial for the survival of all organisms. In many eukaryotes, a feedback control mechanism called the spindle checkpoint ensures chromosome segregation fidelity by delaying cell cycle progression until all chromosomes achieve proper attachment to the mitotic spindle. Kinetochores are the macromolecular complexes that act as the interface between chromosomes and spindle microtubules. While most eukaryotes have canonical kinetochore proteins that are widely conserved, kinetoplastids such as Trypanosoma brucei have a seemingly unique set of kinetochore proteins including KKT1-25. It remains poorly understood how kinetoplastids regulate cell cycle progression or ensure chromosome segregation fidelity. Here, we report a crystal structure of the C-terminal domain of KKT14 from Apiculatamorpha spiralis and uncover that it is a pseudokinase. Its structure is most similar to the kinase domain of a spindle checkpoint protein Bub1. In addition, KKT14 has a putative ABBA motif that is present in Bub1 and its paralogue BubR1. We also find that the N-terminal part of KKT14 interacts with KKT15, whose WD40 repeat beta-propeller is phylogenetically closely related to a direct interactor of Bub1/BubR1 called Bub3. Our findings indicate that KKT14-KKT15 are divergent orthologues of Bub1/BubR1-Bub3, which promote accurate chromosome segregation in trypanosomes.

KKT4 is a multi-domain kinetochore protein specific to kinetoplastids, such as Trypanosoma brucei. It lacks significant sequence similarity to known kinetochore proteins in other eukaryotes. Our recent X-ray structure of the C-terminal region of KKT4 shows that it has a tandem BRCT (BRCA1 C Terminus) domain fold with a sulfate ion bound in a typical binding site for a phosphorylated serine or threonine. Here we present the 1H, 13C and 15N resonance assignments for the BRCT domain of KKT4 (KKT4463-645) from T. brucei. We show that the BRCT domain can bind phosphate ions in solution using residues involved in sulfate ion binding in the X-ray structure. We have used these assignments to characterise the secondary structure and backbone dynamics of the BRCT domain in solution. Mutating the residues involved in phosphate ion binding in T. brucei KKT4 BRCT results in growth defects confirming the importance of the BRCT phosphopeptide-binding activity in vivo. These results may facilitate rational drug design efforts in the future to combat diseases caused by kinetoplastid parasites.

Inhibition of hypoxanthine-guanine-xanthine phosphoribosyltransferase activity decreases the pool of 6-oxo and 6-amino purine nucleoside monophosphates required for DNA and RNA synthesis, resulting in a reduction in cell growth. Therefore, inhibitors of this enzyme have potential to control infections, caused by Plasmodium falciparum and Plasmodium vivax, Trypanosoma brucei, Mycobacterium tuberculosis, and Helicobacter pylori. Five compounds synthesized here that contain a purine base covalently linked by a prolinol group to one or two phosphonate groups have Ki values ranging from 3 nM to >10 μM, depending on the structure of the inhibitor and the biological origin of the enzyme. X-ray crystal structures show that, on binding, these prolinol-containing inhibitors stimulated the movement of active site loops in the enzyme. Against TBr in cell culture, a prodrug exhibited an EC50 of 10 μM. Thus, these compounds are excellent candidates for further development as drug leads against infectious diseases as well as being potential anticancer agents.

In dinoflagellates, a unique and extremely divergent genomic and nuclear organization has evolved. The highly unusual features of dinoflagellate nuclei and genomes include permanently condensed liquid crystalline chromosomes, primarily packaged by proteins other than histones, genes organized in very long unidirectional gene arrays, a general absence of transcriptional regulation, high abundance of the otherwise very rare DNA modification 5-hydroxymethyluracil (5-hmU), and many others. While most of these fascinating properties are originally identified in the 1970s and 1980s, they have not yet been investigated using modern genomic tools.

In this work, we address some of the outstanding questions regarding dinoflagellate genome organization by mapping the genome-wide distribution of 5-hmU (using both immunoprecipitation-based and basepair-resolution chemical mapping approaches) and of chromatin accessibility in the genome of the Symbiodiniaceae dinoflagellate Breviolum minutum. We find that the 5-hmU modification is preferentially enriched over certain classes of repetitive elements, often coincides with the boundaries between gene arrays, and is generally correlated with decreased chromatin accessibility, the latter otherwise being largely uniform along the genome. We discuss the potential roles of 5-hmU in the functional organization of dinoflagellate genomes and its relationship to the transcriptional landscape of gene arrays.

Our results provide the first window into the 5-hmU and chromatin accessibility landscapes in dinoflagellates.

The chromosomal passenger complex (CPC) is an important regulator of cell division, which shows dynamic subcellular localization throughout mitosis, including kinetochores and the spindle midzone. In traditional model eukaryotes such as yeasts and humans, the CPC consists of the catalytic subunit Aurora B kinase, its activator INCENP, and the localization module proteins Borealin and Survivin. Intriguingly, Aurora B and INCENP as well as their localization pattern are conserved in kinetoplastids, an evolutionarily divergent group of eukaryotes that possess unique kinetochore proteins and lack homologs of Borealin or Survivin. It is not understood how the kinetoplastid CPC assembles nor how it is targeted to its subcellular destinations during the cell cycle. Here, we identify two orphan kinesins, KIN-A and KIN-B, as bona fide CPC proteins in Trypanosoma brucei, the kinetoplastid parasite that causes African sleeping sickness. KIN-A and KIN-B form a scaffold for the assembly of the remaining CPC subunits. We show that the C-terminal unstructured tail of KIN-A interacts with the KKT8 complex at kinetochores, while its N-terminal motor domain promotes CPC translocation to spindle microtubules. Thus, the KIN-A:KIN-B complex constitutes a unique 'two-in-one' CPC localization module, which directs the CPC to kinetochores from S phase until metaphase and to the central spindle in anaphase. Our findings highlight the evolutionary diversity of CPC proteins and raise the possibility that kinesins may have served as the original transport vehicles for Aurora kinases in early eukaryotes.

The hemoparasite Trypanosoma equiperdum belongs to the Trypanozoon subgenus and includes several species that are pathogenic to animals and humans in tropical and subtropical areas across the world. As with all eukaryotic organisms, Ca2+ is essential for these parasites to perform cellular processes thus ensuring their survival across their life cycle. Despite the established paradigm to study proteins related to Ca2+ homeostasis as potential drug targets, so far little is known about Ca2+ entry into trypanosomes. Therefore, in the present study, the presence of a plasma membrane Ca2+-channel in T. equiperdum (TeCC), activated by sphingosine and inhibited by verapamil, is described. The TeCC was cloned and analyzed using bioinformatic resources, which confirmed the presence of several domains, motifs, and a topology similar to the Ca2+ channels found in higher eukaryotes. Biochemical and confocal microscopy assays using antibodies raised against an internal region of human L-type Ca2+ channels indicate the presence of a protein with similar predicted molar mass to the sequence analyzed, located at the plasma membrane of T. equiperdum. Physiological assays based on Fura-2 signals and Mn2+ quenching performed on whole parasites showed a unidirectional Ca2+ entry, which is activated by sphingosine and blocked by verapamil, with the distinctive feature of insensitivity to nifedipine and Bay K 8644. This suggests a second Ca2+ entry for T. equiperdum, different from the store-operated Ca2+ entry (SOCE) previously described. Moreover, the evidence presented here for the TeCC indicates molecular and pharmacological differences with their mammal counterparts, which deserve further studies to evaluate the potential of this channel as a drug target.

Diseases caused by Leishmania and Trypanosoma parasites are a major health problem in tropical countries. Because of their complex life cycle involving both vertebrate and insect hosts, and >1 billion years of evolutionarily distance, the cell biology of trypanosomatid parasites exhibits pronounced differences to animal cells. For example, the actin cytoskeleton of trypanosomatids is divergent when compared with other eukaryotes. To understand how actin dynamics are regulated in trypanosomatid parasites, we focused on a central actin-binding protein profilin. Co-crystal structure of Leishmania major actin in complex with L. major profilin revealed that, although the overall folds of actin and profilin are conserved in eukaryotes, Leishmania profilin contains a unique α-helical insertion, which interacts with the target binding cleft of actin monomer. This insertion is conserved across the Trypanosomatidae family and is similar to the structure of WASP homology-2 (WH2) domain, a small actin-binding motif found in many other cytoskeletal regulators. The WH2-like motif contributes to actin monomer binding and enhances the actin nucleotide exchange activity of Leishmania profilin. Moreover, Leishmania profilin inhibited formin-catalyzed actin filament assembly in a mechanism that is dependent on the presence of the WH2-like motif. By generating profilin knockout and knockin Leishmania mexicana strains, we show that profilin is important for efficient endocytic sorting in parasites, and that the ability to bind actin monomers and proline-rich proteins, and the presence of a functional WH2-like motif, are important for the in vivo function of Leishmania profilin. Collectively, this study uncovers molecular principles by which profilin regulates actin dynamics in trypanosomatids.

The genomic region that lies between the telomere and chromosome body, termed the subtelomere, is heterochromatic, repeat-rich, and frequently undergoes rearrangement. Within this region, large-scale structural changes enable gene diversification, and, as such, large multicopy gene families are often found at the subtelomere. In some parasites, genes associated with proliferation, invasion, and survival are often found in these regions, where they benefit from the subtelomere's highly plastic, rapidly changing nature. The increasing availability of complete (or near complete) parasite genomes provides an opportunity to investigate these typically poorly defined and overlooked genomic regions and potentially reveal relevant gene families necessary for the parasite's lifestyle.

Using the latest chromosome-scale genome assembly and hallmark repeat richness observed at chromosome termini, we have identified and characterised the subtelomeres of Schistosoma mansoni, a metazoan parasitic flatworm that infects over 250 million people worldwide. Approximately 12% of the S. mansoni genome is classified as subtelomeric, and, in line with other organisms, we find these regions to be gene-poor but rich in transposable elements. We find that S. mansoni subtelomeres have undergone extensive interchromosomal recombination and that these sites disproportionately contribute to the 2.3% of the genome derived from segmental duplications. This recombination has led to the expansion of subtelomeric gene clusters containing 103 genes, including the immunomodulatory annexins and other gene families with unknown roles. The largest of these is a 49-copy plexin domain-containing protein cluster, exclusively expressed in the tegument-the tissue located at the host-parasite physical interface-of intramolluscan life stages.

We propose that subtelomeric regions act as a genomic playground for trial-and-error of gene duplication and subsequent divergence. Owing to the importance of subtelomeric genes in other parasites, gene families implicated in this subtelomeric expansion within S. mansoni warrant further characterisation for a potential role in parasitism.

Although located at the chromosome end, telomeres are an essential chromosome component that helps maintain genome integrity and chromosome stability from protozoa to mammals. The role of telomere proteins in chromosome end protection is conserved, where they suppress various DNA damage response machineries and block nucleolytic degradation of the natural chromosome ends, although the detailed underlying mechanisms are not identical. In addition, the specialized telomere structure exerts a repressive epigenetic effect on expression of genes located at subtelomeres in a number of eukaryotic organisms. This so-called telomeric silencing also affects virulence of a number of microbial pathogens that undergo antigenic variation/phenotypic switching. Telomere proteins, particularly the RAP1 homologs, have been shown to be a key player for telomeric silencing. RAP1 homologs also suppress the expression of Telomere Repeat-containing RNA (TERRA), which is linked to their roles in telomere stability maintenance. The functions of RAP1s in suppressing telomere recombination are largely conserved from kinetoplastids to mammals. However, the underlying mechanisms of RAP1-mediated telomeric silencing have many species-specific features. In this review, I will focus on Trypanosoma brucei RAP1's functions in suppressing telomeric/subtelomeric DNA recombination and in the regulation of monoallelic expression of subtelomere-located major surface antigen genes. Common and unique mechanisms will be compared among RAP1 homologs, and their implications will be discussed.

Aminoacyl-tRNA synthetases (aaRSs), essential components of the protein synthesizing machinery, have been often chosen for devising therapeutics against parasitic diseases. Due to their relevance in drug development, the current study was designed to explore functional and structural aspects of Leishmania donovani glutamyl-tRNA synthetase (LdGluRS). Hence, LdGluRS was cloned into an expression vector and purified to homogeneity using chromatographic techniques. Purified protein showed maximum enzymatic activity at physiological pH, with more binding capacity towards its cofactor (Adenosine triphosphate, 0.06 ± 0.01 mM) than the cognate substrate (L-glutamate, 9.5 ± 0.5 mM). Remarkably, salicylate inhibited LdGluRS competitively with respect to L-glutamate and exhibited druglikeness with negligible effect on human macrophages. The protein possessed more α-helices (43 %) than β-sheets (12 %), whereas reductions in thermal stability and cofactor-binding affinity, along with variation in mode of inhibition after mutation signified the role of histidine (H60) as a catalytic residue. LdGluRS could also generate a pro-inflammatory milieu in human macrophages by upregulating cytokines. The docking study demonstrated the placement of salicylate into LdGluRS substrate-binding site, and the complex was found to be stable during molecular dynamics (MD) simulation. Altogether, our study highlights the understanding of molecular inhibition and structural features of glutamyl-tRNA synthetase from kinetoplastid parasites.

Metatranscriptomics generates large volumes of sequence data about transcribed genes in natural environments. Taxonomic annotation of these datasets depends on availability of curated reference sequences. For marine microbial eukaryotes, current reference libraries are limited by gaps in sequenced organism diversity and barriers to updating libraries with new sequence data, resulting in taxonomic annotation of about half of eukaryotic environmental transcripts. Here, we introduce Marine Functional EukaRyotic Reference Taxa (MarFERReT), a marine microbial eukaryotic sequence library designed for use with taxonomic annotation of eukaryotic metatranscriptomes. We gathered 902 publicly accessible marine eukaryote genomes and transcriptomes and assessed their sequence quality and cross-contamination issues, selecting 800 validated entries for inclusion in MarFERReT. Version 1.1 of MarFERReT contains reference sequences from 800 marine eukaryotic genomes and transcriptomes, covering 453 species- and strain-level taxa, totaling nearly 28 million protein sequences with associated NCBI and PR2 Taxonomy identifiers and Pfam functional annotations. The MarFERReT project repository hosts containerized build scripts, documentation on installation and use case examples, and information on new versions of MarFERReT.

S-adenosylhomocysteine (SAH) hydrolase is the enzyme responsible for breaking down SAH into adenosine and homocysteine. It has long been believed that a deficiency of this enzyme leads to SAH accumulation, subsequently inhibiting methyltransferases responsible for nucleic acids and proteins, which severely affects cell proliferation. To investigate whether targeting this enzyme could be a viable strategy to combat Trypanosoma brucei, the causative agent of human African trypanosomiasis, we created a null mutant of the SAH hydrolase gene in T. brucei using the Cre/loxP system and conducted a phenotype analysis. Surprisingly, the null mutant, where all five SAH hydrolase gene loci were deleted, exhibited normal proliferation despite the observed SAH accumulation. These findings suggest that inhibiting SAH hydrolase may not be an effective approach to suppressing T. brucei proliferation, making the enzyme a less promising target for antitrypanosome drug development.

Trypanosoma cruzi, the protozoan causative of Chagas disease (ChD), exhibits striking genetic and phenotypic intraspecific diversity, along with ecoepidemiological complexity. Human-pathogen interactions lead to distinct clinical presentations of ChD. In 2009, an international consensus classified T. cruzi strains into six discrete typing units (DTUs), TcI to TcVI, later including TcBat, and proposed reproducible genotyping schemes for DTU identification. This article aims to review the impact of classifying T. cruzi strains into DTUs on our understanding of biological, ecoepidemiological, and pathogenic aspects of T. cruzi. We will explore the likely origin of DTUs and the intrinsic characteristics of each group of strains concerning genome organization, genomics, and susceptibility to drugs used in ChD treatment. We will also provide an overview of the association of DTUs with mammalian reservoirs, and summarize the geographic distribution, and the clinical implications, of prevalent specific DTUs in ChD patients. Throughout this review, we will emphasize the crucial roles of both parasite and human genetics in defining ChD pathogenesis and chemotherapy outcome.

UPF1-like helicases play roles in telomeric heterochromatin formation and X-chromosome inactivation, and also in monogenic variant surface glycoprotein (VSG) expression via VSG exclusion-factor-2 (VEX2), a UPF1-related protein in the African trypanosome. We show that VEX2 associates with chromatin specifically at the single active VSG expression site on chromosome 6, forming an allele-selective connection, via VEX1, to the trans-splicing locus on chromosome 9, physically bridging two chromosomes and the VSG transcription and splicing compartments. We further show that the VEX-complex is multimeric and self-regulates turnover to tightly control its abundance. Using single cell transcriptomics following VEX2-depletion, we observed simultaneous derepression of many other telomeric VSGs and multi-allelic VSG expression in individual cells. Thus, an allele-selective, inter-chromosomal, and self-limiting VEX1-2 bridge supports monogenic VSG expression and multi-allelic VSG exclusion.

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