Detection and molecular identification of leishmania RNA virus (LRV) in Iranian Leishmania species

Clinical polymorphism of zoonotic cutaneous leishmaniasis: combination of the clinical and the parasitological diagnosis

Zoonotic cutaneous leishmaniasis (ZCL) is a neglected tropical disease caused by Leishmania (L.) major. This zoonosis is characterized by a broad-spectrum clinical polymorphism and may be underestimated and poorly treated since it is a simulator of various dermatoses. The aim of our study was to analyze the clinical polymorphism of patients with ZCL. A total of 142 patients with confirmed CL based on the microscopic examination of skin lesion biopsies were included in this study. Molecular typing of Leishmania species revealed that all patients were infected with L. major. In total, 14 clinical forms were observed. Six were typical and eight were atypical. The typical ZCL forms are grouped as follows: papular (26.76%), ulcero-crusted (26.05%), ulcerated (13.38%), impetiginous (9.86%), nodular (9.15%), and papulo-nodular (5.63%) lesions. In atypical ZCL forms, we described erythematous (2.81%), erysipeloid (1.4%), sporotrichoid, (1.4%), keratotic (0.7%) lupoid (0.7%), lichenoid (0.7%), psoriasiform (0.7%), and zosteriform (0.7%) lesions. Here, the lichenoid and the keratotic forms caused by L. major were reported for the first time in Tunisia. These findings will help physicians to be aware of the unusual lesions of ZCL that could be confused with other dermatological diseases. For this reason, it will be necessary to improve the diagnosis of CL especially in endemic areas. Such large clinical polymorphism caused by L. major may be the result of a complex association between the vector microbiota, the parasite, and the host immune state, and further studies should be carried out in order to reveal the mechanisms involved in clinical polymorphism of ZCL.

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.

The effects of Leishmania RNA virus 2 (LRV2) on the virulence factors of L. major and pro-inflammatory biomarkers: an in vitro study on human monocyte cell line (THP-1)

BACKGROUND

Cutaneous Leishmaniasis (CL) is a parasitic disease with diverse outcomes. Clinical diversity is influenced by various factors such as Leishmania species and host genetic background. The role of Leishmania RNA virus (LRV), as an endosymbiont, is suggested to not only affect the pathogenesis of Leishmania, but also impact host immune responses. This study aimed to investigate the influence of LRV2 on the expression of a number of virulence factors (VFs) of Leishmania and pro-inflammatory biomarkers.

MATERIALS AND METHODS

Sample were obtained from CL patients from Golestan province. Leishmania species were identified by PCR (LIN 4, 17), and the presence of LRV2 was checked using the semi-nested PCR (RdRp gene). Human monocyte cell line (THP-1) was treated with three isolates of L. major with LRV2 and one isolate of L. major without LRV2. The treatments with four isolates were administered for the time points: zero, 12, 24, 36, and 48 h after co-infection. The expression levels of Leishmania VFs genes including GP63, HSP83, and MPI, as well as pro-inflammatory biomarkers genes including NLRP3, IL18, and IL1β, were measured using quantitative real-time PCR.

RESULTS

The expression of GP63, HSP83, and MPI revealed up-regulation in LRV2 + isolates compared to LRV2- isolates. The expression of the pro-inflammatory biomarkers including NLRP3, IL1β, and IL18 genes in LRV2- were higher than LRV2 + isolates.

CONCLUSION

This finding suggests that LRV2 + may have a probable effect on the Leishmania VFs and pro-inflammatory biomarkers in the human macrophage model.

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.

Parasite hybridization promotes spreading of endosymbiotic viruses

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 endosymbiotic 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 parasites and their endosymbiotic Leishmania RNA virus. We show that parasite populations circulate in isolated pockets of suitable habitat and are associated with single viral lineages that appear in low prevalence. In contrast, groups of hybrid parasites were geographically and ecologically dispersed, and commonly infected from a pool of genetically diverse viruses. Our results suggest that parasite hybridization, likely due to increased human migration and ecological perturbations, increased the frequency of endosymbiotic interactions known to play a key role in disease severity.

Effect of Leishmania RNA virus 2 on virulence factors and cytokines gene expression in a human macrophage infected with Leishmania major: A preliminary study

Cutaneous leishmaniasis (CL) is one of the most important infectious parasitic diseases in the world caused by the Leishmania parasite. In recent decades, the presence of a virus from the Totiviridae family has been proven in some Leishmania species. Although the existence of LRV2 in the Old world Leishmania species has been confirmed, almost no studies have been done to determine the potential impact of LRV2 on the immunopathogenicity of the Leishmania parasite. In this preliminary study, we measured the expression of target genes, including Glycoprotein 63 (gp63), Heat Shock Protein 70 (hsp70), Cysteine Protease b (cpb), Interleukin 1 beta (IL-1β), IL8 and IL-12 in LRV2 positive Leishmania major strain (LRV2⁺L. major) and LRV2 negative L. major strain (LRV2^-L. major). We exposed THP-1, a human leukemia monocytic cell line, to promastigotes of both strains. After the initial infection, RNA was extracted at different time points, and the relative gene expression was determined using a real-time quantitative reverse transcription-polymerase chain reaction (qRT-PCR). Findings showed that the presence of LRV2 in L. major was able to increase the expression of gp63, hsp70, and cpb genes; also, we observed lower levels of expression in cytokine genes of IL-1β, IL-8, IL-12 in the presence of LRV2⁺, which are critical factors in the host's immune response against leishmaniasis. These changes could suggest that the presence of LRV2 in L. major parasite may change the outcome of the disease and increase the probability of Leishmania survival; nevertheless, further studies are needed to confirm our results.

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.

Identification of Leishmania species and frequency distribution of LRV1 and LRV2 viruses on cutaneous leishmaniasis patients in Isfahan Province, Iran

PURPOSE

