Short tRNA anticodon stem and mutant eRF1 allow stop codon reassignment.

Cognate tRNAs deliver specific amino acids to translating ribosomes according to the standard genetic code, and three codons with no cognate tRNAs serve as stop codons. Some protists have reassigned all stop codons as sense codons, neglecting this fundamental principle1-4. Here we analyse the in-frame stop codons in 7,259 predicted protein-coding genes of a previously undescribed trypanosomatid, Blastocrithidia nonstop. We reveal that in this species in-frame stop codons are underrepresented in genes expressed at high levels and that UAA serves as the only termination codon. Whereas new tRNAsGlu fully cognate to UAG and UAA evolved to reassign these stop codons, the UGA reassignment followed a different path through shortening the anticodon stem of tRNATrpCCA from five to four base pairs (bp). The canonical 5-bp tRNATrp recognizes UGG as dictated by the genetic code, whereas its shortened 4-bp variant incorporates tryptophan also into in-frame UGA. Mimicking this evolutionary twist by engineering both variants from B. nonstop, Trypanosoma brucei and Saccharomyces cerevisiae and expressing them in the last two species, we recorded a significantly higher readthrough for all 4-bp variants. Furthermore, a gene encoding B. nonstop release factor 1 acquired a mutation that specifically restricts UGA recognition, robustly potentiating the UGA reassignment. Virtually the same strategy has been adopted by the ciliate Condylostoma magnum. Hence, we describe a previously unknown, universal mechanism that has been exploited in unrelated eukaryotes with reassigned stop codons.

Blastocrithidia nonstop mitochondrial genome and its expression are remarkably insulated from nuclear codon reassignment.

The canonical stop codons of the nuclear genome of the trypanosomatid Blastocrithidia nonstop are recoded. Here, we investigated the effect of this recoding on the mitochondrial genome and gene expression. Trypanosomatids possess a single mitochondrion and protein-coding transcripts of this genome require RNA editing in order to generate open reading frames of many transcripts encoded as 'cryptogenes'. Small RNAs that can number in the hundreds direct editing and produce a mitochondrial transcriptome of unusual complexity. We find B. nonstop to have a typical trypanosomatid mitochondrial genetic code, which presumably requires the mitochondrion to disable utilization of the two nucleus-encoded suppressor tRNAs, which appear to be imported into the organelle. Alterations of the protein factors responsible for mRNA editing were also documented, but they have likely originated from sources other than B. nonstop nuclear genome recoding. The population of guide RNAs directing editing is minimal, yet virtually all genes for the plethora of known editing factors are still present. Most intriguingly, despite lacking complex I cryptogene guide RNAs, these cryptogene transcripts are stochastically edited to high levels.

In silico prediction of the metabolism of Blastocrithidia nonstop, a trypanosomatid with non-canonical genetic code.

BACKGROUND: Almost all extant organisms use the same, so-called canonical, genetic code with departures from it being very rare. Even more exceptional are the instances when a eukaryote with non-canonical code can be easily cultivated and has its whole genome and transcriptome sequenced. This is the case of Blastocrithidia nonstop, a trypanosomatid flagellate that reassigned all three stop codons to encode amino acids. RESULTS: We in silico predicted the metabolism of B. nonstop and compared it with that of the well-studied human parasites Trypanosoma brucei and Leishmania major. The mapped mitochondrial, glycosomal and cytosolic metabolism contains all typical features of these diverse and important parasites. We also provided experimental validation for some of the predicted observations, concerning, specifically presence of glycosomes, cellular respiration, and assembly of the respiratory complexes. CONCLUSIONS: In an unusual comparison of metabolism between a parasitic protist with a massively altered genetic code and its close relatives that rely on a canonical code we showed that the dramatic differences on the level of nucleic acids do not seem to be reflected in the metabolisms. Moreover, although the genome of B. nonstop is extremely AT-rich, we could not find any alterations of its pyrimidine synthesis pathway when compared to other trypanosomatids. Hence, we conclude that the dramatic alteration of the genetic code of B. nonstop has no significant repercussions on the metabolism of this flagellate.

