The complete coding region of the maxicircle as a superior phylogenetic marker for exploring evolutionary relationships between members of the Leishmaniinae

In silico prediction of the metabolism of Blastocrithidia nonstop, a trypanosomatid with non-canonical genetic code

BACKGROUND

Almost all extant organisms use the same, so-called canonical, genetic code with departures from it being very rare. Even more exceptional are the instances when a eukaryote with non-canonical code can be easily cultivated and has its whole genome and transcriptome sequenced. This is the case of Blastocrithidia nonstop, a trypanosomatid flagellate that reassigned all three stop codons to encode amino acids.

RESULTS

We in silico predicted the metabolism of B. nonstop and compared it with that of the well-studied human parasites Trypanosoma brucei and Leishmania major. The mapped mitochondrial, glycosomal and cytosolic metabolism contains all typical features of these diverse and important parasites. We also provided experimental validation for some of the predicted observations, concerning, specifically presence of glycosomes, cellular respiration, and assembly of the respiratory complexes.

CONCLUSIONS

In an unusual comparison of metabolism between a parasitic protist with a massively altered genetic code and its close relatives that rely on a canonical code we showed that the dramatic differences on the level of nucleic acids do not seem to be reflected in the metabolisms. Moreover, although the genome of B. nonstop is extremely AT-rich, we could not find any alterations of its pyrimidine synthesis pathway when compared to other trypanosomatids. Hence, we conclude that the dramatic alteration of the genetic code of B. nonstop has no significant repercussions on the metabolism of this flagellate.

Genetic characterization and description of Leishmania (Leishmania) ellisi sp. nov.: a new human-infecting species from the USA

In a 2018 report, an unusual case of cutaneous leishmaniasis was described in a 72-year-old female patient residing in Arizona, United States of America (USA). Preliminary analysis of the 18S rDNA and glyceraldehyde-3-phosphate dehydrogenase genes supported the conclusion that the Leishmania strain (strain 218-L139) isolated from this case was a novel species, though a complete taxonomic description was not provided. Identification of Leishmania at the species level is critical for clinical management and epidemiologic investigations so it is important that novel human-infecting species are characterized taxonomically and assigned a unique scientific name compliant with the ICZN code. Therefore, we sought to provide a complete taxonomic description of Leishmania strain 218-L139. Phylogenetic analysis of several nuclear loci and partial maxicircle genome sequences supported its position within the subgenus Leishmania and further clarified the distinctness of this new species. Morphological characterization of cultured promastigotes and amastigotes from the original case material is also provided. Thus, we conclude that Leishmania (Leishmania) ellisi is a new cause of autochthonous cutaneous leishmaniasis in the USA.

RNA viruses of Crithidia bombi, a parasite of bumblebees

Leishbuviridae (Bunyavirales) are a diverse monophyletic group of negative-sense single-stranded RNA virus infecting parasitic flagellates of the family Trypanosomatidae. The presence of RNA viruses in trypanosomatids can influence the virulence of the latter. Here, we performed a screening for viruses in Crithidia bombi - a common parasite of important pollinators Bombus spp. (bumblebees) that negatively affects its host in stressful conditions. The majority (8/10) of C. bombi isolates collected in Europe and North America were positive for a virus that we named Crithidia bombi leishbuvirus 1 with high conservation of amino acid sequences between isolates. The results of our comparative phylogenetic analyses of the trypanosomatids and their viruses suggest that the high mobility of bumblebees and frequent coinfections by different strains of C. bombi determine an extensive viral exchange between the latter.

Cyclospora cayetanensis comprises at least 3 species that cause human cyclosporiasis

The apicomplexan parasite Cyclospora cayetanensis causes seasonal foodborne outbreaks of the gastrointestinal illness cyclosporiasis. Prior to the coronavirus disease-2019 pandemic, annually reported cases were increasing in the USA, leading the US Centers for Disease Control and Prevention to develop a genotyping tool to complement cyclosporiasis outbreak investigations. Thousands of US isolates and 1 from China (strain CHN_HEN01) were genotyped by Illumina amplicon sequencing, revealing 2 lineages (A and B). The allelic composition of isolates was examined at each locus. Two nuclear loci (CDS3 and 360i2) distinguished lineages A and B. CDS3 had 2 major alleles: 1 almost exclusive to lineage A and the other to lineage B. Six 360i2 alleles were observed – 2 exclusive to lineage A (alleles A1 and A2), 2 to lineage B (B1 and B2) and 1 (B4) was exclusive to CHN_HEN01 which shared allele B3 with lineage B. Examination of heterozygous genotypes revealed that mixtures of A- and B-type 360i2 alleles occurred rarely, suggesting a lack of gene flow between lineages. Phylogenetic analysis of loci from whole-genome shotgun sequences, mitochondrial and apicoplast genomes, revealed that CHN_HEN01 represents a distinct lineage (C). Retrospective examination of epidemiologic data revealed associations between lineage and the geographical distribution of US infections plus strong temporal associations. Given the multiple lines of evidence for speciation within human-infecting Cyclospora, we provide an updated taxonomic description of C. cayetanensis, and describe 2 novel species as aetiological agents of human cyclosporiasis: Cyclospora ashfordi sp. nov. and Cyclospora henanensis sp. nov. (Apicomplexa: Eimeriidae).

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.

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.

