Mitochondrial genomics of human pathogenic parasite Leishmania (Viannia) panamensis

High-throughput analysis of the Trypanosoma cruzi minicirculome (mcDNA) unveils structural variation and functional diversity

Trypanosoma cruzi causes Chagas disease and has a unique extranuclear genome enclosed in a structure called the kinetoplast, which contains circular genomes known as maxi- and minicircles. While the structure and function of maxicircles are well-understood, many aspects of minicircles remain to be discovered. Here, we performed a high-throughput analysis of the minicirculome (mcDNA) in 50 clones isolated from Colombia's diverse T. cruzi I populations. Results indicate that mcDNA comprises four diverse subpopulations with different structures, lengths, and numbers of interspersed semi-conserved (previously termed ultra-conserved regions mHCV) and hypervariable (mHVPs) regions. Analysis of mcDNA ancestry and inter-clone differentiation indicates the interbreeding of minicircle sequence classes is placed along diverse strains and hosts. These results support evidence of the multiclonal dynamics and random bi-parental segregation. Finally, we disclosed the guide RNA repertoire encoded by mcDNA at a clonal scale, and several attributes of its abundance and function are discussed.

A phased genome assembly of a Colombian Trypanosoma cruzi TcI strain and the evolution of gene families

Chagas is an endemic disease in tropical regions of Latin America, caused by the parasite Trypanosoma cruzi. High intraspecies variability and genome complexity have been challenges to assemble high quality genomes needed for studies in evolution, population genomics, diagnosis and drug development. Here we present a chromosome-level phased assembly of a TcI T. cruzi strain (Dm25). While 29 chromosomes show a large collinearity with the assembly of the Brazil A4 strain, three chromosomes show both large heterozygosity and large divergence, compared to previous assemblies of TcI T. cruzi strains. Nucleotide and protein evolution statistics indicate that T. cruzi Marinkellei separated before the diversification of T. cruzi in the known DTUs. Interchromosomal paralogs of dispersed gene families and histones appeared before but at the same time have a more strict purifying selection, compared to other repeat families. Previously unreported large tandem arrays of protein kinases and histones were identified in this assembly. Over one million variants obtained from Illumina reads aligned to the primary assembly clearly separate the main DTUs. We expect that this new assembly will be a valuable resource for further studies on evolution and functional genomics of Trypanosomatids.

Identification of a unique conserved region from a kinetoplastid genome of Leishmania orientalis (formerly named Leishmania siamensis) strain PCM2 in Thailand

Mitochondrial DNAs (mtDNAs) appear in almost all eukaryotic species and are useful molecular markers for phylogenetic studies and species identification. Kinetoplast DNAs (kDNAs) are structurally complex circular mtDNA networks in kinetoplastids, divided into maxicircles and minicircles. Despite several kDNAs of many Leishmania species being examined, the kDNAs of the new species, Leishmania orientalis (formerly named Leishmania siamensis) strain PCM2, have not been explored. This study aimed to investigate the maxicircle and minicircle DNAs of L. orientalis strain PCM2 using hybrid genome sequencing technologies and bioinformatic analyses. The kDNA sequences were isolated and assembled using the SPAdes hybrid assembler from the Illumina short-read and PacBio long-read data. Circular contigs of the maxicircle and minicircle DNAs were reconstructed and confirmed by BLASTn and rKOMICs programs. The kDNA genome was annotated by BLASTn before the genome comparison and phylogenetic analysis by progressiveMauve, MAFFT, and MEGA programs. The maxicircle of L. orientalis strain PCM2 (18,215 bp) showed 99.92% similarity and gene arrangement to Leishmania enriettii strain LEM3045 maxicircle with variation in the 12s rRNA gene and divergent region. Phylogenetics of the whole sequence, coding regions, divergent regions, and 12s rRNA gene also confirmed this relationship and subgenera separation. The identified 105 classes of minicircles (402-1177 bp) were clustered monophyletically and related to the Leishmania donovani minicircles. The kinetoplast maxicircle and minicircle DNAs of L. orientalis strain PCM2 contained a unique conserved region potentially useful for specific diagnosis of L. orientalis and further exploration of this parasite population genetics in Thailand and related regions.

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.