Leishmaniasis is one of the most serious health problems in developing countries. Iran is one of the endemic regions of cutaneous leishmaniasis. Leishmania RNA virus (LRV) is a dsRNA virus member of the Totiviridae family, which was first detected in the promastigotes of Leishmania braziliensis guyanensis. Our study aimed to investigate possible changes in the predominant and causative strains of CL and screening the LRV1 and LRV2 species genome from Leishmania species isolated from the lesions of patients.

MATERIALS AND METHODS

Direct smear samples obtained from 62 patients with leishmaniasis referring to the Skin Diseases and Leishmaniasis Research Center in Isfahan province during 2021-2022 were examined. Total DNA extraction procedures and conservation of site-specific multiplex PCR and nested PCR were performed for detecting Leishmania species. The molecular identification of LRV1 and LRV2 viruses, samples were used for total RNA extraction and real-time (RT)-PCR analysis, followed by conducting a restriction enzyme assay to confirm the PCR products.

RESULTS

Of the total Leishmania isolates, 54 and 8 isolates were identified as L. major and L. tropica, respectively. LRV2 was identified in 18 samples affected by L. major, while LRV1 was only detected in one of the samples with L. tropica. No LRV2 was found in any samples with L. tropica. The results showed that there was a significant relationship between LRV1 and the type of leishmaniasis (Sig. ​= ​0.009, P ​≤ ​0.05), while this relationship was not observed between LRV2 and the type of leishmaniasis.

CONCLUSIONS

The presence of a significant number of LRV2 in isolated samples, as well as the recognition of LRV1 in one of the Old World leishmaniasis species, which is a new result, could pave the way for investigating further aspects of this disease and successful treatment strategies in future studies.

First report of Leishmania RNA virus 1 in Leishmania (Viannia) braziliensis clinical isolates from Rio de Janeiro State - Brazil

BACKGROUND

Leishmania parasites carry a double-stranded RNA virus (Leishmania RNA virus - LRV) that has been divided in LRV1 and LRV2.

OBJECTIVES

Leishmania (Viannia) braziliensis clinical isolates were assessed in order to determine LRV presence.

METHODS

Two-round polymerase chain reaction (PCR and nested PCR) was performed to detect LRV1 or LRV2 in L. (V.) braziliensis clinical isolates (n = 12).

FINDINGS

LRV1 was detected in three clinical isolates which was phylogenetically related to other sequences reported from other American tegumentary leishmaniasis (ATL) endemic areas of Brazil. Patients infected with L. (V.) braziliensis LRV-negative showed only cutaneous lesions while LRV-positive reported different manifestations.

MAIN CONCLUSION

Data presented here show for the first time that LRV1 is circulating in L. (V.) braziliensis clinical isolates from Rio de Janeiro State in Brazil.

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.

Leishmania RNA virus 2 (LRV2) exacerbates dermal lesions caused by Leishmania major and comparatively unresponsive to meglumine antimoniate treatment

PURPOSE

The present study investigated the possible role of Leishmania RNA virus 2 (LRV2) in the severity of dermal lesions and treatment failure due to Leishmania major.

METHODS

The drug susceptibility of 14 clinical isolates of L.major, including resistant (n = 7) and sensitive (n = 7) isolates, was checked in the J774A.1 macrophage cell line. The presence of LRV2 among isolates was investigated by the RdRp gene and semi-nested PCR. Moreover, 1 × 10⁶ sensitive L. major LRV2⁺ and LRV2^- promastigotes were inoculated subcutaneously into the base tails of the 40 BALB/c mice divided into 4 groups (n = 10 in each group), including clinical LRV2⁺, clinical LRV2^-, positive control LRV2⁺ and negative control LRV2^-. The groups were infected with a unique isolate. The lesion size and parasite burden were evaluated.

RESULTS

Sensitive and resistant isolates were determined by the drug susceptibility method. A higher presence of LRV2 was observed among MA-resistant isolates (6/7) compared with susceptible isolates (4/7), which was not statistically significant (P = 0.237). On the other hand, a comparison of the lesion sizes between the LRV2⁺ and LRV2^- BALB/c mice groups revealed that the mean size of the lesion in the LRV2⁺ groups was significantly higher than the LRV2^- (P = 0.034). In the same direction, there was an increased parasite burden in mice inoculated with LRV2⁺ groups compared with the LRV2- BALB/c mice groups (P = 0.002).

CONCLUSIONS

Our findings showed that the presence of LRV2 could be one of the factors contributing to exacerbating CL. Although we found a higher presence of LRV2 in the resistant isolates, it seems that further investigations are recommended to determine the detailed association between lesions' aggravation and being comparatively unresponsive to treatment.

Successful Isolation of Leishmania RNA Virus (LRV) from Leishmania major in a Cutaneous Leishmaniasis Focus in Central Iran: An Update on Cases

Purpose

Cutaneous leishmaniasis (CL) is a major vector-borne disease that affects people globally, including Iran. Different factors are associated with leishmaniasis pathogenicity; recently, a link of the possible relationship between Leishmania RNA Virus (LRV) and disease severity was proposed, especially in the New World leishmaniasis (NWL). This study was aimed to investigate the presence of LRV2 in Leishmania isolates in Aran o Bidgol, Isfahan province.

Methods

Samples were collected from 110 CL-suspected patients referred to the health center. In this study, we aimed to investigate CL cases (parasitologically and clinically), identify Leishmania species (by ITS1-PCR–RFLP), and finally detection of LRV2 (by RdRp-semi-nested PCR).

Results

Parasitological methods showed 60 positive cases, based on the HaeIII enzyme restriction profile, 59 cases were caused by L. major and 1 case by L. tropica. Our project is the first study on LRV2 isolation in Aran o Bidgol city and the LRV was successfully detected from a single L. major isolated in a women’s hand lesion. Using BLAST, 94.8–100% similarity was observed in the RdRp sequence of current LRV isolate with those available in GenBank from Iran or overseas.

