Analyses of Leishmania-LRV Co-Phylogenetic Patterns and Evolutionary Variability of Viral Proteins

Exploring the significant genetic diversity of Iranian isolates of Leishmania RNA virus 2 using whole genome sequence analysis

BACKGROUND

Our work presents the whole genome sequence and phylogenetic analysis of five Leishmania RNA virus 2 (LRV2) isolates obtained from patients with cutaneous leishmaniasis (CL) in Iran.

METHODS

The whole genome sequencing of LRV2 was performed using a primer walking approach. The resulting sequences were analyzed for genetic and haplotype diversity, highlighting their independent evolution and significant genetic divergence.

RESULTS

The whole genome sequence of the current LRV2 showed high genetic and haplotype diversity. The study also revealed the existence of three distinct clades of LRV2, with the LRV2 sequences infecting L. major, L. aethiopica, and sauroleishmania belonging to separate lineages. These lineages have seemingly evolved independently, as the geographic distribution of their flagellate hosts does not overlap with the Leishmania species. The divergence between these three clades is attributed to considerable antiquity, leading to genetic modifications within the viruses residing in them and resulting in structural differences in their genome.

CONCLUSIONS

These findings contribute to our understanding of the genetic diversity and evolution of LRVs, providing valuable insights into their role in Leishmania infections. Further investigations are needed to understand the significance of these polymorphic sites and their potential impact on viral characteristics and disease outcomes.

Detection of Leishmania RNA Virus 2 (LRV2) among Clinical Isolates of Leishmania Major in Four Endemic Regions of Iran

PURPOSE

Leishmania RNA viruses (LRV) are double-stranded RNA viruses (dsRNA viruses) that play a role in the pathogenesis of Leishmania parasites. Cutaneous leishmaniasis (CL) is endemic in various parts of Iran. Our aimed was to investigate presence of LRV among the Leishmania major isolates in four endemic regions of Iran.

METHODS

In a cross-sectional study, we assessed the presence of LRV1 and LRV2 in 181 clinical isolates of L. major from four endemic cities in Iran using reverse transcription polymerase chain reaction (RT-PCR). After RNA extraction and cDNA synthesis, RT-PCR tests were conducted with LRV1 and LRV2 specific primers. Human beta-actin and kmp genes served as internal and external controls, respectively, and the Allele ID software was used to optimize melting curves.

RESULTS

LRV2 was detected in 27.6% (50 out of 181) of L. major isolates, while no LRV1 was found. We did not observe a statistically significant difference in the presence of LRV2 based on age group, number, or location of lesions.

CONCLUSION

This study confirms the presence of LRV2 in clinical isolates of L. major from endemic regions of Iran. Further researches with larger sample sizes is recommended to explore the association between LRV and clinical symptoms as well as treatment response.

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

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

Epidemiology and Genetic Characterization of Leishmania RNA Virus in Leishmania (Viannia) spp. Isolates from Cutaneous Leishmaniasis Endemic Areas in Panama

Leishmania (Viannia) spp. can harbor a double-stranded RNA virus known as Leishmania RNA virus 1 (LRV-1), whose presence has been reported in nine countries across the Americas and seven Leishmania species. Here, we studied 100 Leishmania (Viannia) isolates from patients with cutaneous leishmaniasis collected from different endemic areas in Panama from 2016 to 2022. We identified L. (V.) panamensis, L. (V.) guyanensis, L. (V.) braziliensis/guyanensis hybrid, and L. (V.) panamensis sp.1. (genetic variant). LRV-1 was detected by RT-PCR in 9% of L. (Viannia) isolates (eight cases in L. (V.) panamensis, and one in L. (V.) guyanensis). Phylogenetic analysis based on sequencing data classified all LRV-1 isolates within genotype A, suggesting that LRV phylogenetic proximity is closely aligned with geographical distribution or to the phylogenetic proximity of the Leishmania host in the case of the L. (V.) panamensis and L. (V.) guyanensis in Panama.

Identification of diverse RNA viruses in Obscuromonas flagellates (Euglenozoa: Trypanosomatidae: Blastocrithidiinae)