Trypanosoma cruzi strain and starvation-driven mitochondrial RNA editing and transcriptome variability.

Trypanosoma cruzi is a unicellular protistan parasitic species that is comprised of strains and isolates exhibiting high levels of genetic and metabolic variability. In the insect vector, it is known to be highly responsive to starvation, a signal for progression to a life stage in which it can infect mammalian cells. Most mRNAs encoded in its mitochondrion require the targeted insertion and deletion of uridines to become translatable transcripts. This study defined differences in uridine-insertion/deletion RNA editing among three strains and established the mechanism whereby abundances of edited (and, thus, translatable) mitochondrial gene products increase during starvation. Our approach utilized our custom T-Aligner toolkit to describe transcriptome-wide editing events and reconstruct editing products from high-throughput sequencing data. We found that the relative abundance of mitochondrial transcripts and the proportion of mRNAs that are edited varies greatly between analyzed strains, a characteristic that could potentially impact metabolic capacity. Starvation typically led to an increase in overall editing activity rather than affecting a specific step in the process. We also determined that transcripts CR3, CR4, and ND3 produce multiple open reading frames that, if translated, would generate different proteins. Finally, we quantitated the inherent flexibility of editing in T. cruzi and found it to be higher relative to that in a related trypanosomatid lineage. Over time, new editing domains or patterns could prove advantageous to the organism and become more widespread within individual transcriptomes or among strains.

Phlebotomine sand fly (Diptera: Phlebotominae) diversity in the foci of cutaneous leishmaniasis in the Surxondaryo Region of Uzbekistan: 50 years on.

In Uzbekistan, the number of reported leishmaniasis cases is rising at the alarming rate. In this work, we studied the phlebotomine sand fly (Diptera: Phlebotominae) diversity in the foci of cutaneous leishmaniasis in the Surxondaryo Region of Uzbekistan and compared it with the data obtained for the same area 50 years ago, when infection prevalence was reportedly low. We found that the implicated vector for zoonotic leishmaniasis, P. papatasi, remained eudominant; the proportion of implicated anthroponotic leishmaniasis vector, P. sergenti, rose significantly from averaged 5.4 to 41.4%; Phlebotomus alexandri, a suspected visceral leishmaniasis vector, was eudominant at two sites, and a second suspected vector for this disease, P. longiductus, was newly recorded in the region. We conclude that the increase in the documented cases of cutaneous leishmaniasis in the Surxondaryo Region of Uzbekistan may be connected to the changes in fauna of sand flies vectoring Leishmania spp.

Diversity of RNA viruses in the cosmopolitan monoxenous trypanosomatid Leptomonas pyrrhocoris.

