To gain insight into the microbiological safety of food products routinely traded across international borders in Southeast Asian countries, beef imported from Malaysia to southern Thailand was examined for contamination with Escherichia coli O157 and its subsequent spread into the imported areas. We screened 31 samples exported from Malaysia and 36 domestic Thai samples. Isolation methods including an O157 antigen-targeted immunomagnetic separation technique, screening on CHROMagar O157 medium, and serotype confirmation of E. coli isolates by specific agglutination tests were employed. Fourteen strains of E. coli O157:H7 were isolated from eight Malaysian samples (25.8%) and six strains from four Thai samples (11.1%). These strains were of the stx1-stx2+eae+ genotype except one Malaysian strain which was of the stx1- stx2- eae+ genotype. All 19 O157:H7 strains possessing the stx2 gene produced little or no Stx2 (reversed passive latex agglutination titer ≤ 4). Of the 19 strains, five Malaysian (38.5%) and two Thai (33.3%) strains exhibited resistance to a set of antibiotics. Finally, the results of two DNA fingerprinting
analyses (O157 IS-printing targeted to IS629 and pulsed-field gel electrophoresis, PFGE) of the O157:H7 strains possessing the stx2 gene, indicated that the Malaysian and Thai strains are closely related. Therefore, E. coli O157:H7 might be transferred from Malaysia to southern Thailand through beef trade.
Massively parallel sequencing (MPS) can identify sequence variation within short tandem repeat (STR) alleles as well as their nominal allele lengths that traditionally have been obtained by capillary electrophoresis. Using the MiSeq FGx Forensic Genomics System (Illumina), STRait Razor, and in-house excel workbooks, genetic variation was characterized within STR repeat and flanking regions of 27 autosomal, 7 X-chromosome and 24 Y-chromosome STR markers in 777 unrelated individuals from four population groups. Seven hundred and forty six autosomal, 227 X-chromosome, and 324 Y-chromosome STR alleles were identified by sequence compared with 357 autosomal, 107 X-chromosome, and 189 Y-chromosome STR alleles that were identified by length. Within the observed sequence variation, 227 autosomal, 156 X-chromosome, and 112 Y-chromosome novel alleles were identified and described. One hundred and seventy six autosomal, 123 X-chromosome, and 93 Y-chromosome sequence variants resided within STR repeat regions, and 86 autosomal, 39 X-chromosome, and 20 Y-chromosome variants were located in STR flanking regions. Three markers, D18S51, DXS10135, and DYS385a-b had 1, 4, and 1 alleles, respectively, which contained both a novel repeat region variant and a flanking sequence variant in the same nucleotide sequence. There were 50 markers that demonstrated a relative increase in diversity with the variant sequence alleles compared with those of traditional nominal length alleles. These population data illustrate the genetic variation that exists in the commonly used STR markers in the selected population samples and provide allele frequencies for statistical calculations related to STR profiling with MPS data.
On December 1st, 2014, the sinking of Oryong 501 occurred in the Bering Sea off the east coast of Russia. A total of 60 crew members, including 35 Indonesians, 13 Filipinos, 11 South Koreans and 1 Russian inspector were on board out of which only seven survived. Through an international rescue operation, the dead bodies of 27 were found and the remaining 26 crew are still missing. After transferring the dead bodies to the Busan Harbor in South Korea, the operation to identify the deceased began involving DVI teams from three countries: Korea, Indonesia and the Philippines. When a deep sea fishing boat sinks, it is very difficult to obtain antemortem data of the crew who had been on board for a long time. This is especially so if the crews are multinational. Further, the accuracy of the antemortem data provided by the families may be questionable, and the provided data is often not standardized. Despite the fact that the antemortem data were received in different formats, the identification process for the bodies of the 27 crew from the Oryong sinking was quickly completed through the cooperation among the three DVI teams. This case is an excellent example of how efficiently a DVI operation can be conducted in the Asia Pacific region. Issues raised during this operation should enable even better preparation for similar events in the future.
Enterobacteriaceae is a large family within the Gram-negative bacteria that primarily inhabits in the gastrointestinal tract of human and animals. The bacteria within this group are readily survived in the environment with some species found living free in the water where energy sources are scarce, making them ideal indicators for faecal contamination of the river water. Some species within the family have been used as indicator for the presence of pathogenic bacteria whilst on the other hand some species have been directly associated with various diseases in human and animals. The main aim of this research study was to determine the distribution and characteristics of the Enterobacteriaceae in water samples collected from river and waterfalls within a community resort. The health risk associated with the bacteria was analysed with regard to their susceptibility to antibiotics. Samples were collected from surface water and water falling down directly from waterfalls of river within the community resort. The samples collected were plated onto Eosine Methylene Blue agar (EMBA) for the isolation of the Enterobacteriaceae. Bacterial colonies growing on the agar were randomly picked, purified, stocked and then identified using API 20E identification kit. DNA fingerprinting using (GTG)5-PCR was utilised to determine their genetic profiles before the isolates were grouped into a dendrogram using RAPDistance software package. The level of antibiotic susceptibility of the bacteria isolates was analysed using disc diffusion technique. This study confirmed the presence of Enterobacter, Klebsiella, Citrobacter, Pantoea and Serratia in the water samples with their single and multiple antibiotic resistance and susceptible characteristics. The dendrogram presented in this study shows genetic similarities and differences among the strains, suggesting while there is a potential for single distribution of a clone, there is also possibility of the distribution of different strains within species in the water environment. Therefore, awareness on the potential risk associated with genetically diverse intermediate and resistant enteric bacteria in the recreational water should be communicated to the public especially communities within the study area.
Serratia marcescens is a well-known cause of nosocomial infections and outbreaks, particularly in immunocompromised patients with severe underlying disease. An outbreak due to S. marcescens infection was detected from 13 to 22 February 2001 at the intensive care unit (ICU) of our institution. We used pulsed-field gel electrophoresis (PFGE) typing to analyse the outbreak strains involved.
As a new crop in Malaysia, forty-four Bambara groundnut (Vigna subterranea L. verdc.) genotypes were sampled from eleven distinct populations of different origins to explore the genetic structure, genetic inconsistency, and fixation index. The Bambara groundnut, an African underutilized legume, has the capacity to boost food and nutrition security while simultaneously addressing environmental sustainability, food availability, and economic inequalities. A set of 32 ISSRs were screened out of 96 primers based on very sharp, clear, and reproducible bands which detected a total of 510 loci with an average of 97.64% polymorphism. The average calculated value of PIC = 0.243, RP = 5.30, H = 0.285, and MI = 0.675 representing the efficiency of primer set for genetic differentiation among the genotypes. The ISSR primers revealed the number of alleles (Na = 1.97), the effective number of alleles (Ne = 1.38), Nei's genetic diversity (h = 0.248), and a moderate level of gene flow (Nm = 2.26) across the genotypes studied. The estimated Shannon's information index (I = 0.395) indicates a high level of genetic variation exists among the accessions. Based on Nei's genetic dissimilarity a UPMGA phylogenetic tree was constructed and grouped the entire genotypes into 3 major clusters and 6 subclusters. PCA analysis revealed that first principal component extracted maximum variation (PC1 = 13.92%) than second principal component (PC2 = 12.59%). Bayesian model-based STRUCTURE analysis assembled the genotypes into 3 (best ΔK = 3) genetic groups. The fixation-index (Fst) analysis narrated a very great genetic diversity (Fst = 0.19 to 0.40) exists within the accessions of these 3 clusters. This investigation specifies the effectiveness of the ISSR primers system for the molecular portrayal of V. subterranea genotypes that could be used for genetic diversity valuation, detection, and tagging of potential genotypes with quick, precise, and authentic measures for this crop improvement through effective breeding schemes.
A total of 422 Mycobacterium tuberculosis isolates from eight countries were subjected to IS6110 and IS1081 DNA fingerprinting by means of restriction fragment analysis to characterize M. tuberculosis strains from each country. Chinese, Mongolian, Hong Kong, Filipino, and Korean isolates had comparatively more copies of IS6110 (proportion with eight or more copies; 95% +/- 5%), while Thai, Malaysian, and Vietnamese isolates had fewer copies (proportion with eight or more copies, 60% +/- 4%). We found a number of novel IS1081 types in this study. One IS1081 type was present in 56% of Filipino isolates, had a specific 6.6-kb PvuII fragment in its IS6110 DNA fingerprint, and was termed the "Filipino family." The IS1081 types of Thai isolates had interposing characteristics between the characteristics of northeastern Asian and southeastern Asian IS1081 types. A 1.3-kb single-copy IS6110 fragment was found only in Vietnamese M. tuberculosis isolates. Although M. tuberculosis isolates from each country had comparatively similar characteristics depending on the classification factor, each country's isolates showed characteristic DNA fingerprints and differed slightly from the isolates from the other countries in either the mode number of IS6110 copies or the distribution of IS1081 types.
In this study, PCR-RFLP analysis (PRA) targeting hsp65 and rpoB gene regions was evaluated for the identification of mycobacterial species isolated from Malaysian patients. Overall, the hsp65 PRA identified 92.2 % of 90 isolates compared to 85.6 % by the rpoB PRA. With 47 rapidly growing species, the hsp65 PRA identified fewer (89.4 %) species than the rpoB PRA (95.7 %), but with 23 slow-growing species the reverse was true (91.3 % identification by the hsp65 PRA but only 52.5 % by the rpoB PRA). There were 16 isolates with discordant PRA results, which were resolved by 16S rRNA and hsp65 gene sequence analysis. The findings in this study suggest that the hsp65 PRA is more useful than the rpoB PRA for the identification of Mycobacterium species, particularly with the slow-growing members of the genus. In addition, this study reports 5 and 12 novel restriction patterns for inclusion in the hsp65 and rpoB PRA algorithms, respectively.
The worldwide-distributed aquatic fungus Articulospora tetracladia Ingold is a dominant sporulating species in streams of the Northwest Iberian Peninsula. To elucidate the genetic diversity of A. tetracladia, we analyzed isolates collected from various types of plant litter or foam in streams from North and Central Portugal and North Spain, between 2000 and 2010. Genetic diversity of these fungal populations was assessed by denaturing gradient gel electrophoresis (DGGE) fingerprints and by using ITS1-5.8S-ITS2 barcodes. Moreover, ITS1-5.8S-ITS2 barcodes of A. tetracladia reported in other parts of the world (Central Europe, United Kingdom, Canada, Japan and Malaysia) were retrieved from the National Center for Biotechnology (NCBI) and the National Institute of Technology and Evaluation Biological Resource Center (NBRC) to probe into genetic diversity of A. tetracladia. PCR-DGGE of ITS2 region of 50 Iberian fungal isolates distinguished eight operational taxonomic units (OTUs), which were similar to those obtained from neighboring trees based on ITS2 gene sequences. On the other hand, ITS1-5.8S-ITS2 barcodes of 68 fungal isolates yielded nine OTUs, but five fungal isolates were not assigned to any of these OTUs. Molecular diversity was highest for OTU-8, which included only European isolates. Two haplotypes were observed within OTU-8 and OTU-9, while only one haplotype was found within each of the remaining OTUs. Malaysia did not share haplotypes with other countries. Overall results indicate that, apart from the Malaysian genotypes, A. tetracladia genotypes were geographically widespread irrespective of sampling time, sites or substrates. Furthermore, PCR-DGGE appeared to be a rapid tool for assessing intraspecific diversity of aquatic hyphomycetes.
We barcoded 25 in vitro isolates (representing 92 samples) of Giardia duodenalis from humans and other animals, which have been assembled by the Upcroft team at the Queensland Institute of Medical Research over a period of almost three decades. We used mutation scanning-coupled sequencing of loci in the triosephosphate isomerase, glutamate dehydrogenase and β-giardin genes, combined with phylogenetic analysis, to genetically characterise them. Specifically, the isolates (n514) of G. duodenalis from humans from Australia (AD113; BRIS/83/HEPU/106; BRIS/87/HEPU/713; BRIS/89/HEPU/1003; BRIS/92/HEPU/1541; BRIS/92/HEPU/1590; BRIS/92/HEPU/2443; BRIS/93/HEPU/1706), Malaysia (KL/92/IMR/1106) and Afghanistan (WB), a cat from Australia (BAC2), a sheep from Canada (OAS1) and a sulphur-crested cockatoo from Australia (BRIS/95/HEPU/2041) represented assemblage A (sub-assemblage AI-1, AI-2 or AII-2); isolates (n510) from humans from Australia (BRIS/91/HEPU/1279; BRIS/92/HEPU/2342; BRIS/92/HEPU/2348; BRIS/93/HEPU/1638; BRIS/93/HEPU/1653; BRIS/93/HEPU/1705; BRIS/93/HEPU/1718; BRIS/93/HEPU/1727), Papua New Guinea (BRIS/92/HEPU/1487) and Canada (H7) represented assemblage B (sub-assemblage BIV) and an isolate from cattle from Australia (BRIS/92/HEPU/1709) had a match to assemblage E. Isolate BRIS/90/HEPU/1229 from a human from Australia was shown to represent a mixed population of assemblages A and B. These barcoded isolates (including stocks and derived lines) now allow direct comparisons of experimental data among laboratories and represent a massive resource for transcriptomic, proteomic, metabolic and functional genomic studies using advanced molecular technologies.
17 Y-STRs (DYS456, DYS389I, DYS390, DYS389II, DYS458, DYS19, DYS385a/b, DYS393, DYS391, DYS439, DYS635 or Y-GATA C4, DYS392, Y-GATA H4, DYS437, DYS438 and DYS448) have been analyzed in 320 male individuals from Sarawak, an eastern state of Malaysia on the Borneo island using the AmpFlSTR Y-filer (Applied Biosystems, Foster City, CA). These individuals were from three indigenous ethnic groups in Sarawak comprising of 103 Ibans, 113 Bidayuhs and 104 Melanaus. The observed 17-loci haplotypes and the individual allele frequencies for each locus were estimated, whilst the locus diversity, haplotype diversity and discrimination capacity were calculated in the three groups. Analysis of molecular variance (AMOVA) indicated that 87.6% of the haplotypic variation was found within population and 12.4% between populations (fixation index F(ST)=0.124, p=0.000). This study has revealed that the indigenous populations in Sarawak are distinctly different to each other, and to the three major ethnic groups in Malaysia (Malays, Chinese and Indians), with the Melanaus having a strikingly high degree of shared haplotypes within. There are rare unusual variants and microvariants that were not present in Malaysian Malay, Chinese or Indian groups. In addition, occurrences of DYS385 duplications which were only noticeably present in Chinese group previously was also observed in the Iban group whilst null alleles were detected at several Y-loci (namely DYS19, DYS392, DYS389II and DYS448) in the Iban and Melanau groups.
Twenty-one Vibrio parahaemolyticus isolates representing 21 samples of coastal seawater from three beaches in peninsular Malaysia were found to be sensitive to streptomycin, norfloxacin and chloramphenicol. Resistance was observed to penicillin (100%), ampicillin (95.2%), carbenicilin (95.2%), erythromycin (95.2%), bacitracin (71.4%), cephalothin (28.6%), moxalactam (28.6%), kanamycin (19.1%), tetracycline (14.3%), nalidixic acid (9.5%) and gentamicin (9.5%). Plasmids of 2.6 to 35.8 mDa were detected among plasmid-containing isolates. All isolates carried the Vp-toxR gene specific to V. parahaemolyticus and were negative for the tdh gene, but only one isolate was positive for the trh gene. DNA fingerprinting of the isolates using ERIC-PCR and PFGE showed that the isolates belong to two major clonal groups, with several isolates from different locations in the same group, indicating the presence of similar strains in the different locations.
Seven isolates of Burkholderia pseudomallei from cases of melioidosis in human (2 isolates) and animal (2 isolates), cat (one isolate) and from soil samples (2 isolates) were examined for in vitro sensitivity to 14 antimicrobial agents and for presence of plasmid DNA. Randomly amplified polymorphic DNA (RAPD) analysis was used to type the isolates, using two arbitrary primers. All isolates were sensitive to chloramphenicol, kanamycin, carbenicillin, rifampicin, enrofloxacin, tetracycline and sulfamethoxazole-trimethoprim. No plasmid was detected in all the isolates tested. RADP fingerprinting demonstrated genomic relationship between isolates, which provides an effective method to study the epidemiology of the isolates examined.
Staphylococcus aureus is a persistent human pathogen responsible for a variety of infections ranging from soft-tissue infections to bacteremia. The objective of this study was to determine genetic relatedness between methicillin-resistant S. aureus (MRSA) and methicillin-susceptible S. aureus (MSSA) strains. We isolated 35 MRSA and 21 MSSA strains from sporadic cases at the main tertiary hospital in Terengganu, Malaysia, screening them for the presence of virulence genes. Their genetic relatedness was determined by accessory gene regulator (agr) types, PCR-restriction fragment length polymorphism (RFLP) of the coa gene, pulsed-field gel electrophoresis (PFGE), S. aureus protein A (spa), and multilocus-sequence typing (MLST). We found that 57% of MRSA and 43% of MSSA strains harbored enterotoxin genes. The majority (87.5%) of the strains were agr type I. PCR-RFLP and PFGE genotyping of the coa gene revealed that MRSA strains were genetically related, whereas MSSA strains had higher heterogeneity. The combined genotype, MLST-spa type ST239-t037, was shared among MRSA and MSSA strains, indicating that MRSA strains could have evolved from MSSA strains. Two combined MLST-spa types were present in MRSA strains, whereas 7 different MLST-spa types were detected in MSSA strains, including 2 combined types (ST779-t878 and ST1179-t267) that have not been reported in Malaysia. In conclusion, enterotoxin genes were more prevalent in MRSA than in MSSA strains in the Terengganu hospital. The MSSA strains were genetically more diverse than the MRSA strains.
The isolation of therapeutic and functional protease inhibitors in vitro via combinatorial chemistry and phage display technology has been described previously. Here we report the construction of a combinatorial mouse-human chimeric antibody fragment (Fab) antibody library targeted against the protease of the tropical pathogen, Burkholderia pseudomallei. The resulting library was biopanned against the protease, and selected clones were analyzed for their ability to function as protease inhibitors. Three families of Fabs were identified by restriction fingerprinting, all of which demonstrated high specificity towards the protease of B. pseudomallei. Purified Fabs also demonstrated the capacity to inhibit B. pseudomallei protease activity in vitro, and this inhibitory property was exclusive to the pathogenic protease. Thus these recombinant antibodies are candidates for immunotherapy and tools to aid in further elucidation of the mechanism of action of the B. pseudomallei protease.
The first ever mass identification operation of skeletal remains conducted for the clandestine graves in Malaysia consisted of 165 individuals unearthed from 28 human trafficking transit camps located in Wang Kelian, along the Thai-Malaysia border. A DVI response was triggered in which expert teams comprising of pathologists, anthropologists, odontologists, radiologists and DNA experts were gathered at the identified operation centre. The Department of Forensic Medicine, Hospital Sultanah Bahiyah, Alor Star, Kedah, located approximately 75km away from Wang Kelian, was temporarily converted into a victim identification centre (VIC) as it is the nearest available forensic facility to the mass grave site. The mortuary operation was conducted over a period of 3 months from June to September 2015, and was divided into two phases; phase 1 involving the postmortem examination of the remains of 116 suspected individuals and for phase 2 the remains of 49 suspected individuals. The fact that the graves were of unknown individuals afforded the mass identification operation a sufficient duration of 2 weeks as preparatory phase enabling procedurals and daily victim identification workflow to be established, and the setting up of a temporary body storage for the designated mortuary. The temporary body storage has proven to be a significant factor in enabling the successful conclusion of the VIC operation to the final phase of temporary controlled burials. Recognition from two international observers, Mr. Andréas Patiño Umaña, from the International Committee of Red Cross (ICRC) and Prof. Noel Woodford from Victoria Institute of Forensic Medicine (VIFM) had proven the mortuary operation was in compliance to the international quality and standards. The overall victim identification and mortuary operation identified a number of significant challenges, in particular the management of commingled human remains as well as the compilation of postmortem data in the absence of antemortem data for future reconciliation. CF Index in this DVI operation is 9%, indicating the primary identifications in this operation were effective. Limitations and further improvements of the mass identification operation will be discussed. This paper details the planning, preparations and management of the mass identification operation on the exhumed human remains which was also a forensic humanitarian service for the dead with dignity and respect.
The use of single nucleotide polymorphisms (SNPs) in forensic genetics has been limited to challenged samples with low template and/or degraded DNA. The recent introduction of massively parallel sequencing (MPS) technologies has expanded the potential applications of these markers and increased the discrimination power of well-established loci by considering variation in the flanking regions of target loci. The ForenSeq Signature Preparation Kit contains 165 SNP amplicons for ancestry- (aiSNPs), identity- (iiSNPs), and phenotype-inference (piSNPs). In this study, 714 individuals from four major populations (African American, AFA; East Asian, ASN; US Caucasian, CAU; and Southwest US Hispanic, HIS) previously reported by Churchill et al. [Forensic Sci Int Genet. 30 (2017) 81-92; DOI: https://doi.org/10.1016/j.fsigen.2017.06.004] were assessed using STRait Razor v2s to determine the level of diversity in the flanking regions of these amplicons. The results show that nearly 70% of loci showed some level of flanking region variation with 22 iiSNPs and 8 aiSNPs categorized as microhaplotypes in this study. The heterozygosities of these microhaplotypes approached, and in one instance surpassed, those of some core STR loci. Also, the impact of the flanking region on other forensic parameters (e.g., power of exclusion and power of discrimination) was examined. Sixteen of the 94 iiSNPs had an effective allele number greater than 2.00 across the four populations. To assess what effect the flanking region information had on the ancestry inference, genotype probabilities and likelihood ratios were determined. Additionally, concordance with the ForenSeq UAS and Nextera Rapid Capture was evaluated, and patterns of heterozygote imbalance were identified. Pairwise comparison of the iiSNP diplotypes determined the probability of detecting a mixture (i.e., observing ≥ 3 haplotypes) using these loci alone was 0.9952. The improvement in random match probabilities for the full regions over the target iiSNPs was found to be significant. When combining the iiSNPs with the autosomal STRs, the combined match probabilities ranged from 6.40 × 10-73 (ASN) to 1.02 × 10-79 (AFA).
Of 97 strains of Vibrio cholerae isolated from various seafoods in Malaysia in 1998 and 1999, 20 strains carried the ctx gene and produced cholera toxin. Fourteen, one, and five of these toxigenic strains belonged to the O139, O1 Ogawa, and rough serotypes, respectively. The rough strains had the rfb gene of the O1 serotype. The toxigenic strains varied in their biochemical characteristics, the amount of cholera toxin produced, their antibiograms, and the presence or absence of the pTLC plasmid sequence. DNA fingerprinting analysis by arbitrarily primed PCR, ribotyping, and a pulsed-field gel electrophoresis method classified the toxigenic strains into 3, 7, and 10 types, respectively. The relatedness of these toxigenic strains to clinical strains isolated in other countries and from international travelers was examined by using a dendrogram constructed from the pulsed-field gel electrophoresis profiles. The results of the examination of the antibiogram and the possession of the toxin-linked cryptic plasmid were consistent with the dendrogram-based relatedness: the O139 strains isolated from Malaysian seafoods could be separated into two groups that appear to have been introduced from the Bengal area independently. The rough strains of Malaysian seafood origin formed one group and belonged to a cluster unique to the Thailand-Malaysia-Laos region, and this group may have persisted in this area for a long period. The single O1 Ogawa strain detected in Malaysian seafood appears to have an origin and route of introduction different from those of the O139 and the rough strains.
Burkholderia pseudomallei, the causative agent of melioidosis, is intrinsically resistant to many antibiotics. Ceftazidime (CAZ), the synthetic β-lactam, is normally used as the first-line antibiotic therapy for treatment of melioidosis. However, acquired CAZ resistance can develop in vivo during treatment with CAZ, leading to mortality if therapy is not switched to a different antibiotic(s) in a timely manner. In this study, susceptibilities of 81 B. pseudomallei isolates to nine different antimicrobial agents were determined using the disk diffusion method, broth microdilution test and Etest. Highest percentage of susceptibility was demonstrated to CAZ, amoxicillin/clavulanic acid, meropenem, imipenem, and trimethoprim/sulfamethoxazole. Although these drugs demonstrated the highest percentage of susceptibility in B. pseudomallei, the overall results underline the importance of the emergence of resistance in this organism. PCR results showed that, of the 81 B. pseudomallei, six multidrug resistant (MDR) isolates carried bpeB, amrB, and BPSS1119 and penA genes. Genotyping of the isolates using random amplified polymorphic DNA analysis showed six different PCR fingerprinting patterns generated from the six MDR isolates clusters (A) and eight PCR fingerprinting patterns generated for the remaining 75 non-MDR isolates clusters (B).
Anisakis spp. (Nematoda: Anisakidae) parasitize a wide range of marine animals, mammals serving as the definitive host and different fish species as intermediate or paratenic hosts. In this study, 18 fish species were investigated for Anisakis infection. Katsuwonus pelamis, Euthynnus affinis, Caranx sp., and Auxis thazard were infected with high prevalence of Anisakis type I, while Cephalopholis cyanostigma and Rastrelliger kanagurta revealed low prevalence. The mean intensity of Anisakis larvae in K. pelamis and A. thazard was 49.7 and 5.6, respectively. A total of 73 Anisakis type I larvae collected from K. pelamis and A. thazard were all identified as Anisakis typica by PCR-RFLP analysis. Five specimens of Anisakis from K. pelamis and 15 specimens from A. thazard were sequenced using ITS1-5.8S-ITS2 region and 6 specimens from A. thazard and 4 specimens from K. pelamis were sequenced in mtDNA cox2 region. Alignments of the samples in the ITS region showed 2 patterns of nucleotides. The first pattern (genotype) of Anisakis from A. thazard had 100% similarity with adult A. typica from dolphins from USA, whereas the second genotype from A. thazard and K. pelamis had 4 base pairs different in ITS1 region with adult A. typica from USA. In the mtDNA cox2 regions, Anisakis type I specimens from A. thazard and K. pelamis showed similarity range from 94% to 99% with A. typica AB517571/DQ116427. The difference of 4 bp nucleotides in ITS1 regions and divergence into 2 subgroups in mtDNA cox2 indicating the existence of A. typica sibling species in the Makassar Strait.