Novel detection of Leishmania RNA virus-1 (LRV-1) in clinical isolates of Leishmania Viannia panamensis.

American tegumentary leishmaniasis comprises a discrete set of clinical presentations endemic to Latin America. Leishmania RNA virus-1 (LRV-1) is a double-stranded RNA virus identified in 20­25% of the Leishmania Viannia braziliensis and L. V. guyanensis, however not in L. V. panamensis. This is the first report of LRV-1 in L. V. panamensis and its associations with clinical phenotypes of ATL. Unique surplus discard clinical isolates of L. V. panamensis were identified from the Public Health Ontario Laboratory (PHOL) and the Leishmania Clinic of the Instituto de Medicina Tropical 'Alexander von Humboldt' between 2012 and 2019 and screened for LRV-1 by real-time polymerase chain reaction. Patient isolates were stratified according to clinical phenotype. Of 30 patients with L. V. panamensis, 14 (47%) and 16 (53%) patients had severe and non-severe ATL, respectively. Five (36%) of 14 severe cases and 2 (12%) of 16 non-severe cases were positive for LRV-1, respectively. No differences in sex were observed for clinical phenotype and LRV-1 status. Although an association between LRV-1 status and clinical phenotype was not demonstrated, this is the first description of the novel detection of LRV-1 in L. V. panamensis, a species that has been documented predominantly in Central America.

A fine mapping of single nucleotide variants and haplotype analysis of IL13 gene in patients with Leishmania guyanensis-cutaneous leishmaniasis and plasma cytokines IL-4, IL-5, and IL-13.

Introduction: Leishmaniasis continues to pose a substantial health burden in 97 countries worldwide. The progression and outcome of Leishmania infection are influenced by various factors, including the cytokine milieu, the skin microbiota at the infection site, the specific Leishmania species involved, the genetic background of the host, and the parasite load. In endemic regions to leishmaniasis, only a fraction of individuals infected actually develops the disease. Overexpression of IL-13 in naturally resistant C57BL/6 mice renders them susceptible to L. major infection. Haplotypes constructed from several single nucleotide variant (SNV) along a chromosome fragment may provide insight into any SNV near the fragment that may be genuinely associated with a phenotype in genetic association studies. Methods: We investigated nine SNVs (SNV1rs1881457A>C, SNV2rs1295687C>G, SNV3rs2069744C>T, SNV4rs2069747C>T, SNV5rs20541A>G, SNV6rs1295685A>G, SNV7rs848A>C, SNV8rs2069750G >C, and SNV9rs847T>C) spanning the entire IL13 gene in patients with L. guyanensis cutaneous leishmaniasis (Lg-CL). Results: Our analysis did not reveal any significant association between the SNVs and susceptibility/protection against Lg-CL development. However, haplotype analysis, excluding SNV4rs2069747 and SNV8rs2069750 due to low minor allele frequency, revealed that carriers of the haplotype CCCTAAC had a 93% reduced likelihood developing Lg-CL. Similarly, the haplotypes ACCCGCT (ORadj=0.02 [95% CI 0.00-0.07]; p-value, 6.0×10-19) and AGCTAAC (ORadj=0.00[95% CI 0.00-0.00]; p-value 2.7×10-12) appeared to provide protection against the development of Lg-CL. Conversely, carriers of haplotype ACCTGCC have 190% increased likelihood of developing Lg-CL (ORadj=2.9 [95%CI 1.68-5.2]; p-value, 2.5×10-6). Similarly, haplotype ACCCAAT (ORadj=2.7 [95%CI 1.5-4.7]; p-value, 3.2×10-5) and haplotype AGCCGCC are associated with susceptibility to the development of Lg-CL (ORadj=1.7[95%CI 1.04-2.8]; p-value, 0.01). In our investigation, we also found a correlation between the genotypes of rs2069744, rs20541, rs1295685, rs847, and rs848 and plasma IL-5 levels among Lg-Cl patients. Furthermore, rs20541 showed a correlation with plasma IL-13 levels among Lg-Cl patients, while rs2069744 and rs848 showed a correlation with plasma IL-4 levels among the same group. Conclusions: Overall, our study identifies three haplotypes of IL13 associated with resistance to disease development and three haplotypes linked to susceptibility. These findings suggest the possibility of a variant outside the gene region that may contribute, in conjunction with other genes, to differences in susceptibility and partially to the pathology.

