MtDNA control region sequence polymorphisms and phylogenetic analysis of Malay population living in or around Kuala Lumpur in Malaysia

Autosomal STR Profiling and Databanking in Malaysia: Current Status and Future Prospects

Science and technology are extensively used in criminal investigation. From the mid- to late-1980s, one of the scientific discoveries that has had a particularly remarkable impact on this field has been the use of highly variable DNA sequence regions (minisatellites) in the human genome for individual identification. The technique was initially referred to as DNA fingerprinting, but is now more widely referred to as DNA profiling. Since then, many new developments have occurred within this area of science. These include the introduction of new genetic markers (microsatellites also known as short tandem repeats/STRs), the use of the polymerase chain reaction for target amplification, the development of DNA databases (databanking), and the advancement and/or improvement of genotyping protocols and technologies. In 2019, we described the progress of DNA profiling and DNA databanking in Malaysia for the first time. This report included information on DNA analysis regulations and legislation, STR genotyping protocols, database management, and accreditation status. Here, we provide an update on the performance of our DNA databank (numbers of DNA profiles and hits) plus the technical issues associated with correctly assigning the weight of evidence for DNA profiles in an ethnically diverse population, and the potential application of rapid DNA testing in the country. A total of 116,534 DNA profiles were obtained and stored in the Forensic DNA Databank of Malaysia (FDDM) by 2019, having increased from 70,570 in 2017. The number of hits increased by more than three-fold in just two years, where 17 and 69 hits between the DNA profiles stored in the FDDM and those from crime scenes, suspects, detainees, drug users, convicts, missing persons, or volunteers were recorded in 2017 and 2019, respectively. Forensic DNA analysis and databanking are thus progressing well in Malaysia and have already contributed to many criminal investigations. However, several other issues are discussed here, including the need for STR population data for uncharacterized population groups, and pilot trials for adopting rapid DNA profiling technology. These aspects should be considered by policy makers and law enforcement agencies in order to increase the reliability and efficiency of DNA profiling in criminal cases and in kinship analysis in Malaysia.

A genetic chronology for the Indian Subcontinent points to heavily sex-biased dispersals

BACKGROUND

India is a patchwork of tribal and non-tribal populations that speak many different languages from various language families. Indo-European, spoken across northern and central India, and also in Pakistan and Bangladesh, has been frequently connected to the so-called "Indo-Aryan invasions" from Central Asia ~3.5 ka and the establishment of the caste system, but the extent of immigration at this time remains extremely controversial. South India, on the other hand, is dominated by Dravidian languages. India displays a high level of endogamy due to its strict social boundaries, and high genetic drift as a result of long-term isolation which, together with a very complex history, makes the genetic study of Indian populations challenging.

RESULTS

We have combined a detailed, high-resolution mitogenome analysis with summaries of autosomal data and Y-chromosome lineages to establish a settlement chronology for the Indian Subcontinent. Maternal lineages document the earliest settlement ~55-65 ka (thousand years ago), and major population shifts in the later Pleistocene that explain previous dating discrepancies and neutrality violation. Whilst current genome-wide analyses conflate all dispersals from Southwest and Central Asia, we were able to tease out from the mitogenome data distinct dispersal episodes dating from between the Last Glacial Maximum to the Bronze Age. Moreover, we found an extremely marked sex bias by comparing the different genetic systems.

CONCLUSIONS

Maternal lineages primarily reflect earlier, pre-Holocene processes, and paternal lineages predominantly episodes within the last 10 ka. In particular, genetic influx from Central Asia in the Bronze Age was strongly male-driven, consistent with the patriarchal, patrilocal and patrilineal social structure attributed to the inferred pastoralist early Indo-European society. This was part of a much wider process of Indo-European expansion, with an ultimate source in the Pontic-Caspian region, which carried closely related Y-chromosome lineages, a smaller fraction of autosomal genome-wide variation and an even smaller fraction of mitogenomes across a vast swathe of Eurasia between 5 and 3.5 ka.

Carriers of human mitochondrial DNA macrohaplogroup M colonized India from southeastern Asia

Background

