The Leishmania major maxicircle divergent region is variable in different isolates and cell types

Mitochondrial DNA Structure in Trypanosoma cruzi

Kinetoplastids display a single, large mitochondrion per cell, with their mitochondrial DNA referred to as the kinetoplast. This kinetoplast is a network of concatenated circular molecules comprising a maxicircle (20-64 kb) and up to thousands of minicircles varying in size depending on the species (0.5-10 kb). In Trypanosoma cruzi, maxicircles contain typical mitochondrial genes found in other eukaryotes. They consist of coding and divergent/variable regions, complicating their assembly due to repetitive elements. However, next-generation sequencing (NGS) methods have resolved these issues, enabling the complete sequencing of maxicircles from different strains. Furthermore, several insertions and deletions in the maxicircle sequences have been identified among strains, affecting specific genes. Unique to kinetoplastids, minicircles play a crucial role in a particular U-insertion/deletion RNA editing system by encoding guide RNAs (gRNAs). These gRNAs are essential for editing and maturing maxicircle mRNAs. In Trypanosoma cruzi, although only a few studies have utilized NGS methods to date, the structure of these molecules suggests a classification into four main groups of minicircles. This classification is based on their size and the number of highly conserved regions (mHCRs) and hypervariable regions (mHVRs).

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.

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.

Laboratory diagnostics for human Leishmania infections: a polymerase chain reaction-focussed review of detection and identification methods

Leishmania infections span a range of clinical syndromes and impact humans from many geographic foci, but primarily the world's poorest regions. Transmitted by the bite of a female sand fly, Leishmania infections are increasing with human movement (due to international travel and war) as well as with shifts in vector habitat (due to climate change). Accurate diagnosis of the 20 or so species of Leishmania that infect humans can lead to the successful treatment of infections and, importantly, their prevention through modelling and intervention programs. A multitude of laboratory techniques for the detection of Leishmania have been developed over the past few decades, and although many have drawbacks, several of them show promise, particularly molecular methods like polymerase chain reaction. This review provides an overview of the methods available to diagnostic laboratories, from traditional techniques to the now-preferred molecular techniques, with an emphasis on polymerase chain reaction-based detection and typing methods.

Comparative genomics of Leishmania isolates from Brazil confirms the presence of Leishmania major in the Americas

Leishmania (Leishmania) major is an important agent of cutaneous leishmaniasis, having as a vector sandflies belonging to the genus Phlebotomus. Although this species has been described as restricted to the Old World, parasites similar to L. major have been isolated from South American patients who have never travelled abroad. These parasites were named "L. major-like", and several studies have been carried out to characterise them biochemically, molecularly, and biologically. However, the phylogenetic origin of these isolates is still unknown. In the present study we characterised three L. major-like isolates, named BH49, BH121 and BH129, using comparative genomics approaches. We evaluated the presence of gene and segmental duplications/deletions and the presence of aneuploidies that could explain the differences in infectivity observed in the BH49 and BH121 isolates. All isolates presented a pattern of mosaic aneuploidy and gene copy number variation, which are common in the genus Leishmania. Virulence factors such as phosphatases and peptidases were found to have increased gene copy numbers in the infective isolate, which could explain the difference in infectivity previously observed between BH121 and BH49. Phylogenetic analyses revealed that BH49, BH121 and BH129 L. major-like grouped with L. major isolates, and suggest they were imported from the Old World in at least two independent events. We suggest that new epidemiological inquiries should also evaluate L. major infections in South America, to assess the epidemiological importance of this species in the New World.

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.

Separating the Wheat from the Chaff: RNA Editing and Selection of Translatable mRNA in Trypanosome Mitochondria