Revisiting the heterogeneous global genomic population structure of Leishmania infantum

Leishmania infantum is the main causative agent responsible for visceral leishmaniasis (VL), a disease with global distribution. The genomic structure and genetic variation of this species have been widely studied in different parts of the world. However, in some countries, this information is still yet unknown, as is the genomic behaviour of the main antigens used in VL diagnosis (rK39 and rK28), which have demonstrated variable sensitivity and specificity in a manner dependent on the geographic region analysed. The objective of this study was to explore the genomic architecture and diversity of four Colombian L. infantum isolates obtained in this study and to compare these results with the genetic analysis of 183 L. infantum isolates from across the world (obtained from public databases), as well as to analyse the whole rK39 and rK28 antigen sequences in our dataset. The results showed that, at the global level, L. infantum has high genetic homogeneity and extensive aneuploidy. Furthermore, we demonstrated that there are distinct populations of L. infantum circulating in various countries throughout the globe and that populations of distant countries have close genomic relationships. Additionally, this study demonstrated the high genetic variability of the rK28 antigen worldwide. In conclusion, our study allowed us to (i) expand our knowledge of the genomic structure of global L. infantum; (ii) describe the intra-specific genomic variability of this species; and (iii) understand the genomic characteristics of the main antigens used in the diagnosis of VL. Additionally, this is the first study to report whole-genome sequences of Colombian L. infantum isolates.

Maxicircle architecture and evolutionary insights into Trypanosoma cruzi complex

We sequenced maxicircles from T. cruzi strains representative of the species evolutionary diversity by using long-read sequencing, which allowed us to uncollapse their repetitive regions, finding that their real lengths range from 35 to 50 kb. T. cruzi maxicircles have a common architecture composed of four regions: coding region (CR), AT-rich region, short (SR) and long repeats (LR). Distribution of genes, both in order and in strand orientation are conserved, being the main differences the presence of deletions affecting genes coding for NADH dehydrogenase subunits, reinforcing biochemical findings that indicate that complex I is not functional in T. cruzi. Moreover, the presence of complete minicircles into maxicircles of some strains lead us to think about the origin of minicircles. Finally, a careful phylogenetic analysis was conducted using coding regions of maxicircles from up to 29 strains, and 1108 single copy nuclear genes from all of the DTUs, clearly establishing that taxonomically T. cruzi is a complex of species composed by group 1 that contains clades A (TcI), B (TcIII) and D (TcIV), and group 2 (1 and 2 do not coincide with groups I and II described decades ago) containing clade C (TcII), being all hybrid strains of the BC type. Three variants of maxicircles exist in T. cruzi: a, b and c, in correspondence with clades A, B, and C from mitochondrial phylogenies. While A and C carry maxicircles a and c respectively, both clades B and D carry b maxicircle variant; hybrid strains also carry the b- variant. We then propose a new nomenclature that is self-descriptive and makes use of both the phylogenetic relationships and the maxicircle variants present in T. cruzi.

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.

Total Ortholog Median Matrix as an alternative unsupervised approach for phylogenomics based on evolutionary distance between protein coding genes

The increasing number of available genomic data allowed the development of phylogenomic analytical tools. Current methods compile information from single gene phylogenies, whether based on topologies or multiple sequence alignments. Generally, phylogenomic analyses elect gene families or genomic regions to construct phylogenomic trees. Here, we presented an alternative approach for Phylogenomics, named TOMM (Total Ortholog Median Matrix), to construct a representative phylogram composed by amino acid distance measures of all pairwise ortholog protein sequence pairs from desired species inside a group of organisms. The procedure is divided two main steps, (1) ortholog detection and (2) creation of a matrix with the median amino acid distance measures of all pairwise orthologous sequences. We tested this approach within three different group of organisms: Kinetoplastida protozoa, hematophagous Diptera vectors and Primates. Our approach was robust and efficacious to reconstruct the phylogenetic relationships for the three groups. Moreover, novel branch topologies could be achieved, providing insights about some phylogenetic relationships between some taxa.

Mitochondrial DNAs provide insight into trypanosome phylogeny and molecular evolution

Background

Trypanosomes are single-celled eukaryotic parasites characterised by the unique biology of their mitochondrial DNA. African livestock trypanosomes impose a major burden on agriculture across sub-Saharan Africa, but are poorly understood compared to those that cause sleeping sickness and Chagas disease in humans. Here we explore the potential of the maxicircle, a component of trypanosome mitochondrial DNA to study the evolutionary history of trypanosomes.

Results