Trypanosomatids (Euglenozoa) are a diverse group of unicellular flagellates predominately infecting insects (monoxenous species) or circulating between insects and vertebrates or plants (dixenous species). Monoxenous trypanosomatids harbor a wide range of RNA viruses belonging to the families Narnaviridae, Totiviridae, Qinviridae, Leishbuviridae, and a putative group of tombus-like viruses. Here, we focus on the subfamily Blastocrithidiinae, a previously unexplored divergent group of monoxenous trypanosomatids comprising two related genera: Obscuromonas and Blastocrithidia. Members of the genus Blastocrithidia employ a unique genetic code, in which all three stop codons are repurposed to encode amino acids, with TAA also used to terminate translation. Obscuromonas isolates studied here bear viruses of three families: Narnaviridae, Qinviridae, and Mitoviridae. The latter viral group is documented in trypanosomatid flagellates for the first time. While other known mitoviruses replicate in the mitochondria, those of trypanosomatids appear to reside in the cytoplasm. Although no RNA viruses were detected in Blastocrithidia spp., we identified an endogenous viral element in the genome of B. triatomae indicating its past encounter(s) with tombus-like viruses.

Diversity and dissemination of viruses in pathogenic protozoa

Viruses are the most abundant biological entities on Earth and play a significant role in the evolution of many organisms and ecosystems. In pathogenic protozoa, the presence of viruses has been linked to an increased risk of treatment failure and severe clinical outcome. Here, we studied the molecular epidemiology of the zoonotic disease cutaneous leishmaniasis in Peru and Bolivia through a joint evolutionary analysis of Leishmania braziliensis and their dsRNA Leishmania virus 1. We show that parasite populations circulate in tropical rainforests and are associated with single viral lineages that appear in low prevalence. In contrast, groups of hybrid parasites are geographically and ecologically more dispersed and associated with an increased prevalence, diversity and spread of viruses. Our results suggest that parasite gene flow and hybridization increased the frequency of parasite-virus symbioses, a process that may change the epidemiology of leishmaniasis in the region.

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.

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.

Exploring Host-Specificity: Untangling the Relationship between Leishmania (Viannia) Species and Its Endosymbiont Leishmania RNA Virus 1

A relevant aspect in the epidemiology of Tegumentary Leishmaniasis (TL) are the Leishmania parasites carrying a viral endosymbiont, Leishmania RNA Virus 1 (LRV1), a dsRNA virus. Leishmania parasites carrying LRV1 are prone to causing more severe TL symptoms, increasing the likelihood of unfavorable clinical outcomes. LRV1 has been observed in the cultured strains of five L. (Viannia) species, and host specificity was suggested when studying the LRV1 from L. braziliensis and L. guyanensis strains. The coevolution hypothesis of LRV1 and Leishmania was based on phylogenetic analyses, implying an association between LRV1 genotypes, Leishmania species, and their geographic origins. This study aimed to investigate LRV1 specificity relative to Leishmania (Viannia) species hosts by analyzing LRV1 from L. (Viannia) species. To this end, LRV1 was screened in L. (Viannia) species other than L. braziliensis or L. guyanensis, and it was detected in 11 out of 15 L. naiffi and two out of four L. shawi. Phylogenetic analyses based on partial LRV1 genomic sequencing supported the hypothesis of host specificity, as LRV1 clustered according to their respective Leishmania species' hosts. These findings underscore the importance of investigating Leishmania and LRV1 coevolution and its impact on Leishmania (Viannia) species dispersion and pathogenesis in the American Continent.

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.

Kinetoplast Genome of Leishmania spp. Is under Strong Purifying Selection

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

Oviedo M, K. Quintero Y, L. Fernández O, Gómez MA. Leishmania RNA Virus Is Not Detected in All Species of the Leishmania Viannia Subgenus: The Case of L. (V.) panamensis in Colombia

The endosymbiotic Leishmania RNA virus 1 (LRV1) has been associated with severity and clinical manifestations of American tegumentary leishmaniasis caused by species of the Leishmania (Viannia) subgenus. Between and within Leishmania species, and among endemic countries, the prevalence of LRV is highly variable. The LRV virus has not been detected in L. (V.) panamensis, the second-most prevalent species in Central America and Colombia. However, no systematic screening of LRV has been conducted in L. (V.) panamensis, and thus it is still controversial whether this virus is truly absent from the species. We sought to determine the prevalence of LRV1 in L. (V.) panamensis clinical strains isolated from patients with cutaneous leishmaniasis (CL), from different geographic areas of Colombia. We analyzed 219 clinical strains; 78% were L. (V.) panamensis, 18% were L. (V.) braziliensis, and 4% were L. (V.) guyanensis. Screening for LRV1 was performed by quantitative reverse transcription-polymerase chain reaction. The LRV1 was detected in 18% (7 of 40) of L. (V) braziliensis strains, and was not detected in any of the L. (V.) guyanensis or L. (V.) panamensis strains. The LRV1-positive L. (V). braziliensis strains came from the Amazon Basin. Of the seven LRV1-positive strains, two were isolated from patients with mucocutaneous leishmaniasis, and the remaining from patients with CL. Our results confirm the absence of LRV1 in L. (V.) panamensis in Colombia.

