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.
We classified diversity in eight new complete mitochondrial genome sequences and 41 partial sequences from living Aboriginal Australians into five haplogroups. Haplogroup AuB belongs to global lineage M, and AuA, AuC, AuD, and AuE to N. Within N, we recognize subdivisions, assigning AuA to haplogroup S, AuD to haplogroup O, AuC to P4, and AuE to P8. On available evidence, (S)AuA and (M)AuB are widespread in Australia. (P4)AuC is found in the Riverine region of western New South Wales, and was identified by others in northern Australia. (O)AuD and (P8)AuE were clearly identified only from central Australia. Our eight Australian full mt genome sequences, combined with 20 others (Ingman and Gyllensten 2003 Genome Res. 13:1600-1606) and compared with full mt genome sequences from regions to the north that include Papua New Guinea, Malaya, and Andaman and Nicobar Islands, show that ancestral connections between regions are deep and limited to clustering at the level of the N and M macrohaplogroups. The Australian-specific distribution of the five haplogroups identified indicates genetic isolation over a long period. Ancestral connections within Australia are deeper than those reflected by known linguistic or culturally based affinities. Applying a coalescence analysis to a gene tree for the coding regions of the eight genomic sequences, we made estimates of time depth that support a continuity of presence for the descendants of a founding population already established by 40,000 years ago.
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.
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.
Mitochondrial DNA variation was surveyed in nine populations of the pigtail macaque (Macaca nemestrina), covering all three recognized subspecies in Southeast Asia. To do this, a 2,300 base pair fragment spanning the mitochondrial NAD 3 and NAD 4 genes and flanking tRNA subunits leucine and glycine was targeted for amplification and digested with a battery of 16 restriction endonucleases. Out of a total of 107 individuals, 32 unique haplotypes could be distinguished. Parsimony and neighbor-joining analyses grouped the haplotypes into five strongly supported assemblages representing China/Thailand, Malaysia, Sumatra, Borneo, and Siberut. These results indicate that the mainland and island mtDNA haplotypes are strictly and uniquely limited to the geographic ranges of the recognized morphological subspecies. Cladistic and neighbor-joining analyses indicate that inferred phylogenies of mtDNA haplotypes are congruent with subspecies designations. Furthermore, in support of morphological studies, results indicate that the Mentawai macaque is most likely not a distinct species but a subspecies of M. nemestrina.
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.
This study was performed to establish the genetic variability of Aedes albopictus within Subang Jaya, Selangor, Malaysia, by using the nicotinamide adenine dinucleotide dehydrogenase 5 subunit (ND5) mitochondrial DNA (mtDNA) marker. A total of 90 samples were collected from 9 localities within an area of the Subang Jaya Municipality. Genetic variability was determined through the amplification and sequencing of a fragment of the ND5 gene. Eight distinct mtDNA haplotypes were identified. The evolutionary relationship of the local haplotypes alongside 28 reference strains was used to construct a phylogram, the analysis of which revealed low genetic differentiation in terms of both nucleotide and haplotype diversity. Bayesian method was used to infer the phylogenetic tree, revealing a unique relationship between local isolates. The study corroborates the reliability of ND5 to identify distinct lineages for polymorphism-based studies and supplements the existing body of knowledge regarding its genetic diversity. This in turn could potentially aid existing vector control strategies to help mitigate the risk and spread of the dengue virus.
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.
The Blueline Rasbora (Rasbora sarawakensis) is a small ray-finned fish categorized under the genus Rasbora in the Cyprinidae family. In this study, the complete mitogenome sequence of R. sarawakensis was sequenced using four primers targeting overlapping regions. The mitogenome is 16,709 bp in size, accommodating 22 transfer RNA genes, 13 protein-coding genes, two ribosomal RNA genes and a putative control region. Identical gene organisation was detected between this species and other genus counterparts. The heavy strand houses 28 genes while the light strand stores the other nine genes. Most protein-coding genes employ ATG as start codon, excluding COI gene, which utilizes GTG instead. The central conserved sequence blocks (CSB-F, CSB-E and CSB-D), variable sequence blocks (CSB-3, CSB-2 and CSB-1) as well as the terminal associated sequence (TAS) are conserved in the control region. The maximum likelihood phylogenetic tree revealed the divergence of R. sarawakensis from the basal region of the Rasbora clade, where its evolutionary relationships with R. maculatus and R. pauciperforata are poorly resolved as indicated by the low bootstrap values. This work acts as steppingstone towards further molecular evolution and population genetics studies of Rasbora genus in future.
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.
Mitochondria have multiple functions, including synthesis of adenine triphosphate, production of reactive oxygen species, calcium signalling, thermogenesis and apoptosis. Mitochondria have a significant contribution in regulating the various physiological aspects of reproductive function, from spermatogenesis up to fertilisation. Mitochondrial functionality and intact mitochondrial membrane potential are a pre-requisite for sperm motility, hyperactivation, capacitation, acrosin activity, acrosome reaction and DNA integrity. Optimal mitochondrial activity is therefore crucial for human sperm function and semen quality. However, the precise role of mitochondria in spermatozoa remains to be fully explored. Defects in sperm mitochondrial function severely impair the maintenance of energy production required for sperm motility and may be an underlying cause of asthenozoospermia. Sperm mtDNA is susceptible to oxidative damage and mutations that could compromise sperm function leading to infertility. Males with abnormal semen parameters have increased mtDNA copy number and reduced mtDNA integrity. This review discusses the role of mitochondria in sperm function, along with the causes and impact of its dysfunction on male fertility. Greater understanding of sperm mitochondrial function and its correlation with sperm quality could provide further insights into their contribution in the assessment of the infertile male.
Most new cryptic species are described using conventional tree- and distance-based species delimitation methods (SDMs), which rely on phylogenetic arrangements and measures of genetic divergence. However, although numerous factors such as population structure and gene flow are known to confound phylogenetic inference and species delimitation, the influence of these processes is not frequently evaluated. Using large numbers of exons, introns, and ultraconserved elements obtained using the FrogCap sequence-capture protocol, we compared conventional SDMs with more robust genomic analyses that assess population structure and gene flow to characterize species boundaries in a Southeast Asian frog complex (Pulchrana picturata). Our results showed that gene flow and introgression can produce phylogenetic patterns and levels of divergence that resemble distinct species (up to 10% divergence in mitochondrial DNA). Hybrid populations were inferred as independent (singleton) clades that were highly divergent from adjacent populations (7%-10%) and unusually similar (<3%) to allopatric populations. Such anomalous patterns are not uncommon in Southeast Asian amphibians, which brings into question whether the high levels of cryptic diversity observed in other amphibian groups reflect distinct cryptic species-or, instead, highly admixed and structured metapopulation lineages. Our results also provide an alternative explanation to the conundrum of divergent (sometimes nonsister) sympatric lineages-a pattern that has been celebrated as indicative of true cryptic speciation. Based on these findings, we recommend that species delimitation of continuously distributed "cryptic" groups should not rely solely on conventional SDMs, but should necessarily examine population structure and gene flow to avoid taxonomic inflation.
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.
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.
Environmental DNA detection has emerged as a powerful tool to monitor aquatic species without the need for capture or visual identification and is particularly useful for rare or elusive species. Our objective was to develop an eDNA approach for detecting the southern river terrapin (Batagur affinis) in Malaysia. We designed species-specific primers for a fragment of B. affinis mtDNA and evaluated their effectiveness in silico, in vitro and in situ. The primers amplified 110 bp of the cytochrome b mtDNA sequence of B. affinis from aquarium water samples housing nine juvenile B. affinis. We also successfully detected B. affinis eDNA from river samples taken from a site where turtles were known to be in the vicinity. Prospects and challenges of using an eDNA approach to help determine the distribution of B. affinis, essential information for an effective conservation plan, are discussed.
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.
Dispersal of soil-dwelling organisms via the repeatedly exposed Sunda shelf through much of the Pleistocene in Southeast Asia has not been studied extensively, especially for invertebrates. Here we investigated the phylogeography of an endemic termite species, Macrotermes gilvus (Hagen), to elucidate the spatiotemporal dynamics of dispersal routes of terrestrial fauna in Pleistocene Southeast Asia. We sampled 213 termite colonies from 66 localities throughout the region. Independently inherited microsatellites and mtDNA markers were used to infer the phylogeographic framework of M. gilvus. Discrete phylogeographic analysis and molecular dating based on fossil calibration were used to infer the dynamics of M. gilvus dispersal in time and space across Southeast Asia. We found that the termite dispersal events were consistently dated within the Pleistocene time frame. The dispersal pattern was multidirectional, radiating eastwards and southwards out of Indochina, which was identified as the origin for dispersal events. We found no direct dispersal events between Sumatra and Borneo despite the presence of a terrestrial connection between them during the Pleistocene. Instead, central Java served as an important link allowing termite colonies to be established in Borneo and Sumatra. Our findings support the hypothesis of a north-south dispersal corridor in Southeast Asia and suggest the presence of alternative dispersal routes across Sundaland during the Pleistocene. For the first time, we also propose that a west-east dispersal through over-water rafting likely occurred across the Pleistocene South China Sea. We found at least two independent entry routes for terrestrial species to infiltrate Sumatra and Borneo at different times.