Conclusion

L. major was the main cause of CL in Aran o Bidgol, although L. tropica is also present in a much lower proportion in the area. This is the first report on the presence of LRV2 in Aran o Bidgol and the fifth in Iran.

Increasing the Sensitivity of Leishmania RNA Virus 2 (LRV2) Detection with a Modification in cDNA Synthesis

OBJECTIVE

Leishmania RNA virus was detected the first time in the New World Leishmania species. Recent studies were also showed the presence of Leishmania RNA virus 2 (LRV2) in Old Word Leishmania species including Turkish L. major and L. tropica isolates. This study aimed to increase the sensitivity of qPCR with a modification in the denaturation step of cDNA preparation protocol.

METHODS

In this study, LRV2+ three L. major, two L. tropica strains and L. major control strain (MHOM/SU/73/5-ASKH) were included. Total RNA isolation was done using different numbers of Leishmania promastigotes (10⁸, 10⁵ and 10³). Before cDNA synthesis, samples were denatured at 95 °C for 2 min, as a modification of the kit procedure. qPCR was undertaken using 0.5 mM primers (LRV F-HR/LRV R-HR) diluted in SYBR Green Master mix.

RESULTS

We observed lower Ct values in amplicons with the modified version than with the classical kit protocol for cDNA synthesis, in all of the strains used in the study. The addition of pre-denaturation step at 95 °C showed lower Ct values meaning the sensitivity increased. Different parasite dilutions showed similar results.

CONCLUSION

It is important to increase the sensitivity especially with the aim for detecting LRV in clinical samples obtained from patients probably have less number of parasites. The presence and burden of the virus can help to understand the relationship between the clinical findings and the pathogenicity of the parasite which may lead to changes in the course of treatment.

The Dangerous Liaisons in the Oxidative Stress Response to Leishmania Infection

Leishmania parasites preferentially invade macrophages, the professional phagocytic cells, at the site of infection. Macrophages play conflicting roles in Leishmania infection either by the destruction of internalized parasites or by providing a safe shelter for parasite replication. In response to invading pathogens, however, macrophages induce an oxidative burst as a mechanism of defense to promote pathogen removal and contribute to signaling pathways involving inflammation and the immune response. Thus, oxidative stress plays a dual role in infection whereby free radicals protect against invading pathogens but can also cause inflammation resulting in tissue damage. The induced oxidative stress in parasitic infections triggers the activation in the host of the antioxidant response to counteract the damaging oxidative burst. Consequently, macrophages are crucial for disease progression or control. The ultimate outcome depends on dangerous liaisons between the infecting Leishmania spp. and the type and strength of the host immune response.

Molecular identification of Leishmania RNA virus in cutaneous leishmaniasis patients and rodent reservoirs in Isfahan province, Iran

Leishmania RNA virus (LRV) is a double-strand RNA virus that was first detected in members of the Leishmania viannia in the New World. The present study aimed to investigate the presence of LRV in the Leishmania species isolated from cutaneous leishmaniasis (CL) patients and rodents as reservoirs in Isfahan province an old zoonotic CL focus, center of Iran. Totally, 85 samples were collected from CL patients (n = 80) and rodent reservoirs (n = 5) from different regions of Isfahan province. Species identification was determined using the PCR-RFLP method. Viral dsRNA was extracted and for observation of 5.3 kb dsRNA on an agarose gel. The presence of LRV was surveyed using the Semi-nested PCR method. For phylogenetic analyzes, 6 samples of 13 isolates were sequenced and a phylogenetic tree was drawn by MEGA7 version 7.0.26. Of 80 Leishmania isolates recovered from the patients with CL, 79 and only one were identified as L. major and L. tropica, respectively. Also, the PCR assays detected four L. major and one L. turanica in five assessed Rhombomys opimus as the rodent reservoirs. LRV was detected only in Leishmania species isolated from 13 species of 85 (15.3%) CL including (L. major, n = 12) and (L. tropica, n = 1). Phylogenetic analysis showed that they were belonged to LRV2 and had the highest similarity with Iranian reference LRV2 in GenBank. Our results showed that the LRV2 was present in cutaneous Leishmania species in Isfahan province is the most historical and touristic province of Iran. In the study LRV was not reported from rodent reservoirs, it may be due to the small sample size. Phylogenetic analysis of current sequences demonstrated that these isolates belong to the registered LRV2 of the Old World.

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

Leishmania spp. are important pathogens causing a vector-borne disease with a broad range of clinical manifestations from self-healing ulcers to the life-threatening visceral forms. Presence of Leishmania RNA virus (LRV) confers survival advantage to these parasites by suppressing anti-leishmanial immunity in the vertebrate host. The two viral species, LRV1 and LRV2 infect species of the subgenera Viannia and Leishmania, respectively. In this work we investigated co-phylogenetic patterns of leishmaniae and their viruses on a small scale (LRV2 in L. major) and demonstrated their predominant coevolution, occasionally broken by intraspecific host switches. Our analysis of the two viral genes, encoding the capsid and RNA-dependent RNA polymerase (RDRP), revealed them to be under the pressure of purifying selection, which was considerably stronger for the former gene across the whole tree. The selective pressure also differs between the LRV clades and correlates with the frequency of interspecific host switches. In addition, using experimental (capsid) and predicted (RDRP) models we demonstrated that the evolutionary variability across the structure is strikingly different in these two viral proteins.

Molecular characterization and genetic diversity of cutaneous leishmaniasis from North Eastern Pakistan