Multiparametric approach to assess the disease severity and progression of cutaneous leishmaniasis infection

The pathophysiology of Cutaneous Leishmaniasis (CL), an infection caused by Leishmania tropica (L. tropica) and Leishmania major (L. major) is primarily determined by inflammation-mediated immune cells. The immune response mainly depends on cells and molecules related to T-cells that influence susceptibility and disease development. Understanding the immunological mechanisms that cause tissue injury or lesion healing is critical for developing appropriate treatment strategies. In the present study, T-cells profile and cell-free mitochondrial DNA (CF mt-DNA) were investigated in CL patients infected with L. tropica (n = 34) and L. major (n = 2) and compared with non-infected healthy controls (n = 20). There was a significant (p<0.0001) difference between CD4+ T-cells among L. tropica and L. major CL-infected groups as compared to control while no significant difference (p = 0.8597) was found in the percentages of CD8+ T-cells. When L. tropica and L. major CL-infected individuals were compared to controls, the levels of IL-4 and expression of CF mt-DNA were significantly higher (p<0.0001). Higher levels of CF mt-DNA were detected in CL patients, irrespective of the infective Leishmania species. We proposed that the levels of CF mt-DNA and IL-4 in CL-infected individuals can be used to determine the disease progression. A better understanding of these biomarkers and evaluation of the immune responses in CL patients might benefit the development of vaccines and immunotherapies.

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

Parasites of medical importance, such as Leishmania and Trypanosoma, are characterized by the presence of thousands of circular DNA molecules forming a structure known as kinetoplast, within the mitochondria. The maxicircles, which are equivalent to the mitochondrial genome in other eukaryotes, have been proposed as a promising phylogenetic marker. Using whole-DNA sequencing data, it is also possible to assemble maxicircle sequences as shown here and in previous works. In this study, based on data available in public databases and using a bioinformatics workflow previously reported by our group, we assembled the complete coding region of the maxicircles for 26 prototypical strains of trypanosomatid species. Phylogenetic analysis based on this dataset resulted in a robust tree showing an accurate taxonomy of kinetoplastids, which was also able to discern between closely related Leishmania species that are usually difficult to discriminate by classical methodologies. In addition, we provide a dataset of the maxicircle sequences of 60 Leishmania infantum field isolates from America, Western Europe, North Africa, and Eastern Europe. In agreement with previous studies, our data indicate that L. infantum parasites from Brazil are highly homogeneous and closely related to European strains, which were transferred there during the discovery of America. However, this study showed the existence of different L. infantum populations/clades within the Mediterranean region. A maxicircle signature for each clade has been established. Interestingly, two L. infantum clades were found coexisting in the same region of Spain, one similar to the American strains, represented by the Spanish JPCM5 reference strain, and the other, named “non-JPC like”, may be related to an important leishmaniasis outbreak that occurred in Madrid a few years ago. In conclusion, the maxicircle sequence emerges as a robust molecular marker for phylogenetic analysis and species typing within the kinetoplastids, which also has the potential to discriminate intraspecific variability.

Experimental Hybridization in Leishmania: Tools for the Study of Genetic Exchange

Despite major advances over the last decade in our understanding of Leishmania reproductive strategies, the sexual cycle in Leishmania has defied direct observation and remains poorly investigated due to experimental constraints. Here, we summarize the findings and conclusions drawn from genetic analysis of experimental hybrids generated in sand flies and highlight the recent advances in generating hybrids in vitro. The ability to hybridize between culture forms of different species and strains of Leishmania should invite more intensive investigation of the mechanisms underlying genetic exchange and provide a rich source of recombinant parasites for future genetic analyses.

Molecular Diagnosis of Leishmaniasis: Quantification of Parasite Load by a Real-Time PCR Assay with High Sensitivity