From a mtDNA dominant perspective, the exit from Africa of modern humans to colonize Eurasia occurred once, around 60 kya, following a southern coastal route across Arabia and India to reach Australia short after. These pioneers carried with them the currently dominant Eurasian lineages M and N. Based also on mtDNA phylogenetic and phylogeographic grounds, some authors have proposed the coeval existence of a northern route across the Levant that brought mtDNA macrohaplogroup N to Australia. To contrast both hypothesis, here we reanalyzed the phylogeography and respective ages of mtDNA haplogroups belonging to macrohaplogroup M in different regions of Eurasia and Australasia.

Results

The macrohaplogroup M has a historical implantation in West Eurasia, including the Arabian Peninsula. Founder ages of M lineages in India are significantly younger than those in East Asia, Southeast Asia and Near Oceania. Moreover, there is a significant positive correlation between the age of the M haplogroups and its longitudinal geographical distribution. These results point to a colonization of the Indian subcontinent by modern humans carrying M lineages from the east instead the west side.

Conclusions

The existence of a northern route, previously proposed for the mtDNA macrohaplogroup N, is confirmed here for the macrohaplogroup M. Both mtDNA macrolineages seem to have differentiated in South East Asia from ancestral L3 lineages. Taking this genetic evidence and those reported by other disciplines we have constructed a new and more conciliatory model to explain the history of modern humans out of Africa.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0816-8) contains supplementary material, which is available to authorized users.

Multiplex APLP System for High-Resolution Haplogrouping of Extremely Degraded East-Asian Mitochondrial DNAs

Mitochondrial DNA (mtDNA) serves as a powerful tool for exploring matrilineal phylogeographic ancestry, as well as for analyzing highly degraded samples, because of its polymorphic nature and high copy numbers per cell. The recent advent of complete mitochondrial genome sequencing has led to improved techniques for phylogenetic analyses based on mtDNA, and many multiplex genotyping methods have been developed for the hierarchical analysis of phylogenetically important mutations. However, few high-resolution multiplex genotyping systems for analyzing East-Asian mtDNA can be applied to extremely degraded samples. Here, we present a multiplex system for analyzing mitochondrial single nucleotide polymorphisms (mtSNPs), which relies on a novel amplified product-length polymorphisms (APLP) method that uses inosine-flapped primers and is specifically designed for the detailed haplogrouping of extremely degraded East-Asian mtDNAs. We used fourteen 6-plex polymerase chain reactions (PCRs) and subsequent electrophoresis to examine 81 haplogroup-defining SNPs and 3 insertion/deletion sites, and we were able to securely assign the studied mtDNAs to relevant haplogroups. Our system requires only 1×10-13 g (100 fg) of crude DNA to obtain a full profile. Owing to its small amplicon size (<110 bp), this new APLP system was successfully applied to extremely degraded samples for which direct sequencing of hypervariable segments using mini-primer sets was unsuccessful, and proved to be more robust than conventional APLP analysis. Thus, our new APLP system is effective for retrieving reliable data from extremely degraded East-Asian mtDNAs.

Genetic analysis of mitochondrial DNA control region variations in four tribes of Khyber Pakhtunkhwa, Pakistan

Due to its geo strategic position at the crossroad of Asia, Pakistan has gained crucial importance of playing its pivotal role in subsequent human migratory events, both prehistoric and historic. This human movement became possible through an ancient overland network of trails called "The Silk Route" linking Asia Minor, Middle East China, Central Asia and Southeast Asia. This study was conducted to analyze complete mitochondrial control region samples of 100 individuals of four major Pashtun tribes namely, Bangash, Khattak, Mahsuds and Orakzai in the province of Khyber Pakhtunkhwa, Pakistan. All Pashtun tribes revealed high genetic diversity which is comparable to the other Central Asian, Southeast Asian and European populations. The configuration of genetic variation and heterogeneity further unveiled through Multidimensional Scaling, Principal Component Analysis and phylogenetic analysis. The results revealed that Pashtun are the composite mosaic of West Eurasian ancestry of numerous geographic origin. They received substantial gene flow during different invasive movements and have a high element of the Western provenance. The most common haplogroups reported in this study are: South Asian haplogroups M (28%) and R (8%); whereas, West Asians haplogroups are present, albeit in high frequencies (67%) and widespread over all; HV (15%), U (17%), H (9%), J (8%), K (8%), W (4%), N (3%) and T (3%). Moreover, we linked the unexplored genetic connection between Ashkenazi Jews and Pashtun. The presence of specific haplotypes J1b (4%) and K1a1b1a (5%) pointed to a genetic connection of Jewish conglomeration in Khattak tribe. This was a result of an ancient genetic influx in the early Neolithic period that led to the formation of a diverse genetic substratum in present day Pashtun.

