Phylogenetic analysis of Trypanosomatina (Protozoa: Kinetoplastida) based on minicircle conserved regions

Kinetoplast Genome of Leishmania spp. Is under Strong Purifying Selection

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

Low genetic heterogeneity of Leishmania major in different geographical regions of Iran

To examine the genetic diversity of Leishmania major, 100 Giemsa-stained positive slides were collected from endemic foci of Iran (Northeast, Central, and Southwest provinces) over two consecutive years during 2019-2021. The Leishmania ITS-rDNA gene was amplified and Leishmania sp. was recognized by PCR-RFLP and sequencing. In addition, 178 registered ITS-rDNA sequences from other geographical regions of Iran were retrieved from GenBank, including different host species (human, sandfly and rodent). A total of 40 new haplotypes were discovered using the ITS-rDNA sequence analysis. IR29 (20.6%) and IR34 (61%) were the two most common haplotypes, represented by a star-like feature in the overall population. Analysis of the molecular variance test revealed low genetic diversity of L. major in human cases (Haplotype diversity; 0.341), rodent (Hd; 0.387) and sandfly (Hd; 0.390) sequences. The lowest genetic diversity of L. major was observed in Southwest/Southeast Iran (Hd: 0.104-0.286). The statistically Fst value indicated that L. major is not genetically differentiated between geographic regions of Iran, except for the Northeast-Southwest (Fst: 0.29055) and Central-Southwest (Fst: 0.30294) population pairs. The current study as the first investigation discloses new perspectives for further evaluation in the identification local transmission paradigms and initiating effective prevention strategies.

Complete minicircle genome of Leptomonas pyrrhocoris reveals sources of its non-canonical mitochondrial RNA editing events

Uridine insertion/deletion (U-indel) editing of mitochondrial mRNA, unique to the protistan class Kinetoplastea, generates canonical as well as potentially non-productive editing events. While the molecular machinery and the role of the guide (g) RNAs that provide required information for U-indel editing are well understood, little is known about the forces underlying its apparently error-prone nature. Analysis of a gRNA:mRNA pair allows the dissection of editing events in a given position of a given mitochondrial transcript. A complete gRNA dataset, paired with a fully characterized mRNA population that includes non-canonically edited transcripts, would allow such an analysis to be performed globally across the mitochondrial transcriptome. To achieve this, we have assembled 67 minicircles of the insect parasite Leptomonas pyrrhocoris, with each minicircle typically encoding one gRNA located in one of two similar-sized units of different origin. From this relatively narrow set of annotated gRNAs, we have dissected all identified mitochondrial editing events in L. pyrrhocoris, the strains of which dramatically differ in the abundance of individual minicircle classes. Our results support a model in which a multitude of editing events are driven by a limited set of gRNAs, with individual gRNAs possessing an inherent ability to guide canonical and non-canonical editing.

Novel organization of mitochondrial minicircles and guide RNAs in the zoonotic pathogen Trypanosoma lewisi

Kinetoplastid flagellates are known for several unusual features, one of which is their complex mitochondrial genome, known as kinetoplast (k) DNA, composed of mutually catenated maxi- and minicircles. Trypanosoma lewisi is a member of the Stercorarian group of trypanosomes which is, based on human infections and experimental data, now considered a zoonotic pathogen. By assembling a total of 58 minicircle classes, which fall into two distinct categories, we describe a novel type of kDNA organization in T. lewisi. RNA-seq approaches allowed us to map the details of uridine insertion and deletion editing events upon the kDNA transcriptome. Moreover, sequencing of small RNA molecules enabled the identification of 169 unique guide (g) RNA genes, with two differently organized minicircle categories both encoding essential gRNAs. The unprecedented organization of minicircles and gRNAs in T. lewisi broadens our knowledge of the structure and expression of the mitochondrial genomes of these human and animal pathogens. Finally, a scenario describing the evolution of minicircles is presented.

Molecular detection of Leishmania spp. in Brazilian cross-border south region mammalian hosts