BACKGROUND: Trypanosomatids are parasitic flagellates well known because of some representatives infecting humans, domestic animals, and cultural plants. Many trypanosomatid species bear RNA viruses, which, in the case of human pathogens Leishmania spp., influence the course of the disease. One of the close relatives of leishmaniae, Leptomonas pyrrhocoris, has been previously shown to harbor viruses of the groups not documented in other trypanosomatids. At the same time, this species has a worldwide distribution and high prevalence in the natural populations of its cosmopolitan firebug host. It therefore represents an attractive model to study the diversity of RNA viruses. RESULTS: We surveyed 106 axenic cultures of L. pyrrhocoris and found that 64 (60%) of these displayed 2-12 double-stranded RNA fragments. The analysis of next-generation sequencing data revealed four viral groups with seven species, of which up to five were simultaneously detected in a single trypanosomatid isolate. Only two of these species, a tombus-like virus and an Ostravirus, were earlier documented in L. pyrrhocoris. In addition, there were four new species of Leishbuviridae, the family encompassing trypanosomatid-specific viruses, and a new species of Qinviridae, the family previously known only from metatranscriptomes of invertebrates. Currently, this is the only qinvirus with an unambiguously determined host. Our phylogenetic inferences suggest reassortment in the tombus-like virus owing to the interaction of different trypanosomatid strains. Two of the new Leishbuviridae members branch early on the phylogenetic tree of this family and display intermediate stages of genomic segment reduction between insect Phenuiviridae and crown Leishbuviridae. CONCLUSIONS: The unprecedented wide range of viruses in one protist species and the simultaneous presence of up to five viral species in a single Leptomonas pyrrhocoris isolate indicate the uniqueness of this flagellate. This is likely determined by the peculiarity of its firebug host, a highly abundant cosmopolitan species with several habits ensuring wide distribution and profuseness of L. pyrrhocoris, as well as its exposure to a wider spectrum of viruses compared to other trypanosomatids combined with a limited ability to transmit these viruses to its relatives. Thus, L. pyrrhocoris represents a suitable model to study the adoption of new viruses and their relationships with a protist host.

Properties and Mechanisms of Deletions, Insertions, and Substitutions in the Evolutionary History of SARS-CoV-2.

SARS-CoV-2 has accumulated many mutations since its emergence in late 2019. Nucleotide substitutions leading to amino acid replacements constitute the primary material for natural selection. Insertions, deletions, and substitutions appear to be critical for coronavirus's macro- and microevolution. Understanding the molecular mechanisms of mutations in the mutational hotspots (positions, loci with recurrent mutations, and nucleotide context) is important for disentangling roles of mutagenesis and selection. In the SARS-CoV-2 genome, deletions and insertions are frequently associated with repetitive sequences, whereas C>U substitutions are often surrounded by nucleotides resembling the APOBEC mutable motifs. We describe various approaches to mutation spectra analyses, including the context features of RNAs that are likely to be involved in the generation of recurrent mutations. We also discuss the interplay between mutations and natural selection as a complex evolutionary trend. The substantial variability and complexity of pipelines for the reconstruction of mutations and the huge number of genomic sequences are major problems for the analyses of mutations in the SARS-CoV-2 genome. As a solution, we advocate for the development of a centralized database of predicted mutations, which needs to be updated on a regular basis.

The 29-nucleotide deletion in SARS-CoV: truncated versions of ORF8 are under purifying selection.

BACKGROUND: Accessory proteins have diverse roles in coronavirus pathobiology. One of them in SARS-CoV (the causative agent of the severe acute respiratory syndrome outbreak in 2002-2003) is encoded by the open reading frame 8 (ORF8). Among the most dramatic genomic changes observed in SARS-CoV isolated from patients during the peak of the pandemic in 2003 was the acquisition of a characteristic 29-nucleotide deletion in ORF8. This deletion cause splitting of ORF8 into two smaller ORFs, namely ORF8a and ORF8b. Functional consequences of this event are not entirely clear. RESULTS: Here, we performed evolutionary analyses of ORF8a and ORF8b genes and documented that in both cases the frequency of synonymous mutations was greater than that of nonsynonymous ones. These results suggest that ORF8a and ORF8b are under purifying selection, thus proteins translated from these ORFs are likely to be functionally important. Comparisons with several other SARS-CoV genes revealed that another accessory gene, ORF7a, has a similar ratio of nonsynonymous to synonymous mutations suggesting that ORF8a, ORF8b, and ORF7a are under similar selection pressure. CONCLUSIONS: Our results for SARS-CoV echo the known excess of deletions in the ORF7a-ORF7b-ORF8 complex of accessory genes in SARS-CoV-2. A high frequency of deletions in this gene complex might reflect recurrent searches in "functional space" of various accessory protein combinations that may eventually produce more advantageous configurations of accessory proteins similar to the fixed deletion in the SARS-CoV ORF8 gene.