Variants of CARD8 in Leishmania guyanensis-cutaneous leishmaniasis and influence of the variants genotypes on circulating plasma cytokines IL-1ß, TNFα and IL-8.

Nucleotide-binding oligomerization domain, leucine-rich repeat-containing protein family (NLR) are intracellular pathogen recognition receptors mediating innate immunity, releasing proinflammatory cytokines IL-1ß and IL-18, and promoting pyroptotic cell death, upon sensing pathogenic or endogenous danger signals. In animal models, NLRP3 inflammasome has a dual role, pathogenic or protective in Leishmania-infection, depending on the Leishmania species and mice strain. Caspase recruitment containing domain 8 (CARD8) is a negative regulator of NLRP3 inflammasome and also an inhibitor of transcription factor NFĸB, a major transcription factor of proinflammatory cytokines. We investigated whether single nucleotide variants in CARD8 may partially explain why only a proportion of individuals coming from the same area of endemicity of leishmaniasis develop cutaneous leishmaniasis caused by Leishmania guyanensis. We genotyped four single nucleotide variants of the CARD8 gene by direct nucleotide sequencing in 1741 individuals from an endemic area of leishmaniasis, constituting 850 patients with CL and 891 healthy controls. The frequencies of the genotypes of the variants rs2288877 T>C, rs73944113 C>T, and rs2043211 A>T are similar among the patients with CL and HC, while the variant rs2288876 A>G) reveals an excess of the genotype AA among the patients with CL (44%) compared to 37% in the HC group. Allele A of the variant rs2288876 A>G) is associated with susceptibility to CL (OR = 1.2 [95%CI 1.03-1.4]; P = 0.01). Haplotype analysis showed that individuals harboring the haplotype CCAA have 280% odds of developing CL caused by L. guyanensis (OR = 3.8 [95% CI 2.0-7.7]; p = 0.00004). The variants rs2288877 T>C and rs2288876 A>G correlate with the plasma level of IL-8. Spearman correlation showed a significant positive correlation between the rs2288876 A>G allele A and the level of IL-8 (ρ = 0.22; p = 0.0002). CARD8 may partially contribute to the development of CL caused by L. guyanensis.

Detection of Leishmania RNA Virus 1 in Leishmania (Viannia) panamensis Isolates, Panama.

We detected Leishmania RNA virus 1 (LRV1) in 11 isolates of Leishmania (Viannia) panamensis collected during 2014-2019 from patients from different geographic areas in Panama. The distribution suggested a spread of LRV1 in L. (V.) panamensis parasites. We found no association between LRV1 and an increase in clinical pathology.

Diversity and geographical distribution of Leishmania species and the emergence of Leishmania (Leishmania) infantum and L. (Viannia) panamensis in Central-Western Venezuela.

Transmission of cutaneous leishmaniasis in Venezuela reveals diverse and changing epidemiological landscapes, as well as a spectrum of clinical phenotypes presumed to be linked to a variety of Leishmania species. Central-western Venezuela constitutes one of the highest endemic epicenters in the country, and updated molecular epidemiological information is still lacking. Therefore, in this study we aimed to characterize the landscape of circulating Leishmania species across central-western Venezuela through the last two decades, performed comparisons of haplotype and nucleotide diversity, and built a geospatial map of parasite species distribution. A total of 120 clinical samples were collected from patients across the cutaneous disease spectrum, retrieving parasitic DNA, and further characterizing by PCR and sequencing of the HSP70 gene fragment. This data was later collated with further genetic, geospatial and epidemiological analyses. A peculiar pattern of species occurrence including Leishmania (Leishmania) amazonensis (77.63% N=59), Leishmania (Leishmania) infantum (14.47% N=11), Leishmania (Viannia) panamensis (5.26% N=4) and Leishmania (Viannia) braziliensis (2.63% N=2) was revealed, also highlighting a very low genetic diversity amongst all analyzed sequences. Geographical distribution showed that most cases are widely distributed across the greater urban-sub urban area of the Irribaren municipality. L.(L.) amazonensis appears to be widely dispersed throughout Lara state. Statistical analyses failed to reveal significance for any comparisons, leading to conclude a lack of association between the infective Leishmania species and clinical phenotypes. To the best of our knowledge, this is an unprecedented study which addresses comprehensively the geographical distribution of Leishmania species in central-western Venezuela throughout the last two decades, and the first to incriminate L. (L.) infantum as an etiologic agent of cutaneous leishmaniasis in this region. Our findings support that Leishmania endemism in central-western Venezuela is caused mainly by L.(L.) amazonensis. Future studies are needed to unveil additional details on the ecological intricacies and transmission aspects of leishmaniasis (i.e. sampling phlebotomines and mammals) and to adopt adequate public health prevention and control strategies and mitigate disease impact in this endemic region.

