The knowledge of the diversity of haemosporidian parasites is of primary importance as their representatives include agents of bird malaria. We investigated the occurrence of Haemoproteus spp. and Plasmodium spp. in bird populations from a single locality in the State of Selangor, Peninsular Malaysia, and report on the parasite prevalence of the two genera. A combination of methods (molecular and morphological) was used for detecting these parasites. Seventy-nine bird individuals were caught using mist-nets in July and August 2010 at Gombak Field Station of the University of Malaya, Kuala Lumpur. In total, 23 birds were identified as positive for Haemoproteus or Plasmodium infection and one individual was recognized as carrying mixed infection. The total prevalence of haemosporidians in the collected samples was 30.3%. Infections with parasites of the genus Haemoproteus were predominant compared to those of the genus Plasmodium. In total, 10 new cyt b lineages of Haemoproteus spp. and 3 new cyt b lineages of Plasmodium spp. were recorded in this study. From all recorded haemosporidian lineages (16 in total), 3 were known from previous studies - hCOLL2, hYWT2 and pNILSUN1. Two of them are linked with their corresponding morphospecies - Haemoproteus pallidus (COLL2) and Haemoproteus motacillae (YWT2). The morphological analysis in the present study confirmed the results obtained by the PCR method relative to prevalence, with 25.3% total prevalence of Haemoproteus and Plasmodium parasites. The intensities of infection varied between 0.01% and 19%. Most infections were light, with intensities below 0.1%. The present study is the first molecular survey of the protozoan blood parasites of the order Haemosporida recorded in Malaysia.
It has been discovered that Plasmodium knowlesi (P. knowlesi) is transmitted from macaque to man. Thus, the aim of the present study was to determine P. knowlesi genetic diversity in both human (n = 147) and long-tailed macaque (n = 26) samples from high- and low-endemicity localities. Genotyping was performed using seven neutral microsatellite loci markers. The size of the alleles, multiplicity of infection (MOI), mean number of alleles (Na), expected heterozygosity (HE), linkage disequilibrium (LD), and genetic differentiation (FST) were determined. In highly endemic P. knowlesi localities, the MOI for human and long-tailed macaque isolates was 1.04 and 1.15, respectively, while the Na was 11.14 and 7.86, respectively. Based on the allele frequency distribution for all loci, and with FST
There have been three major rabies epidemics in China since the 1950s. To gain more insights into the molecular epidemiology of rabies viruses (RVs) for the third (the current) epidemic, we isolated RV from dogs and humans in major endemic areas, and characterized these isolates genetically by sequencing the entire glycoprotein (G) gene and the G-L non-coding region. These sequences were also compared phylogenetically with RVs isolated in China during previous epidemics and those around the world. Comparison of the entire G genes among the Chinese isolates revealed up to 21.8% divergence at the nucleotide level and 17.8% at the amino acid level. The available Chinese isolates could be divided into two distinct clades, each of which could be further divided into six lineages. Viruses in clade I include most of the Chinese viruses as well as viruses from southeast Asian countries including Indonesia, Malaysia, the Philippines, Thailand, and Vietnam. The viruses in the other clade were found infrequently in China, but are closely related to viruses distributed worldwide among terrestrial animals. Interestingly, most of the viruses isolated during the past 10 years belong to lineage A viruses within clade I whereas most of the viruses isolated before 1996 belong to other lineages within clades I and II. Our results indicated that lineages A viruses have been predominant during the past 10 years and thus are largely responsible for the third and the current epidemic in China. Our results also suggested that the Chinese RV isolates in clade I share a common recent ancestor with those circulating in southeast Asia.
Plasmodium knowlesi, a malaria parasite of macaques, has emerged as an important parasite of humans. Despite the significance of P. knowlesi malaria in parts of Southeast Asia, very little is known about the genetic variation in this parasite. Our aim here was to explore sequence variation in a molecule called the 42kDa merozoite surface protein-1 (MSP-1), which is found on the surface of blood stages of Plasmodium spp. and plays a key role in erythrocyte invasion. Several studies of P. falciparum have reported that the C-terminus (a 42kDa fragment) of merozoite surface protein-1 (MSP-142; consisting of MSP-119 and MSP-133) is a potential candidate for a malaria vaccine. However, to date, no study has yet investigated the sequence diversity of the gene encoding P. knowlesi MSP-142 (comprising Pk-msp-119 and Pk-msp-133) among isolates in Malaysia. The present study explored this aspect. Twelve P. knowlesi isolates were collected from patients from hospitals in Selangor and Sabah Borneo, Malaysia, between 2012 and 2014. The Pk-msp-142 gene was amplified by PCR and directly sequenced. Haplotype diversity (Hd) and nucleotide diversity (л) were studied among the isolates. There was relatively high genetic variation among P. knowlesi isolates; overall Hd and л were 1±0.034 and 0.01132±0.00124, respectively. A total of nine different haplotypes related to amino acid alterations at 13 positions, and the Pk-MSP-119 sequence was found to be more conserved than Pk-msp-133. We have found evidence for negative selection in Pk-msp-42 as well as the 33kDa and 19kDa fragments by comparing the rate of non-synonymous versus synonymous substitutions. Future investigations should study large numbers of samples from disparate geographical locations to critically assess whether this molecule might be a potential vaccine target for P. knowlesi.
Helicobacter pylori is the most predominant bacterium in almost 50% of the world's population and colonization causes a persistent inflammatory response leading to chronic gastritis. It shows high genetic diversity and individuals generally harbour a distinct bacterial population. With the advancement of whole-genome sequencing technology, new H. pylori subpopulations have been identified that show admixture between various H. pylori strains. Genotypic variation of H. pylori may be related to the presence of virulence factors among strains and is associated with different outcomes of infection in different individuals. This review summarizes the genetic diversity in H. pylori strain populations and its virulence characteristics responsible for variable outcomes in different ethnic groups.
The species diversity and genetic structure of mosquitoes belonging to the Anopheles maculatus group in Southeast Asia were investigated using the internal transcribed spacer 2 (ITS2) of ribosomal DNA (rDNA). A molecular phylogeny indicates the presence of at least one hitherto unrecognised species. Mosquitoes of chromosomal form K from eastern Thailand have a unique ITS2 sequence that is 3.7% divergent from the next most closely related taxon (An. sawadwongporni) in the group. In the context of negligible intraspecific variation at ITS2, this suggests that chromosomal form K is most probably a distinct species. Although An. maculatus sensu stricto from northern Thailand and southern Thailand/peninsular Malaysia differ from each other in chromosomal banding pattern and vectorial capacity, no intraspecific variation was observed in the ITS2 sequences of this species over this entire geographic area despite an extensive survey. A PCR-based identification method was developed to distinguish five species of the group (An. maculatus, An. dravidicus, An. pseudowillmori, An. sawadwongporni and chromosomal form K) to assist field-based studies in northwestern Thailand. Sequences from 187 mosquitoes (mostly An. maculatus and An. sawadwongporni) revealed no intraspecific variation in specimens from Thailand, Cambodia, mainland China, Malaysia, Taiwan and Vietnam, suggesting that this identification method will be widely applicable in Southeast Asia. The lack of detectable genetic structure also suggests that populations of these species are either connected by gene flow and/or share a recent common history.
Plasmodium knowlesi is the leading cause of malaria in Malaysia. Serine Repeat Antigens (SERAs) have an essential role in the parasite life cycle. However, genetic characterization on P. knowlesi SERA3 Ag2 (PkSERA3 Ag2) is lacking. In the present study, nucleotide diversity, natural selection, and haplotypes of PkSERA3 Ag2 in clinical samples from Peninsular Malaysia and Malaysian Borneo were investigated. A total of 50 P. knowlesi clinical samples were collected from Peninsular Malaysia and Malaysian Borneo. The PkSERA3 Ag2 gene was amplified using PCR, and subsequently cloned and sequenced. Genetic diversity, haplotype, natural selection as well as genetic structure and differentiation of PkSERA3 Ag2 were analysed. In addition, in silico analyses were performed to identify repeat motifs, B-cell epitopes, and antigenicity indices of the protein. Analysis of 114 PkSERA3 Ag2 sequences revealed high nucleotide diversity of the gene in Malaysia. A codon-based Z-test indicated that the gene underwent purifying selection. Haplotype and population structure analyses identified two distinct PkSERA3 Ag2 clusters (K = 2, ΔK = 721.14) but no clear genetic distinction between PkSERA3 Ag2 from Peninsular Malaysia and Malaysian Borneo. FST index indicated moderate differentiation of the gene. In silico analyses revealed unique repeat motifs among PkSERA3 Ag2 isolates. Moreover, the amino acid sequence of PkSERA3 Ag2 exhibited potential B-cell epitopes and possessed high antigenicity indices. These findings enhance the understanding of PkSERA3 Ag2 gene as well as its antigenic properties. Further validation is necessary to ascertain the utility of PkSERA3 Ag2 as a serological marker for P. knowlesi infection.
Genetic characterization of Theileria species infecting bovines in India was attempted targeting the 18S ribosomal RNA region of the parasite. Blood samples of bovines (n = 452), suspected for haemoprotozoan infections, from 9 different states of the country were microscopically examined for Theileria species infection. Four Theileria spp. positive blood samples from each state were randomly utilized for PCR amplification of the 18S rRNA gene (approx. 1529 bp) followed by cloning and sequencing. The sequence data analysis of all the 36 isolates revealed that 33 isolates had high sequence similarity with published sequences of T. annulata, whereas 3 isolates (MF287917, MF287924 and MF287928) showed close similarity with published sequences of T. orientalis. Sequence homology within the isolates ranged between 95.8 and 100% and variation in the length of targeted region was also noticed in different isolates (1527-1538 nt). Phylogenetic tree created for T. annulata sequences revealed that a total of 24 Indian isolates formed a major clade and grouped together with isolates originating from countries like China, Spain, Turkey and USA. Remaining 09 isolates clustered in a separate group and were closely related to the TA5 isolate of T. annulata (a new genotype) originating from India and also with the isolates from East Asian countries like Japan and Malaysia. All the three T. orientalis isolates had minimal intraspecific variation (99-100% homology) amongst themselves. Further, in the phylogenetic analysis T. orientalis Indian isolates were found to cluster away from other 14 isolates of T. buffeli/sergenti/orientalis originating from different countries (Australia, China, Indonesia and Spain). However, these 3 isolates clustered together with the T. buffeli Indian isolate (EF126184). Present study confirmed the circulation of different genotypes of T. annulata in India, along with T. orientalis isolates.
Getah virus (GETV) was first isolated in Malaysia in 1955. Since then, epidemics in horses and pigs caused by GETV have resulted in huge economic losses. At present, GETV has spread across Eurasia and Southeast Asia, including mainland China, Korea, Japan, Mongolia, and Russia. Data show that the Most Recent Common Ancestor (MRCA) of GETV existed about 145years ago (95% HPD: 75-244) and gradually evolved into four distinct evolutionary populations: Groups I-IV. The MRCA of GETVs in Group III, which includes all GETVs isolated from mosquitoes, pigs, horses, and other animals since the 1960s (from latitude 19°N to 60°N), existed about 51years ago (95% HPD: 51-72). Group III is responsible for most viral epidemics among domestic animals. An analysis of the GETV E2 protein sequence and structure revealed seven common amino acid mutation sites. These sites are responsible for the structural and electrostatic differences detected between widespread Group III isolates and the prototype strain MM2021. These differences may account for the recent geographical radiation of the virus. Considering the economic significance of GETV infection in pigs and horses, we recommend the implementation of strict viral screening and monitoring programs.
Given the HIV epidemic in Malaysia, genetic information on opportunistic pathogens, such as Cryptosporidium and Giardia, in HIV/AIDS patients is pivotal to enhance our understanding of epidemiology, patient care, management and disease surveillance. In the present study, 122 faecal samples from HIV/AIDS patients were examined for the presence of Cryptosporidium oocysts and Giardia cysts using a conventional coproscopic approach. Such oocysts and cysts were detected in 22.1% and 5.7% of the 122 faecal samples, respectively. Genomic DNAs from selected samples were tested in a nested-PCR, targeting regions of the small subunit (SSU) of nuclear ribosomal RNA and the 60kDa glycoprotein (gp60) genes (for Cryptosporidium), and the triose-phosphate isomerase (tpi) gene (for Giardia), followed by direct sequencing. The sequencing of amplicons derived from SSU revealed that Cryptosporidium parvum was the most frequently detected species (64% of 25 samples tested), followed by C. hominis (24%), C. meleagridis (8%) and C. felis (4%). Sequencing of a region of gp60 identified C. parvum subgenotype IIdA15G2R1 and C. hominis subgenotypes IaA14R1, IbA10G2R2, IdA15R2, IeA11G2T3R1 and IfA11G1R2. Sequencing of amplicons derived from tpi revealed G. duodenalis assemblage A, which is of zoonotic importance. This is the first report of C. hominis, C. meleagridis and C. felis from Malaysian HIV/AIDS patients. Future work should focus on an extensive analysis of Cryptosporidium and Giardia in such patients as well as in domestic and wild animals, in order to improve the understanding of transmission patterns and dynamics in Malaysia. It would also be particularly interesting to establish the relationship among clinical manifestation, CD4 cell counts and genotypes/subgenotypes of Cryptosporidium and Giardia in HIV/AIDS patients. Such insights would assist in a better management of clinical disease in immuno-deficient patients as well as improved preventive and control strategies.
The performances of a commonly used Plasmodium falciparum-detecting rapid diagnostic test (RDT) were determined in symptomatic individuals living in Cameroon. Discrepancies between RDT and light microscopy (LM) results were further investigated, with a focus on non-falciparum malaria (NFM) which are still largely understudied in sub-Saharan Africa (sSA) countries. In the present study, a total of 355 individuals aged 1-65 years were enrolled in the study. Their signs/symptoms and sociodemographic characteristics were documented. The RDT reliability was evaluated using LM as gold standard method. Polymerase chain reaction (PCR) of Plasmodium 18S gene was performed for samples with discordant results between LM and RDT (i.e., RDT-/LM+, and RDT+/LM-). The PCR amplicons of NFM species were sequenced and BLASTed. The prevalence of malaria infection by LM was 95.7% (95% CI: 93.1-97.4%). The sensitivity and specificity of RDT for P. falciparum detection was 94.0% and 66.7%, respectively. By PCR assay, P. ovale curtisi (PoC) was found in 5 of the 30 discordant samples, and on sequence analysis these isolates were found to be phylogenetically closer to sequences reported from China-Myanmar border and Malaysia. This is the first report on molecular characterization of P. ovale subspecies in Cameroon. The study also outlines the good diagnostic performances of the RDT for detection of P. falciparum. Though, the presence of PoC indicated the importance of having RDTs targeting the NFM species in malaria diagnosis and treatment, which is presently limited in the country.
Dengue virus type 4 (DENV-4) circulates in tropical and subtropical countries from Asia and the Americas. Despite the importance of dengue virus distribution, little is known about the worldwide viral spread. Following a Bayesian phylogenetic approach we inferred the evolutionary history of 310 isolates sampled from 37 countries during the time period 1956-2008 and the spreading dynamics for genotypes I and II. The region (tropical rainforest biome) comprised by Malaysia-Thailand was the most likely ancestral area from which the serotype has originated and spread. Interestingly, cross-correlation analysis on demographic time series with the Asian sequences showed a statistically significant negative correlation that could be suggestive of competition among genotypes within the same serotype.
The emergence of a new coronavirus, in around late December 2019 which had first been reported in Wuhan, China has now developed into a massive threat to global public health. The World Health Organization (WHO) has named the disease caused by the virus as COVID-19 and the virus which is the culprit was renamed from the initial novel respiratory 2019 coronavirus to SARS-CoV-2. The person-to-person transmission of this virus is ongoing despite drastic public health mitigation measures such as social distancing and movement restrictions implemented in most countries. Understanding the source of such an infectious pathogen is crucial to develop a means of avoiding transmission and further to develop therapeutic drugs and vaccines. To identify the etiological source of a novel human pathogen is a dynamic process that needs comprehensive and extensive scientific validations, such as observed in the Middle East respiratory syndrome (MERS), severe acute respiratory syndrome (SARS), and human immunodeficiency virus (HIV) cases. In this context, this review is devoted to understanding the taxonomic characteristics of SARS-CoV-2 and HIV. Herein, we discuss the emergence and molecular mechanisms of both viral infections. Nevertheless, no vaccine or therapeutic drug is yet to be approved for the treatment of SARS-CoV-2, although it is highly likely that new effective medications that target the virus specifically will take years to establish. Therefore, this review reflects the latest repurpose of existing antiviral therapeutic drug choices available to combat SARS-CoV-2.
Indonesia is one of the biggest dengue endemic countries, and, thus, is an important place to investigate the evolution of dengue virus (DENV). We have continuously isolated DENV in Surabaya, the second biggest city in Indonesia, since 2008. We previously reported sequential changes in the predominant serotype from DENV type 2 (DENV-2) to DENV type 1 (DENV-1) in November 2008 and from DENV-1 to DENV-2 in July 2013. The predominance of DENV-2 continued in 2014, but not in 2015. We herein phylogenetically investigated DENV-2 transitions in Surabaya between 2008 and 2014 to analyze the divergence and evolution of DENV-2 concomitant with serotype shifts. All DENV-2 isolated in Surabaya were classified into the Cosmopolitan genotype, and further divided into 6 clusters. Clusters 1-3, dominated by Surabaya strains, were defined as the "Surabaya lineage". Clusters 4-6, dominated by strains from Singapore, Malaysia, and many parts of Indonesia, were the "South East Asian lineage". The most recent common ancestor of these strains existed in 1988, coinciding with the time that an Indonesian dengue outbreak took place. Cluster 1 appeared to be unique because no other DENV-2 isolate was included in this cluster. The predominance of DENV-2 in 2008 and 2013-14 were caused by cluster 1, whereas clusters 2 and 3 sporadically emerged in 2011 and 2012. The characteristic amino acids of cluster 1, E-170V and E-282Y, may be responsible for its prevalence in Surabaya. No amino acid difference was observed in the envelope region between strains in 2008 and 2013-14, suggesting that the re-emergence of DENV-2 in Surabaya was due to the loss or decrease of herd immunity in the 5-year period when DENV-2 subsided. The South East Asian lineage primarily emerged in Surabaya in 2014, probably imported from other parts of Indonesia or foreign countries.
Dengue is hyperendemic in most of Southeast Asia. In this region, all four dengue virus serotypes are persistently present. Major dengue outbreak cycle occurs in a cyclical pattern involving the different dengue virus serotypes. In Malaysia, since the 1980s, the major outbreak cycles have involved dengue virus type 3 (DENV3), dengue virus type 1 (DENV1) and dengue virus type 2 (DENV2), occurring in that order (DENV3/DENV1/DENV2). Only limited information on the DENV3 cycles, however, have been described. In the current study, we examined the major outbreak cycle involving DENV3 using data from 1985 to 2016. We examined the genetic diversity of DENV3 isolates obtained during the period when DENV3 was the dominant serotype and during the inter-dominant transmission period. Results obtained suggest that the typical DENV3/DENV1/DENV2 cyclical outbreak cycle in Malaysia has recently been disrupted. The last recorded major outbreak cycle involving DENV3 occurred in 2002, and the expected major outbreak cycle involving DENV3 in 2006-2012 did not materialize. DENV genome analyses revealed that DENV3 genotype II (DENV3/II) was the predominant DENV3 genotype (67%-100%) recovered between 1987 and 2002. DENV3 genotype I (DENV3/I) emerged in 2002 followed by the introduction of DENV3 genotype III (DENV3/III) in 2008. These newly emerged DENV3 genotypes replaced DENV3/II, but there was no major upsurge of DENV3 cases that accompanied the emergence of these viruses. DENV3 remained in the background of DENV1 and DENV2 until now. Virus genome sequence analysis suggested that intrinsic differences within the different dengue virus genotypes could have influenced the transmission efficiency of DENV3. Further studies and continuous monitoring of the virus are needed for better understanding of the DENV transmission dynamics in hyperendemic regions.
From 1989 to 2011 in Kuala Lumpur, Malaysia, multiple genotypes from both respiratory syncytial virus (RSV) subgroups were found co-circulating each year. RSV-A subgroup predominated in 12 out of 17years with the remaining years predominated by RSV-B subgroup. Local RSV strains exhibited temporal clustering with RSV strains reported in previous epidemiological studies. Every few years, the existing predominant genotype was replaced by a new genotype. The RSV-A genotypes GA2, GA5 and GA7 were replaced by NA1 and NA2, while BA became the predominant RSV-B genotype. A unique local cluster, BA12, was seen in 2009, and the recently-described ON1 genotype with 72-nt duplication emerged in 2011. Our findings will have important implications for future vaccine intervention.
Fasciola gigantica and hybrid Fasciola flukes, responsible for the disease fasciolosis, are found in Southeast Asian countries. In the present study, we performed molecular species identification of Fasciola flukes distributed in Terengganu, Malaysia using multiplex PCR for phosphoenolpyruvate carboxykinase (pepck) and PCR-restriction fragment length polymorphism (RFLP) for DNA polymerase delta (pold). Simultaneously, phylogenetic analysis based on mitochondrial NADH dehydrogenase subunit 1 (nad1) was performed for the first time on Malaysian Fasciola flukes to infer the dispersal direction among neighboring countries. A total of 40 flukes used in this study were identified as F. gigantica. Eight nad1 haplotypes were identified in the F. gigantica population of Terengganu. Median-joining network analysis revealed that the Malaysian population was related to those obtained from bordering countries such as Thailand and Indonesia. However, genetic differentiation was detected using population genetics analyses. Nevertheless, the nucleotide diversity (π) value suggested that F. gigantica with the predominant haplotypes was introduced into Malaysia from Thailand and Indonesia. The dispersal direction suggested by population genetics in the present study may not be fully reliable since Fasciola flukes were collected from a single location in one state of Malaysia. Further studies analyzing more samples from many locations are required to validate the dispersal direction proposed herein.
Shigella is an intracellular bacterial pathogen that causes bacterial dysentery called shigellosis. The assessment of pro- and anti-inflammatory mediators produced by immune cells against this bacteria are vital in identifying the effectiveness of the immune reaction in protecting the host. In Malaysia, Shigella is ranked as the third most common bacteria causing diarrheal disease among children below 5 years old. In the present study, we aim to examine the differential cytokine gene expressions of macrophages in response to two types of clinical strains of Shigella flexneri 2a (S. flexneri 2a) isolated from patients admitted in Hospital Universiti Sains Malaysia, Kelantan, Malaysia. THP-1-derived macrophages, as the model of human macrophages, were infected separately with S. flexneri 2a mild (SH062) and virulence (SH057) strains for 6, 12, and 24 h, respectively. The gene expression level of inflammatory mediators was identified using real-time quantitative polymerase chain reaction (RT-qPCR). The production of nitric oxide (NO) by the macrophages was measured by using a commercialized NO assay kit. The ability of macrophages to kill the intracellular bacteria was assessed by intracellular killing assay. Induction of tumor necrosis factor-alpha (TNFα), interleukin (IL)-1β, IL-6, IL-12, inducible NO synthase (iNOS), and NO, confirmed the pro-inflammatory reaction of the THP-1-derived macrophages in response to S. flexneri 2a, especially against the SH507 strain. The SH057 also induced a marked increase in the expression levels of the anti-inflammatory cytokine mRNAs at 12 h and 24 h post-infection. In the intracellular killing assay, both strains showed less viable, indicating the generation of pro-inflammatory cytokines in the presence of iNOS and NO was crucial in the stimulation of macrophages for the host defense against shigellosis. Transcription analysis of THP-1-derived macrophages in this study identifies differentially expressed cytokine genes that correlated with the virulence factor of S. flexneri 2a.
Shigellosis is one of the most common diseases found in the developing countries, especially those countries that are prone flood. The causative agent for this disease is the Shigella species. This organism is one of the third most common enteropathogens responsible for childhood diarrhea. Since Shigella can survive gastric acidity and is an intracellular pathogen, it becomes difficult to treat. Also, uncontrolled use of antibiotics has led to development of resistant strains which poses a threat to public health. Therefore, there is a need for long term control of Shigella infection which can be achieved by designing a proper and effective vaccine. In this study, emphasis was made on designing a candidate that could elicit both B-cell and T-cell immune response. Hence B- and T-cell epitopes of outer membrane channel protein (OM) and putative lipoprotein (PL) from S. flexneri 2a were computationally predicted using immunoinformatics approach and a chimeric construct (chimeric-OP) containing the immunogenic epitopes selected from OM and PL was designed, cloned and expressed in E. coli system. The immunogenicity of the recombinant chimeric-OP was assessed using Shigella antigen infected rabbit antibody. The result showed that the chimeric-OP was a synthetic peptide candidate suitable for the development of vaccine and immunodiagnostics against Shigella infection.
Acute hepatopancreatic necrosis disease (AHPND) is an important shrimp disease of economic importance which causes mass mortality of cultivated penaeid shrimps in Southeast Asian countries, Mexico and South America. This disease was originally caused by Vibrio parahaemolyticus (VPAHPND) which is reported to harbour a transferable plasmid carrying the virulent PirAB-like toxin genes (pirABvp). However, little is known about the pathogenicity of VPAHPND. To extend our understanding, comparative genomic analyses was performed in this study to identify the genetic differences and to understand the phylogenetic relationship of VPAHPND strains. Seven Vibrio parahaemolyticus strains (five VPAHPND strains and two non-VPAHPND strains) were sequenced and 31 draft genomes of V. parahaemolyticus were retrieved from NCBI database and incorporated into the genomic comparison to elucidate their genomic diversity. The study showed that the genome sizes of the VPAHPND strains were approximately 5 Mbp. Ten sequence types (STs) were identified among the VPAHPND strains using in silico-Multilocus Sequence Typing analysis (MLST) and ST 970 was the predominant ST. Phylogenetic analysis based on MLST and single nucleotide polymorphisms (SNP) showed that the VPAHPND strains were genetically diverse. Based on the comparative genomic analysis, several functional proteins were identified from diiferent categories associated with virulence-related proteins, secretory proteins, conserved domain proteins, transporter proteins, and phage proteins. The CRISPR analysis showed that VPAHPND strains contained less number of CRISPRs elements than non-VPAHPND strains while six prophages regions were identified in the genomes, suggested the lack of CRISPR might promote prophage insertion. The genomic information in this study provide improved understanding of the virulence of these VPAHPND strains.