Northward genetic penetration across the Himalayas viewed from Sherpa people

The Himalayas have been suggested as a natural barrier for human migrations, especially the northward dispersals from the Indian Subcontinent to Tibetan Plateau. However, although the majority of Sherpa have a Tibeto-Burman origin, considerable genetic components from Indian Subcontinent have been observed in Sherpa people living in Tibet. The western Y chromosomal haplogroups R1a1a-M17, J-M304, and F*-M89 comprise almost 17% of Sherpa paternal gene pool. In the maternal side, M5c2, M21d, and U from the west also count up to 8% of Sherpa people. Those lineages with South Asian origin indicate that the Himalayas have been permeable to bidirectional gene flow.

Population data of mitochondrial DNA HVS-I and HVS-II sequences for 208 Henan Han Chinese

The two hypervariable segments (HVS-I and HVS-II) of mitochondrial DNA (mtDNA) control region were sequenced for a population of 208 unrelated healthy individuals sampled from Suiping County, Henan Province, China. A total of 192 different haplotypes were identified, of which 179 haplotypes were unique (93.23%). The variation of the mtDNA HVS-I and HVS-II was confined to 166 nucleotide positions, of which 115 were observed in the HVS-I and 51 in the HVS-II. The haplotype diversity and random match probability were 0.9991 and 0.0061, respectively. Following the principle of the updated East Asian mtDNA phylogeny tree, individual samples were assigned to the specific haplogroups based on the information both from control region and coding-region obtained. Haplogroup D was the most common haplogroup (25.96%). The northern China-prevalent haplogroups (A, C, D, G, M8, Y, and Z) and the southern China-prevalent haplogroups (B, F, M7, N9, and R9) accounted for 48.56% and 46.63%, respectively, of the Henan Han mtDNA gene pool. The mtDNA hypervariable region was highly polymorphic in Henan Han population. These sequences could serve as mtDNA reference data for forensic casework in Henan population as well as for population genetic study.

Genetic structure of Qiangic populations residing in the western Sichuan corridor

The Qiangic languages in western Sichuan (WSC) are believed to be the oldest branch of the Sino-Tibetan linguistic family, and therefore, all Sino-Tibetan populations might have originated in WSC. However, very few genetic investigations have been done on Qiangic populations and no genetic evidences for the origin of Sino-Tibetan populations have been provided. By using the informative Y chromosome and mitochondrial DNA (mtDNA) markers, we analyzed the genetic structure of Qiangic populations. Our results revealed a predominantly Northern Asian-specific component in Qiangic populations, especially in maternal lineages. The Qiangic populations are an admixture of the northward migrations of East Asian initial settlers with Y chromosome haplogroup D (D1-M15 and the later originated D3a-P47) in the late Paleolithic age, and the southward Di-Qiang people with dominant haplogroup O3a2c1*-M134 and O3a2c1a-M117 in the Neolithic Age.

Analysis of human mitochondrial DNA polymorphisms in the Japanese population

The highly polymorphic nature and high amplification efficiency of mitochondrial DNA (mtDNA) is valuable for the analysis of biological evidence in forensic casework, such as the identification of individuals and assignment of race/ethnicity. To be useful, a mtDNA polymorphism database for the Japanese population requires an understanding of the range of haplotype variation and phylogenies of mtDNA sequences. To extend current knowledge on the haplotypes in the Japanese population, this study defines new lineages and provides more detail about some of those previously described. We compared the hypervariable regions (HVRs) of 270 healthy, unrelated Japanese individuals and demonstrated 192 haplotypes. Combining HVR1 and HVR2, the genetic diversity was 0.9935, thus providing a high level of identification capability. Haplogroup status was defined for 160 individuals using HVR1, HVR2, and particular coding region polymorphisms; these individuals belonged to 94 haplotypes, four of which were new lineages. The complete mtDNA sequence was also determined from seven individuals.

Large-scale mtDNA screening reveals a surprising matrilineal complexity in east Asia and its implications to the peopling of the region