Visceral leishmaniasis is an endemic zoonotic disease identified especially in developing territories. Brazil's northeast, southeast and midwest have been endemic for several years; currently, the infection is spreading to the south. Dogs are the main reservoirs; however, other mammal species have also been infected. Herein, we have identified the infecting Leishmania species in dogs and horses from the south of Brazil, a new outbreak of the infection. Blood samples were collected in the urban area of Uruguaiana city. DNA was extracted from peripheral blood, kinetoplast DNA (kDNA) and ribosomal DNA (rDNA) fragments were obtained by polymerase chain reaction (PCR) and sequenced. Out of 123 samples, 25 of them (14 dogs and 11 horses) were positive for Leishmania spp. Sequence alignment and phylogenetic analysis revealed that the kDNA in positive samples was similar to four species previously reported: L. infantum/L. chagasi, L. donovani, L. major. Despite kDNA minicircles regions are very useful due to high sensitivity to Leishmania spp. DNA detection, the sequence polymorphism among minicircles can be an obstacle to interspecific differentiation. Our results suggest that these strains are circulating in Brazil south region cross-border and indicate the susceptibility of new outbreak for visceral leishmaniasis infection in horses domiciled in endemic region for canine and human visceral leishmaniasis.

Encapsulation of the HSP-90 Chaperone Inhibitor 17-AAG in Stable Liposome Allow Increasing the Therapeutic Index as Assessed, in vitro, on Leishmania (L) amazonensis Amastigotes-Hosted in Mouse CBA Macrophages

The current long-term treatment for leishmaniasis causes severe side effects and resistance in some cases. An evaluation of the anti-leishmanial potential of an HSP90-inhibitor, 17-allylamino-17-demethoxygeldanamycin (17-AAG), demonstrated its potent effect against Leishmania spp. in vitro and in vivo. We have previously shown that 17-AAG can kill L. (L) amazonensis promastigotes with an IC₅₀ of 65 nM and intracellular amastigote at concentrations as low as 125 nM. As this compound presents low solubility and high toxicity in human clinical trials, we prepared an inclusion complex containing hydroxypropyl-β-cyclodextrin and 17-AAG (17-AAG:HPβCD) to improve its solubility. This complex was characterized by scanning electron microscopy, and X-ray diffraction. Liposomes-containing 17-AAG:HPβCD was prepared and evaluated for encapsulation efficiency (EE%), particle size, polydispersity index (PDI), pH, and zeta potential, before and after accelerated and long-term stability testing. An evaluation of leishmanicidal activity against promastigotes and intracellular amastigotes of L. (L) amazonensis was also performed. The characterization techniques utilized confirmed the formation of the inclusion complex, HPβCD:17-AAG, with a resulting 33-fold-enhancement in compound water solubility. Stability studies revealed that 17-AAG:HPβCD-loaded liposomes were smaller than 200 nm, with 99% EE. Stability testing detected no alterations in PDI that was 0.295, pH 7.63, and zeta potential +22.6, suggesting liposome stability, and suitability for evaluating leishmanicidal activity. Treatment of infected macrophages with 0.006 nM of 17-AAG:HPβCD or 17-AAG:HPβCD-loaded liposomes resulted in almost complete amastigote clearance inside macrophages after 48 h. This reduction is similar to the one observed in infected macrophages treated with 2 μM amphotericin B. Our results showed that nanotechnology and drug delivery systems could be used to increase the antileishmanial efficacy and potency of 17-AAG in vitro, while also resulting in reduced toxicity that indicates these formulations may represent a potential therapeutic strategy against leishmaniasis.

Minicircle Kinetoplast Genome of Insect Trypanosomatid Leptomonas pyrrhocoris

We present here the structure of a minicircle population based on transcriptome sequencing of Leptomonas pyrrhocoris. We show that minicircle DNA molecules are dimeric. As in dixenous species, the entire molecule of minicircle DNA is transcribed. This is the first minicircle transcriptome of monoxenous trypanosomatid species determined using NGS technology.

New Approaches to Systematics of Trypanosomatidae: Criteria for Taxonomic (Re)description

While dixenous trypanosomatids represent one of the most dangerous pathogens for humans and domestic animals, their monoxenous relatives have frequently become model organisms for studies of diversity of parasitic protists and host-parasite associations. Yet, the classification of the family Trypanosomatidae is not finalized and often confusing. Here we attempt to make a blueprint for future studies in this field. We would like to elicit a discussion about an updated procedure, as traditional taxonomy was not primarily designed to be used for protists, nor can molecular phylogenetics solve all the problems alone. The current status, specific cases, and examples of generalized solutions are presented under conditions where practicality is openly favored over rigid taxonomic codes or blind phylogenetic approach.

Exploring the environmental diversity of kinetoplastid flagellates in the high-throughput DNA sequencing era

