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

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

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

OmpA protein sequence-based typing and virulence-associated gene profiles of Pasteurella multocida isolates associated with bovine haemorrhagic septicaemia and porcine pneumonic pasteurellosis in Thailand

Abstract Background Pasteurella multocida is a Gram-negative bacterium that causes economically significant infections of a broad range of animal species. Pneumonic and septicaemic pasteurellosis caused by this bacterium remain important problems in pigs, cattle, and water buffaloes in Thailand. The aim of this study was to characterise the virulence-associated gene profiles and to develop an OmpA molecular typing scheme for classifying 191 bovine and porcine isolates of P. multocida collected between 1989 and 2012 in Thailand using polymerase chain reactions (PCRs), nucleotide sequencing, and sequence and structural bioinformatics analyses. Results PCR screening successfully characterised the profiles of 25 virulence-associated genes in all isolates. The gene profiles separated these isolates into bovine and porcine clusters based on eight genes (hgbB, hsf1, tadD, nanH, pfhA, plpE, pmHAS, and tbpA). Phylogenetic analyses of the nucleotide and protein sequences corresponding to the ompA gene, which encodes a major outer membrane surface protein, showed two major bovine and porcine clusters. Structural prediction and analysis of the dN/dS ratio revealed four hypervariable extracellular loops of the OmpA transmembrane domains. These four loops were used to develop an OmpA typing scheme. This scheme classified 186 isolates into five major loop sequence types (LST8, LST12, LST15, LST18, and LST19), consistent with the phylogenetic results. The loop regions of the bovine isolates were predicted to be more antigenic than those of the porcine isolates. Thus, molecular evolution of the OmpA proteins could be used to classify P. multocida isolates into different capsular types, host types, and, possibly, pathogenicity levels. Conclusions Together with the virulence-associated gene profiles, the typing reported in this work provides a better understanding of P. multocida virulence. Effective monitoring and potential strain-specific subunit vaccines could be developed based on these loop oligopeptides

Comparative draft genomes of Leishmania orientalis isolate PCM2 (formerly named Leishmania siamensis ) and Leishmania martiniquensis isolate PCM3 from the southern province of Thailand

(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

Additional file 2: Table S2. of OmpA protein sequence-based typing and virulence-associated gene profiles of Pasteurella multocida isolates associated with bovine haemorrhagic septicaemia and porcine pneumonic pasteurellosis in Thailand

Information of 191 P. multocida isolates associated with diseased pigs, cattle and water buffaloes in Thailand collected from 1989 to 2012. These isolates were collected and maintained by the National Institute of Animal Health, Department of Livestock, Ministry of Agriculture and Cooperatives. Each isolate was named with two-letter code. (DOCX 38 kb

Additional file 1: Table S1. of OmpA protein sequence-based typing and virulence-associated gene profiles of Pasteurella multocida isolates associated with bovine haemorrhagic septicaemia and porcine pneumonic pasteurellosis in Thailand

Clustering summary of 186 bovine and porcine isolates of P. multocida collected from Thailand and 41 strains from the NCBI database. Three clustering methods were shown including nucleotide and protein phylogenetic analyses, and the OmpA protein sequence-based typing using four extracellular loop sequence types (LSTs) of the transmembrane domain. The LST types in the last column were created by the combination of different types of the four extracellular loops 1, 2, 3, and 4. DNA and protein cluster codes correlated with those in Fig. 2a and b. (DOCX 65 kb