In order to achieve a thorough coverage of the basal lineages in the Chinese matrilineal pool, we have sequenced the mitochondrial DNA (mtDNA) control region and partial coding region segments of 6,093 mtDNAs sampled from 84 populations across China. By comparing with the available complete mtDNA sequences, 194 of those mtDNAs could not be firmly assigned into the available haplogroups. Completely sequencing 51 representatives selected from these unclassified mtDNAs identified a number of novel lineages, including five novel basal haplogroups that directly emanate from the Eurasian founder nodes (M and N). No matrilineal contribution from the archaic hominid was observed. Subsequent analyses suggested that these newly identified basal lineages likely represent the genetic relics of modern humans initially peopling East Asia instead of being the results of gene flow from the neighboring regions. The observation that most of the newly recognized mtDNA lineages have already differentiated and show the highest genetic diversity in southern China provided additional evidence in support of the Southern Route peopling hypothesis of East Asians. Specifically, the enrichment of most of the basal lineages in southern China and their rather ancient ages in Late Pleistocene further suggested that this region was likely the genetic reservoir of modern humans after they entered East Asia.

Tracing the Austronesian footprint in Mainland Southeast Asia: a perspective from mitochondrial DNA

As the relic of the ancient Champa Kingdom, the Cham people represent the major Austronesian speakers in Mainland Southeast Asia (MSEA) and their origin is evidently associated with the Austronesian diffusion in MSEA. Hitherto, hypotheses stemming mainly from linguistic and cultural viewpoints on the origin of the Cham people remain a welter of controversies. Among the points of dissension is the muddled issue of whether the Cham people arose from demic or cultural diffusion from the Austronesians. Addressing this issue also helps elucidate the dispersal mode of the Austronesian language. In the present study, we have analyzed mitochondrial DNA (mtDNA) control-region and coding-region sequence variations in 168 Cham and 139 Kinh individuals from Vietnam. Around 77% and 95% matrilineal components in the Chams and the Kinhs, respectively, could be assigned into the defined mtDNA haplogroups. Additionally, three common East Eurasian haplogroups B, R9, and M7 account for the majority (>60%) of maternal components in both populations. Entire sequencing of 20 representative mtDNAs selected from the thus far unclassified lineages, together with four new mtDNA genome sequences from Thailand, led to the identification of one new haplogroup M77 and helped to re-evaluate several haplogroups determined previously. Comparing the Chams with other Southeast Asian populations reveals that the Chams had a closer affinity with the Mon-Khmer populations in MSEA than with the Austronesian populations from Island Southeast Asia (ISEA). Further analyses failed to detect the potential homelands of the Chams in ISEA. Therefore, our results suggested that the origin of the Cham was likely a process of assimilation of massive local Mon-Khmer populations accompanied with language shift, thus indicating that the Austronesian diffusion in MSEA was mainly mediated by cultural diffusion, at least from the matrilineal genetic perspective, an observation in agreement with the hypothesis of the Nusantao Maritime Trading and Communication Networks.

Sequence polymorphisms of mtDNA HV1, HV2, and HV3 regions in the Malay population of Peninsular Malaysia

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.

Mitochondrial DNA diversity in a population from Santa Catarina (Brazil): predominance of the European input

The state of Santa Catarina (Brazil) is known to have represented a cultural crossroads in South America due to several historic migrations mainly from Europe and Africa. We set out to scrutinize whether the genetic imprint of these migrations could be traced through analysis of the matrilineal gene pool of the Catarinenses. The entire control region of the mitochondrial DNA was studied in 80 healthy and maternally unrelated individuals. The analysis of haplogroup distribution revealed that this population is extremely heterogeneous, showing the coexistence of matrilineal lineages with three different phylogeographic origins. European lineages are the most frequent due mainly to the impact of relatively recent migratory waves from Europe. In spite of this, Native American lineages and African lineages incorporated with the slave trade are also present in noticeable proportions. The strikingly high variability generated by intense gene flow is mirrored in a high sequence diversity (0.9930) and power of discrimination (0.9806). Thus, analysis of the entire mitochondrial DNA control region emerges as a valuable tool for forensic genetic purposes in this highly admixed population, an attribute common to several present-day Latin American populations.