Chromosome-scale genome sequencing, assembly and annotation of six genomes from subfamily Leishmaniinae

TriTrypDB: An integrated functional genomics resource for kinetoplastida

Parasitic diseases caused by kinetoplastid parasites are a burden to public health throughout tropical and subtropical regions of the world. TriTrypDB (https://tritrypdb.org) is a free online resource for data mining of genomic and functional data from these kinetoplastid parasites and is part of the VEuPathDB Bioinformatics Resource Center (https://veupathdb.org). As of release 59, TriTrypDB hosts 83 kinetoplastid genomes, nine of which, including Trypanosoma brucei brucei TREU927, Trypanosoma cruzi CL Brener and Leishmania major Friedlin, undergo manual curation by integrating information from scientific publications, high-throughput assays and user submitted comments. TriTrypDB also integrates transcriptomic, proteomic, epigenomic, population-level and isolate data, functional information from genome-wide RNAi knock-down and fluorescent tagging, and results from automated bioinformatics analysis pipelines. TriTrypDB offers a user-friendly web interface embedded with a genome browser, search strategy system and bioinformatics tools to support custom in silico experiments that leverage integrated data. A Galaxy workspace enables users to analyze their private data (e.g., RNA-sequencing, variant calling, etc.) and explore their results privately in the context of publicly available information in the database. The recent addition of an annotation platform based on Apollo enables users to provide both functional and structural changes that will appear as 'community annotations' immediately and, pending curatorial review, will be integrated into the official genome annotation.

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.

Making the Most of Its Short Reads: A Bioinformatics Workflow for Analysing the Short-Read-Only Data of Leishmania orientalis (Formerly Named Leishmania siamensis) Isolate PCM2 in Thailand

Simple Summary

Leishmaniasis is a parasitic disease caused by flagellated protozoa of the genus Leishmania. Multiple genome sequencing platforms have been employed to complete Leishmania genomes at the expense of high cost. This study proposes an integrative bioinformatic workflow for assembling only the short-read data of Leishmania orientalis isolate PCM2 from Thailand and produce an acceptable-quality genome for further genomic analysis. This workflow gives extensive information required for identifying strain-specific markers and virulence-associated genes useful for drug and vaccine development before a more exhaustive and expensive investigation.

Abstract

Background: Leishmania orientalis (formerly named Leishmania siamensis) has been neglected for years in Thailand. The genomic study of L. orientalis has gained much attention recently after the release of the first high-quality reference genome of the isolate LSCM4. The integrative approach of multiple sequencing platforms for whole-genome sequencing has proven effective at the expense of considerably expensive costs. This study presents a preliminary bioinformatic workflow including the use of multi-step de novo assembly coupled with the reference-based assembly method to produce high-quality genomic drafts from the short-read Illumina sequence data of L. orientalis isolate PCM2. Results: The integrating multi-step de novo assembly by MEGAHIT and SPAdes with the reference-based method using the L. enriettii genome and salvaging the unmapped reads resulted in the 30.27 Mb genomic draft of L. orientalis isolate PCM2 with 3367 contigs and 8887 predicted genes. The results from the integrated approach showed the best integrity, coverage, and contig alignment when compared to the genome of L. orientalis isolate LSCM4 collected from the northern province of Thailand. Similar patterns of gene ratios and frequency were observed from the GO biological process annotation. Fifty GO terms were assigned to the assembled genomes, and 23 of these (accounting for 61.6% of the annotated genes) showed higher gene counts and ratios when results from our workflow were compared to those of the LSCM4 isolate. Conclusions: These results indicated that our proposed bioinformatic workflow produced an acceptable-quality genome of L. orientalis strain PCM2 for functional genomic analysis, maximising the usage of the short-read data. This workflow would give extensive information required for identifying strain-specific markers and virulence-associated genes useful for drug and vaccine development before a more exhaustive and expensive investigation.

Leishmania enriettii visceralises in the trachea, lungs, and spleen of Cavia porcellus

BACKGROUND

Leishmania (Mundinia) enriettii is a species commonly found in the guinea pig, Cavia porcellus. Although it is a dermotropic species, there is still an uncertainty regarding its ability to visceralise during Leishmania life cycle.

OBJECTIVE

Here, we investigated the ability of L. enriettii (strain L88) to visceralise in lungs, trachea, spleen, and liver of C. porcellus, its natural vertebrate host.

METHODS

Animals were infected sub-cutaneously in the nose and followed for 12 weeks using histological (hematoxilin-eosin) and molecular tools (polymerase chain reaction-restriction fragment length polymorphism - PCR-RFLP). To isolate parasite from C. porcellus, animals were experimentally infected for viscera removal and PCR typing targeting hsp70 gene.

FINDINGS

Histological analysis revealed intense and diffuse inflammation with the presence of amastigotes in the trachea, lung, and spleen up to 12 weeks post-infection (PI). Molecular analysis of paraffin-embedded tissues detected parasite DNA in the trachea and spleen between the 4th and 8th weeks PI. At the 12th PI, no parasite DNA was detected in any of the organs. To confirm that the spleen could serve as a temporary site for L. enriettii, we performed additional in vivo experiments. During 6th week PI, the parasite was isolated from the spleen confirming previous histopathological and PCR observations.

MAIN CONCLUSION

Leishmania enriettii (strain L88) was able to visceralise in the trachea, lung, and spleen of C. porcellus.

Comparative Draft Genomes of Leishmania orientalis Isolate PCM2 (Formerly Named Leishmania siamensis) and Leishmania martiniquensis Isolate PCM3 from the Southern Province of Thailand

Simple Summary

This study successfully sequenced the draft genomes of the southern isolates of Leishmania orientalis and Leishmania martiniquensis in Thailand. The comparison with the genomes of the northern isolates revealed species-level similarity with a level of genome and proteome variation, suggesting the different emerging strains. Comparing the proteins of these southern strains to those of the northern ones and 14 other Leishmania species showed six protein groups with numerous unique proteins: 53 for the southern strain PCM2 of L. orientalis and 97 for the strain PCM3 of L. martiniquensis. Some of these proteins were related to virulence, drug resistance, drug target, and stress response, which could be targeted for further experimental characterization. Therefore, the findings could initiate further genetic and population genomic investigation, and the close monitoring of L. orientalis and L. martiniquensis in Thailand and neighboring regions.

Abstract

(1) Background: Autochthonous leishmaniasis, a sandfly-borne disease caused by the protozoan parasites Leishmania orientalis (formerly named Leishmania siamensis) and Leishmania martiniquensis, has been reported for immunocompromised and immunocompetent patients in the southern province of Thailand. Apart from the recent genomes of the northern isolates, limited information is known on the emergence and genetics of these parasites. (2) Methods: This study sequenced and compared the genomes of L. orientalis isolate PCM2 and L. martiniquensis isolate PCM3 with those of the northern isolates and other 14 Leishmania species using short-read whole-genome sequencing methods and comparative bioinformatic analyses. (3) Results: The genomes of the southern isolates of L. orientalis and L. martiniquensis were 30.01 Mbp and 32.39 Mbp, and the comparison with the genomes of the northern isolates revealed species-level similarity with a level of genome and proteome variation, suggesting the different strains. Comparative proteome analysis showed six protein groups with 53 unique proteins for the strain PCM2 and 97 for the strain PCM3. Certain proteins were related to virulence, drug resistance, and stress response. (4) Conclusion: Therefore, the findings could indicate the need for more genetic and population genomic investigation, and the close monitoring of L. orientalis and L. martiniquensis in Thailand and neighboring regions.