This study utilized genetic and morphometric approaches to assess the molecular and morphometric differentiation among commercially important species of mud crab. Molecular investigations were derived from 542 bp mitochondrial DNA COI on 249 individuals within genus Scylla from nine states in Malaysia represents four marine regions; South China Sea, Sulu Sea, Straits of Singapore and Straits of Malacca. Four specimens were obtained from Indonesia to give a robust analysis in this study. For species delimitation, Automatic Barcode Gap Discovery (ABGD) method on a web interface was employed. Analysis on phylogenetics was implemented utilizing Neighbour joining (NJ) and Maximum Parsimony (MP) methods. The inter- and intraspecies genetic distances (Ds ) was computed using Kimura 2-parameter distance and executed in MEGA version 5.05. All samples were genetically and morphologically identified and clustered into four distinct species. Among the species, S. olivacea was the most abundant (n = 111), on the other hand the occurrence of S. paramamosain in Malaysia was very low (n = 29). No single individual of S. serrata from Malaysia was recorded in this study. Both genetic distance and phylogenetic approaches exhibited a correlative monophyletic association among all specimens analysed. This present study is crucial as it reports the reassessment of all species within genus Scylla in Malaysia, eventually could be employed as a reference source for subsequent research mainly on mariculture and other conservation efforts for the species.
We analyzed the complete mitochondrial cytochrome b (cytb) gene and fragments of four nuclear loci: ApoB, RAG2, IRBP1 and BRCA1. These data allowed us to provide new insights into the diversity of the Asiatic water shrews of Indochina. A new, highly divergent genetic lineage of Chimarrogale was found in southern Vietnam, and this lineage included specimens from the provinces of Kon Tum, Dak Lak, and Lam Dong. Such finding represents the newest and southernmost records of Chimarrogale in Indochina. Morphological analysis classified the specimens from southern Vietnam as C. varennei proper, which is restricted to that region, whereas the polymorphic C. himalayica, which contained at least four cytochrome b haplogroups, occurred in central and northern Vietnam and southern China. This distinct C. varennei lineage closely related to the C. platycephalus + C. leander clade suggests the existence of an unknown glacial refuge in Tay Nguyen Plateau, southern Vietnam. Because the Bornean C. phaeura (i) was sister-group of the rest of Chimarrogale sensu lato and (ii) had a high genetic divergence (~15% for cytochrome b) and geographical isolation, we suggest that C. phaeura be placed into a separate genus, Crossogale Thomas, 1921. This genus should also include C. sumatrana (Sumatra) and C. hantu (Peninsular Malaysia). On those grounds, we propose a new classification system for Asiatic water shrews.
Extrachromosomal (ec) DNA in eukaryotic cells has been known for decades. The structures described range from linear double stranded (ds) DNA to circular dsDNA, distinct from mitochondrial (mt) DNA. The sizes of circular forms are described from some hundred base pairs (bp) up to more than 150 kbp. The number of molecules per cell ranges from several hundred to a thousand. Semi-quantitative determinations of circular dsDNA show proportions as high as several percentages of the total DNA per cell. These ecDNA fractions harbor sequences that are known to be present in chromosomal DNA (chrDNA) too. Sequencing projects on, for example the human genome, have to take into account the ecDNA sequences which are simultaneously ascertained; corrections cannot be performed retrospectively. Concerning the results of sequencings derived from extracted whole DNA: if the ecDNA fractions contained therein are not taken into account, erroneous conclusions at the chromosomal level may result.
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.
A clear understanding of population structure is essential for assessing conservation status and implementing management strategies. A small, non-migratory population of humpback whales in the Arabian Sea is classified as "Endangered" on the IUCN Red List of Threatened Species, an assessment constrained by a lack of data, including limited understanding of its relationship to other populations. We analysed 11 microsatellite markers and mitochondrial DNA sequences extracted from 67 Arabian Sea humpback whale tissue samples and compared them to equivalent datasets from the Southern Hemisphere and North Pacific. Results show that the Arabian Sea population is highly distinct; estimates of gene flow and divergence times suggest a Southern Indian Ocean origin but indicate that it has been isolated for approximately 70,000 years, remarkable for a species that is typically highly migratory. Genetic diversity values are significantly lower than those obtained for Southern Hemisphere populations and signatures of ancient and recent genetic bottlenecks were identified. Our findings suggest this is the world's most isolated humpback whale population, which, when combined with low population abundance estimates and anthropogenic threats, raises concern for its survival. We recommend an amendment of the status of the population to "Critically Endangered" on the IUCN Red List.
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 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.
It has recently been shown that ancestors of New Guineans and Bougainville Islanders have inherited a proportion of their ancestry from Denisovans, an archaic hominin group from Siberia. However, only a sparse sampling of populations from Southeast Asia and Oceania were analyzed. Here, we quantify Denisova admixture in 33 additional populations from Asia and Oceania. Aboriginal Australians, Near Oceanians, Polynesians, Fijians, east Indonesians, and Mamanwa (a "Negrito" group from the Philippines) have all inherited genetic material from Denisovans, but mainland East Asians, western Indonesians, Jehai (a Negrito group from Malaysia), and Onge (a Negrito group from the Andaman Islands) have not. These results indicate that Denisova gene flow occurred into the common ancestors of New Guineans, Australians, and Mamanwa but not into the ancestors of the Jehai and Onge and suggest that relatives of present-day East Asians were not in Southeast Asia when the Denisova gene flow occurred. Our finding that descendants of the earliest inhabitants of Southeast Asia do not all harbor Denisova admixture is inconsistent with a history in which the Denisova interbreeding occurred in mainland Asia and then spread over Southeast Asia, leading to all its earliest modern human inhabitants. Instead, the data can be most parsimoniously explained if the Denisova gene flow occurred in Southeast Asia itself. Thus, archaic Denisovans must have lived over an extraordinarily broad geographic and ecological range, from Siberia to tropical Asia.
Haruan (Channa striatus) is in great demand in the Malaysian domestic fish market. In the present study, mtDNA cyt b was used to investigate genetic variation of C. striatus among populations in Peninsular Malaysia. The overall population of C. striatus demonstrated a high level of haplotype diversity (h) and a low-to-moderate level of nucleotide diversity (π). Analysis of molecular variance (AMOVA) results showed a significantly different genetic differentiation among 6 populations (F(ST)=0.37566, P=0.01). Gene flow (Nm) was high and ranged from 0.32469 to infinity (∞). No significant relationship between genetic distance and geographic distance was detected. A UPGMA tree based on the distance matrix of net interpopulation nucleotide divergence (d(A)) and haplotype network of mtDNA cyt b revealed that C. striatus is divided into 2 major clades. The neutrality and mismatch distribution tests for all populations suggested that C. striatus in the study areas had undergone population expansion. The estimated time of population expansion in the mtDNA cyt b of C. striatus populations occurred 0.72-6.19 million years ago. Genetic diversity of mtDNA cyt b and population structure among Haruan populations in Peninsular Malaysia will be useful in fisheries management for standardization for Good Agriculture Practices (GAP) in fish-farming technology, as well as providing the basis for Good Manufacturing Practices (GMP).
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.