Assembly of a Large Collection of Maxicircle Sequences and Their Usefulness for Leishmania Taxonomy and Strain Typing

Advances and Challenges in the Diagnosis of Leishmaniasis

Leishmaniasis remains a significant public health challenge, particularly in endemic regions with limited resources. Traditional diagnostic methods, including microscopy, culture, and serology, though widely utilized, often suffer from limitations such as variable  sensitivity, time delays, and the need for specialized infrastructure. Some of these limitations have been addressed with the emergence of molecular diagnostic techniques. Quantitative PCR (q-PCR), loop-mediated isothermal amplification (LAMP), and recombinase polymerase amplification (RPA) assays have improved  the diagnostic sensitivity and specificity, enabling species identification and detection of asymptomatic infections. Further, nanodiagnostics and portable sequencing technologies such as the MinION™, along with lab-on-chip platforms, are revolutionizing the diagnostic landscape of leishmaniasis by offering point-of-care (POC) options for remote settings and field-based diagnosis. This review provides an in-depth analysis of these cutting-edge advances, discusses their application in resource-constrained settings, and evaluates their potential to reshape the future of leishmaniasis diagnosis and management.

Mucosal leishmaniasis is associated with the Leishmania RNA virus and inappropriate cutaneous leishmaniasis treatment

BACKGROUND

Mucosal leishmaniasis (ML) is a severe clinical form of leishmaniasis that is characterized by the destruction of the nasal and/or the oral mucosae and appears as a late complication in 5% to 10% of cutaneous leishmaniasis (CL) cases produced by species belonging to Leishmania (Viannia) subgenus. Some strains of Leishmania spp. carry an RNA virus known as Leishmania RNA virus (LRV) that may contribute to the appearance of ML.

METHODS

To examine the role of LRV type 1 (LRV1) as a risk factor associated with ML, a retrospective case-control study involving 103 patients was conducted. Cases were defined as patients with ML (n = 33), and controls corresponded to patients with CL and without mucosal lesions (n = 70). Clinical data were recorded from the patient's medical records. Cryopreserved biopsies were used to detect LRV1 and identify Leishmania species.

RESULTS

