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 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.
Chiari malformation type 1 (CM 1) is a rare and complex neurological condition. This congenital condition affects the lower posterior fossa, where the brain connects to the spinal cord. Although the exact cause of CM 1 remains unclear, genetic predisposition plays a considerable role in structural defects of the cerebellum. Here, we report on a 15-year-old female patient with CM 1 who exhibited both nuclear and mitochondrial genetic variants, a combination that has not been previously described. We identified a silent mutation in exon 2 (c. 556 G>A, p. Lys185=) of PAX1 and a DKK1 variant in intron 3 (548-3 t>C) in the nuclear DNA. We also screened the D-loop region of mitochondrial DNA as it exhibits a higher susceptibility to mutations than other mitochondrial DNA regions. Several hotspot variants were revealed, including those in positions 303-309 and 16519 (t>C), as well as some variants that had not been documented in MITOMAP. Our findings highlight the potential role of genetic alterations in D-loop in CM 1.
The increased rate at which complete mitogenomes are being sequenced and their increasing use for phylogenetic studies have resulted in a bioinformatic bottleneck in preparing and utilising such data for phylogenetic analysis. Hence, we present MitoPhAST, an automated tool that (1) identifies annotated protein-coding gene features and generates a standardised, concatenated and partitioned amino acid alignment directly from complete/partial GenBank/EMBL-format mitogenome flat files, (2) generates a maximum likelihood phylogenetic tree using optimised protein models and (3) reports various mitochondrial genes and sequence information in a table format. To demonstrate the capacity of MitoPhAST in handling a large dataset, we used 81 publicly available decapod mitogenomes, together with eight new complete mitogenomes of Australian freshwater crayfishes, including the first for the genus Gramastacus, to undertake an updated test of the monophyly of the major groups of the order Decapoda and their phylogenetic relationships. The recovered phylogenetic trees using both Bayesian and ML methods support the results of studies using fragments of mtDNA and nuclear markers and other smaller-scale studies using whole mitogenomes. In comparison to the fragment-based phylogenies, nodal support values are generally higher despite reduced taxon sampling suggesting there is value in utilising more fully mitogenomic data. Additionally, the simple table output from MitoPhAST provides an efficient summary and statistical overview of the mitogenomes under study at the gene level, allowing the identification of missing or duplicated genes and gene rearrangements. The finding of new mtDNA gene rearrangements in several genera of Australian freshwater crayfishes indicates that this group has undergone an unusually high rate of evolutionary change for this organelle compared to other major families of decapod crustaceans. As a result, freshwater crayfishes are likely to be a useful model for studies designed to understand the evolution of mtDNA rearrangements. We anticipate that our bioinformatics pipeline will substantially help mitogenome-based studies increase the speed, accuracy and efficiency of phylogenetic studies utilising mitogenome information. MitoPhAST is available for download at https://github.com/mht85/MitoPhAST.
Edentistoma octosulcatum Tömösváry, 1882, is a rare, superficially millipede-like centipede known only from Borneo and the Philippines. It is unique within the order Scolopendromorpha for its slow gait, robust tergites, and highly modified gizzard and mandible morphology. Not much is known about the biology of the species but it has been speculated to be arboreal with a possibly vegetarian diet. Until now its phylogenetic position within the subfamily Otostigminae has been based only on morphological characters, being variably ranked as a monotypic tribe (Arrhabdotini) or classified with the Southeast Asian genus Sterropristes Attems, 1934. The first molecular data for E. octosulcatum sourced from a newly collected specimen from Sarawak were analysed with and without morphology. Parsimony analysis of 122 morphological characters together with two nuclear and two mitochondrial loci resolves Edentistoma as sister group to three Indo-Australian species of Rhysida, this clade in turn grouping with Ethmostigmus, whereas maximum likelihood and parsimony analyses of the molecular data on their own ally Edentistoma with species of Otostigmus. A position of Edentistoma within Otostigmini (rather than being its sister group as predicted by the Arrhabdotini hypothesis) is consistently retrieved under different analytical conditions, but support values within the subfamily remain low for most nodes. The species exhibits strong pushing behaviour, suggestive of burrowing habits. Evidence against a suggested vegetarian diet is provided by observation of E. octosulcatum feeding on millipedes in the genus Trachelomegalus.
The phylogenetic relationships of long-tailed macaque (Macaca fascicularis fascicularis) populations distributed in Peninsular Malaysia in relation to other regions remain unknown. The aim of this study was to reveal the phylogeography and population genetics of Peninsular Malaysia's M. f. fascicularis based on the D-loop region of mitochondrial DNA. Sixty-five haplotypes were detected in all populations, with only Vietnam and Cambodia sharing four haplotypes. The minimum-spanning network projected a distant relationship between Peninsular Malaysian and insular populations. Genetic differentiation (F(ST), Nst) results suggested that the gene flow among Peninsular Malaysian and the other populations is very low. Phylogenetic tree reconstructions indicated a monophyletic clade of Malaysia's population with continental populations (NJ = 97%, MP = 76%, and Bayesian = 1.00 posterior probabilities). The results demonstrate that Peninsular Malaysia's M. f. fascicularis belonged to Indochinese populations as opposed to the previously claimed Sundaic populations. M. f. fascicularis groups are estimated to have colonized Peninsular Malaysia ~0.47 million years ago (MYA) directly from Indochina through seaways, by means of natural sea rafting, or through terrestrial radiation during continental shelf emersion. Here, the Isthmus of Kra played a central part as biogeographical barriers that then separated it from the remaining continental populations.
Samples of 100 random healthy unrelated Iraqi male persons from the Arab ethnic group of Iraqi population were collected for mtDNA coding region sequencing by using the Sanger technique and to establish the degree of variation characteristic of a fragment. Portion of coding region encompassing positions 11,719-12,184 was amplified in accordance with the Anderson reference sequence. PCR products were purified by EZ-10 spin column then sequenced and detected by using the ABI 3130xL DNA Analyzer. This is to intend the detection of polymorphisms of mtDNA. Four new polymorphic positions 11,741, 11,756, 11,878, and 12,133 are described which may be suitable in the future to be the sources for human identification purpose in Iraq. The obtained data can be used to identify variable nucleotide positions characterized by frequent occurrence most promising for identification variants. The calculated value GD = 0.95 and RMP = 0.048 of the genetic diversity should be understood as high in the context of coding function of the analysed DNA fragment. The relatively high gene diversity and a relatively low random match probability were observed in this study.
The complete mitochondrial genome of a highland freshwater crayfish, Cherax monticola, was recovered by shotgun sequencing. The mitogenome consists of 15,917 base pairs containing 13 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs and a non-coding AT-rich region. The base composition of C. monticola is 33.46% for T, 21.48% for C, 33.71% for A and 11.35% for G, with an AT bias of 67.17%.
The commercial freshwater crayfish Cherax quadricarinatus complete mitochondrial genome was recovered from partial genome sequencing using the MiSeq Personal Sequencer. The mitogenome has 15,869 base pairs consisting of 13 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs, and a non-coding AT-rich region. The base composition of C. quadricarinatus is 32.16% for T, 23.39% for C, 33.26% for A, and 11.19% for G, with an AT bias of 65.42%.
The mitogenome of the black yabby, Geocharax gracilis, was sequenced using the MiSeq Personal Sequencer. It has 15,924 base pairs consisting of 13 protein-coding genes, 2 ribosomal subunit genes, 23 transfer RNAs, and a non-coding AT-rich region. The base composition of G. gracilis mitogenome is 32.18% for T, 22.32% for C, 34.83% for A, and 10.68% for G, with an AT bias of 67.01%. The mitogenome gene order is typical for that of parastacid crayfish with the exception of some minor rearrangements involving tRNA genes.
This study provides a first description of the phylogeographic patterns and evolutionary history of two species of the mudskipper genus Periophthalmus. These amphibious gobies are distributed throughout the whole Indo-Pacific region and Atlantic coast of Africa, in peritidal habitats of soft-bottom coastal ecosystems. Three sequence datasets of two widely distributed species, Periophthalmus argentilineatus and P. kalolo, were obtained by amplifying and sequencing two mtDNA markers (D-loop and 16S rDNA) and the nDNA rag1 region. The three datasets were then used to perform phylogeographic, demographic and population genetic analyses. Our results indicate that tectonic events and past climatic oscillations strongly contributed to shape present genetic differentiation, phylogeographic and demographic patterns. We found support for the monophyly of P. kalolo, and only shallow genetic differentiation between East-African and Indo-Malayan populations of this species. However, our collections of the morphospecies P. argentilineatus include three molecularly distinct lineages, one of them more closely related to P. kalolo. The presence of Miocenic timings for the most recent common ancestors of some of these morphologically similar clades, suggests the presence of strong stabilising selection in mudskippers' habitats. At population level, demographic analyses and palaeoecological records of mangrove ecosystems suggest that Pleistocene bottlenecks and expansion plus secondary contact events of the studied species were associated with recurrent sea transgressions during interglacials, and sea regressions or stable regimes during glacials, respectively.
A primary factor in population management and wildlife conservation is the delineation of population units derived from descriptions of population genetic structure. Yet, predicting factors that influence the patterns of gene flow in a population particularly at landscape scales remains a major challenge in evolutionary biology. Here we report a population genetic study of the mud crab Scylla olivacea examined based on a 542 bp segment of the mitochondrial DNA cytochrome c oxidase I gene among 91 individuals from six localities in the west and east coast of Peninsular Malaysia. In total 55 unique haplotypes were distinguished with 45 private haplotypes and a single common haplotype shared among all populations studied. The other ten haplotypes were shared among various populations. The sharing of this haplotype reflects the connection of the mangrove areas between east and west coast of Peninsular Malaysia. High haplotype diversity (h = 0.968 ± 0.021; mean ± SD) and low nucleotide diversity (π = 0.120 ± 0.015; mean ± SD) were displayed, which may be indicative of genetic bottleneck events. No significant phylogenetic lineages were recognized using neighbour-joining and maximum parsimony methods. Hierarchical AMOVA analysis indicated that 99.33 % of the genetic variation was contained within populations and 0.67 % occurred among populations, suggesting no geographical patterning among populations studied, supported by F st test. Mismatch distribution analysis showed that the observed distribution of the pairwise mutation differences among haplotypes was multimodal, which is not concordant with a sudden range expansion scenario. However, neutrality tests showed non-significant negative values suggesting that the populations studied may have experienced past population growth, but the expansion may have been restricted to separate local areas that resulted in the non-significant negative Fu's Fs and Tajima's D value. Overall, this present preliminary study was able to be a reference on the phylogenetic relationships and assessment of genetic structure of Scylla sp. in Malaysia.
We sequenced the complete mitochondrial genome of the Tioman Island rock gecko, Cnemaspis limi, which is known as an endemic species to Malaysia. The complete mitogenome is 16,680 bp in size, consisting of 37 genes coding for 13 proteins, 22 transfer RNAs, two ribosomal RNAs and one control region. The A + T content of the overall base composition of H-strand is 53.09% (T: 23.20%, C: 32.48%, A: 29.89% and G: 14.43%). The major non-coding region (control region) is 1254 bp in length with the A + T content of 55.09% and four replicates of a 76-bp repeat within this region.
Malaysia remains as a crossroad of different cultures and peoples, and it has long been recognized that studying its population history can provide crucial insight into the prehistory of Southeast Asia as a whole. The earliest inhabitants were the Orang Asli in Peninsular Malaysia and the indigenous groups in Sabah and Sarawak. Although they were the earliest migrants in this region, these tribes are divided geographically by the South China Sea. We analyzed DNA sequences of 18 Orang Asli using mitochondrial DNA extracted from blood samples, each representing one sub-tribe, and from five Sarawakian Iban. Mitochondrial DNA was extracted from hair samples in order to examine relationships with the main ethnic groups in Malaysia. The D-loop region and cytochrome b genes were used as the candidate loci. Phylogenetic relationships were investigated using maximum parsimony and neighbor joining algorithms, and each tree was subjected to bootstrap analysis with 1000 replicates. Analyses of the HVS I region showed that the Iban are not a distinct group from the Orang Asli; they form a sub-clade within the Orang Asli. Based on the cytochrome b gene, the Iban clustered with the Orang Asli in the same clade. We found evidence for considerable gene flow between Orang Asli and Iban. We concluded that the Orang Asli, Iban and the main ethnic groups of Malaysia are probably derived from a common ancestor. This is in agreement with a single-route migration theory, but it does not dismiss a two-route migration theory.
Termites of the genus Odontotermes are important decomposers in the Old World tropics and are sometimes important pests of crops, timber and trees. The species within the genus often have overlapping size ranges and are difficult to differentiate based on morphology. As a result, the taxonomy of Odontotermes in Peninsular Malaysia has not been adequately worked out. In this study, we examined the phylogeny of 40 samples of Odontotermes from Peninsular Malaysia using two mitochondrial DNA regions, that is, the 16S ribosomal RNA and cytochrome oxidase subunit I genes, to aid in elucidating the number of species in the peninsula. Phylogenies were reconstructed from the individual gene and combined gene data sets using parsimony and likelihood criteria. The phylogenies supported the presence of up to eleven species in Peninsular Malaysia, which were identified as O. escherichi, O. hainanensis, O. javanicus, O. longignathus, O. malaccensis, O. oblongatus, O. paraoblongatus, O. sarawakensis, and three possibly new species. Additionally, some of our taxa are thought to comprise a complex of two or more species. The number of species found in this study using DNA methods was more than the initial nine species thought to occur in Peninsular Malaysia. The support values for the clades and morphology of the soldiers provided further evidence for the existence of eleven or more species. Higher resolution genetic markers such as microsatellites would be required to confirm the presence of cryptic species in some taxa.
PCR has been extensively used for amplification of DNA sequences. We conducted a study to obtain the best amplification conditions for cytochrome b (Cyt b), cytochrome c oxidase I (COI) and 12S rRNA (12S) gene fragments of Malayan gaur mtDNA. DNA from seven Malayan gaur samples were extracted for PCR amplification. Various trials and combinations were tested to determine the best conditions of PCR mixture and profile to obtain the best PCR products for sequencing purposes. Four selected target factors for enhancing PCR, annealing temperature, concentration of primer pairs, amount of Taq polymerase, and PCR cycle duration, were optimized by keeping the amount of DNA template (50 ng/μL) and concentration of PCR buffer (1X), MgCl(2) (2.5 mM) and dNTP mixture (200 μM each) constant. All genes were successfully amplified, giving the correct fragment lengths, as assigned for both forward and reverse primers. The optimal conditions were determined to be: 0.1 μM primers for Cyt b and COI, 0.3 μM primers for 12S, 1 U Taq polymerase for all genes, 30 s of both denaturation and annealing cycles for Cyt b, 1 min of both stages for 12S and COI and annealing temperature of 58.4 ° C for Cyt b, 56.1 ° C for 12S and 51.3 ° C for COI. PCR products obtained under these conditions produced excellent DNA sequences.
We examined genetic variation in blood cockles in an effort to obtain information useful for the sustainability, management, and the stability of this species as a major commodity in the fisheries sector. Ten populations of cockles were sampled from the north to the south of the west coast of peninsular Malaysia. The cockles were collected in collaboration with the Fisheries Research Institute, Penang. The population genetic analysis of the cockles were studied via RAPD-PCR and mtDNA sequencing. Three hundred individuals were analyzed with RAPD-PCR experiments. High gene diversity over all loci was observed (Shannon index = 0.549 ± 0.056 and Nei's gene diversity = 0.4852 ± 0.0430 among 35 loci). The second method, mtDNA sequencing, was employed to complement the information obtained from RAPD-PCR. The gene selected for mtDNA sequencing was cytochrome c oxidase I (COI). One hundred and fifty individuals were sequenced, yielding a partial gene of 585 bp. Statistical analysis showed homogeneity in general but did reveal some degree of variability between the populations in Johor and the rest of the populations. The Mantel test showed a positive but nonsignificant correlation between geographic and genetic distances (r = 0.2710, P = 0.622), as in the RAPD analysis. We propose that the homogeneity between distant populations is caused by two factors: 1) the translocation of the spats; 2) larvae are carried by current movement from the north of the peninsula to the south. The different genetic composition found in Johor could be due to pollution, mutagenic substances or physical factors such as the depth of the water column. This population genetic study is the first for this species in peninsular Malaysia. The data from this study have important implications for fishery management, conservation of blood cockles and translocation policies for aquaculture and stock enhancement programs.
Recent morphological and molecular studies led to the recognition of two extant species of clouded leopards; Neofelis nebulosa from mainland southeast Asia and Neofelis diardi from the Sunda Islands of Borneo and Sumatra, including the Batu Islands. In addition to these new species-level distinctions, preliminary molecular data suggested a genetic substructure that separates Bornean and Sumatran clouded leopards, indicating the possibility of two subspecies of N. diardi. This suggestion was based on an analysis of only three Sumatran and seven Bornean individuals. Accordingly, in this study we re-evaluated this proposed subspecies differentiation using additional molecular (mainly historical) samples of eight Bornean and 13 Sumatran clouded leopards; a craniometric analysis of 28 specimens; and examination of pelage morphology of 20 museum specimens and of photographs of 12 wild camera-trapped animals. Molecular (mtDNA and microsatellite loci), craniomandibular and dental analyses strongly support the differentiation of Bornean and Sumatran clouded leopards, but pelage characteristics fail to separate them completely, most probably owing to small sample sizes, but it may also reflect habitat similarities between the two islands and their recent divergence. However, some provisional discriminating pelage characters are presented that need further testing. According to our estimates both populations diverged from each other during the Middle to Late Pleistocene (between 400 and 120 kyr). We present a discussion on the evolutionary history of Neofelis diardi sspp. on the Sunda Shelf, a revised taxonomy for the Sunda clouded leopard, N. diardi, and formally describe the Bornean subspecies, Neofelis diardi borneensis, including the designation of a holotype (BM.3.4.9.2 from Baram, Sarawak) in accordance with the rules of the International Code of Zoological Nomenclature.
The Asian subterranean termite, Coptotermes gestroi, originally from northeast India through Burma, Thailand, Malaysia, and the Indonesian archipelago, is a major termite pest introduced in several countries around the world, including Brazil. We sequenced the mitochondrial COII gene from individuals representing 23 populations. Phylogenetic analysis of COII gene sequences from this and other studies resulted in two main groups: (1) populations of Cleveland (USA) and four populations of Malaysia and (2) populations of Brazil, four populations of Malaysia, and one population from each of Thailand, Puerto Rico, and Key West (USA). Three new localities are reported here, considerably enlarging the distribution of C. gestroi in Brazil: Campo Grande (state of Mato Grosso do Sul), Itajaí (state of Santa Catarina), and Porto Alegre (state of Rio Grande do Sul).
By investigating genealogical relationships, we estimated the phylogenetic history and biogeography in the megophryid genus Leptobrachium (sensu lato, including Vibrissaphora) from southern China, Indochina, Thailand and the Sundaland. The genealogical relationships among the 30 named and unnamed taxa were estimated using 2009 bp of sequences from the mitochondrial DNA genes 12S rRNA, tRNA(val), and 16S rRNA using maximum parsimony, maximum likelihood, and Bayesian inference methods. The genus Leptobrachium was a well-supported monophyletic group that contained two major clades. One clade had three subclades primarily from disjunct regions including Borneo, Peninsular Malaysia and Java, and Thailand. The Bornean subclade included one species each from the Philippines and Sumatra. The other major clade consisted of two subclades, one from Indochina and the other from southern China (Vibrissaphora). Divergence times estimated an old evolutionary history of each subclade, one that could not be explained by the geohistory of Southeast Asian major landmasses.