The class Kinetoplastea encompasses both free-living and parasitic species from a wide range of hosts. Several representatives of this group are responsible for severe human diseases and for economic losses in agriculture and livestock. While this group encompasses over 30 genera, most of the available information has been derived from the vertebrate pathogenic genera Leishmaniaand Trypanosoma. Recent studies of the previously neglected groups of Kinetoplastea indicated that the actual diversity is much higher than previously thought. This article discusses the known segment of kinetoplastid diversity and how gene-directed Sanger sequencing and next-generation sequencing methods can help to deepen our knowledge of these interesting protists.

The associations of Leishmania major and Leishmania tropica aspects by focusing their morphological and molecular features on clinical appearances in Khuzestan Province, Iran

Cutaneous leishmaniasis has various phenotypic aspects consisting of polymorphic amastigotes with different genetic ranges. Samples were collected from suspected patients of Khuzestan province. Prepared smears were stained, scaled, and measured using ocular micrometer. The Cyt b, ITS-rDNA, and microsatellite genes of Leishmania were amplified and Leishmania species were identified by molecular analyses. Of 150 examined suspected patients, 102 were identified to Leishmania species (90 L. major, nine L. tropica, and three unidentified). The amastigotes of 90 L. major had regular and different irregular shapes within three clinical lesions with no and/or low genetic diversity. Three haplotypes of Cyt b of L. major were found but no variation was observed using ITS-rDNA gene. Interesting findings were that all nine L. tropica had regular amastigote shapes with more genetic variations, also a patient which had coinfection of L. major, L. tropica, and Crithidia. At least two L. major and L. tropica were identified in suspected patients of the regions. Different irregular amastigotes' shapes of L. major can be explained by various reservoir hosts and vectors. In contrast, more molecular variations in L. tropica could be justified by genetic characters. Unidentified Leishmania could be mixed pathogens or nonpathogens with mammals' Leishmania or Crithidia.

Diversity of trypanosomatids (Kinetoplastea: Trypanosomatidae) parasitizing fleas (Insecta: Siphonaptera) and description of a new genus Blechomonas gen. n

To further investigate the diversity of Trypanosomatidae we have examined the species present within the flea (Siphonaptera) population in the Czech Republic. Out of 1549 fleas, 239 were found to be infected by trypanosomatids. Axenic cultures were established from 90 infected specimens and 29 of them were further characterized. Molecular phylogenetic analysis of the SL RNA, gGAPDH, and SSU rRNA genes revealed a striking diversity within this group and analyzed isolates were classified into 16 Typing units (TUs) of which 15 typified new species. In addition to one Trypanosoma species, two TUs grouped within the sub-family Leishmaniinae, two clustered together with Herpetomonas, wheras 11 TUs formed a novel clade branching off between Trypanosoma spp. and remaining trypanosomatids. We propose to recognize this clade as a new genus Blechomonas and a new subfamily Blechomonadinae, and provide molecular and morphological description of 11 TUs representing this genus. Our finding of such an ancient host-specific group sheds new light at the origin of Trypanosomatidae and the roots of dixenous parasitism. The strict host restriction of Blechomonas to Siphonaptera with adult fleas' dependence on blood meal may reflect passing of parasites from larvae through pupae to adults and implies potential transmission to the warm-blooded vertebrates.

Analysis of kinetoplast cytochrome b gene of 16 Leishmania isolates from different foci of China: different species of Leishmania in China and their phylogenetic inference

BACKGROUND

Leishmania species belong to the family Trypanosomatidae and cause leishmaniasis, a geographically widespread disease that infects humans and other vertebrates. This disease remains endemic in China. Due to the large geographic area and complex ecological environment, the taxonomic position and phylogenetic relationship of Chinese Leishmania isolates remain uncertain. A recent internal transcribed spacer 1 and cytochrome oxidase II phylogeny of Chinese Leishmania isolates has challenged some aspects of their traditional taxonomy as well as cladistics hypotheses of their phylogeny. The current study was designed to provide further disease background and sequence analysis.

METHODS

We systematically analyzed 50 cytochrome b (cyt b) gene sequences of 19 isolates (16 from China, 3 from other countries) sequenced after polymerase chain reaction (PCR) using a special primer for cyt b as well as 31 sequences downloaded from GenBank. After alignment, the data were analyzed using the maximum parsimony, Bayesian and netwok methods.

RESULTS