The genetic diversity of Leishmania spp. in North Eastern Pakistan remains undetermined despite increased cases of cutaneous leishmaniasis (CL). This study was designed to decipher the molecular characterization and genetic diversity of Leishmania spp. in North Eastern Pakistan. Out of 13761 CL suspected cases, 567 cases were microscopically positive and confirmed as Leishmania spp. by internal transcribed spacer (ITS) gene amplification through the PCR- RFLP technique. Further, isolates were directly sequenced to conduct phylogenetic analysis for genetic diversity. Among suspected CL cases, Mirpur showed the highest proportion of CL infection with 4.85% (118/2431) of the cases, while the Neelum district showed the lowest percentage at 3.29% (9/273). The slide positivity rate, annual blood examination rate, and annual parasite incidence rate were 3.84, 0.27, and 0.01% respectively, and the incidence of CL in the age group 1-20 years old was higher in males (50.92%) than females (25.75%). The RFLP analysis and sequencing confirmed the occurrence of Leishmania tropica, Leishmania major, and Leishmania infantum. Leishmania tropica (p = 0.02) confirmed significantly higher nucleotides variation than L. major (p = 0.05). Current findings confirmed the prior assumption that anthroponotic CL is the primary CL form present in North Eastern Pakistan. Moreover, this is the first report based on molecular identification of L. major, and L. infantum from North Eastern Pakistan. This remarkable heterogeneity in the Leishmania spp. is the leading cause of treatment failure and emergence of new haplotypes. Therefore more extensive investigations are recommended from all geographical regions of North Eastern Pakistan, especially those using a large sample size.

The Geographical Distribution of Human Cutaneous and Visceral Leishmania Species Identified by Molecular Methods in Iran: A Systematic Review With Meta-Analysis

Leishmaniasis is one of the most common vector-borne parasitic diseases in Iran. Leishmania species identification is necessary for epidemiological aspects, precise prognosis, control and treatment of the disease. We systematically searched all the studies, reports, and documentation related to species identification and geographical distribution of causative agents of cutaneous (CL), mucosal (ML), and visceral leishmaniasis (VL) using DNA-based molecular diagnostic techniques in Iran. International databases including PubMed, ScienceDirect, Embase, Google Scholar, Scopus, and Web of Science were systemically searched for English articles and Iran's databases including SID, IranMedex and Magiran were searched for Persian reports and articles. Searches were performed from 1999 to 2019 (20 years). The current review was conducted using the keywords: cutaneous leishmaniasis, visceral leishmaniasis, Leishmania species, Human, Molecular, PCR, and Iran. The study quality was evaluated using the NOS checklist. This meta-analysis procedure was accomplished using STATA, version 2.7.9. Of the 3,426 records identified in the initial search, 154 articles met inclusion criteria and qualified for the systematic review and meta-analysis. In subgroup analysis, the pooled frequency of causative agents of CL isolates was 67.3% (95% CI: 59.51–74.67%) for L. major and 32.1% (95% CI: 24.72–39.87%) for L. tropica. In addition, the pooled frequency of causative agents of VL isolates was 97.1% (95% CI: 94.6–98.8%) for L. infantum and 2.9% (95% CI: 1.12–5.37%) for L. tropica. The findings of this study showed that the main causative agents of CL and VL in Iran are L. major and L. infantum, respectively. Moreover, kinetoplast DNA (kDNA) and internal transcriber spacer (ITS) were the most used markers for identifying Leishmania species. The current study provides valuable data to encourage and direct researchers as well as public health managers in the comprehensive leishmaniasis control and prevention planning in Iran.

Determinants of Unresponsiveness to Treatment in Cutaneous Leishmaniasis: A Focus on Anthroponotic Form Due to Leishmania tropica

Cutaneous leishmaniasis (CL) is a curable disease; however, due to various risk factors, unresponsiveness to CL treatments is inevitable. The treatment of CL has been firmly correlated with multiple determinants, such as demographical, clinical, and environmental factors, the host’s immune response, poor treatment adherence, the parasite’s genetic make-up, and Leishmania RNA virus. This study primarily focuses on the risk factors associated with different therapeutic outcomes following meglumine antimoniate (MA; Glucantime^®) treatment and policy approaches to prevent unresponsiveness in CL patients with a focus on anthroponotic form (ACL). Findings suggest that effective preventive and therapeutic measures should be more vigorously implemented, particularly in endemic areas. Accordingly, extensive training is essential to monitor drug unresponsiveness regularly, especially in tropical regions where the disease is prevalent. Since humans are the fundamental reservoir host of ACL due to L. tropica, prompt detection, early diagnosis, and timely and effective treatment could help control this disease. Furthermore, major challenges and gaps remain: efficacious vaccine, new tools, and expert staff are crucial before CL can be definitively controlled.

The Maze Pathway of Coevolution: A Critical Review over the Leishmania and Its Endosymbiotic History

The description of the genus Leishmania as the causative agent of leishmaniasis occurred in the modern age. However, evolutionary studies suggest that the origin of Leishmania can be traced back to the Mesozoic era. Subsequently, during its evolutionary process, it achieved worldwide dispersion predating the breakup of the Gondwana supercontinent. It is assumed that this parasite evolved from monoxenic Trypanosomatidae. Phylogenetic studies locate dixenous Leishmania in a well-supported clade, in the recently named subfamily Leishmaniinae, which also includes monoxenous trypanosomatids. Virus-like particles have been reported in many species of this family. To date, several Leishmania species have been reported to be infected by Leishmania RNA virus (LRV) and Leishbunyavirus (LBV). Since the first descriptions of LRVs decades ago, differences in their genomic structures have been highlighted, leading to the designation of LRV1 in L. (Viannia) species and LRV2 in L. (Leishmania) species. There are strong indications that viruses that infect Leishmania spp. have the ability to enhance parasitic survival in humans as well as in experimental infections, through highly complex and specialized mechanisms. Phylogenetic analyses of these viruses have shown that their genomic differences correlate with the parasite species infected, suggesting a coevolutionary process. Herein, we will explore what has been described in the literature regarding the relationship between Leishmania and endosymbiotic Leishmania viruses and what is known about this association that could contribute to discussions about the worldwide dispersion of Leishmania.