Genomic analysis of Leishmania turanica strains from different regions of Central Asia

The evolution in Leishmania is governed by the opposite forces of clonality and sexual reproduction, with vicariance being an important factor. As such, Leishmania spp. populations may be monospecific or mixed. Leishmania turanica in Central Asia is a good model to compare these two types. In most areas, populations of L. turanica are mixed with L. gerbilli and L. major. Notably, co-infection with L. turanica in great gerbils helps L. major to withstand a break in the transmission cycle. Conversely, the populations of L. turanica in Mongolia are monospecific and geographically isolated. In this work, we compare genomes of several well-characterized strains of L. turanica originated from monospecific and mixed populations in Central Asia in order to shed light on genetic factors, which may drive evolution of these parasites in different settings. Our results illustrate that evolutionary differences between mixed and monospecific populations of L. turanica are not dramatic. On the level of large-scale genomic rearrangements, we confirmed that different genomic loci and different types of rearrangements may differentiate strains originated from mixed and monospecific populations, with genome translocations being the most prominent example. Our data suggests that L. turanica has a significantly higher level of chromosomal copy number variation between the strains compared to its sister species L. major with only one supernumerary chromosome. This suggests that L. turanica (in contrast to L. major) is in the active phase of evolutionary adaptation.

Multiple Regulations of Parasitic Protozoan Viruses: A Double-Edged Sword for Protozoa

Parasite infections affect human and animal health significantly and contribute to a major burden on the global economy. Parasitic protozoan viruses (PPVs) affect the protozoan parasites' morphology, phenotypes, pathogenicity, and growth rates. This discovery provides an opportunity to develop a novel preventive and therapeutic strategy for parasitic protozoan diseases (PPDs). Currently, there is greater awareness regarding PPVs; however, knowledge of viruses and their associations with host diseases remains limited. Parasite-host interactions become more complex owing to PPVs; however, few studies have investigated underlying viral regulatory mechanisms in parasites. In this study, we reviewed relevant studies to identify studies that investigated PPV development and life cycles, the triangular association between viruses, parasites, and hosts, and the effects of viruses on protozoan pathogenicity. This study highlights that viruses can alter parasite biology, and viral infection of parasites may exacerbate the adverse effects of virus-containing parasites on hosts or reduce parasite virulence. PPVs should be considered in the prevention of parasitic epidemics and outbreaks, although their effects on the host and the complexity of the triangular association between PPVs, protozoans, and hosts remain unclear. IMPORTANCE PPVs-based regulation of parasitic protozoa can provide a theoretical basis and direction for PPD prevention and control, although PPVs and PPV regulatory mechanisms remain unclear. In this review, we investigated the differences between PPVs and the unique properties of each virus regarding virus discovery, structures, and life cycles, focused on the Trichomonas vaginalis virus, Giardia lamblia virus, Leishmania RNA virus, and the Cryptosporidium parvum virus 1. The triangular association between PPVs, parasitic protozoa, and hosts reveals the "double-edged sword" property of PPVs, which maintains a balance between parasitic protozoa and hosts in both positive and negative respects. These studies discuss the complexity of parasitic protozoa and their co-existence with hosts and suggest novel pathways for using PPVs as tools to gain a deeper understanding of protozoal infection and treatment.

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.

Virion structure of Leishmania RNA virus 1