In the mitochondria of trypanosomes and related kinetoplastid protists, most mRNAs undergo a long and sophisticated maturation pathway before they can be productively translated by mitochondrial ribosomes. Some of the aspects of this pathway (identity of the promotors, transcription initiation, and termination signals) remain obscure, and some (post-transcriptional modification by U-insertion/deletion, RNA editing, 3'-end maturation) have been illuminated by research during the last decades. The RNA editing creates an open reading frame for a productive translation, but the fully edited mRNA often represents a minor fraction in the pool of pre-edited and partially edited precursors. Therefore, it has been expected that the final stages of the mRNA processing generate molecular hallmarks, which allow for the efficient and selective recognition of translation-competent templates. The general contours and several important details of this process have become known only recently and represent the subject of this review.

Molecular Diagnosis of Visceral Leishmaniasis

Visceral leishmaniasis (VL), a deadly parasitic disease, is a major public health concern globally. Countries affected by VL have signed the London Declaration on Neglected Tropical Diseases and committed to eliminate VL as a public health problem by 2020. To achieve and sustain VL elimination, it will become progressively important not to miss any remaining cases in the community who can maintain transmission. This requires accurate identification of symptomatic and asymptomatic carriers using highly sensitive diagnostic tools at the primary health service setting. The rK39 rapid diagnostic test (RDT) is the most widely used tool and with its good sensitivity and specificity is the first choice for decentralized diagnosis of VL in endemic areas. However, this test cannot discriminate between current, subclinical, or past infections and is useless for diagnosis of relapses and as a prognostic (cure) test. Importantly, as the goal of elimination of VL as a public health problem is approaching, the number of people susceptible to infection will increase. Therefore, correct diagnosis using a highly sensitive diagnostic test is crucial for applying appropriate treatment and management of cases. Recent advances in molecular techniques have improved Leishmania detection and quantification, and therefore this technology has  become increasingly relevant due to its possible application in a variety of clinical sample types. Most importantly, given current problems in identifying asymptomatic individuals because of poor correlation between the main methods of detection, molecular tests are valuable for VL elimination programs, especially to monitor changes in burden of infection in specific communities. This review provides a comprehensive overview of the available VL diagnostics and discusses the usefulness of molecular methods in the diagnosis, quantification, and species differentiation as well as their clinical applications.

Recent advances in trypanosomatid research: genome organization, expression, metabolism, taxonomy and evolution

Unicellular flagellates of the family Trypanosomatidae are obligatory parasites of invertebrates, vertebrates and plants. Dixenous species are aetiological agents of a number of diseases in humans, domestic animals and plants. Their monoxenous relatives are restricted to insects. Because of the high biological diversity, adaptability to dramatically different environmental conditions, and omnipresence, these protists have major impact on all biotic communities that still needs to be fully elucidated. In addition, as these organisms represent a highly divergent evolutionary lineage, they are strikingly different from the common 'model system' eukaryotes, such as some mammals, plants or fungi. A number of excellent reviews, published over the past decade, were dedicated to specialized topics from the areas of trypanosomatid molecular and cell biology, biochemistry, host-parasite relationships or other aspects of these fascinating organisms. However, there is a need for a more comprehensive review that summarizing recent advances in the studies of trypanosomatids in the last 30 years, a task, which we tried to accomplish with the current paper.

Leishmania infections: Molecular targets and diagnosis

Progress in the diagnosis of leishmaniases depends on the development of effective methods and the discovery of suitable biomarkers. We propose firstly an update classification of Leishmania species and their synonymies. We demonstrate a global map highlighting the geography of known endemic Leishmania species pathogenic to humans. We summarize a complete list of techniques currently in use and discuss their advantages and limitations. The available data highlights the benefits of molecular markers in terms of their sensitivity and specificity to quantify variation from the subgeneric level to species complexes, (sub) species within complexes, and individual populations and infection foci. Each DNA-based detection method is supplied with a comprehensive description of markers and primers and proposal for a classification based on the role of each target and primer in the detection, identification and quantification of leishmaniasis infection. We outline a genome-wide map of genes informative for diagnosis that have been used for Leishmania genotyping. Furthermore, we propose a classification method based on the suitability of well-studied molecular markers for typing the 21 known Leishmania species pathogenic to humans. This can be applied to newly discovered species and to hybrid strains originating from inter-species crosses. Developing more effective and sensitive diagnostic methods and biomarkers is vital for enhancing Leishmania infection control programs.

Analysis of the mitochondrial maxicircle of Trypanosoma lewisi, a neglected human pathogen

BACKGROUND

The haemoflagellate Trypanosoma lewisi is a kinetoplastid parasite which, as it has been recently reported to cause human disease, deserves increased attention. Characteristic features of all kinetoplastid flagellates are a uniquely structured mitochondrial DNA or kinetoplast, comprised of a network of catenated DNA circles, and RNA editing of mitochondrial transcripts. The aim of this study was to describe the kinetoplast DNA of T. lewisi.

METHODS/RESULTS

In this study, purified kinetoplast DNA from T. lewisi was sequenced using high-throughput sequencing in combination with sequencing of PCR amplicons. This allowed the assembly of the T. lewisi kinetoplast maxicircle DNA, which is a homologue of the mitochondrial genome in other eukaryotes. The assembly of 23,745 bp comprises the non-coding and coding regions. Comparative analysis of the maxicircle sequence of T. lewisi with Trypanosoma cruzi, Trypanosoma rangeli, Trypanosoma brucei and Leishmania tarentolae revealed that it shares 78%, 77%, 74% and 66% sequence identity with these parasites, respectively. The high GC content in at least 9 maxicircle genes of T. lewisi (ATPase6; NADH dehydrogenase subunits ND3, ND7, ND8 and ND9; G-rich regions GR3 and GR4; cytochrome oxidase subunit COIII and ribosomal protein RPS12) implies that their products may be extensively edited. A detailed analysis of the non-coding region revealed that it contains numerous repeat motifs and palindromes.

CONCLUSIONS

We have sequenced and comprehensively annotated the kinetoplast maxicircle of T. lewisi. Our analysis reveals that T. lewisi is closely related to T. cruzi and T. brucei, and may share similar RNA editing patterns with them rather than with L. tarentolae. These findings provide novel insight into the biological features of this emerging human pathogen.

Multiple mitochondrial introgression events and heteroplasmy in trypanosoma cruzi revealed by maxicircle MLST and next generation sequencing

Background

Mitochondrial DNA is a valuable taxonomic marker due to its relatively fast rate of evolution. In Trypanosoma cruzi, the causative agent of Chagas disease, the mitochondrial genome has a unique structural organization consisting of 20–50 maxicircles (∼20 kb) and thousands of minicircles (0.5–10 kb). T. cruzi is an early diverging protist displaying remarkable genetic heterogeneity and is recognized as a complex of six discrete typing units (DTUs). The majority of infected humans are asymptomatic for life while 30–35% develop potentially fatal cardiac and/or digestive syndromes. However, the relationship between specific clinical outcomes and T. cruzi genotype remains elusive. The availability of whole genome sequences has driven advances in high resolution genotyping techniques and re-invigorated interest in exploring the diversity present within the various DTUs.

Methodology/Principal Findings

To describe intra-DTU diversity, we developed a highly resolutive maxicircle multilocus sequence typing (mtMLST) scheme based on ten gene fragments. A panel of 32 TcI isolates was genotyped using the mtMLST scheme, GPI, mini-exon and 25 microsatellite loci. Comparison of nuclear and mitochondrial data revealed clearly incongruent phylogenetic histories among different geographical populations as well as major DTUs. In parallel, we exploited read depth data, generated by Illumina sequencing of the maxicircle genome from the TcI reference strain Sylvio X10/1, to provide the first evidence of mitochondrial heteroplasmy (heterogeneous mitochondrial genomes in an individual cell) in T. cruzi.

Conclusions/Significance

mtMLST provides a powerful approach to genotyping at the sub-DTU level. This strategy will facilitate attempts to resolve phenotypic variation in T. cruzi and to address epidemiologically important hypotheses in conjunction with intensive spatio-temporal sampling. The observations of both general and specific incidences of nuclear-mitochondrial phylogenetic incongruence indicate that genetic recombination is geographically widespread and continues to influence the natural population structure of TcI, a conclusion which challenges the traditional paradigm of clonality in T. cruzi.

Chagas disease, caused by the protozoan parasite Trypanosoma cruzi, is an important public health problem in Latin America. While molecular techniques can differentiate the major T. cruzi genetic lineages, few have sufficient resolution to describe diversity among closely related strains. The online availability of three mitochondrial genomes allowed us to design a multilocus sequence typing (mtMLST) scheme to exploit these rapidly evolving markers. We compared mtMLST with current nuclear typing tools using isolates belonging to the oldest and most widely occurring lineage TcI. T. cruzi is generally believed to reproduce clonally. However, in this study, distinct branching patterns between mitochondrial and nuclear phylogenetic trees revealed multiple incidences of genetic exchange within different geographical populations and major lineages. We also examined Illumina sequencing data from the TcI genome strain which revealed multiple different mitochondrial genomes within an individual parasite (heteroplasmy) that were, however, not sufficiently divergent to represent a major source of typing error. We strongly recommend this combined nuclear and mitochondrial genotyping methodology to reveal cryptic diversity and genetic exchange in T. cruzi. The level of resolution that this mtMLST provides should greatly assist attempts to elucidate the complex interactions between parasite genotype, clinical outcome and disease distribution.

Maxicircle (mitochondrial) genome sequence (partial) of Leishmania major: gene content, arrangement and composition compared with Leishmania tarentolae

We report 8420 bp of DNA sequence data from the maxicircle (mitochondrial) genome of Leishmania major (MHOM/SU/73/5ASKH), a much larger portion of this genome than has been reported previously from any Leishmania species infecting humans. This region contains 10 partial and complete genes: 5 protein-encoding genes (COII, COIII, ND1, ND7 and Cyt b); two ribosomal RNA subunits (12S and 9S) and three unidentified open reading frames (MURF1, MURF4 (ATPase6) and MURF5), as in the lizard-infecting species L. tarentolae. The genes from L. major exhibit 85-87% identity with those of L. tarentolae at the nucleotide level and 71-94% identity at the amino acid level. Most differences between sequences from the two species are transversions. The gene order and arrangement within the maxicircle of L. major are similar to those in L. tarentolae, but base composition and codon usage differ between the species. Codons assigned for initiation for protein-coding genes available for comparison are similar in five genes in the two species. Pre-editing was identified in some of the protein-coding genes. Short intergenic non-coding regions are also present in L. major as they are in L. tarentolae. Intergenic regions between 9S rRNA and MURF5, MURF1 and ND1 genes are G+C rich and considered to be extensive RNA editing regions. The RNA editing process is likely to be conserved in similar pattern in L. major as in L. tarentolae.

Identification of geographically distributed sub-populations of Leishmania (Leishmania) major by microsatellite analysis

BACKGROUND

Leishmania (Leishmania) major, one of the agents causing cutaneous leishmaniasis (CL) in humans, is widely distributed in the Old World where different species of wild rodent and phlebotomine sand fly serve as animal reservoir hosts and vectors, respectively. Despite this, strains of L. (L.) major isolated from many different sources over many years have proved to be relatively uniform. To investigate the population structure of the species highly polymorphic microsatellite markers were employed for greater discrimination among it's otherwise closely related strains, an approach applied successfully to other species of Leishmania.

RESULTS

Multilocus Microsatellite Typing (MLMT) based on 10 different microsatellite markers was applied to 106 strains of L. (L.) major from different regions where it is endemic. On applying a Bayesian model-based approach, three main populations were identified, corresponding to three separate geographical regions: Central Asia (CA); the Middle East (ME); and Africa (AF). This was congruent with phylogenetic reconstructions based on genetic distances. Re-analysis separated each of the populations into two sub-populations. The two African sub-populations did not correlate well with strains' geographical origin. Strains falling into the sub-populations CA and ME did mostly group according to their place of isolation although some anomalies were seen, probably, owing to human migration.

CONCLUSION

The model- and distance-based analyses of the microsatellite data exposed three main populations of L. (L.) major, Central Asia, the Middle East and Africa, each of which separated into two sub-populations. This probably correlates with the different species of rodent host.

Conserved repeats in the kinetoplast maxicircle divergent region of Leishmania sp. and Leptomonas seymouri

The maxicircle control region [also termed divergent region (DR)] composed of various repeat elements remains the most poorly studied part of the kinetoplast genome. Only three extensive DR sequences demonstrating no significant similarity were available for trypanosomatids (Leishmania tarentolae, Crithidia oncopelti, Trypanosoma brucei). Recently, extensive DR sequences have been obtained for Leishmania major and Trypanosoma cruzi. In this work we have sequenced DR fragments of Leishmania turanica, Leishmania mexicana, Leishmania chagasi and two monogenetic trypanosomatids Leptomonas seymouri and Leptomonas collosoma. With the emergence of the additional extensive sequences some conserved features of DR structure become evident. A conserved palindromic sequence has been revealed in the DRs of the studied Leishmania species, L. seymouri, and T. cruzi. The overall DR structure appears to be similar in all the Leishmania species, their relative L. seymouri, and T. brucei: long relatively GC-rich repeats are interspersed with clusters of short AT-rich repeats. C. oncopelti, L. collosoma, and T. cruzi have a completely different DR structure. Identification of conserved sequences and invariable structural features of the DR may further our understanding of the functioning of this important genome fragment.

Radically different maxicircle classes within the same kinetoplast: an artefact or a novel feature of the kinetoplast genome?

We discuss here some results which suggest that radically different maxicircle classes coexist within the same kinetoplast. These data, although tentative and incomplete, may provide a new outlook on the kinetoplast genome structure and expression.