Sequences of six haplotypes representing 10 Chinese isolates formed a monophyletic group and clustered with Leishmania tarentolae. The isolates GS1, GS7, XJ771 of this study from China clustered with other isolates of Leishmania donovani complex. The isolate JS1 was a sister to Leishmania tropica, which represented an L. tropica complex instead of clustering with L. donovani complex or with the other 10 Chinese isolates. The isolates KXG-2 and GS-GER20 formed a monophyletic group with Leishmania turanica from central Asia. In the different phylogenetic trees, all of the Chinese isolates occurred in at least four groups regardless of geographic distribution.

CONCLUSIONS

The undescribed Leishmania species of China, which are clearly causative agents of canine leishmaniasis and human visceral leishmaniasis and are related to Sauroleishmania, may have evolved from a common ancestral parasite that came from the Americas and may have split off earlier than the other old world Leishmania. Our results also suggest the following: the isolates GS7, GS1 and XJ771 occur as part of the L. donovani complex; the JS1 isolate is L. tropica; and the isolate GS-GER20 identified as Leishmania gerbilli is close to KXG-2 which is L. turanica.

Diversity and phylogeny of insect trypanosomatids: all that is hidden shall be revealed

Monoxenous trypanosomatids, which are usually regarded as benign dwellers of the insect alimentary tract, represent a relatively obscure group within the family Trypanosomatidae. This field of study has long been in disarray with the genus level taxonomy of this group remaining artificial, species criteria elusive, host specificity and occurrence poorly known, and their diversity mostly unexplored. The time has arrived to remedy this situation: a phylogenetic approach has been applied to taxa recognition and description, and a culture-independent (PCR-based) approach for detection and identification of organisms in nature has made it feasible to study the diversity of the group. Although more than 100 typing units have been discovered recently, these appear to represent a small segment of trypanosomatid biodiversity, which still remains to be uncovered.

Molecular approaches for a better understanding of the epidemiology and population genetics ofLeishmania

Molecular approaches are being used increasingly for epidemiological studies of visceral and cutaneous leishmaniases. Several molecular markers resolving genetic differences betweenLeishmaniaparasites at species and strain levels have been developed to address key epidemiological and population genetic questions. The current gold standard, multilocus enzyme typing (MLEE), needs cultured parasites and lacks discriminatory power. PCR assays identifying species directly with clinical samples have proven useful in numerous field studies. Multilocus sequence typing (MLST) is potentially the most powerful phylogenetic approach and will, most probably, replace MLEE in the future. Multilocus microsatellite typing (MLMT) is able to discriminate below the zymodeme level and seems to be the best candidate for becoming the gold standard for distinction of strains. Population genetic studies by MLMT revealed geographical and hierarchic population structure inL. tropica, L. majorand theL. donovanicomplex. The existence of hybrids and gene flow betweenLeishmaniapopulations suggests that sexual recombination is more frequent than previously thought. However, typing and analytical tools need to be further improved. Accessible databases should be created and sustained for integrating data obtained by different researchers. This would allow for global analyses and help to avoid biases in analyses due to small sample sizes.

Rational sub-division of plant trypanosomes (Phytomonas spp.) based on minicircle conserved region analysis

The sequences of minicircle conserved regions from various plant trypanosomatids have been determined and analyzed. The goal of this study was to add another tool to the arsenal of molecular probes for distinguishing between the different trypanosomatids occurring in plants: systemic trypanosomatids multiplying in the sap, those from the laticiferous tubes, and those developing in fruits, seeds or flowers but not in the plant itself and that are frequently considered as opportunistic insect trypanosomatids. As some plant intraphloemic trypanosomatids are the causative agents of important diseases, a clear definition of the different types of trypanosomatids is critical. The conserved region of the mitochondrial minicircle provides several specific features in a small sequence region containing three functionally elements required for minicircle replication. Trees generated from the analysis recapitulated trees drawn from analyses of isoenzymes, RAPD, and particular gene sequences, supporting the validity of the small region used in this work. Three groups of isolates were significant and in accordance with previous work. The peculiarity of phloem-restricted trypanosomatids associated with wilts of coconut and oil palm in Latin America - group H - is confirmed. In agreement with previous studies on their biological and serological properties the results highlighted this group called 'phloemicola'. It always differentiated from all other latex and fruit isolates or opportunistic trypanosomatids, like insect trypanosomatids. We can assert that phloemicola is the only well-defined taxon among all plant trypanosomatids. A group of non-pathogenic latex isolates from South American euphorbs (G), and a heterogenous group (A) including one fruit, one possible latex and one insect isolate are clearly distinct groups. The group of Mediterranean isolates from latex (D), even with a low boostrap, stood out well from other groups. The remainder of the isolates fell into a heterogeneous cluster. At least eight different groups in the plant trypanosomatids were identified.