Real-time PCR was developed to quantify Leishmania infantum kinetoplast DNA and optimized to achieve a sensitivity of 1 parasite/mL. For this purpose, we cloned the conserved kDNA fragment of 120 bp into competent cells and correlated them with serial dilutions of DNA extracted from reference parasite cultures calculating that a parasite cell contains approximately 36 molecules of kDNA. This assay was applied to estimate parasite load in clinical samples from visceral, cutaneous leishmaniasis patients and infected dogs and cats comparing with conventional diagnosis. The study aimed to propose a real-time PCR for the detection of Leishmania DNA from clinical samples trying to solve the diagnostic problems due to the low sensitivity of microscopic examination or the low predictive values of serology and resolve problems related to in vitro culture. The quantitative PCR assay in this study allowed detection of Leishmania DNA and quantification of considerably low parasite loads in samples that had been diagnosed negative by conventional techniques. In conclusion, this quantitative PCR can be used for the diagnosis of both human, canine and feline Leishmaniasis with high sensitivity and specificity, but also for evaluating treatment and the endpoint determination of leishmaniasis.

In silico and in vitro evaluation of primers for molecular differentiation of Leishmania species

Leishmaniasis is a zoonotic disease caused by over 20 species of protozoan parasites of the genus Leishmania. Infection is commonly spread by sandflies and produces a wide spectrum of clinical signs and symptoms. Therefore, from an epidemiological and therapeutic standpoint, it is important to detect and differentiate Leishmania spp. The objective of this study was to combinate in silico and in vitro strategies to evaluate the analytical specificity of primers previously described in the literature. According to electronic PCR (e-PCR) analysis, 23 out of 141 pairs of primers selected through literature search matched their previously reported analytical specificity. In vitro evaluation of nine of these primer pairs by quantitative PCR (qPCR) confirmed the analytical specificity of five of them at the level of Leishmania spp., L. mexicana complex or Leishmania and Viannia subgenera. Based on these findings, the combination of e-PCR and qPCR is suggested to be a valuable approach to maximize the specificity of new primer pairs for the laboratory diagnosis of infections with Leishmania spp.

Common Structural Patterns in the Maxicircle Divergent Region of Trypanosomatidae

Maxicircles of all kinetoplastid flagellates are functional analogs of mitochondrial genome of other eukaryotes. They consist of two distinct parts, called the coding region and the divergent region (DR). The DR is composed of highly repetitive sequences and, as such, remains the least explored segment of a trypanosomatid genome. It is extremely difficult to sequence and assemble, that is why very few full length maxicircle sequences were available until now. Using PacBio data, we assembled 17 complete maxicircles from different species of trypanosomatids. Here we present their large-scale comparative analysis and describe common patterns of DR organization in trypanosomatids.

Leishmania Mitochondrial Genomes: Maxicircle Structure and Heterogeneity of Minicircles

The mitochondrial DNA (mtDNA), which is present in almost all eukaryotic organisms, is a useful marker for phylogenetic studies due to its relative high conservation and its inheritance manner. In Leishmania and other trypanosomatids, the mtDNA (also referred to as kinetoplast DNA or kDNA) is composed of thousands of minicircles and a few maxicircles, catenated together into a complex network. Maxicircles are functionally similar to other eukaryotic mtDNAs, whereas minicircles are involved in RNA editing of some maxicircle-encoded transcripts. Next-generation sequencing (NGS) is increasingly used for assembling nuclear genomes and, currently, a large number of genomic sequences are available. However, most of the time, the mitochondrial genome is ignored in the genome assembly processes. The aim of this study was to develop a pipeline to assemble Leishmania minicircles and maxicircle DNA molecules, exploiting the raw data generated in the NGS projects. As a result, the maxicircle molecules and the plethora of minicircle classes for Leishmania major, Leishmania infantum and Leishmania braziliensis have been characterized. We have observed that whereas the heterogeneity of minicircle sequences existing in a single cell hampers their use for Leishmania typing and classification, maxicircles emerge as an extremely robust genetic marker for taxonomic studies within the clade of kinetoplastids.