Characterization of a new Leishmania major strain for use in a controlled human infection model

Leishmaniasis is widely regarded as a vaccine-preventable disease, but the costs required to reach pivotal Phase 3 studies and uncertainty about which candidate vaccines should be progressed into human studies significantly limits progress in vaccine development for this neglected tropical disease. Controlled human infection models (CHIMs) provide a pathway for accelerating vaccine development and to more fully understand disease pathogenesis and correlates of protection. Here, we describe the isolation, characterization and GMP manufacture of a new clinical strain of Leishmania major. Two fresh strains of L. major from Israel were initially compared by genome sequencing, in vivo infectivity and drug sensitivity in mice, and development and transmission competence in sand flies, allowing one to be selected for GMP production. This study addresses a major roadblock in the development of vaccines for leishmaniasis, providing a key resource for CHIM studies of sand fly transmitted cutaneous leishmaniasis.

Detection of Leishmania RNA Virus in Clinical Samples from Cutaneous Leishmaniasis Patients Varies according to the Type of Sample

Leishmania RNA virus (LRV) is a double-stranded RNA virus belonging to the Totiviridae family detected as cytoplasmic inclusions in some strains of the human parasite Leishmania spp. Experimental evidence supports the hypothesis that human coinfection with Leishmania spp.-LRV triggers an exacerbated immune response in the host that can be responsible for the observed complicated outcomes in cutaneous leishmaniasis (CL), such as mucosal leishmaniasis (ML) and treatment failure of CL. However, the reported frequencies of LRV associated with complicated outcomes in patient's series are highly variable, diminishing the relevance on the virus presence in the pathogenesis of the disease. To assess whether or not the inconsistent information about the frequency of LRV associated with CL complicated outcomes could be related to the virus detection approach, the present study evaluated the LRV presence in clinical samples using a diagnostic algorithm according to the type of the sample. In 36 samples with diagnosis of complicated forms of CL (15 of ML and 21 of CL antimony treatment failure) and six samples with non-Leishmania spp. infection, the LRV presence was assessed by RT-PCR, RT-qPCR, and nested RT-PCR. Viral load was estimated in parasite clinical isolates. By combining the methods, LRV1 presence was confirmed in 45% (9/20) of isolates and 37.5% (6/16) of the incisional biopsies. Remarkably, in some cases (4/8), LRV1 was undetectable in the isolates but present in their respective biopsies, and less frequently, the opposite was observed (1/8), suggesting the possibility of loss of parasites harboring LRV1 during the in vitro growth.

Presence and diversity of Leishmania RNA virus in an old zoonotic cutaneous leishmaniasis focus, northeastern Iran: haplotype and phylogenetic based approach

OBJECTIVE

Leishmania RNA virus (LRV) is a double-stranded RNA (dsRNA) virus that circulates within many species of the Leishmania parasite. In this study, we aimed to investigate the presence of LRV2 circulating in Leishmania isolates in an old focus of ZCL located in northeastern of Iran.

METHODS

Leishmania isolates were collected from 85 patients that confirmed to have cutaneous leishmaniasis (CL) based on parasitological examination. To identify the Leishmania isolates, species-specific primer sets were applied for molecular identification. The presence of LRV2 was performed by RdRp-semi nested-PCR. The genetic diversity were calculated using MEGA and DnaSP. To assess haplotype diversity, 31 LRV2 strains in different regions were surveyed using analysis a 292-bp section of the RdRp sequences.

RESULTS

Out of 85 patients, 83 (97.6 %) were diagnosed with L. major and 2 (2.4 %) with L. tropica. LRV2 virus was detected in 59 (69.4%) of the CL cases. For the first time, LRV2 was reported in one L. tropica strain in Iran. The current LRV2 sequences indicated the highest similarities to an Old World LRV2. Moreover, 10 unique haplotypes were identified based on the analyzed sequences of the RdRp gene.

CONCLUSIONS

Our results indicated the highest occurrence of Leishmania/LRV2 co-circulation in this known ZCL focus from northeastern Iran. Phylogenetic analyses of LRV2 sequences confirmed that these isolates belong to the order of LRV2 from the Old World. This study offered an insight into LRV2 haplotype that the informative issue can be used for genetic research of LRV2 in other regions.

Interferon Alpha Favors Macrophage Infection by Visceral Leishmania Species Through Upregulation of Sialoadhesin Expression

Type I interferons (IFNs) induced by an endogenous Leishmania RNA virus or exogenous viral infections have been shown to exacerbate infections with New World Cutaneous Leishmania parasites, however, the impact of type I IFNs in visceral Leishmania infections and implicated mechanisms remain to be unraveled. This study assessed the impact of type I IFN on macrophage infection with L. infantum and L. donovani and the implication of sialoadhesin (Siglec-1/CD169, Sn) as an IFN-inducible surface receptor. Stimulation of bone marrow-derived macrophages with type I IFN (IFN-α) significantly enhanced susceptibility to infection of reference laboratory strains and a set of recent clinical isolates. IFN-α particularly enhanced promastigote uptake. Enhanced macrophage susceptibility was linked to upregulated Sn surface expression as a major contributing factor to the infection exacerbating effect of IFN-α. Stimulation experiments in Sn-deficient macrophages, macrophage pretreatment with a monoclonal anti-Sn antibody or a novel bivalent anti-Sn nanobody and blocking of parasites with soluble Sn restored normal susceptibility levels. Infection of Sn-deficient mice with bioluminescent L. infantum promastigotes revealed a moderate, strain-dependent role for Sn during visceral infection under the used experimental conditions. These data indicate that IFN-responsive Sn expression can enhance the susceptibility of macrophages to infection with visceral Leishmania promastigotes and that targeting of Sn may have some protective effects in early infection.

