The complex climatic and geological history of Southeast Asia has shaped this region's high biodiversity. In particular, sea level fluctuations associated with repeated glacial cycles during the Pleistocene both facilitated, and limited, connectivity between populations. In this study, we used data from two mitochondrial and three anonymous nuclear markers to determine whether a fresh/brackish water killifish, Aplocheilus panchax, Hamilton, 1822, could be used to further understand how climatic oscillations and associated sea level fluctuations have shaped the distribution of biota within this region, and whether such patterns show evidence of isolation within palaeodrainage basins. Our analyses revealed three major mitochondrial clades within A. panchax. The basal divergence of A. panchax mitochondrial lineages was approximately 3.5 Ma, whilst the subsequent divergence timings of these clades occurred early Pleistocene (~2.6 Ma), proceeding through the Pleistocene. Continuous phylogeographic analysis showed a clear west-east dispersal followed by rapid radiation across Southeast Asia. Individuals from Krabi, just north of the Isthmus of Kra, were more closely related to the Indian lineages, providing further evidence for a freshwater faunal disjunction at the Isthmus of Kra biogeographic barrier. Our results suggest that Sulawesi, across the Wallace Line, was colonised relatively recently (~30 ka). Nuclear DNA is less geographically structured, although Mantel tests indicated that nuclear genetic distances were correlated with geographic proximity. Overall, these results imply that recent gene flow, as opposed to historical isolation, has been the key factor determining patterns of nuclear genetic variation in A. panchax, however, some evidence of historical isolation is retained within the mitochondrial genome. Our study further validates the existence of a major biogeographic boundary at the Kra Isthmus, and also demonstrates the use of widely distributed fresh/brackishwater species in phylogeographic studies, and their ability to disperse across major marine barriers in relatively recent time periods.
The multi-racial, multi-lingual, multi-religious, and multi-cultural Malaysia is situated at the crossroads of Southeast Asia. It has a total population of about 23 million, comprising of more than fifty ethnic groups: Malays, Chinese, Indian, and the minorities of Orang Asli in Peninsular; Iban, Bidayuh, Melanau, Kenyah Kayan, Ukit, Penan, Sekapan, Lahanan, Lun Bawang, Kelabit, Berawan, Punan Bah etc. in Sarawak; and Kadazan, Bajau, Murut, Paitan, Suluk Bonggi, Illanun, Bengkahak Tidung etc. in Sabah. The origin of the multi-ethnic character of this country traces back to decades of human migration from various regions of Southeast Asia. With her rich human biodiversity, a study of the human population genetics is imperative, either for forensic database purposes, or as the stepping stone for medical applications.
Ancestry-informative markers (AIMs) can be used to infer the ancestry of an individual to minimize the inaccuracy of self-reported ethnicity in biomedical research. In this study, we describe three methods for selecting AIM SNPs for the Malay population (Malay AIM panel) using different approaches based on pairwise FST, informativeness for assignment (In), and PCA-correlated SNPs (PCAIMs). These Malay AIM panels were extracted from genotype data stored in SNP arrays hosted by the Malaysian node of the Human Variome Project (MyHVP) and the Singapore Genome Variation Project (SGVP). In particular, genotype data from a total of 165 Malay individuals were analyzed, comprising data on 117 individual genotypes from the Affymetrix SNP-6 SNP array platform and data on 48 individual genotypes from the OMNI 2.5 Illumina SNP array platform. The HapMap phase 3 database (1397 individuals from 11 populations) was used as a reference for comparison with the Malay genotype data. The accuracy of each resulting Malay AIM panel was evaluated using a machine learning "ancestry-predictive model" constructed by using WEKA, a comprehensive machine learning platform written in Java. A total of 1250 SNPs were finally selected, which successfully identified Malay individuals from other world populations with an accuracy of 90%, but the accuracy decreased to 80% using 157 SNPs according to the pairwise FST method, while a panel of 200 SNPs selected using In and PCAIMs could be used to identify Malay individuals with an accuracy of approximately 80%.
Malay, the main ethnic group in Peninsular Malaysia, is represented by various sub-ethnic groups such as Melayu Banjar, Melayu Bugis, Melayu Champa, Melayu Java, Melayu Kedah Melayu Kelantan, Melayu Minang and Melayu Patani. Using data retrieved from the MyHVP (Malaysian Human Variome Project) database, a total of 135 individuals from these sub-ethnic groups were profiled using the Affymetrix GeneChip Mapping Xba 50-K single nucleotide polymorphism (SNP) array to identify SNPs that were ancestry-informative markers (AIMs) for Malays of Peninsular Malaysia. Prior to selecting the AIMs, the genetic structure of Malays was explored with reference to 11 other populations obtained from the Pan-Asian SNP Consortium database using principal component analysis (PCA) and ADMIXTURE. Iterative pruning principal component analysis (ipPCA) was further used to identify sub-groups of Malays. Subsequently, we constructed an AIMs panel for Malays using the informativeness for assignment (In) of genetic markers, and the K-nearest neighbor classifier (KNN) was used to teach the classification models. A model of 250 SNPs ranked by In, correctly classified Malay individuals with an accuracy of up to 90%. The identified panel of SNPs could be utilized as a panel of AIMs to ascertain the specific ancestry of Malays, which may be useful in disease association studies, biomedical research or forensic investigation purposes.
Many insect pests, including the brown planthopper (BPH), undergo windborne migration that is challenging to observe and track. It remains controversial about their migration patterns and largely unknown regarding the underlying genetic basis. By analyzing 360 whole genomes from around the globe, we clarify the genetic sources of worldwide BPHs and illuminate a landscape of BPH migration showing that East Asian populations perform closed-circuit journeys between Indochina and the Far East, while populations of Malay Archipelago and South Asia undergo one-way migration to Indochina. We further find round-trip migration accelerates population differentiation, with highly diverged regions enriching in a gene desert chromosome that is simultaneously the speciation hotspot between BPH and related species. This study not only shows the power of applying genomic approaches to demystify the migration in windborne migrants but also enhances our understanding of how seasonal movements affect speciation and evolution in insects.
Domestication to captive rearing conditions, along with targeted selective breeding have genetic consequences that vary from those in wild environments. Nile tilapia (Oreochromis niloticus) is one of the most translocated and farmed aquaculture species globally, farmed throughout Asia, North and South America, and its African native range. In Egypt, a breeding program established the Abbassa Strain of Nile tilapia (AS) in 2002 based on local broodstock sourced from the Nile River. The AS has been intensively selected for growth and has gone through genetic bottlenecks which have likely shifted levels and composition of genetic diversity within the strain. Consequently, there are questions on the possible genetic impact AS escapees may have on endemic populations of Nile tilapia. However, to date there have been no genetic studies comparing genetic changes in the domesticated AS to local wild populations. This study used 9,827 genome-wide SNPs to investigate population genetic structure and signatures of selection in the AS (generations 9-11) and eight wild Nile tilapia populations from Egypt. SNP analyses identified two major genetic clusters (captive and wild populations), with wild populations showing evidence of isolation-by-distance among the Nile Delta and upstream riverine populations. Between genetic clusters, approximately 6.9% of SNPs were identified as outliers with outliers identified on all 22 O. niloticus chromosomes. A lack of localized outlier clustering on the genome suggests that no genes of major effect were presently detected. The AS has retained high levels of genetic diversity (Ho_All = 0.21 ± 0.01; He_All = 0.23 ± 0.01) when compared to wild populations (Ho_All = 0.18 ± 0.01; He_All = 0.17 ± 0.01) after 11 years of domestication and selective breeding. Additionally, 565 SNPs were unique within the AS line. While these private SNPs may be due to domestication signals or founder effects, it is suspected that introgression with blue tilapia (Oreochromis aureus) has occurred. This study highlights the importance of understanding the effects of domestication in addition to wild population structure to inform future management and dissemination decisions. Furthermore, by conducting a baseline genetic study of wild populations prior to the dissemination of a domestic line, the effects of aquaculture on these populations can be monitored over time.
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.
Peninsular Malaysia is populated by the Malays, Chinese, Indians, and Orang Asli. We have analyzed 17 Y-STRs loci for 243 randomly unrelated individuals, which include 153 Malays (7 Acheh, 13 Champa, 11 Rawa, 9 Kedah, 23 Minang, 15 Bugis, 43 Kelantan, 14 Jawa, and 18 Bugis) and 90 Orang Asli [54 Semang (16 Kensiu, 13 Lanoh, 25 Bateq); 30 Senoi (21 Semai, 9 Che Wong); and 6 Proto-Malay (6 Orang Kanaq)] from selected settlements in Peninsular Malaysia using the AmpFlSTR Yfiler™ kit (Applied Biosystems™). The overall haplotype diversity is 0.9966, i.e., 0.9984 for the Malays and 0.9793 for the Orang Asli. A total of 158 haplotypes (65.02%) were individually unique. The p value and pairwise Rst analysis was calculated to show the genetic structure of the samples with other world populations (from YHRD website). Based on the Y-STR data, Champa, Acheh, Kedah, Minang, and Kelantan are clustered together. Lanoh and Kensiu (Semang) are very closely related, suggesting similar paternal ancestry. Jawa Malays and Indonesian Java, plus the Bugis Malays and Australian Aborigines shared high degree of paternal lineage affinity. This study presents data for very precious relict groups, who are the earliest inhabitants of Peninsular Malaysia.
"Bumiputra" or "son of the soil" is a term used to represent the Malays and other indigenous populations of Malaysia. The Malays are Austronesian speaking population and originated from different parts of the Indo-Malay Archipelago. The migration of Malay population from different parts of Indo-Malay Archipelago were mainly due to trading purposes which shaped the current Malay sub-ethnic groups with unique culture and with distinctive dialects. In this study, HLA typing was carried out using Sequence-based Typing (SBT) method on 109 individuals comprising of four Malay sub-ethnic groups namely Kelantan (n = 28), Champa (n = 29), Patani (n = 25) and Mandailing (n = 27) Malays. The HLA data is available in the Allele Frequencies Net Database (AFND).
Prevalence of hypertension (HTN) varies substantially across different populations. HTN is not only common - affecting at least one third of the world's adult population - but is also the most important driver for cardiovascular diseases. Yet up to a third of hypertensive patients are resistant to therapy, contributed by secondary hypertension but more commonly the hitherto inability to precisely predict response to specific antihypertensive agents. Population and individual genomics information could be useful in guiding the selection and predicting the response to treatment - an approach known as precision medicine. However this cannot be achieved without the knowledge of genetic variations that influence blood pressure (BP). A number of evolutionary factors including population demographics and forces of natural selection may be involved. This article explores some ideas on how natural selection influences BP regulation in ethnically and geographically diverse populations that could lead to them being susceptible to HTN. We explore how such evolutionary factors could impact the implementation of precision medicine in HTN. Finally, in order to ensure the success of precision medicine in HTN, we call for more initiatives to understand the genetic architecture within and between diverse populations with ancestry from different parts of the world, and to precisely classify the intermediate phenotypes of HTN.
Microsatellites are the most popular markers for parentage assignment and population genetic studies. To meet the demand for international comparability for genetic studies of Asian seabass, a standard panel of 28 microsatellites has been selected and characterized using the DNA of 24 individuals from Thailand, Malaysia, Indonesia and Australia. The average allele number of these markers was 10.82 +/- 0.71 (range: 6-19), and the expected heterozygosity averaged 0.76 +/- 0.02 (range: 0.63-1.00). All microsatellites showed Mendelian inheritance. In addition, eight standard size controls have been developed by cloning a set of microsatellite alleles into a pGEM-T vector to calibrate allele sizes determined by different laboratories, and are available upon request. Seven multiplex PCRs, each amplifying 3-5 markers, were optimized to accurately and rapidly genotype microsatellites. Parentage assignment using 10 microsatellites in two crosses (10 x 10 and 20 x 20) demonstrated a high power of these markers for revealing parent-sibling connections. This standard set of microsatellites will standardize genetic diversity studies of Asian seabass, and the multiplex PCR sets will facilitate parentage assignment.
Invasive apple snails, Pomacea canaliculata and P. maculata, have a widespread distribution globally and are regarded as devastating pests of agricultural wetlands. The two species are morphologically similar, which hinders species identification via morphological approaches and species-specific management efforts. Advances in molecular genetics may contribute effective diagnostic tools to potentially resolve morphological ambiguity. DNA barcoding has revolutionized the field of taxonomy by providing an alternative, simple approach for species discrimination, where short sections of DNA, the cytochrome c oxidase subunit I (COI) gene in particular, are used as 'barcodes' to delineate species boundaries. In our study, we aimed to assess the effectiveness of two mitochondrial markers, the COI and 16S ribosomal deoxyribonucleic acid (16S rDNA) markers for DNA barcoding of P. canaliculata and P. maculata. The COI and 16S rDNA sequences of 40 Pomacea specimens collected from six localities in Peninsular Malaysia were analyzed to assess their barcoding performance using phylogenetic methods and distance-based assessments. The results confirmed both markers were suitable for barcoding P. canaliculata and P. maculata. The phylogenies of the COI and 16S rDNA markers demonstrated species-specific monophyly and were largely congruent with the exception of one individual. The COI marker exhibited a larger barcoding gap (6.06-6.58%) than the 16S rDNA marker (1.54%); however, the magnitude of barcoding gap generated within the barcoding region of the 16S rDNA marker (12-fold) was bigger than the COI counterpart (approximately 9-fold). Both markers were generally successful in identifying P. canaliculata and P. maculata in the similarity-based DNA identifications. The COI + 16S rDNA concatenated dataset successfully recovered monophylies of P. canaliculata and P. maculata but concatenation did not improve individual datasets in distance-based analyses. Overall, although both markers were successful for the identification of apple snails, the COI molecular marker is a better barcoding marker and could be utilized in various population genetic studies of P. canaliculata and P. maculata.
Macaca fascicularis, also known as the cynomolgus macaque, is an important non-human primate animal model used in biomedical research. It is an Old-World primate widely distributed in Southeast Asia and is one of the most abundant macaque species in Malaysia. However, the genetic structure of wild cynomolgus macaque populations in Malaysia has not been thoroughly elucidated. In this study, we developed genic-simple sequence repeat (genic-SSR) markers from an in-house transcriptome dataset generated from the Malaysian cynomolgus macaque via RNA sequencing, and applied these markers on 26 cynomolgus macaque individuals. A collection of 14,751 genic-SSRs were identified, where 13,709 were perfect SSRs. Dinucleotide repeats were the most common repeat motifs with a frequency of 65.05%, followed by trinucleotide repeats (20.55%). Subsequently, we designed 300 pairs of primers based on perfect di- and trinucleotide SSRs, in which 105 SSRs were associated with functional genes. A subset of 30 SSR markers were randomly selected and validated, yielding 19 polymorphic markers with an average polymorphism information content value of 0.431. The development of genic-SSR markers in this study is indeed timely to provide useful markers for functional and population genetic studies of the cynomolgus macaque and other related non-human primate species.
BACKGROUND: Apolipoprotein E (apoE) is encoded by a polymorphic gene located on chromosome 19. The three common apoE alleles are epsilon2, epsilon3 and epsilon4. We studied the frequencies of the apoE alleles and genotypes in the three ethnic groups-Malay, Chinese and Indian-in Malaysia using DNA amplification followed by agarose gel electrophoresis.
METHODS: EDTA blood was collected and DNA was extracted using proteinase K-SDS digestion and purified by phenol-chloroform extraction. The apoE gene sequence was amplified using the PCR and apoE genotyping was performed by restriction enzyme digestion with HhaI.
RESULTS: Genotyping of the apoE gene produces six genotypes-E2/E2, E2/E3, E3/E3, E2/E4, E3/E4 and E4/E4. The most common apoE genotype in the Malays, Chinese and Indians studied was E3/E3, thus the most common apoE allele was epsilon3. The three common apoE genotypes were E3/E3 followed by E3/E4 and E2/E3, except in the Indians where E2/E3 was not detected. The three apoE alleles were confirmed in the Malays, Chinese and Indians except for the epsilon2 allele which was absent in the Indians.
CONCLUSION: The combined frequency of the apoE alleles in the Malays, Chinese and Indians was 0.058, 0.829 and 0.114 for epsilon2, epsilon3 and epsilon4, respectively.
MiniSTR loci has demonstrated to be an effective approach to recover genetic information from degraded sample, due to the improved PCR efficiency of their reduced PCR amplicon sizes. This study constructed a partial miniSGM panel and investigated the performance of four miniSTR loci, D2S1338, D16S539, D18S51 and FGA, in three ethnic populations residing in Singapore. The suitability of the miniSTR primers for Singapore populations was assessed for loci D16S539, D18S51 and FGA.
Laboratories intending to adopt cycle sequencing of PCR products in their routine analysis often face a confusing range of methods and kits. Through the study of mitochondrial cytochrome b, we have shown that clean and highly reproducible sequences could be obtained by using a combination of existing simple and economical methods in the preparation of DNA templates, PCR, purification of PCR products and sequencing. Our protocol is useful in itself or as a standard in typing other PCR-amplified DNA at the population level.
The Malay people are an important ethnic composition in Southeast Asia, but their genetic make-up and population structure remain poorly studied. Here we conducted a genome-wide study of four geographical Malay populations: Peninsular Malaysian Malay (PMM), Singaporean Malay (SGM), Indonesian Malay (IDM) and Sri Lankan Malay (SLM). All the four Malay populations showed substantial admixture with multiple ancestries. We identified four major ancestral components in Malay populations: Austronesian (17%-62%), Proto-Malay (15%-31%), East Asian (4%-16%) and South Asian (3%-34%). Approximately 34% of the genetic makeup of SLM is of South Asian ancestry, resulting in its distinct genetic pattern compared with the other three Malay populations. Besides, substantial differentiation was observed between the Malay populations from the north and the south, and between those from the west and the east. In summary, this study revealed that the genetic identity of the Malays comprises a mixed entity of multiple ancestries represented by Austronesian, Proto-Malay, East Asian and South Asian, with most of the admixture events estimated to have occurred 175 to 1,500 years ago, which in turn suggests that geographical isolation and independent admixture have significantly shaped the genetic architectures and the diversity of the Malay populations.
Tibetan high-altitude adaptation (HAA) has been studied extensively, and many candidate genes have been reported. Subsequent efforts targeting HAA functional variants, however, have not been that successful (e.g., no functional variant has been suggested for the top candidate HAA gene, EPAS1). With WinXPCNVer, a method developed in this study, we detected in microarray data a Tibetan-enriched deletion (TED) carried by 90% of Tibetans; 50% were homozygous for the deletion, whereas only 3% carried the TED and 0% carried the homozygous deletion in 2,792 worldwide samples (p < 10(-15)). We employed long PCR and Sanger sequencing technologies to determine the exact copy number and breakpoints of the TED in 70 additional Tibetan and 182 diverse samples. The TED had identical boundaries (chr2: 46,694,276-46,697,683; hg19) and was 80 kb downstream of EPAS1. Notably, the TED was in strong linkage disequilibrium (LD; r(2) = 0.8) with EPAS1 variants associated with reduced blood concentrations of hemoglobin. It was also in complete LD with the 5-SNP motif, which was suspected to be introgressed from Denisovans, but the deletion itself was absent from the Denisovan sequence. Correspondingly, we detected that footprints of positive selection for the TED occurred 12,803 (95% confidence interval = 12,075-14,725) years ago. We further whole-genome deep sequenced (>60×) seven Tibetans and verified the TED but failed to identify any other copy-number variations with comparable patterns, giving this TED top priority for further study. We speculate that the specific patterns of the TED resulted from its own functionality in HAA of Tibetans or LD with a functional variant of EPAS1.