Sequence variation of the cytochrome b gene of various human infecting members of the genus Leishmania and their phylogeny

The Cytochrome b (Cyt b) gene has proved to be useful for identification and classification of many mammals and plants. In order to evaluate the utility of this gene for discrimination of Leishmania parasites as well as for exploring their phylogenetic relationships, we determined the nucleotide sequences of the Cyt b gene from 13 human-infecting Leishmania species (14 strains) from the New and Old Worlds. The Cyt b genes, approximately 1080 base pairs, were found to be A/T rich, and their 5' terminal-editing regions were highly conserved. The nucleotide sequence variation among them was enough to discriminate parasite species; 245 nucleotide positions were polymorphic and 190 positions were parsimony informative. The phylogenetic relationships based on this gene, showed good agreement with the classification of Lainson & Shaw (1987) except for the inclusion of L. (L.) major in the L. (L.) tropica complex and the placement of L. tarentolae in another genus. These data show that the Cyt b gene is useful for phylogenetic study of Leishmania parasites.

Phylogenetic position of the giant anuran trypanosomes Trypanosoma chattoni, Trypanosoma fallisi, Trypanosoma mega, Trypanosoma neveulemairei, and Trypanosoma ranarum inferred from 18S rRNA gene sequences

Phylogenetic relationships within the kinetoplastid flagellates were inferred from comparisons of small-subunit ribosomal RNA gene sequences. These included 5 new gene sequences, Trypanosoma fallisi (2,239 bp), Trypanosoma chattoni (2,180 bp), Trypanosoma mega (2,211 bp), Trypanosoma neveulemairei (2,197 bp), and Trypanosoma ranarum (2,203 bp). Trees produced using maximum-parsimony and distance-matrix methods (least-squares, neighbor-joining, and maximum-likelihood), supported by strong bootstrap and quartet-puzzle analyses, indicated that the trypanosomes are a monophyletic group that divides into 2 major lineages, the salivarian trypanosomes and the nonsalivarian trypanosomes. The nonsalivarian trypanosomes further divide into 2 lineages, 1 containing trypanosomes of birds, mammals, and reptiles and the other containing trypanosomes of fish, reptiles, and anurans. Among the giant trypanosomes, T. chattoni is clearly shown to be distantly related to all the other anuran trypanosome species. Trypanosoma mega is closely associated with T. fallisi and T. ranarum, whereas T. neveulemairei and Trypanosoma rotatorium are sister taxa. The branching order of the anuran trypanosomes suggests that some toad trypanosomes may have evolved by host switching from frogs to toads.

Kinetoplast DNA minicircles of phloem-restricted Phytomonas associated with wilt diseases of coconut and oil palms have a two-domain structure

We report the cloning and sequencing of the first minicircle from a phloem-restricted, pathogenic Phytomonas sp. (Hart 1) isolated from a coconut palm with hartrot disease. The minicircle possessed a two-domain structure of two conserved regions, each containing three conserved sequence blocks (CSB). Based on the sequence around CSB 3 from Hart 1, PCR primers were designed to allow specific amplification of Phytomonas minicircles. This primer pair demonstrated specificity for at least six groups of plant trypanosomatids and did not amplify from insect trypanosomatids. The PCR results were consistent with a two-domain structure for other plant trypanosomatids.

Some current problems in the systematics of Trypanosomatids

Trypanosomatids have been traditionally allocated to a number of genera that were described based on morphological features and host range. Recently molecular studies have provided new data that has allowed a reexamination of the genera. While in some cases the molecular data has been in agreement with the morphological characters they have also reinforced existing doubts about some current generic divisions as well as raising new concerns. A revision of the trypanosomatid genera is required. Suggested features of such a revision would include: (1) The possible division of Trypanosoma into new genera to reflect the wide genetic diversity of this group; (2) The inclusion of Leishmania, Sauroleishmania and Endotrypanum within a single genus given their high genetic affinity; (3) The complete revision of the monogenetic typanosomatid genera to reflect monophyletic groups; (4) A more precise redescription of Phytomonas so as to only include the monophyletic plant flagellates.