Shining the spotlight on the neglected: new high-quality genome assemblies as a gateway to understanding the evolution of Trypanosomatidae.

BACKGROUND: Protists of the family Trypanosomatidae (phylum Euglenozoa) have gained notoriety as parasites affecting humans, domestic animals, and agricultural plants. However, the true extent of the group's diversity spreads far beyond the medically and veterinary relevant species. We address several knowledge gaps in trypanosomatid research by undertaking sequencing, assembly, and analysis of genomes from previously overlooked representatives of this protistan group. RESULTS: We assembled genomes for twenty-one trypanosomatid species, with a primary focus on insect parasites and Trypanosoma spp. parasitizing non-human hosts. The assemblies exhibit sizes consistent with previously sequenced trypanosomatid genomes, ranging from approximately 18 Mb for Obscuromonas modryi to 35 Mb for Crithidia brevicula and Zelonia costaricensis. Despite being the smallest, the genome of O. modryi has the highest content of repetitive elements, contributing nearly half of its total size. Conversely, the highest proportion of unique DNA is found in the genomes of Wallacemonas spp., with repeats accounting for less than 8% of the assembly length. The majority of examined species exhibit varying degrees of aneuploidy, with trisomy being the most frequently observed condition after disomy. CONCLUSIONS: The genome of Obscuromonas modryi represents a very unusual, if not unique, example of evolution driven by two antidromous forces: i) increasing dependence on the host leading to genomic shrinkage and ii) expansion of repeats causing genome enlargement. The observed variation in somy within and between trypanosomatid genera suggests that these flagellates are largely predisposed to aneuploidy and, apparently, exploit it to gain a fitness advantage. High heterogeneity in the genome size, repeat content, and variation in chromosome copy numbers in the newly-sequenced species highlight the remarkable genome plasticity exhibited by trypanosomatid flagellates. These new genome assemblies are a robust foundation for future research on the genetic basis of life cycle changes and adaptation to different hosts in the family Trypanosomatidae.

Revisiting epidemiology of leishmaniasis in central Asia: lessons learnt.

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

Complete minicircle genome of Leptomonas pyrrhocoris reveals sources of its non-canonical mitochondrial RNA editing events.

Uridine insertion/deletion (U-indel) editing of mitochondrial mRNA, unique to the protistan class Kinetoplastea, generates canonical as well as potentially non-productive editing events. While the molecular machinery and the role of the guide (g) RNAs that provide required information for U-indel editing are well understood, little is known about the forces underlying its apparently error-prone nature. Analysis of a gRNA:mRNA pair allows the dissection of editing events in a given position of a given mitochondrial transcript. A complete gRNA dataset, paired with a fully characterized mRNA population that includes non-canonically edited transcripts, would allow such an analysis to be performed globally across the mitochondrial transcriptome. To achieve this, we have assembled 67 minicircles of the insect parasite Leptomonas pyrrhocoris, with each minicircle typically encoding one gRNA located in one of two similar-sized units of different origin. From this relatively narrow set of annotated gRNAs, we have dissected all identified mitochondrial editing events in L. pyrrhocoris, the strains of which dramatically differ in the abundance of individual minicircle classes. Our results support a model in which a multitude of editing events are driven by a limited set of gRNAs, with individual gRNAs possessing an inherent ability to guide canonical and non-canonical editing.

RNA viruses of Crithidia bombi, a parasite of bumblebees.

Leishbuviridae (Bunyavirales) are a diverse monophyletic group of negative-sense single-stranded RNA virus infecting parasitic flagellates of the family Trypanosomatidae. The presence of RNA viruses in trypanosomatids can influence the virulence of the latter. Here, we performed a screening for viruses in Crithidia bombi - a common parasite of important pollinators Bombus spp. (bumblebees) that negatively affects its host in stressful conditions. The majority (8/10) of C. bombi isolates collected in Europe and North America were positive for a virus that we named Crithidia bombi leishbuvirus 1 with high conservation of amino acid sequences between isolates. The results of our comparative phylogenetic analyses of the trypanosomatids and their viruses suggest that the high mobility of bumblebees and frequent coinfections by different strains of C. bombi determine an extensive viral exchange between the latter.