Raising the suspicion of a non-autochthonous infection: identification of Leishmania guyanensis from Costa Rica exhibits a Leishmaniavirus related to Brazilian north-east and French Guiana viral genotypes.

BACKGROUND: Costa Rica has a history of neglecting prevention, control and research of leishmaniasis, including limited understanding on Leishmania species causing human disease across the country and a complete lack of knowledge on the Leishmania RNA virus, described as a factor linked to the worsening and metastasis of leishmanial lesions. OBJECTIVES: The aim of this work was to describe a case of cutaneous leishmaniasis by Leishmania (Viannia) guyanensis, bearing infection with Leishmaniavirus 1 (LRV1) in Costa Rica, raising the suspicion of imported parasites in the region. METHODS: The Leishmania strain was previously identified by routine hsp70 polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) in Costa Rica and subsequently characterised by isoenzyme electrophoresis and Sanger sequencing in Brazil. Screening for LRV1 was conducted with a dual RT-PCR approach and sequencing of the fragment obtained. FINDINGS: Since 2016 Costa Rica performs Leishmania isolation and typing as part of its epidemiological surveillance activities. Amongst 113 strains typed until 2019, only one was characterised as a L. (V.) guyanensis, corresponding to the first confirmed report of this species in the country. Interestingly, the same strain tested positive for LRV1. Sequencing of the viral orf1 and 2, clustered this sample with other LRV1 genotypes of South American origin, from the Northeast of Brazil and French Guiana. MAIN CONCLUSION: The unique characteristics of this finding raised the suspicion that it was not an autochthonous strain. Notwithstanding its presumed origin, this report points to the occurrence of said endosymbiont in Central American Leishmania strains. The possibility of its local dispersion represents one more challenge faced by regional health authorities in preventing and controlling leishmaniasis.

Sequencing of hsp70 for discernment of species from the Leishmania (Viannia) guyanensis complex from endemic areas in Colombia.

BACKGROUND: Colombia is ranked very high among countries with the highest numbers of endemic Leishmania species (n = 9) causing human disease. Although much effort has been devoted to generating simple and specific tools for Leishmania species identification, challenges remain in the discrimination of species belonging to the Leishmania (Viannia) guyanensis complex: L. (V.) guyanensis and L. (V.) panamensis. METHODS: A set of seven reference strains of species belonging to the L. (Leishmania) and L. (Viannia) subgenera, clinical strains from human cases of cutaneous leishmaniasis (CL; n = 26) and samples collected from sylvatic mammals and sand flies (n = 7) from endemic areas in Colombia were analyzed in this study. The heat-shock protein 70 gene (hsp70) was amplified by PCR from DNA extracted from logarithmic-phase promastigotes or tissue samples, and the PCR products were sequenced. Sequence alignment was performed against a set of previously published and curated sequences, and phylogenetic analysis based on the maximum-likelihood and Bayesian inference approaches was conducted. Haplotype diversity among strains and species of the L. (V.) guyanensis complex was explored using a median-joining network. RESULTS: Sequencing of the hsp70 gene for L. (Viannia) spp. typing was comparable to species identification using isoenzyme electrophoresis or monoclonal antibodies. Complete species matching was found, except for one sylvatic sample with an identity yet unsolved. Among the L. (V.) panamensis clinical strains, two distinctive phylogenetic clusters were found to correlate with two different zymodemes: L. (V.) panamensis Z2.2 and Z2.3. Analysis of samples from sylvatic environments identified novel records of naturally infected wild mammal and sand fly species. CONCLUSIONS: Our results support the adequacy of hsp70 gene sequencing as a single-locus approach for discrimination of L. (Viannia) spp., as well as for exploring the genetic diversity within the L. (V.) guyanensis complex.