Relationship of Leishmania RNA Virus (LRV) and treatment failure in clinical isolates of Leishmania major

OBJECTIVE

Leishmaniasis is caused by different Leishmania spp. Treatment failure (TF) of cutaneous leishmaniasis (CL) is a serious issue that may be due to various reasons, previous studies suggested Leishmania RNA virus (LRV) as a potential cause of TF. Two variant groups of LRV1 and LRV2 are reported. In this study, the presence of LRV1/LRV2 was compared in TF with treatment response (TR) isolates of L. major. Clinical isolates of 15 TF and 15 TR were collected from CL patients referred to the Health Centers of Isfahan. Genomic DNA was extracted to identify Leishmania spp. using ITS1-PCR-RFLP. Identification of LRV1/LRV2 was performed using SYBR Green Real-Time PCR. The statistical analysis to test relationship between the treatment response with Glucantime and the presence of LRV were performed using SPSS 16.0 with Fisher's Exact test. P value of less than 0.05 was considered significant.

RESULTS

ITS1-PCR-RFLP results showed that every isolate was identified as L. major. The results showed no LRV1 in any of the samples but 7 TR isolates and 2 TF isolates showed positive for LRV2. Statistical analysis showed no significant difference between the presence of LRV2 and response to Glucantime (p-value = 0.1086). Therefore, other mechanisms might be responsible for TF.

The First Non-LRV RNA Virus in Leishmania

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.

Molecular Characterization of Leishmania RNA virus 2 in Leishmania major from Uzbekistan

Here we report sequence and phylogenetic analysis of two new isolates of Leishmania RNA virus 2 (LRV2) found in Leishmania major isolated from human patients with cutaneous leishmaniasis in south Uzbekistan. These new virus-infected flagellates were isolated in the same region of Uzbekistan and the viral sequences differed by only nineteen SNPs, all except one being silent mutations. Therefore, we concluded that they belong to a single LRV2 species. New viruses are closely related to the LRV2-Lmj-ASKH documented in Turkmenistan in 1995, which is congruent with their shared host (L. major) and common geographical origin.

Transcriptional Analysis of Human Skin Lesions Identifies Tryptophan-2,3-Deoxygenase as a Restriction Factor for Cutaneous Leishmania

Disease manifestation after infection with cutaneous Leishmania species is the result of a complex interplay of diverse factors, including the immune status of the host, the infecting parasite species, or the parasite load at the lesion site. Understanding how these factors impact on the pathology of cutaneous leishmaniasis (CL) may provide new targets to manage the infection and improve clinical outcome. We quantified the relative expression of 170 genes involved in a diverse range of biological processes, in the skin biopsies from patients afflicted with CL caused by infection with either L. major or L. tropica. As compared to healthy skin, CL lesions bear elevated levels of transcripts involved in the immune response, and conversely, present a significant downregulation in the expression of genes involved in epidermal integrity and arginine or fatty acid metabolism. The expression of transcripts encoding for cytotoxic mediators and chemokines in lesions was inversely correlated with the expression of genes involved in epidermal integrity, suggesting that cytotoxicity is a major mediator of CL pathology. When comparing the transcriptional profiles of lesions caused by either L. major or L. tropica, we found them to be very similar, the later presenting an aggravated inflammatory/cytotoxic profile. Finally, we identified genes positively correlated with the parasite load in lesions. Among others, these included Th2 or regulatory cytokines, such as IL4 or IL10. Remarkably, a single gene among our dataset, encoding for tryptophan-2,3-deoxygenase (TDO), presented a negative correlation with the parasite load, suggesting that its expression may restrict parasite numbers in lesions. In agreement, treatment of macrophages infected with L. major in vitro with a TDO inhibitor led to an increase in parasite transcripts. Our work provides new insights into the factors that impact CL pathology and identifies TDO as a restriction factor for cutaneous Leishmania.

Global status of synchronizing Leishmania RNA virus in Leishmania parasites: A systematic review with meta-analysis

Leishmaniasis is one of the most neglected tropical diseases caused by protozoan parasites belonging to the genus Leishmania. There is much evidence regarding prevalence of Leishmania RNAvirus (LRV) causing Old World leishmaniasis (OWL) and New World leishmaniasis (NWL); however, a combined evidence-based knowledge on this topic is not still available. The purpose of this systematic review and meta-analysis was to address the global status of synchronizing LRV in Leishmania in the available literature. The data were systematically collected from the English electronic databases up to May 2018. Then, the studies were screened based on the inclusion and exclusion criteria. The random-effect model was used by forest plot with 95% confidence interval (CI). Overall, 877 samples from 17 articles were included in this study. Given species of Leishmania, the highest prevalence of LRV belonged to Leishmania (L.) Viannia (V.) guyanensis and L. V. braziliensis. Additionally, the virus was detected also in L. V. amazonensis, L. V. panamanensis, L. V. lainsoni, L. aethiopica, L. major and L. infantum. By random-effect model, the global prevalence of LRV was estimated to be 26.2% (95% CI: 14.4% - 40.1%). The high prevalence of LRV among causative agents of NWLisolated from the metastatic clinical forms suggests potential association of LRV with metastatic clinical forms in New World endemic regions. A comprehensive investigation on experimental and clinical aspects of LRV is needed to fully appraise the role of these viruses in pathogenicity of Leishmania parasites and their drug resistance.

RNA Viruses in Blechomonas (Trypanosomatidae) and Evolution of Leishmaniavirus

