The genus Myotis includes the largest number of species in the family Vespertilionidae (Chiroptera), and its members are distributed throughout most of the world. To re-evaluate the phylogenetic position of East Asian Myotis species with respect to Myotis species worldwide, we analyzed mitochondrial gene sequences of NADH dehydrogenase subunit 1 and cytochrome b from 24 East Asian individuals as well as 42 vespertilionid bats determined previously. The results suggest that: (1) some individuals having the same species name in Europe and Japan do not form a monophyletic clade, indicating that some bat species exhibit morphological convergence, (2) Japanese Myotis mystacinus forms a sister relationship with Myotis brandtii (Palaearctic), and both species are included in the American clade implying that an ancestor of these species originated in North America, and (3) the Black whiskered bat, Myotis pruinosus, is endemic to Japan and forms sister relationships with Myotis yanbarensis and Myotis montivagus collected from Okinawa (Japan) and Selangor (Malaysia), respectively, implying that M. pruinosus originated from the south. The systematics of Japanese and East Asian Myotis bats were revisited by considering their phylogenetic relationships. Our study provides the first extensive phylogenetic hypothesis of the genus Myotis that includes East Asian and Japanese species.
Vivipary with precocious seedlings in mangrove plants was thought to be a hindrance to long-range dispersal. To examine the extent of seedling dispersal across oceans, we investigated the phylogeny and genetic structure among East Asiatic populations of Kandelia candel based on organelle DNAs. In total, three, 28 and seven haplotypes of the chloroplast DNA (cpDNA) atpB-rbcL spacer, cpDNA trnL-trnF spacer, and mitochondrial DNA (mtDNA) internal transcribed spacer (ITS) were identified, respectively, from 202 individuals. Three data sets suggested consistent phylogenies recovering two differentiated lineages corresponding to geographical regions, i.e. northern South-China-Sea + East-China-Sea region and southern South-China-Sea region (Sarawak). Phylogenetically, the Sarawak population was closely related to the Ranong population of western Peninsula Malaysia instead of other South-China-Sea populations, indicating its possible origin from the Indian Ocean Rim. No geographical subdivision was detected within the northern geographical region. An analysis of molecular variance (AMOVA) revealed low levels of genetic differentiation between and within mainland and island populations (phiCT = 0.015, phiSC = 0.037), indicating conspicuous long-distance seedling dispersal across oceans. Significant linkage disequilibrium excluded the possibility of recurrent homoplasious mutations as the major force causing phylogenetic discrepancy between mtDNA and the trnL-trnF spacer within the northern region. Instead, relative ages of alleles contributed to non-random chlorotype-mitotype associations and tree inconsistency. Widespread distribution and random associations (chi2 = 0.822, P = 0.189) of eight hypothetical ancestral cytotypes indicated the panmixis of populations of the northern geographical region as a whole. In contrast, rare and recently evolved alleles were restricted to marginal populations, revealing some preferential directional migration.
Recent electrophoretic data have indicated that Schistosoma japonicum in mainland China may be a species complex, with the existence of a cryptic species being predicted from the analysis of schistosome populations from Sichuan province. To investigate the Sichuan form of S. japonicum, 4.9 kbp of mitochondrial DNA from each of three samples of the parasite from China (two from Sichuan and one from Hunan) and one from Sorsogon in the Philippines were amplified, sequenced and characterized. The sequence data were compared with those from the related South-east Asian species of S. mekongi (Khong Island, Laos) and S. mlayensis (Baling, Malaysia) and that from S. japonicm from Anhui (China). At both the nucleotide and amino-acid levels, the variation among the five S. japonicum samples was limited (< 1%). This was consistent with the conclusions drawn from previous molecular studies, in which minimal variation among S. japonicum populations was also detected. In contrast, S. mekongi and S. malayensis, species recognized as separate but closely related, differ from each other by about 10%, and each differs by 25%-26% from S. japonicum. Phylogenetic trees provided a graphic representation of these differences, showing all S. japonicum sequences to be very tightly clustered and distant from S. mekongi and S. malayensis, the last two being clearly distinct from each other. The results thus indicate no significant intra-specific genetic variation among S. japonicum samples collected from different geographical areas and do not support the idea of a distinct form in Sichuan.
Populations of the Asian elephant (Elephas maximus) have been reduced in size and become highly fragmented during the past 3,000 to 4,000 years. Historical records reveal elephant dispersal by humans via trade and war. How have these anthropogenic impacts affected genetic variation and structure of Asian elephant populations? We sequenced mitochondrial DNA (mtDNA) to assay genetic variation and phylogeography across much of the Asian elephant's range. Initially we compare cytochrome b sequences (cyt b) between nine Asian and five African elephants and use the fossil-based age of their separation (approximately 5 million years ago) to obtain a rate of about 0.013 (95% CI = 0.011-0.018) corrected sequence divergence per million years. We also assess variation in part of the mtDNA control region (CR) and adjacent tRNA genes in 57 Asian elephants from seven countries (Sri Lanka, India, Nepal, Myanmar, Thailand, Malaysia, and Indonesia). Asian elephants have typical levels of mtDNA variation, and coalescence analyses suggest their populations were growing in the late Pleistocene. Reconstructed phylogenies reveal two major clades (A and B) differing on average by HKY85/gamma-corrected distances of 0.020 for cyt b and 0.050 for the CR segment (corresponding to a coalescence time based on our cyt b rate of approximately 1.2 million years). Individuals of both major clades exist in all locations but Indonesia and Malaysia. Most elephants from Malaysia and all from Indonesia are in well-supported, basal clades within clade A. thus supporting their status as evolutionarily significant units (ESUs). The proportion of clade A individuals decreases to the north, which could result from retention and subsequent loss of ancient lineages in long-term stable populations or, perhaps more likely, via recent mixing of two expanding populations that were isolated in the mid-Pleistocene. The distribution of clade A individuals appears to have been impacted by human trade in elephants among Myanmar, Sri Lanka, and India, and the subspecies and ESU statuses of Sri Lankan elephants are not supported by molecular data.
For elucidation of the taxonomic status of the Japanese Fasciola species, whole mitochondrial DNA of Fasciola hepatica from Australia, F. gigantica from Malaysia, and Fasciola sp. from Japan was digested with three four-base-cutting endonucleases: HinfI, MspI, and RsaI. The resulting digestion patterns showed that for each enzyme there were some bands specific for each geographical isolate and that the Japanese Fasciola sp. shared more bands with F. gigantica than with F. hepatica. Nucleotide sequences of two regions, the second internal transcribed spacer (ITS2) of the nuclear ribosomal RNA cluster and mitochondrial cytochrome c oxidase subunit I (COI), were also compared among them. The ITS2 sequence was highly conserved among the three isolates. F. gigantica and the Japanese Fasciola sp. were identical, but they differed from the Australian F. hepatica at six sites, one of which was a deletion. The COI sequence was less conserved but implied a similar relationship between the isolates. There seems no reason to regard the Japanese Fasciola sp. as anything other than a strain of F. gigantica.
Illegal hunting is a major threat to the elephants of Africa, with more elephants killed by poachers than die from natural causes. DNA from tusks has been used to infer the source populations for confiscated ivory, relying on nuclear genetic markers. However, mitochondrial DNA (mtDNA) sequences can also provide information on the geographic origins of elephants due to female elephant philopatry. Here, we introduce the Loxodonta Localizer (LL; www.loxodontalocalizer.org), an interactive software tool that uses a database of mtDNA sequences compiled from previously published studies to provide information on the potential provenance of confiscated ivory. A 316 bp control region sequence, which can be readily generated from DNA extracted from ivory, is used as a query. The software generates a listing of haplotypes reported among 1917 African elephants in 24 range countries, sorted in order of similarity to the query sequence. The African locations from which haplotype sequences have been previously reported are shown on a map. We demonstrate examples of haplotypes reported from only a single locality or country, examine the utility of the program in identifying elephants from countries with varying degrees of sampling, and analyze batches of confiscated ivory. The LL allows for the source of confiscated ivory to be assessed within days, using widely available molecular methods that do not depend on a particular platform or laboratory. The program enables identification of potential regions or localities from which elephants are being poached, with capacity for rapid identification of populations newly or consistently targeted by poachers.
The uniparentally inherited mitochondrial DNA (mtDNA) is in the limelight for the past two decades, in studies relating to demographic history of mankind and in forensic kinship testing. In this study, human mtDNA hypervariable segments 1, 2, and 3 (HV1, HV2, and HV3) were analyzed in 248 unrelated Malay individuals in Peninsular Malaysia. Combined analyses of HV1, HV2, and HV3 revealed a total of 180 mtDNA haplotypes with 149 unique haplotypes and 31 haplotypes occurring in more than one individual. The genetic diversity was estimated to be 99.47%, and the probability of any two individuals sharing the same mtDNA haplotype was 0.93%. The most frequent mtDNA haplotype (73, 146, 150, 195, 263, 315.1C, 16140, 16182C, 16183C, 16189, 16217, 16274, and 16335) was shared by 11 (4.44%) individuals. The nucleotide diversity and mean of pair-wise differences were found to be 0.036063 ± 0.020101 and 12.544022 ± 6.230486, respectively.
Control region polymorphisms in the mitochondrial DNA of 124 unrelated individuals from the Malay population living in or around Kuala Lumpur in Malaysia were investigated and phylogenetic haplogroup lineages were determined. The intergenic COII/tRNALys 9-bp deletion, 3010 and 5178 mutations, and several coding region polymorphisms were examined to discriminate some phylogenetic haplogroups. Sequence comparison of the control regions led to the identification of 117 mitochondrial haplotypes, in which 103 types were observed in only one individual and the other nine types were shared by more than two individuals. Gene diversity was estimated to be 0.997. Phylogenetic haplogroup determination revealed that the gene pool of the modern Malay population in Malaysia consisted mainly of southeast Asian, east Asian, unidentified and unique, and aboriginal southeast-specific haplogroups. These results suggest a multi-original nature for the modern Malay population. The present database may help not only in personal identification but also in determining geographic origin in forensic casework in Malaysian, Southeast Asian and East Asian populations.
The emergence of Massively Parallel Sequencing technologies enabled the analysis of full mitochondrial (mt)DNA sequences from forensically relevant samples that have, so far, only been typed in the control region or its hypervariable segments. In this study, we evaluated the performance of a commercially available multiplex-PCR-based assay, the Precision ID mtDNA Whole Genome Panel (Thermo Fisher Scientific), for the amplification and sequencing of the entire mitochondrial genome (mitogenome) from even degraded forensic specimens. For this purpose, more than 500 samples from 24 different populations were selected to cover the vast majority of established superhaplogroups. These are known to harbor different signature sequence motifs corresponding to their phylogenetic background that could have an effect on primer binding and, thus, could limit a broad application of this molecular genetic tool. The selected samples derived from various forensically relevant tissue sources and were DNA extracted using different methods. We evaluated sequence concordance and heteroplasmy detection and compared the findings to conventional Sanger sequencing as well as an orthogonal MPS platform. We discuss advantages and limitations of this approach with respect to forensic genetic workflow and analytical requirements.
The Americas were the last inhabitable continents to be occupied by humans, with a growing multidisciplinary consensus for entry 15-25 thousand years ago (kya) from northeast Asia via the former Beringia land bridge [1-4]. Autosomal DNA analyses have dated the separation of Native American ancestors from the Asian gene pool to 23 kya or later [5, 6] and mtDNA analyses to ∼25 kya [7], followed by isolation ("Beringian Standstill" [8, 9]) for 2.4-9 ky and then a rapid expansion throughout the Americas. Here, we present a calibrated sequence-based analysis of 222 Native American and relevant Eurasian Y chromosomes (24 new) from haplogroups Q and C [10], with four major conclusions. First, we identify three to four independent lineages as autochthonous and likely founders: the major Q-M3 and rarer Q-CTS1780 present throughout the Americas, the very rare C3-MPB373 in South America, and possibly the C3-P39/Z30536 in North America. Second, from the divergence times and Eurasian/American distribution of lineages, we estimate a Beringian Standstill duration of 2.7 ky or 4.6 ky, according to alternative models, and entry south of the ice sheet after 19.5 kya. Third, we describe the star-like expansion of Q-M848 (within Q-M3) starting at 15 kya [11] in the Americas, followed by establishment of substantial spatial structure in South America by 12 kya. Fourth, the deep branches of the Q-CTS1780 lineage present at low frequencies throughout the Americas today [12] may reflect a separate out-of-Beringia dispersal after the melting of the glaciers at the end of the Pleistocene.
Confiscated slow lorises (Nycticebus spp.) at Bangpra Water-Bird Breeding Center (BWBC) in Thailand provided an opportunity to demonstrate the application of noninvasive genetic approaches for species identification when morphology of the animals was ambiguous. The slow lorises at BWBC had been assigned to either N. bengalensis or N. pygmaeus, based on body size. However, the morphology of N. bengalensis is highly variable and overlaps with that of N. coucang (sensu stricto). Phylogenetic analysis of cytochrome b and d-loop mitochondrial regions placed all confiscated N. pygmaeus with the published sequences of N. pygmaeus and distinguished them from other Nycticebus. All other confiscated individuals formed a monophyletic clade, most individuals grouping with published N. bengalensis sequences from wild populations in Vietnam and distinct from Peninsular Malaysian and Sumatran N. coucang, Javan N. javanicus and Bornean N. menagensis. Six individuals within the N. bengalensis clade formed a separate subgroup that did not group with any reference material as indicated by phylogenetic and haplotype network analyses. Whether these trafficked individuals are undiscovered wild populations will require further investigation. Additional genetic studies of wild slow loris populations in different regions are therefore urgently required for reference to aid the protection and conservation of these threatened species.
Southeast Asia (SEA) is enriched with a complex history of peopling. Malaysia, which is located at the crossroads of SEA, has been recognized as one of the hubs for early human migration. To unravel the genomic complexity of the native inhabitants of Malaysia, we sequenced 12 samples from 3 indigenous populations from Peninsular Malaysia and 4 native populations from North Borneo to a high coverage of 28-37×. We showed that the Negritos from Peninsular Malaysia shared a common ancestor with the East Asians, but exhibited some level of gene flow from South Asia, while the North Borneo populations exhibited closer genetic affinity towards East Asians than the Malays. The analysis of time of divergence suggested that ancestors of Negrito were the earliest settlers in the Malay Peninsula, whom first separated from the Papuans ~ 50-33 thousand years ago (kya), followed by East Asian (~ 40-15 kya), while the divergence time frame between North Borneo and East Asia populations predates the Austronesian expansion period implies a possible pre-Neolithic colonization. Substantial Neanderthal ancestry was confirmed in our genomes, as was observed in other East Asians. However, no significant difference was observed, in terms of the proportion of Denisovan gene flow into these native inhabitants from Malaysia. Judging from the similar amount of introgression in the Southeast Asians and East Asians, our findings suggest that the Denisovan gene flow may have occurred before the divergence of these populations and that the shared similarities are likely an ancestral component.
The present investigation was carried out in an attempt to study the phylogenetic analysis of different breeds of domestic chickens in Peninsular Malaysia inferred from partial cytochrome b gene information and random amplified polymorphic DNA (RAPD) markers. Phylogenetic analysis using both neighbor-joining (NJ) and maximum parsimony (MP) methods produced three clusters that encompassed Type-I village chickens, the red jungle fowl subspecies and the Japanese Chunky broilers. The phylogenetic analysis also revealed that majority of the Malaysian commercial chickens were randomly assembled with the Type-II village chickens. In RAPD assay, phylogenetic analysis using neighbor-joining produced six clusters that were completely distinguished based on the locality of chickens. High levels of genetic variations were observed among the village chickens, the commercial broilers, and between the commercial broilers and layer chickens. In this study, it was found that Type-I village chickens could be distinguished from the commercial chickens and Type-II village chickens at the position of the 27th nucleotide of the 351 bp cytochrome b gene. This study also revealed that RAPD markers were unable to differentiate the type of chickens, but it showed the effectiveness of RAPD in evaluating the genetic variation and the genetic relationships between chicken lines and populations.
Surrounded by speakers of Indo-European, Dravidian and Tibeto-Burman languages, around 11 million Munda (a branch of Austroasiatic language family) speakers live in the densely populated and genetically diverse South Asia. Their genetic makeup holds components characteristic of South Asians as well as Southeast Asians. The admixture time between these components has been previously estimated on the basis of archaeology, linguistics and uniparental markers. Using genome-wide genotype data of 102 Munda speakers and contextual data from South and Southeast Asia, we retrieved admixture dates between 2000-3800 years ago for different populations of Munda. The best modern proxies for the source populations for the admixture with proportions 0.29/0.71 are Lao people from Laos and Dravidian speakers from Kerala in India. The South Asian population(s), with whom the incoming Southeast Asians intermixed, had a smaller proportion of West Eurasian genetic component than contemporary proxies. Somewhat surprisingly Malaysian Peninsular tribes rather than the geographically closer Austroasiatic languages speakers like Vietnamese and Cambodians show highest sharing of IBD segments with the Munda. In addition, we affirmed that the grouping of the Munda speakers into North and South Munda based on linguistics is in concordance with genome-wide data.
The proboscis monkey, Nasalis larvatus, is an endemic species to the island of Borneo. It is listed in the IUCN Red List as Endangered with a decreasing population trend. Nevertheless, biological information, especially on the genetic diversity of the species, is still incomplete. Its fragmented distribution poses difficulties in gathering genetic samples along with its widespread distribution across Borneo. This study aims to determine the genetic variation and structure of N. larvatus with an emphasis on Malaysian Borneo populations to elucidate its gene flow. The genetic variation and structure of N. larvatus were examined using 50 sequences of the 1,434-bp cytochrome oxidase subunit I (COI) gene region of mitochondrial DNA. The COI sequences revealed low genetic variation among N. larvatus populations in Malaysian Borneo. This low genetic variability could be the result of inbreeding pressure that may have occurred due to the absence of population expansion in this species over the last 30,000 years. This is supported in our analysis of molecular variance, which showed that groups of N. larvatus are significantly differentiated possibly due to natural geographic barriers. This study provides baseline information on the genetic diversity among proboscis monkey populations in Borneo for the future genetic assessment of the species.
A total of 74 shrimp specimens were sequenced at a 584 bp segment of the cytochrome oxidase subunit I (COI) gene to examine patterns of DNA barcode variation in a mangrove biodiversity hotspot. The Maximum Likelihood tree, barcode gap analysis, Automatic Barcode Gap Discovery analysis and sequence comparisons with data available from Barcode of Life Data System and GenBank recovered 18 taxa of which 15 were identified to species level, 2 at genus level and a single taxon at order level. Two deep mitochondrial DNA lineage divergences were found in the giant tiger prawn, Penaeus monodon. It is suggested that one of the lineages is a consequence of an introduction from aquaculture activity. These results have provided a reliable barcode library for cataloguing shrimps in this area.
The genetic diversity of the endangered crocodile Tomistoma schlegelii was characterized using the protein coding ND 6-tRNA(glu)-cyt b and the cytochrome b-control region (cyt b-CR) markers. Concatenate data revealed six haplotypes with an overall haplotype diversity of 0.769 ± 0.039; nucleotide diversity was 0.00535 ± 0.00172. A nearest-neighbor analysis showed that all individuals clustered with four geographic regions (Sumatra, Peninsular Malaysia, Sarawak, and East Kalimantan) and were genetically differentiated. With the exception of the individuals from haplotype H2, which occurred in both Peninsular Malaysia and Sarawak, all other haplotypes were geographically distinct. The H4 lineage, which was found to be the most divergent, clustered exclusively in the basal clade in all phylogenetic trees, and the haplotype network was unconnected at the 95% reconnection limit, suggesting further investigation to establish its possible status as a distinct evolutionary significant unit or a cryptic species.
In this study the genetic diversity of local freshwater leeches (Hirudinaria spp.) was inferred using mtDNA COI gene analysis and compared with the gross external variations of 26 freshwater leech specimens obtained from the wild and leech farms. Based on a neighbor-joining tree generated from 516 COI base sequences, four distinct clades of Hirudinaria were seen with interspecific genetic divergence in the range of 7.6-14.5%. The external morphological variations based on the presence of stripes, location of gonopores, and anus separated the samples into four morphologically distinct groups matching the four clades obtained from the molecular data. Two black stripes at the ventral region were observed only in specimens found clustered with clades that contained the GenBank-reported H. manillensis, whereas the brown or dark green coloration without stripes on the ventral region was seen in samples that clustered with H. javanica and H. bpling clades.
The true mahseer (Tor spp.) is one of the highest valued fish in the world due to its high nutritional value and great unique taste. Nevertheless, its morphological characterization and single mitochondrial gene phylogeny in the past had yet to resolve the ambiguity in its taxonomical classification. In this study, we sequenced and assembled 11 complete mahseer mitogenomes collected from Java of Indonesia, Pahang and Terengganu of Peninsular Malaysia as well as Sarawak of East Malaysia. The mitogenome evolutionary relationships among closely related Tor spp. samples were investigated based on maximum likelihood phylogenetic tree construction. Compared to the commonly used COX1 gene fragment, the complete COX1, Cytb, ND2, ND4 and ND5 genes appear to be better phylogenetic markers for genetic differentiation at the population level. In addition, a total of six population-specific mitolineage haplotypes were identified among the mahseer samples analyzed, which this offers hints towards its taxonomical landscape.
Studying the genetic history of the Orang Asli of Peninsular Malaysia can provide crucial clues to the peopling of Southeast Asia as a whole. We have analyzed mitochondrial DNA (mtDNAs) control-region and coding-region markers in 447 mtDNAs from the region, including 260 Orang Asli, representative of each of the traditional groupings, the Semang, the Senoi, and the Aboriginal Malays, allowing us to test hypotheses about their origins. All of the Orang Asli groups have undergone high levels of genetic drift, but phylogeographic traces nevertheless remain of the ancestry of their maternal lineages. The Semang have a deep ancestry within the Malay Peninsula, dating to the initial settlement from Africa >50,000 years ago. The Senoi appear to be a composite group, with approximately half of the maternal lineages tracing back to the ancestors of the Semang and about half to Indochina. This is in agreement with the suggestion that they represent the descendants of early Austroasiatic speaking agriculturalists, who brought both their language and their technology to the southern part of the peninsula approximately 4,000 years ago and coalesced with the indigenous population. The Aboriginal Malays are more diverse, and although they show some connections with island Southeast Asia, as expected, they also harbor haplogroups that are either novel or rare elsewhere. Contrary to expectations, complete mtDNA genome sequences from one of these, R9b, suggest an ancestry in Indochina around the time of the Last Glacial Maximum, followed by an early-Holocene dispersal through the Malay Peninsula into island Southeast Asia.