In this work, we describe the first Leishmania-infecting leishbunyavirus—the first virus other than Leishmania RNA virus (LRV) found in trypanosomatid parasites. Its host is Leishmania martiniquensis, a human pathogen causing infections with a wide range of manifestations from asymptomatic to severe visceral disease. This virus (LmarLBV1) possesses many characteristic features of leishbunyaviruses, such as tripartite organization of its RNA genome, with ORFs encoding RNA-dependent RNA polymerase, surface glycoprotein, and nucleoprotein on L, M, and S segments, respectively. Our phylogenetic analyses suggest that LmarLBV1 originated from leishbunyaviruses of monoxenous trypanosomatids and, probably, is a result of genomic re-assortment. The LmarLBV1 facilitates parasites’ infectivity in vitro in primary murine macrophages model. The discovery of a virus in L. martiniquensis poses the question of whether it influences pathogenicity of this parasite in vivo, similarly to the LRV in other Leishmania species.

Maxicircles of all kinetoplastid flagellates are functional analogs of mitochondrial genome of other eukaryotes. They consist of two distinct parts, called the coding region and the divergent region (DR). The DR is composed of highly repetitive sequences and, as such, remains the least explored segment of a trypanosomatid genome. It is extremely difficult to sequence and assemble, that is why very few full length maxicircle sequences were available until now. Using PacBio data, we assembled 17 complete maxicircles from different species of trypanosomatids. Here we present their large-scale comparative analysis and describe common patterns of DR organization in trypanosomatids.

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.

Protein phosphorylation/dephosphorylation is an important regulatory mechanism that controls many key physiological processes. Numerous pathogens successfully use kinases and phosphatases to internalize, replicate, and survive, modifying the host′s phosphorylation profile or signal transduction pathways. Multiple phosphatases and kinases from diverse bacterial pathogens have been implicated in human infections before. In this work, we have identified and characterized the dual specificity protein/lipid phosphatase LmDUSP1 as a novel virulence factor governing Leishmania mexicana infection. The LmDUSP1-encoding gene (LmxM.22.0250 in L. mexicana) has been acquired from bacteria via horizontal gene transfer. Importantly, its orthologues have been associated with virulence in several bacterial species, such as Mycobacterium tuberculosis and Listeria monocytogenes. Leishmania mexicana with ablated LmxM.22.0250 demonstrated severely attenuated virulence in the experimental infection of primary mouse macrophages, suggesting that this gene facilitates Leishmania pathogenicity in vertebrates. Despite significant upregulation of LmxM.22.0250 expression in metacyclic promastigotes, its ablation did not affect the ability of mutant cells to differentiate into virulent stages in insects. It remains to be further investigated which specific biochemical pathways involve LmDUSP1 and how this facilitates the parasite′s survival in the host. One of the interesting possibilities is that LmDUSP1 may target host′s substrate(s), thereby affecting its signal transduction pathways.

Most trypanosomatid flagellates do not have catalase. In the evolution of this group, the gene encoding catalase has been independently acquired at least three times from three different bacterial groups. Here, we demonstrate that the catalase of Vickermania was obtained by horizontal gene transfer from Gammaproteobacteria, extending the list of known bacterial sources of this gene. Comparative biochemical analyses revealed that the enzymes of V. ingenoplastis, Leptomonas pyrrhocoris, and Blastocrithidia sp., representing the three independent catalase-bearing trypanosomatid lineages, have similar properties, except for the unique cyanide resistance in the catalase of the latter species.

Trypanosomatids are among the most extensively studied protists due to their parasitic interactions with insects, vertebrates, and plants. Recently, Blastocrithidia nonstop was found to depart from the canonical genetic code, with all three stop codons reassigned to encode amino acids (UAR for glutamate and UGA for tryptophan), and UAA having dual meaning also as a termination signal (glutamate and stop). To explore features linked to this phenomenon, we analyzed the genomes of four Blastocrithidia and four Obscuromonas species, the latter representing a sister group employing the canonical genetic code. We found that all Blastocrithidia species encode cognate tRNAs for UAR codons, possess a distinct 4 bp anticodon stem tRNA ^Trp _CCA decoding UGA, and utilize UAA as the only stop codon. The distribution of in-frame reassigned codons is consistently non-random, suggesting a translational burden avoided in highly expressed genes. Frame-specific enrichment of UAA codons immediately following the genuine UAA stop codon, not observed in Obscuromonas , points to a specific mode of termination. All Blastocrithidia species possess specific mutations in eukaryotic release factor 1 and a unique acidic region following the prion-like N-terminus of eukaryotic release factor 3 that may be associated with stop codon readthrough. We infer that the common ancestor of the genus Blastocrithidia already exhibited a GC-poor genome with the non-canonical genetic code. Our comparative analysis highlights features associated with this extensive stop codon reassignment. This cascade of mutually dependent adaptations, driven by increasing AU-richness in transcripts and frequent emergence of in-frame stops, underscores the dynamic interplay between genome composition and genetic code plasticity to maintain vital functionality. IMPORTANCE: The genetic code, assigning amino acids to codons, is almost universal, yet an increasing number of its alterations keep emerging, mostly in organelles and unicellular eukaryotes. One such case is the trypanosomatid genus Blastocrithidia , where all three stop codons were reassigned to amino acids, with UAA also serving as a sole termination signal. We conducted a comparative analysis of four Blastocrithidia species, all with the same non-canonical genetic code, and their close relatives of the genus Obscuromonas , which retain the canonical code. This across-genome comparison allowed the identification of key traits associated with genetic code reassignment in Blastocrithidia . This work provides insight into the evolutionary steps, facilitating an extensive departure from the canonical genetic code that occurred independently in several eukaryotic lineages.