In this work, we analyzed viral prevalence in trypanosomatid parasites (Blechomonas spp.) infecting Siphonaptera and discovered nine species of viruses from three different groups (leishbunyaviruses, narnaviruses, and leishmaniaviruses). Most of the flagellate isolates bore two or three viral types (mixed infections). Although no new viral groups were documented in Blechomonas spp., our findings are important for the comprehension of viral evolution. The discovery of bunyaviruses in blechomonads was anticipated, since these viruses have envelopes facilitating their interspecific transmission and have already been found in various trypanosomatids and metatranscriptomes with trypanosomatid signatures. In this work, we also provided evidence that even representatives of the family Narnaviridae are capable of host switching and evidently have accomplished switches multiple times in the course of their evolution. The most unexpected finding was the presence of leishmaniaviruses, a group previously solely confined to the human pathogens Leishmania spp. From phylogenetic inferences and analyses of the life cycles of Leishmania and Blechomonas, we concluded that a common ancestor of leishmaniaviruses most likely infected Leishmania first and was acquired by Blechomonas by horizontal transfer. Our findings demonstrate that evolution of leishmaniaviruses is more complex than previously thought and includes occasional host switching.IMPORTANCE Flagellates belonging to the genus Leishmania are important human parasites. Some strains of different Leishmania species harbor viruses (leishmaniaviruses), which facilitate metastatic spread of the parasites, thus aggravating the disease. Up until now, these viruses were known to be hosted only by Leishmania Here, we analyzed viral distribution in Blechomonas, a related group of flagellates parasitizing fleas, and revealed that they also bear leishmaniaviruses. Our findings shed light on the entangled evolution of these viruses. In addition, we documented that Blechomonas can be also infected by leishbunyaviruses and narnaviruses, viral groups known from other insects' flagellates.

The criminal association of Leishmania parasites and viruses

In nature, humans infected with protozoan parasites can encounter viruses, which could alter their host immune response. The impact of viruses on human parasitic diseases remains largely unexplored due to the highly sterilized environment in experimental studies and the difficulty to draw a correlation between co-infection and pathology. Recent studies show that viral infections exacerbate pathology and promote dissemination of some Leishmania infections, based on a hyper-inflammatory reaction driven by type I interferons. Thus, not only the infecting parasite species, but also bystander viral infections could be a major determinant of the outcome of Leishmania infection. In this review, we focus on the contribution of viral co-infection to the exacerbation of leishmaniasis's pathology and its possible impact on treatment and vaccination strategies.

New insights into the genetic diversity of Leishmania RNA Virus 1 and its species-specific relationship with Leishmania parasites

Cutaneous leishmaniasis is a neglected parasitic disease that manifests in infected individuals under different phenotypes, with a range of factors contributing to its broad clinical spectrum. One factor, Leishmania RNA Virus 1 (LRV1), has been described as an endosymbiont present in different species of Leishmania. LRV1 significantly worsens the lesion, exacerbating the immune response in both experimentally infected animals and infected individuals. Little is known about the composition and genetic diversity of these viruses. Here, we investigated the relationship between the genetic composition of LRV1 detected in strains of Leishmania (Viannia) braziliensis and L. (V.) guyanensis and the interaction between the endosymbiont and the parasitic species, analyzing an approximately 850 base pair region of the viral genome. We also included one LRV1 sequence detected in L. (V.) shawi, representing the first report of LRV1 in a species other than L. braziliensis and L. guyanensis. The results illustrate the genetic diversity of the LRV1 strains analyzed here, with smaller divergences detected among viral sequences from the same parasite species. Phylogenetic analyses showed that the LRV1 sequences are grouped according to the parasite species and possibly according to the population of the parasite in which the virus was detected, corroborating the hypothesis of joint evolution of the viruses with the speciation of Leishmania parasites.

PD-L1 May Mediate T-Cell Exhaustion in a Case of Early Diffuse Leishmaniasis Caused by Leishmania (L.) amazonensis

Introduction

Diffuse cutaneous leishmaniasis (DCL) is a rare disease form associated with Leishmania (L.) amazonensis in South America. It represents the “anergic” pole of American Tegumentary Leishmaniasis, and the explanation for its resistance to treatment remains elusive. We aimed to study some possible immunological mechanisms involved in the poor DCL treatment response by evaluating some cell surface molecules obtained from a patient with DCL by flow cytometry.

Case presentation

A 65-year-old DCL patient who initially failed to respond to the standard treatment for the disease showed vacuolated macrophages filled with amastigotes in lesion biopsy, and L. (L.) amazonensis was identified through ITS1PCR amplification. The Leishmania skin test and indirect immunofluorescence analysis revealed negative results. Peripheral blood from the patient was collected after a few months of treatment, when the patient presented with no lesion. Peripheral blood mononuclear cells were analyzed ex vivo and in vitro after 48 h of stimulation with soluble L. (L.) amazonensis antigen (SLA). Cell death, surface molecules, and intracellular molecules, such as IFN-γ and granzyme B, were analyzed in the cells using flow cytometry. Analysis of the surface markers showed an increased expression of the inhibitory molecule programmed death ligand 1 (PD-L1) in the monocytes restimulated with SLA (approximately 65%), whereas the negative controls were 35% positive for PD-L1. Conversely, compared with the negative controls, we observed a decrease in CD4⁺IFN-γ⁺ T cells (8.32 versus 1.7%) and CD8⁺IFN-γ⁺ T cells (14% versus 1%). We also observed a relevant decrease in the granzyme B levels in the CD8⁺ T cells, from 31% in the negative controls to 5% after SLA restimulation.

Conclusion

The dysfunctional activation of PD-L1 inhibitory pathway after Leishmania antigen stimulation and reduced levels of IFN-gamma and granzyme B-producing cells could be closely related to unresponssiveness to standard drug treatment of DCL patient.

RNA viruses in trypanosomatid parasites: a historical overview

Viruses of trypanosomatids are now being extensively studied because of their diversity and the roles they play in flagellates' biology. Among the most prominent examples are leishmaniaviruses implicated in pathogenesis of Leishmania parasites. Here, we present a historical overview of this field, starting with early reports of virus-like particles on electron microphotographs, and culminating in detailed molecular descriptions of viruses obtained using modern next generation sequencing-based techniques. Because of their diversity, different life cycle strategies and host specificity, we believe that trypanosomatids are a fertile ground for further explorations to better understand viral evolution, routes of transitions, and molecular mechanisms of adaptation to different hosts.