The presence of Leishmania RNA virus 1 (LRV1) enables Leishmania protozoan parasites to cause more severe disease than the virus-free strains. The structure of LRV1 virus-like particles has been determined previously, however, the structure of the LRV1 virion has not been characterized. Here we used cryo-electron microscopy and single-particle reconstruction to determine the structures of the LRV1 virion and empty particle isolated from Leishmania guyanensis to resolutions of 4.0 Å and 3.6 Å, respectively. The capsid of LRV1 is built from sixty dimers of capsid proteins organized with icosahedral symmetry. RNA genomes of totiviruses are replicated inside the virions by RNA polymerases expressed as C-terminal extensions of a sub-population of capsid proteins. Most of the virions probably contain one or two copies of the RNA polymerase, however, the location of the polymerase domains in LRV1 capsid could not be identified, indicating that it varies among particles. Importance. Every year over 200 000 people contract leishmaniasis and more than five hundred people die of the disease. The mucocutaneous form of leishmaniasis produces lesions that can destroy the mucous membranes of the nose, mouth, and throat. Leishmania parasites carrying Leishmania RNA virus 1 (LRV1) are predisposed to cause aggravated symptoms in the mucocutaneous form of leishmaniasis. Here, we present the structure of the LRV1 virion determined using cryo-electron microscopy.

Leishmania guyanensis suppressed inducible nitric oxide synthase provoked by its viral endosymbiont

Inducible nitric oxide synthase (iNOS) is essential to the production of nitric oxide (NO), an efficient effector molecule against intracellular human pathogens such as Leishmania protozoan parasites. Some strains of Leishmania are known to bear a viral endosymbiont termed Leishmania RNA virus 1 (LRV1). Recognition of LRV1 by the innate immune sensor Toll-like receptor-3 (TLR3) leads to conditions worsening the disease severity in mice. This process is governed by type I interferon (type I IFNs) arising downstream of TLR3 stimulation and favoring the formation of secondary metastatic lesions. The formation of these lesions is mediated by the inflammatory cytokine IL-17A and occurs in the absence, or low level of, protective cytokine IFN-γ. Here, we described that the presence of LRV1 led to the initial expression of iNOS and low production of NO that failed to control infection. We subsequently showed that LRV1-triggered type I IFN was essential but insufficient to induce robust iNOS induction, which requires strong activation of nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB). Leishmania guyanensis carrying LRV1 (LgyLRV1+) parasites mitigated strong iNOS production by limiting NF-kB activation via the induction of tumor necrosis factor-alpha-induced protein 3 (TNFAIP3), also known as A20. Moreover, our data suggested that production of LRV1-induced iNOS could be correlated with parasite dissemination and metastasis via elevated secretion of IL-17A in the draining lymph nodes. Our findings support an additional strategy by which LRV1-bearing Leishmania guyanensis evaded killing by nitric oxide and suggest that low levels of LRV1-induced NO might contribute to parasite metastasis.

Elimination of LRVs Elicits Different Responses in Leishmania spp

Leishmaniaviruses (LRVs) have been demonstrated to enhance progression of leishmaniasis, a vector-transmitted disease with a wide range of clinical manifestations that is caused by flagellates of the genus Leishmania. Here, we used two previously proposed strategies of the LRV ablation to shed light on the relationships of two Leishmania spp. with their respective viral species (L. guyanensis, LRV1 and L. major, LRV2) and demonstrated considerable difference between two studied systems. LRV1 could be easily eliminated by the expression of exogenous capsids regardless of their origin (the same or distantly related LRV1 strains, or even LRV2), while LRV2 was only partially depleted in the case of the native capsid overexpression. The striking differences were also observed in the effects of complete viral elimination with 2'C-methyladenosine (2-CMA) on the transcriptional profiles of these two Leishmania spp. While virtually no differentially expressed genes were detected after the LRV1 removal from L. guyanensis, the response of L. major after ablation of LRV2 involved 87 genes, the analysis of which suggested a considerable stress experienced even after several passages following the treatment. This effect on L. major was also reflected in a significant decrease of the proliferation rate, not documented in L. guyanensis and naturally virus-free strain of L. major. Our findings suggest that integration of L. major with LRV2 is deeper compared with that of L. guyanensis with LRV1. We presume this determines different effects of the viral presence on the Leishmania spp. infections.

IMPORTANCELeishmania spp. represent human pathogens that cause leishmaniasis, a widespread parasitic disease with mild to fatal clinical manifestations. Some strains of leishmaniae bear leishmaniaviruses (LRVs), and this has been shown to aggravate disease course. We investigated the relationships of two distally related Leishmania spp. with their respective LRVs using different strategies of virus removal. Our results suggest the South American L. guyanensis easily loses its virus with no important consequences for the parasite in the laboratory culture. Conversely, the Old-World L. major is refractory to virus removal and experiences a prominent stress if this removal is nonetheless completed. The drastically different levels of integration between the studied Leishmania spp. and their viruses suggest distinct effects of the viral presence on infections in these species of parasites.