The frequency of LRV1 in the 103 patients was 16.5% (95% CI,10.4-25.12) being higher in samples from cases [33.33% (95% CI,18.89-51.76) than from controls [8.57% (95% CI, 3.82-18.10)]. L. (V.) braziliensis was identified in 63.6% of cases and 55.7% of the controls. Multivariate logistic regression indicated that infection with Leishmania spp. carrying LRV1 (OR = 6.30; 95% CI,1.52-26.10, p = 0.011) acts as risk factors for ML occurrence, while the completed treatment for the cutaneous event decreases the risk of ML (OR = 0.039; 95% CI, 0.01-0.12, p < 0.0001).

CONCLUSIONS

Our data support the association between LRV1 and ML occurrence and emphasize the effect of completed treatment for CL in preventing ML.

Typing of Leishmania isolates from vectors and leporids of the Madrid (Spain) outbreak

In 2009, a large outbreak of leishmaniasis, associated with environmental changes, was declared near Madrid (Spain), in which Phlebotomus perniciosus was the vector, whereas the main reservoirs were hares and rabbits. Analysis of isolates from humans, vectors and leporids from the focus identified the Leishmania infantum ITS-Lombardi genotype. However, multilocus enzyme electrophoresis (MLEE), the reference technique for Leishmania typing, and sequencing of the hsp70 gene, a commonly used marker, were not performed. In the present study, 19 isolates from P. perniciosus (n = 11), hares (n = 5) and rabbits (n = 3) from the outbreak area, all characterized as ITS-Lombardi in previous studies, were analysed by MLEE and hsp70 sequencing. The hsp70 results confirmed that all the analysed strains are L. infantum. However, by MLEE, 4 different zymodemes of L. infantum were identified based on variable mobilities of the NP1 enzyme: MON-34 (NP1100, n = 11), MON-80 (NP1130, n = 6), MON-24 (NP1140, n = 1) and MON-331 (NP1150, n = 1). The relative frequency of these zymodemes does not correspond to their usual occurrence in Spain. Moreover, MON-34 and MON-80 were found in P. perniciosus, hares and rabbits for the first time. These findings continue to provide insights into the outbreak and call for further studies with a higher number of strains.

Genetic variability of Leishmania (Leishmania) infantum causing human visceral leishmaniasis in the Southeastern Brazil

Leishmania infantum is a protozoan that causes visceral leishmaniasis (VL) in the Americas and some regions of Europe. The disease is mainly characterized by hepatosplenomegaly and fever, and can be fatal. Factors related to the host and parasite can contribute to the transmission of Leishmania and the clinical outcome. The intraspecific genetic variability of L. infantum strains may be one of these factors. In this study, we evaluated the genetic variability of L. infantum obtained from bone marrow smear slides from patients in the Sao Paulo State, Brazil. For this, the minicircle of the kDNA hypervariable region was used as target by Sanger sequencing. By analyzing the similarity of the nucleotides and the maximum likelihood tree (Fasttree), we observed a high similarity (98%) among samples. Moreover, we identified four different profiles of L. infantum. In conclusion, L. infantum strains from Sao Paulo State, Brazil, showed low diversity measured by minicircle of the kDNA hypervariable region.

Identification of a conserved maxicircle and unique minicircles as part of the mitochondrial genome of Leishmania martiniquensis strain PCM3 in Thailand

BACKGROUND

The mitochondrial DNA of trypanosomatids, including Leishmania, is known as kinetoplast DNAs (kDNAs). The kDNAs form networks of hundreds of DNA circles that are evidently interlocked and require complex RNA editing. Previous studies showed that kDNA played a role in drug resistance, adaptation, and survival of Leishmania. Leishmania martiniquensis is one of the most frequently observed species in Thailand, and its kDNAs have not been illustrated.

METHODS

This study aimed to extract the kDNA sequences from Illumina short-read and PacBio long-read whole-genome sequence data of L. martiniquensis strain PCM3 priorly isolated from the southern province of Thailand. A circular maxicircle DNA was reconstructed by de novo assembly using the SPAdes program, while the minicircle sequences were retrieved and assembled by the rKOMIC tool. The kDNA contigs were confirmed by blasting to the NCBI database, followed by comparative genomic and phylogenetic analysis.

RESULTS

We successfully constructed the complete circular sequence of the maxicircle (19,008 bp) and 214 classes of the minicircles from L. martiniquensis strain PCM3. The genome comparison and annotation showed that the maxicircle structure of L. martiniquensis strain PCM3 was similar to those of L. enriettii strain LEM3045 (84.29%), L. arabica strain LEM1108 (82.79%), and L. tarentolae (79.2%). Phylogenetic analysis also showed unique evolution of the minicircles of L. martiniquensis strain PCM3 from other examined Leishmania species.

CONCLUSIONS

This was the first report of the complete maxicircle and 214 minicircles of L. martiniquensis strain PCM3 using integrated whole-genome sequencing data. The information will be helpful for further improvement of diagnosis methods and monitoring genetic diversity changes of this parasite.

Identification of Leishmania infantum Epidemiology, Drug Resistance and Pathogenicity Biomarkers with Nanopore Sequencing

The emergence of drug-resistant strains of the parasite Leishmania infantum infecting dogs and humans represents an increasing threat. L. infantum genomes are complex and unstable with extensive structural variations, ranging from aneuploidies to multiple copy number variations (CNVs). These CNVs have recently been validated as biomarkers of Leishmania concerning virulence, tissue tropism, and drug resistance. As a proof-of-concept to develop a novel diagnosis platform (LeishGenApp), four L. infantum samples from humans and dogs were nanopore sequenced. Samples were epidemiologically typed within the Mediterranean L. infantum group, identifying members of the JCP5 and non-JCP5 subgroups, using the conserved region (CR) of the maxicircle kinetoplast. Aneuploidies were frequent and heterogenous between samples, yet only chromosome 31 tetrasomy was common between all the samples. A high frequency of aneuploidies was observed for samples with long passage history (MHOM/TN/80/IPT-1), whereas fewer were detected for samples maintained in vivo (MCRI/ES/2006/CATB033). Twenty-two genes were studied to generate a genetic pharmacoresistance profile against miltefosine, allopurinol, trivalent antimonials, amphotericin, and paromomycin. MHOM/TN/80/IPT-1 and MCRI/ES/2006/CATB033 displayed a genetic profile with potential resistance against miltefosine and allopurinol. Meanwhile, MHOM/ES/2016/CATB101 and LCAN/ES/2020/CATB102 were identified as potentially resistant against paromomycin. All four samples displayed a genetic profile for resistance against trivalent antimonials. Overall, this proof-of-concept revealed the potential of nanopore sequencing and LeishGenApp for the determination of epidemiological, drug resistance, and pathogenicity biomarkers in L. infantum.