Arginase activity of Leishmania isolated from patients with cutaneous leishmaniasis

Cutaneous leishmaniasis (CL) is one of the most important vector-borne parasitic diseases, highly endemic in Iran, and its prevalence is increasing all over the country. Arginase (ARG) activity in isolated Leishmania parasites from CL patients is yet to be explored. This study aimed to compare the ARG activity of isolated Leishmania promastigotes from CL patients with a standard strain of Leishmania major and its influences on the disease pathogenesis. We recruited 16 confirmed CL patients from Qom Province, in central Iran; after detection of Leishmania species using PCR-RFLP, we assessed the levels of ARG in the isolated promastigotes and determined the parasites' growth rate. Only L. major was identified from CL patients. The level of ARG activity in the isolated Leishmania promastigotes from CL patients was significantly higher than that obtained from the standard strain of L. major. No significant correlations between ARG activity and lesion size, number or duration were observed; in contrast, a significant negative correlation was seen between ARG level and Leishmania' growth rate. The obtained results suggest that increased ARG expression and activity in the isolated Leishmania promastigotes might contribute to the higher parasite infectivity and play a major role in the pathogenicity of the CL.

Leptomonas seymouri narna-like virus 1 and not leishmaniaviruses detected in kala-azar samples from India

The great majority of kala-azar/visceral leishmaniasis (VL) cases, which are caused by Leishmania donovani (LD), are reported in Asia. We investigated whether leishmaniaviruses (LRVs) are present in LD isolates. These dsRNA viruses contribute to hyperpathogenicity, as observed in the case of other members of the genus Leishmania. However, LRVs could not be detected in 22 Indian LD isolates tested in the present study, while 70% of these original LD isolates harboured a virus that was not of LD but instead of Leptomonas seymouri (LS) origin. LS is another protozoon that parasitizes the sandfly vector of LD. Historically, LD clinical isolates from India often showed high incidence of LS coinfection. LS was detected in 20 out of the 22 (91%) above-mentioned LD isolates. Leptomonas seymouri narna-like virus 1 (Lepsey NLV1) was identified by whole-genome sequencing in an LD-LS coinfected sample, and its presence was confirmed by PCR and sequencing in 15 (75%) of the 20 LD-LS coinfected samples. The LS-negative LD samples were also virus negative by PCR. That the human host is exposed to an RNA virus in LS, another coinfecting parasite with LD, i.e., the "LD-LS-Lepsey NLV1 triple pathogen" phenomenon, unveils a new paradigm of research towards revisiting the mysteries of Indian leishmaniasis pathogenesis and management.

Unraveling the genetic diversity and phylogeny of Leishmania RNA virus 1 strains of infected Leishmania isolates circulating in French Guiana

Introduction

Leishmania RNA virus type 1 (LRV1) is an endosymbiont of some Leishmania (Vianna) species in South America. Presence of LRV1 in parasites exacerbates disease severity in animal models and humans, related to a disproportioned innate immune response, and is correlated with drug treatment failures in humans. Although the virus was identified decades ago, its genomic diversity has been overlooked until now.

Methodology/Principles findings

We subjected LRV1 strains from 19 L. (V.) guyanensis and one L. (V.) braziliensis isolates obtained from cutaneous leishmaniasis samples identified throughout French Guiana with next-generation sequencing and de novo sequence assembly. We generated and analyzed 24 unique LRV1 sequences over their full-length coding regions. Multiple alignment of these new sequences revealed variability (0.5%–23.5%) across the entire sequence except for highly conserved motifs within the 5’ untranslated region. Phylogenetic analyses showed that viral genomes of L. (V.) guyanensis grouped into five distinct clusters. They further showed a species-dependent clustering between viral genomes of L. (V.) guyanensis and L. (V.) braziliensis, confirming a long-term co-evolutionary history. Noteworthy, we identified cases of multiple LRV1 infections in three of the 20 Leishmania isolates.

Conclusions/Significance

Here, we present the first-ever estimate of LRV1 genomic diversity that exists in Leishmania (V.) guyanensis parasites. Genetic characterization and phylogenetic analyses of these viruses has shed light on their evolutionary relationships. To our knowledge, this study is also the first to report cases of multiple LRV1 infections in some parasites. Finally, this work has made it possible to develop molecular tools for adequate identification and genotyping of LRV1 strains for diagnostic purposes. Given the suspected worsening role of LRV1 infection in the pathogenesis of human leishmaniasis, these data have a major impact from a clinical viewpoint and for the management of Leishmania-infected patients.

Leishmaniasis is a well-known parasitosis due to an infection by the protozoan Leishmania parasites firmly established in South America. In French Guiana, where leishmaniasis is a public health problem, having an annual incidence of 5.6 cases/10,000 inhabitants, 80% of Leishmania spp. parasites are infected by an endosymbiotic virus, Leishmania RNA virus 1 (LRV1). The purpose of this study was to gain insights on the genetic variability and evolution of LRV1 at the genomic level. We subjected 20 Leishmania isolates obtained from cutaneous lesions with next generation sequencing and de novo sequence assembly. This allowed generating 24 LRV1 full-length coding sequences presenting among themselves a significant genetic diversity. The inferred phylogenetic relationships enabled to identify distinct well-supported genotypes and support the hypothesis of co-evolution between LRV1 strains and their parasite hosts. In addition, we identified multiple LRV1 infections in three parasite isolates. Based on these data, future characterization of new strains from other geographic areas and other parasite species should extend knowledge of LRV1 diversification processes. Finally, these results allowed us to identify genomic regions where best to allocate resources for subsequent analyses of LRV1 viral diversity and genotyping that could serve for accurate routine molecular diagnostic applications.