In a previous study of Southeast Asian genetic variation, we characterized mitochondrial DNAs (mtDNAs) from six populations through high-resolution restriction fragment length polymorphism (RFLP) analysis. Our analysis revealed that these Southeast Asian populations were genetically similar to each other, suggesting they had a common origin. However, other patterns of population associations also emerged. Haplotypes from a major founding haplogroup in Papua New Guinea were present in Malaysia; the Vietnamese and Malaysian aborigines (Orang Asli) had high frequencies of haplogroup F, which was also seen in most other Southeast Asian populations; and haplogroup B, defined by the Region V 9-base-pair deletion, was present throughout the region. In addition, the Malaysian and Sabah (Borneo) aborigine populations exhibited a number of unique mtDNA clusters that were not observed in other populations. Unfortunately, it has been difficult to compare these patterns of genetic diversity with those shown in subsequent studies of mtDNA variation in Southeast Asian populations because the latter have typically sequenced the first hypervariable segment (HVS-I) of the control region (CR) sequencing rather than used RFLP haplotyping to characterize the mtDNAs present in them. For this reason, we sequenced the HVS-I of Southeast Asian mtDNAs that had previously been subjected to RFLP analysis, and compared the resulting data with published information from other Southeast Asian and Oceanic groups. Our findings reveal broad patterns of mtDNA haplogroup distribution in Southeast Asia that may reflect different population expansion events in this region over the past 50,000-5,000 years.
Flea-borne pathogens were screened from 100 individual cat fleas using a PCR approach, of which 38 % were infected with at least one bacterium. Overall, 28 % of the flea samples were positive for Bartonella as inferred from ITS DNA region. Of these, 25 % (7/28) were identified as Bartonella clarridgeiae, 42.9 % (12/28) as Bartonella henselae consisted of two different strains, and 32.1 % (9/28) as Bartonella koehlerae, which was detected for the first time in Malaysia. Sequencing of gltA amplicons detected Rickettsia DNA in 14 % of cat flea samples, all of them identified as Rickettsia asembonensis (100 %). None of the flea samples were positive for Mycoplasma DNA in 16S rRNA gene detection. Four fleas were co-infected with Bartonella and Rickettsia DNAs. Statistical analyses reveal no significant association between bacterial infection and mtDNA diversity of the cat flea. Nevertheless, in all types of pathogen infections, infected populations demonstrated lower nucleotide and haplotype diversities compared to uninfected populations. Moreover, lower haplotype numbers were observed in infected populations.
Spirometra larvae are etiological agents of human sparganosis. However, the systematics of spirometrid cestodes has long been controversial. In order to determine the current knowledge on the evolution and genetic structure of Spirometra, an exhaustive population diversity analysis of spirometrid cestodes using the mitochondrial gene: cytochrome c oxidase subunit 1 (cox1) was performed. All publicly available cox1 sequences available in the GenBank and 127 new sequencing genes from China were used as the dataset. The haplotype identify, network, genetic differentiation and phylogenetic analysis were conducted successively. A total of 488 sequences from 20 host species, representing four spirometrid tapeworms (S. decipiens, S. ranarum, S. erinaceieuropaei and Sparganum proliferum) and several unclassified American and African isolates from 113 geographical locations in 17 countries, identified 45 haplotypes. The genetic analysis revealed that there are four clades of spirometrid cestodes: Clade 1 (Brazil + USA) and Clade 2 (Argentina + Venezuela) included isolates from America, Clade 3 contained African isolates and one Korean sample, and the remainders from Asia and Australia belonged to Clade 4; unclassified Spirometra from America and Africa should be considered the separate species within the genus; and the taxonomy of two Korea isolates (S. erinaceieuropaei KJ599680 and S. decipiens KJ599679) was still ambiguous and needs to be further identified. In addition, the demographical analyses supported population expansion for the total spirometrid population. In summary, four lineages were found in the spirometrid tapeworm, and further investigation with deeper sampling is needed to elucidate the population structure.
The estuarine firefly, Pteroptyx tener, aggregates in the thousands in mangrove trees lining tidal rivers in Southeast Asia where they engage one another in a nocturnal, pre-mating ritual of synchronised courtship flashes. Unfortunately, populations of the species by virtue of being restricted to isolated estuarine rivers systems in the region, are at risk of genetic isolation. Because of this concern we undertook the task of sequencing and characterising the mitochondrial DNA genome of P. tener, as the first step towards helping us to characterise and better understand their genetic diversity. We sequenced and assembled the mitochondrial DNA genome of P. tener from two male and female specimens from the district of Kuala Selangor in Peninsular Malaysia and announce the molecules in this publication. We also reconstructed the phylogenetic trees of all available lampyrids mitogenomes and suggest the need to re-examine our current understanding of their classification which have largely been based on morphological data and the cox1 gene. Separately, our analysis of codon usage patterns among lampyrid mitogenomes showed that the codon usage in a majority of the protein-coding genes were non-neutral. Codon usage patterns between mitogenome sequences of P. tener were, however, largely neutral. Our findings demonstrate the usefulness of mitochondrial genes/mitogenomes for analysing both inter- and intra- specific variation in the Lampyridae to aid in species discovery in this highly variable genus; and elucidate the phylogenetic relationships of Pteroptyx spp. from the region.
Bactrocera dorsalis s.s. is a pestiferous tephritid fruit fly distributed from Pakistan to the Pacific, with the Thai/Malay peninsula its southern limit. Sister pest taxa, B. papayae and B. philippinensis, occur in the southeast Asian archipelago and the Philippines, respectively. The relationship among these species is unclear due to their high molecular and morphological similarity. This study analysed population structure of these three species within a southeast Asian biogeographical context to assess potential dispersal patterns and the validity of their current taxonomic status.
Amphiboloidea is a small but widespread group of snails found exclusively, and often abundantly, in mudflat and associated salt marsh or mangrove habitat. This study uses molecular data from three loci (COI, 16S and 28S) to infer phylogenetic relationships in Amphiboloidea and examine its position in Euthyneura. All but two of the named extant species of Amphiboloidea and additional undescribed taxa from across Southeast Asia and the Arabian Gulf were sampled. In contrast to the current morphology-based classification dividing Amphiboloidea into three families, analysis of molecular data supports revision of the classification to comprise two families. Maningrididae is a monotypic family basal to Amphibolidae, which is revised to comprise three subfamilies: Amphibolinae, Phallomedusinae and Salinatorinae. Sequence divergence between Asian populations of Naranjia is relatively large and possibly indicative of species complexes divergent across the Strait of Malacca. Salinatorrosacea and Salinator burmana do not cluster with other Salinator species, and require generic reassignment. In addition, sequences were obtained from an undescribed species of Lactiforis from the Malay Peninsula. Reconstruction of ancestral distributions indicates a plesiomorphic distribution and centre of origin in Australasia, with two genera subsequently diversifying throughout Asia. Increasing the sampling density of amphiboloid taxa in a phylogenetic analysis of Euthyneura did not resolve the identity of the sister taxon to Amphibolidae, but confirmed its inclusion in Pulmonata/Panpulmonata.
Two subspecies of cynomolgus macaques (Macaca fascicularis) are alleged to co-exist in the Philippines, M. f. philippensis in the north and M. f. fascicularis in the south. However, genetic differences between the cynomolgus macaques in the two regions have never been studied to document the propriety of their subspecies status. We genotyped samples of cynomolgus macaques from Batangas in southwestern Luzon and Zamboanga in southwestern Mindanao for 15 short tandem repeat (STR) loci and sequenced an 835 bp fragment of the mtDNA of these animals. The STR genotypes were compared with those of cynomolgus macaques from southern Sumatra, Singapore, Mauritius and Cambodia, and the mtDNA sequences of both Philippine populations were compared with those of cynomolgus macaques from southern Sumatra, Indonesia and Sarawak, Malaysia. We conducted STRUCTURE and PCA analyses based on the STRs and constructed a median joining network based on the mtDNA sequences. The Philippine population from Batangas exhibited much less genetic diversity and greater genetic divergence from all other populations, including the Philippine population from Zamboanga. Sequences from both Batangas and Zamboanga were most closely related to two different mtDNA haplotypes from Sarawak from which they are apparently derived. Those from Zamboanga were more recently derived than those from Batangas, consistent with their later arrival in the Philippines. However, clustering analyses do not support a sufficient genetic distinction of cynomolgus macaques from Batangas from other regional populations assigned to subspecies M. f. fascicularis to warrant the subspecies distinction M. f. philippensis.
The complete mitochondrial genome of the commercially important snout otter clam Lutraria rhynchaena was obtained from low-coverage shotgun sequencing data on the MiSeq platform. The L. rhynchaena mitogenome has 16,927 base pairs (69% A + T content) and made up of 12 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs, and a 953 bp non-coding AT-rich region. This is the first mitogenome to be sequenced from the genus Lutraria, and the seventh to be reported for the family Mactridae.
To elucidate genetic divergence and evolutionary relationship in Fejervarya cancrivora from Indonesia and other Asian countries, allozyme and molecular analyses were carried out using 131 frogs collected from 24 populations in Indonesia, Thailand, Bangladesh, Malaysia, and the Philippines. In the allozymic survey, seventeen enzymatic loci were examined for 92 frogs from eight representative localities. The results showed that F. cancrivora is subdivided into two main groups, the mangrove type and the large- plus Pelabuhan ratu types. The average Nel's genetic distance between the two groups was 0.535. Molecular phylogenetic trees based on nucleotide sequences of the 16S rRNA and Cyt b genes and constructed with the ML, MP, NJ, and BI methods also showed that the individuals of F. cancrivora analyzed comprised two clades, the mangrove type and the large plus Pelabuhan ratu / Sulawesi types, the latter further split into two subclades, the large type and the Pelabuhan ratu / Sulawesi type. The geographical distribution of individuals of the three F. cancrivora types was examined. Ten Individuals from Bangladesh, Thailand, and the Philippines represented the mangrove type; 34 Individuals from Malaysia and Indonesia represented the large type; and 11 individuals from Indonesia represented the Pelabuhan ratu / Sulawesi type. Average sequence divergences among the three types were 5.78-10.22% for the 16S and 12.88-16.38% for Cyt b. Our results suggest that each of the three types can be regarded as a distinct species.
This study examined 396 base pairs of the mitochondrial cytochrome b gene from 110 individuals belonging to the genus Hampala, a group of freshwater cyprinids that inhabit Southeast Asia. The samples were taken from various locations throughout Sarawak, Sabah, and peninsular Malaysia. The nucleotide sequences were subjected to phylogenetic analyses by using the neighbor-joining, maximum parsimony, and maximum likelihood methods. All three methods revealed the reciprocally monophyletic relationship of Hampala macrolepidota to the other Hampala forms, thus strongly supporting its status as a distinct species. Phylogenetic analysis also discovered the existence of two H. bimaculata lineages endemic to Borneo: (1) a newly identified species from the southern and central part of Sarawak assigned as H. bimaculata Type A and (2) the previously described H. bimaculata from northern Sarawak and the west coast of Sabah assigned as H. bimaculata Type B. However, the status of H. sabana and an intermediate form were not elucidated. The results suggest that the intermediate form from the Tawau population is actually a subpopulation of H. sabana, while the highly divergent intermediate form from Kalabakan could represent a cryptic species. The sharing of H. macrolepidota haplotypes in the southern peninsular Malaysia and southern and central Sarawak samples (Hm1 and Hm2) reflected the recent disconnection of the two regions, during the late Pleistocene. Overall, the partial sequence of the mitochondrial cytochrome b gene was useful for resolving the phylogenetic relationships among Hampala fishes in Malaysia.
Little is known about the classification and phylogenetic relationships of the leaf monkeys (Presbytis). We analyzed mitochondrial DNA sequences of cytochrome b (Cyt b) and 12S rRNA to determine the phylogenetic relationships of the genus Presbytis. Gene fragments of 388 and 371 bp of Cyt b and 12S rRNA, respectively, were sequenced from samples of Presbytis melalophos (subspecies femoralis, siamensis, robinsoni, and chrysomelas), P. rubicunda and P. hosei. The genus Trachypithecus (Cercopithecidae) was used as an outgroup. The Cyt b NJ and MP phylogeny trees showed P. m. chrysomelas to be the most primitive, followed by P. hosei, whereas 12S rRNA tree topology only indicated that these two species have close relationships with the other members of the genus. In our analysis, chrysomelas, previously classified as a subspecies of P. melalophos, was not included in either the P. m. femoralis clade or the P. m. siamensis clade. Whether or not there should be a separation at the species level remains to be clarified. The tree topologies also showed that P. m. siamensis is paraphyletic with P. m. robinsoni, and P. m. femoralis with P. rubicunda, in two different clades. Cyt b and 12S rRNA are good gene candidates for the study of phylogenetic relationships at the species level. However, the systematic relationships of some subspecies in this genus remain unclear.
In this study, we sequenced a partial segment of the mitochondrial control region from 21 proboscis monkeys of the Klias peninsula, the last large population remaining on the west coast of Sabah, Malaysia. Our results showed that this population retains substantial genetic variation, and subpopulations from different river systems in the central and southern portions of the Klias share multiple haplotypes. We also compared our data with previously generated sequences from 2 eastern populations of proboscis monkeys in Sabah and found little evidence of regional genetic structure. Based on these results, we argue that conservation efforts should focus on restoring connectivity between central and southern Klias peninsula proboscis monkeys and discuss future analyses needed to better understand the mitochondrial structure of proboscis monkeys in Sabah.
Taxonomic relationships within the Old World fruit bat genus, Cynopterus, have been equivocal for the better part of a century. While nomenclature has been revised multiple times on the basis of phenotypic characters, evolutionary relationships among taxa representing the entire geographic range of the genus have not been determined. We used mitochondrial DNA sequence data to infer phylogenetic relationships among the three most broadly distributed members of the genus: C. brachyotis, C. horsfieldi, and C. sphinx, and to assess whether C. brachyotis represents a single widespread species, or a complex of distinct lineages. Results clearly indicate that C. brachyotis is a complex of lineages. C. sphinx and C. horsfieldi haplotypes formed monophyletic groups nested within the C. brachyotis species complex. We identified six divergent mitochondrial lineages that are currently referred to C. brachyotis. Lineages from India, Myanmar, Sulawesi, and the Philippines are geographically well-defined, while in Malaysia two lineages, designated Sunda and Forest, are broadly sympatric and may be ecologically distinct. Demographic analyses of the Sunda and Forest lineages suggest strikingly different population histories, including a recent and rapid range expansion in the Sunda lineage, possibly associated with changes in sea levels during the Pleistocene. The resolution of the taxonomic issues raised in this study awaits combined analysis of morphometric characters and molecular data. However, since both the Indian and Malaysian Forest C. brachyotis lineages are apparently ecologically restricted to increasingly fragmented forest habitat, we suggest that reevaluation of the conservation status of populations in these regions should be an immediate goal.
Christmas Island is a remote Australian territory located close to the main Indonesian island of Java. Y-chromosome and mitochondrial DNA (mtDNA) markers were used to investigate the genetic structure of the population, which comprises communities of mixed ethnic origin. Analysis of 12 Y-chromosome biallelic polymorphisms revealed a high level of gene diversity and haplotype frequencies that were consistent with source populations in southern China and Southeast Asia. mtDNA hypervariable segment I (HVS-I) sequences displayed high levels of haplotype diversity and nucleotide diversity that were comparable to various Asian populations. Genetic distances revealed extremely low mtDNA differentiation among Christmas Islanders and Asian populations. This was supported by the relatively high proportion of sequence types shared among these populations. The most common mtDNA haplogroups were M* and B, followed by D and F, which are prevalent in East/Southeast Asia. Christmas Islanders of European descent were characterized by the Eurasian haplogroup R*, and a limited degree of admixture was observed. In general, analysis of the genetic data indicated population affinities to southern Chinese (in particular from the Yunnan Province) and Southeast Asia (Thailand, Malaysia, and Cambodia), which was consistent with historical records of settlement. The combined use of these different marker systems provides a useful and appropriate model for the study of contemporary populations derived from different ethnic origins.
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.
The mosquito Aedes albopictus is indigenous to Southeast Asian and is a vector for arbovirus diseases. Studies examining the population genetics structure of A. albopictus have been conducted worldwide; however, there are no documented reports on the population genetic structure of A. albopictus in Malaysia, particularly in Penang. We examined the population genetics of A. albopictus based on a 445-base pair segment of the mitochondrial DNA cytochrome oxidase 1 gene among 77 individuals from 9 localities representing 4 regions (Seberang Perai Utara, Seberang Perai Tengah, Northeast, and Southwest) of Penang. A total of 37 haplotypes were detected, including 28 unique haplotypes. The other 9 haplotypes were shared among various populations. These shared haplotypes reflect the weak population genetic structure of A. albopictus. The phylogenetic tree showed a low bootstrap value with no genetic structure, which was supported by minimum spanning network analysis. Analysis of mismatch distribution showed poor fit of equilibrium distribution. The genetic distance showed low genetic variation, while pairwise FST values showed no significant difference between all regions in Penang except for some localities. High haplotype diversity and low nucleotide diversity was observed for cytochrome oxidase 1 mtDNA. We conclude that there is no population genetic structure of A. albopictus mosquitoes in the Penang area.
For centuries, morphology-based fish identification has been applied without molecular evaluation. Many studies showed that specimens with a similar morphology are frequently found to be quite genetically distinct. One of the fish species that still remains taxonomically problematic is a commercial snapper species, Lutjanus johnii. Because of morphological ambiguities among local fish taxonomists in Malaysia, we examined the ability of the cytochrome oxidase I (COI) gene to genetically examine the taxonomic status of L. johnii. A 626-base pair COI region was successfully amplified and aligned with conspecific sequences that were retrieved from GenBank. The phylogenetic tree obtained showed two major clusters; the first cluster consists of L. johnii from Straits of Malacca, Thailand, Australia, and China while the second cluster comprises L. johnii from China and India. The latter group showed sequence divergence greater than 3.5%. After observing this, we suspected that there might be a cryptic species between the South China Sea and Indian Ocean. This is the first molecular report concerning the commercial species of snapper, L. johnii, in Malaysia, which had only gained provisional recognition from morphological examination.
The mitochondrial genome sequence of the ghost crab, Ocypode ceratophthalmus, is documented (GenBank accession number: LN611669) in this article. This is the first mitogenome for the family Ocypodidae and the second for the order Ocypodoidea. Ocypode ceratophthalmus has a mitogenome of 15,564 base pairs consisting of 13 protein-coding genes, two ribosomal subunit genes, 22 transfer RNAs and a non-coding AT-rich region. The base composition of the O. ceratophthalmus mitogenome is 35.78% for T, 19.36% for C, 33.73% for A and 11.13% for G, with an AT bias of 69.51% and the gene order is the typical arrangement for brachyuran crabs.
The mitochondrial DNA (mtDNA) cytochrome oxidase I (COI) gene has been universally and successfully utilized as a barcoding gene, mainly because it can be amplified easily, applied across a wide range of taxa, and results can be obtained cheaply and quickly. However, in rare cases, the gene can fail to distinguish between species, particularly when exposed to highly sensitive methods of data analysis, such as the Bayesian method, or when taxa have undergone introgressive hybridization, over-splitting, or incomplete lineage sorting. Such cases require the use of alternative markers, and nuclear DNA markers are commonly used. In this study, a dendrogram produced by Bayesian analysis of an mtDNA COI dataset was compared with that of a nuclear DNA ATPS-α dataset, in order to evaluate the efficiency of COI in barcoding Malaysian nerites (Neritidae). In the COI dendrogram, most of the species were in individual clusters, except for two species: Nerita chamaeleon and N. histrio. These two species were placed in the same subcluster, whereas in the ATPS-α dendrogram they were in their own subclusters. Analysis of the ATPS-α gene also placed the two genera of nerites (Nerita and Neritina) in separate clusters, whereas COI gene analysis placed both genera in the same cluster. Therefore, in the case of the Neritidae, the ATPS-α gene is a better barcoding gene than the COI gene.
Phylogeographic patterns and population structure of the pelagic Indian mackerel, Rastrelliger kanagurta were examined in 23 populations collected from the Indonesian-Malaysian Archipelago (IMA) and the West Indian Ocean (WIO). Despite the vast expanse of the IMA and neighbouring seas, no evidence for geographical structure was evident. An indication that R. kanagurta populations across this region are essentially panmictic. This study also revealed that historical isolation was insufficient for R. kanagurta to attain migration drift equilibrium. Two distinct subpopulations were detected between the WIO and the IMA (and adjacent populations); interpopulation genetic variation was high. A plausible explanation for the genetic differentiation observed between the IMA and WIO regions suggest historical isolation as a result of fluctuations in sea levels during the late Pleistocene. This occurrence resulted in the evolution of a phylogeographic break for this species to the north of the Andaman Sea.