First report of putative Leishmania RNA virus 2 (LRV2) in Leishmania infantum strains from canine and human visceral leishmaniasis cases in the southeast of Brazil.

BACKGROUND: Leishmania RNA virus 1 (LRV1) is commonly found in South American Leishmania parasites belonging to the subgenus Viannia, whereas Leishmania RNA virus 2 (LRV2) was previously thought to be restricted to the Old-World pathogens of the subgenus Leishmania. OBJECTIVES: In this study, we investigated the presence of LRV2 in strains of Leishmania (L.) infantum, the causative agent of visceral leishmaniasis (VL), originating from different hosts, clinical forms, and geographical regions. METHODS: A total of seventy-one isolates were screened for LRV2 using semi-nested reverse transcription-polymerase chain reaction (RT-PCR) targeting the RNA-dependent RNA polymerase (RdRp) gene. FINDINGS: We detected LRV2 in two L. infantum isolates (CUR268 and HP-EMO) from canine and human cases, respectively. MAIN CONCLUSIONS: To the best of our knowledge, this is the first detection of LRV2 in the New World.

Parasites of firebugs in Austria with focus on the "micro"-diversity of the cosmopolitan trypanosomatid Leptomonas pyrrhocoris.

In this work, we investigated parasites of the firebug Pyrrhocoris apterus in Austria and demonstrated that in addition to the extensively studied Leptomonas pyrrhocoris, it can also be infected by Blastocrithidia sp. and by a mermithid, which for the first time has been characterized using molecular methods. This diversity can be explained by the gregarious lifestyle, as well as the coprophagous and cannibalistic behavior of the insect hosts that makes them susceptible to various parasites. In addition, we showed no tight association of the L. pyrrhocoris haplotypes and geographical locations (at least, considering the relatively small scale of locations in Austria) implying that the natural populations of L. pyrrhocoris are mixed due to the mobility of their firebug hosts.

Evolution of RNA viruses in trypanosomatids: new insights from the analysis of Sauroleishmania.

RNA viruses play an important role in Leishmania biology and virulence. Their presence was documented in three (out of four) Leishmania subgenera. Sauroleishmania of reptiles remained the only underinvestigated group. In this work, we analyzed the viral occurrence in Sauroleishmania spp. and detected RNA viruses in three out of seven isolates under study. These viruses were of two families-Narnaviridae and Totiviridae. Phylogenetic inferences demonstrated that totiviruses from L. adleri and L. tarentolae group together within a larger cluster of LRV2s, while a narnavirus of L. gymnodactyli appeared as a phylogenetic relative of narnaviruses of Blechomonas spp. Taken together, our work not only expanded the range of trypanosomatids that can host RNA viruses but also shed new light on the evolution and potential routes of viral transmission in these flagellates.

Capsid Structure of Leishmania RNA Virus 1.