We used long-read sequencing to completely assemble maxicircle mitochondrial DNA from four previously uncharacterized African trypanosomes, and leveraged these assemblies to scaffold and assemble a further 103 trypanosome maxicircle gene coding regions from published short-read data. While synteny was largely conserved, there were repeated, independent losses of Complex I genes. Comparison of pre-edited and non-edited genes revealed the impact of RNA editing on nucleotide composition, with non-edited genes approaching the limits of GC loss. African tsetse-transmitted trypanosomes showed high levels of RNA editing compared to other trypanosomes. The gene coding regions of maxicircle mitochondrial DNAs were used to construct time-resolved phylogenetic trees, revealing deep divergence events among isolates of the pathogens Trypanosoma brucei and T. congolense.

Conclusions

Our data represents a new resource for experimental and evolutionary analyses of trypanosome phylogeny, molecular evolution and function. Molecular clock analyses yielded a timescale for trypanosome evolution congruent with major biogeographical events in Africa and revealed the recent emergence of Trypanosoma brucei gambiense and T. equiperdum, major human and animal pathogens.

A review of the systematics, species identification and diagnostics of the Trypanosomatidae using the maxicircle kinetoplast DNA: from past to present

The Trypanosomatid family are a diverse and widespread group of protozoan parasites that belong to the higher order class Kinetoplastida. Containing predominantly monoxenous species (i.e. those having only a single host) that are confined to invertebrate hosts, this class is primarily known for its pathogenic dixenous species (i.e. those that have two hosts), serving as the aetiological agents of the important neglected tropical diseases including leishmaniasis, American trypanosomiasis (Chagas disease) and human African trypanosomiasis. Over the past few decades, a multitude of studies have investigated the diversity, classification and evolutionary history of the trypanosomatid family using different approaches and molecular targets. The mitochondrial-like DNA of the trypanosomatid parasites, also known as the kinetoplast, has emerged as a unique taxonomic and diagnostic target for exploring the evolution of this diverse group of parasitic eukaryotes. This review discusses recent advancements and important developments that have made a significant impact in the field of trypanosomatid systematics and diagnostics in recent years.

Identification of Clinical Infections of Leishmania Imported into Australia: Revising Speciation with Polymerase Chain Reaction-RFLP of the Kinetoplast Maxicircle

Leishmaniasis is a vector-borne disease caused by protozoan parasites of the Leishmania genus. In Australia, leishmaniasis is an imported disease that is presenting itself at increased rates because of international travel, the influx of immigrants, and deployment of military operations to endemic regions. Although Leishmania species are morphologically indistinguishable, there is a strong correlation between some causative species of leishmaniasis and the subsequent response to the treatments available and patient outcome. Consequently, identification of the infective species is imperative as misidentification can result in the administering of an ineffective drug. The aim of this study was to develop a simple diagnostic tool with high sensitivity and specificity, which is capable of detecting the presence of the parasite and accurately differentiating the causative species in question. Using the advantageous properties of the maxi-circle kinetoplast DNA, a polymerase chain reaction (PCR)-restriction fragment length polymorphism (RFLP) targeting the ND7 gene was developed for the analysis of imported cases of human leishmaniasis in Australia. Designed as a dual analysis, concurrent PCR of Leishmania maxi-circle DNA and digestion with two separate enzymes (NlaIII and HpyCH4IV), this study provides an appraisal on 24 imported cases of leishmaniasis between 2008 and 2017. Five Leishmania species were reported, with members of the Viannia subgenus being the most common. The implementation of novel diagnostic procedures for leishmaniasis such as the one reported here is needed to establish a gold standard practice for the diagnosis and treatment of leishmaniasis.

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.

Evolutionary Insight into the Trypanosomatidae Using Alignment-Free Phylogenomics of the Kinetoplast

Advancements in next-generation sequencing techniques have led to a substantial increase in the genomic information available for analyses in evolutionary biology. As such, this data requires the exponential growth in bioinformatic methods and expertise required to understand such vast quantities of genomic data. Alignment-free phylogenomics offer an alternative approach for large-scale analyses that may have the potential to address these challenges. The evolutionary relationships between various species within the trypanosomatid family, specifically members belonging to the genera Leishmania and Trypanosoma have been extensively studies over the last 30 years. However, there is a need for a more exhaustive analysis of the Trypanosomatidae, summarising the evolutionary patterns amongst the entire family of these important protists. The mitochondrial DNA of the trypanosomatids, better known as the kinetoplast, represents a valuable taxonomic marker given its unique presence across all kinetoplastid protozoans. The aim of this study was to validate the reliability and robustness of alignment-free approaches for phylogenomic analyses and its applicability to reconstruct the evolutionary relationships between the trypanosomatid family. In the present study, alignment-free analyses demonstrated the strength of these methods, particularly when dealing with large datasets compared to the traditional phylogenetic approaches. We present a maxicircle genome phylogeny of 46 species spanning the trypanosomatid family, demonstrating the superiority of the maxicircle for the analysis and taxonomic resolution of the Trypanosomatidae.