Genomic diversity and genetic variation of Leishmania panamensis within its endemic range.

Species belonging to the Leishmania (Viannia) subgenus are important causative agents of cutaneous and mucocutaneous leishmaniasis in Central and South America. These parasites possess several distinctive biological features that are influenced by their genetics, population structure, and genome instability. To date, several studies have revealed varying degrees of genetic diversity within Leishmania species. Particularly, in species of the L. (Viannia) subgenus, a generalized high intraspecific genetic diversity has been reported, although, conflicting conclusions have been drawn using different molecular techniques. Despite being the most common Leishmania species circulating in Panama and Colombia, few studies have analyzed clinical samples of Leishmania panamensis using whole-genome sequencing, and their restricted number of samples has limited the information they can provide to understand the population structure of L. panamensis. Here, we used next generation sequencing (NGS) to explore the genetic diversity of L. panamensis within its endemic range, analyzing data from 43 isolates of Colombian and Panamanian origin. Our results show the occurrence of three well-defined geographically correlated groups, and suggests the possible occurrence of additional phylogeographic groups. Furthermore, these results support the existence of a mixed mode of reproduction in L. panamensis, with varying frequencies of events of genetic recombination occurring primarily within subpopulations of closely related strains. This study offers important insights into the population genetics and reproduction mode of L. panamensis, paving the way to better understand their population structure and the emergence and maintenance of key eco-epidemiological traits.

IL-23R Variant rs11805303 Is Associated With Susceptibility to the Development of Cutaneous Leishmaniasis in Leishmania guyanensis-Infected Individuals.

BACKGROUND: Emerging evidence suggests that the interleukin (IL) 17/ IL-23 axis may play a role in the pathogenesis of leishmaniasis. Our aim was to investigate whether the IL-23R variant rs11805303 is a risk factor for the development of cutaneous leishmaniasis (CL) in Leishmania guyanensis-infected individuals. METHODS: We genotyped by polymerase chain reaction-restriction fragment length polymorphism the rs11805303 C/T in 828 patients with CL and 806 healthy individuals. Plasma tumor necrosis factor-α, IL-6, interferon-γ, IL-1ß, and IL-17 were measured with the Bioplex assay. RESULTS: The distribution of the genotypes differed between patients with CL and healthy controls with a common odds ratio of 1.78 (P = 2.2 × 10-11) for the disease-associated T allele. Leishmania guyanensis-infected individuals homozygous for the T allele show a 200% increased risk of progressing to disease development, with a 95% confidence interval ranging from 81% to 400% (P = 9.9 × 10-6) in comparison to individuals homozygous for the C allele. Males homozygous for the T allele have higher plasma levels of IL-17 compared with heterozygous or homozygous CC individuals. CONCLUSIONS: The present association of the IL-23R variant rs11805303 with the development of CL suggests that the IL-17/IL-23 axis may play an important role in the pathogenesis of CL.

Filling the gaps in Leishmania naiffi and Leishmania guyanensis genome plasticity.

Insufficient and irregular data reports on Leishmaniasis, issuing from the developing world, have left much to be desired in terms of understanding the molecular signatures producing distinct infectious phenotypes of the disease. Herein, we report on the complete genome sequencing of Leishmania naiffi and Leishmania guyanensis, sampled from patients in regions of Colombia and Venezuela. In this study, the isolates of cutaneous lesions from both species presented limited structural variation at the chromosomal level, low gene copy number variation, and high genetic heterogeneity. We compared these sequences to the reference genomes hitherto related from Brazil and French Guyana. Although of the same species, we note a consequential genomic disparity between the Venezuelan and French Guyanese isolates of L. guyanensis. Although less significant on the global schema of cutaneous and mucosal disease, such genomic studies of L. naiffi and L. guyanensis substantiate the gaps in understanding of the molecular architecture and multivariate clinical pictures of Leishmaniasis, on an international scale.