Leishmania parasites cause a variety of symptoms, including mucocutaneous leishmaniasis, which results in the destruction of the mucous membranes of the nose, mouth, and throat. The species of Leishmania carrying Leishmania RNA virus 1 (LRV1), from the family Totiviridae, are more likely to cause severe disease and are less sensitive to treatment than those that do not contain the virus. Although the importance of LRV1 for the severity of leishmaniasis was discovered a long time ago, the structure of the virus remained unknown. Here, we present a cryo-electron microscopy reconstruction of the virus-like particle of LRV1 determined to a resolution of 3.65 Å. The capsid has icosahedral symmetry and is formed by 120 copies of a capsid protein assembled in asymmetric dimers. RNA genomes of viruses from the family Totiviridae are synthetized, but not capped at the 5' end, by virus RNA polymerases. To protect viral RNAs from degradation, capsid proteins of the L-A totivirus cleave the 5' caps of host mRNAs, creating decoys to overload the cellular RNA quality control system. Capsid proteins of LRV1 form positively charged clefts, which may be the cleavage sites for the 5' cap of Leishmania mRNAs. The putative RNA binding site of LRV1 is distinct from that of the related L-A virus. The structure of the LRV1 capsid enables the rational design of compounds targeting the putative decapping site. Such inhibitors may be developed into a treatment for mucocutaneous leishmaniasis caused by LRV1-positive species of LeishmaniaIMPORTANCE Twelve million people worldwide suffer from leishmaniasis, resulting in more than 30 thousand deaths annually. The disease has several variants that differ in their symptoms. The mucocutaneous form, which leads to disintegration of the nasal septum, lips, and palate, is caused predominantly by Leishmania parasites carrying Leishmania RNA virus 1 (LRV1). Here, we present the structure of the LRV1 capsid determined using cryo-electron microscopy. Capsid proteins of a related totivirus, L-A virus, protect viral RNAs from degradation by cleaving the 5' caps of host mRNAs. Capsid proteins of LRV1 may have the same function. We show that the LRV1 capsid contains positively charged clefts that may be sites for the cleavage of mRNAs of Leishmania cells. The structure of the LRV1 capsid enables the rational design of compounds targeting the putative mRNA cleavage site. Such inhibitors may be used as treatments for mucocutaneous leishmaniasis.

Differences in Charge Distribution in Leishmania tarentolae Leishmanolysin Result in a Reduced Enzymatic Activity.

Leishmania tarentolae is a non-pathogenic trypanosomatid isolated from lizards widely used for heterologous protein expression and extensively studied to understand the pathogenic mechanisms of leishmaniasis. The repertoire of leishmanolysin genes was reported to be expanded in L. tarentolae genome, but no proteolytic activity was detected. Here, we analyzed L. tarentolae leishmanolysin proteins from the genome to the structural levels and evaluated the enzymatic activity of the wild-type and overexpressing mutants of leishmanolysin. A total of 61 leishmanolysin sequences were retrieved from the L. tarentolae genome. Five of them were selected for phylogenetic analysis, and for three of them, we built 3D models based on the crystallographic structure of L. major ortholog. Molecular dynamics simulations of these models disclosed a less negative electrostatic potential compared to the template. Subsequently, L. major LmjF.10.0460 and L. tarentolae LtaP10.0650 leishmanolysins were cloned in a pLEXSY expression system into L. tarentolae. Proteins from the wild-type and the overexpressing parasites were submitted to enzymatic analysis. Our results revealed that L. tarentolae leishmanolysins harbor a weak enzymatic activity about three times less abundant than L. major leishmanolysin. Our findings strongly suggest that the less negative electrostatic potential of L. tarentolae leishmanolysin can be the reason for the reduced proteolytic activity detected in this parasite.

Differences in mitochondrial NADH dehydrogenase activities in trypanosomatids.

Complex I (NADH dehydrogenase) is the first enzyme in the respiratory chain. It catalyses the electron transfer from NADH to ubiquinone that is associated with proton pumping out of the matrix. In this study, we characterized NADH dehydrogenase activity in seven monoxenous trypanosomatid species: Blechomonas ayalai, Herpetomonas tarakana, Kentomonas sorsogonicus, Leptomonas seymouri, Novymonas esmeraldas, Sergeia podlipaevi and Wallacemonas raviniae. We also investigated the subunit composition of the complex I in dixenous Phytomonas serpens, in which its presence and activity have been previously documented. In addition to P. serpens, the complex I is functionally active in N. esmeraldas and S. podlipaevi. We also identified 24-32 subunits of the complex I in individual species by using mass spectrometry. Among them, for the first time, we recognized several proteins of the mitochondrial DNA origin.