Forensic and phylogeographic characterization of mtDNA lineages from northern Thailand (Chiang Mai)

Analysis of maternal genetic structure of mitochondrial DNA control region from Tai-Kadai-speaking Buyei population in southwestern China

BACKGROUND

Even though the Buyei are a recognised ethnic group in southwestern China, there hasn't been much work done on forensic population genetics, notably using mitochondrial DNA. The sequences and haplogroups of mitochondrial DNA control regions of the Buyei peoples were studied to provide support for the establishment of a reference database for forensic DNA analysis in East Asia.

METHODS AND RESULTS

The mitochondrial DNA control region sequences of 200 Buyei individuals in Guizhou were investigated. The haplotype frequencies and haplogroup distribution of the Buyei nationality in Guizhou were calculated. At the same time, the paired Fst values of the study population and other populations around the world were computed, to explore their genetic polymorphism and population relationship. A total of 179 haplotypes were detected in the Buyei population, with frequencies of 0.005-0.015. All haplotypes were assigned to 89 different haplogroups. The haplotype diversity and random matching probability were 0.999283 and 0.0063, respectively. The paired Fst genetic distances and correlation p-values among the 54 populations revealed that the Guizhou Buyei was most closely related to the Henan Han and the Guizhou Miao, and closer to the Hazara population in Pakistan and the Chiang Mai population.

CONCLUSIONS

The study of mitochondrial DNA based on the maternal genetic structure of the Buyei nationality in Guizhou will benefit the establishment of an East Asian forensic DNA reference database and provide a reference for anthropological research in the future.

The mitochondrial DNA HVI and HVII sequences and haplogroup distribution in a population sample from Vietnam

BACKGROUND

Mitochondrial DNA (mtDNA) analysis has been used in forensics and requires well-established population databases for statistical interpretations. However, high-quality mtDNA data from Vietnamese population samples have been limited.

AIM

To examine the mtDNA sequences and haplogroup compositions of a Vietnamese population to provide an mtDNA dataset that can further be used to construct a Vietnamese-specific reference database.

SUBJECTS AND METHODS

A total of 173 Vietnamese individuals were analysed for two hypervariable regions (HVI and HVII) of mtDNA. Forensic parameters were calculated and haplogroup assignment was performed based on the resulting mtDNA haplotypes. Genetic relationships between the Vietnamese and other Asian populations were investigated through principal component analysis (PCA) and pairwise F_st.

RESULTS

The Vietnamese population sample consisted of 145 different haplotypes with a random match probability of 0.96%, a power of discrimination of 0.9904, and a haplotype diversity of 0.9962. The samples were assigned to 83 haplogroups that were commonly reported in Asia. PCA and pairwise F_st revealed close relationships of the Vietnamese population with other Asian populations, especially with populations in proximity.

CONCLUSION

The results from this study can contribute to the current genetic information content as a supplementary mtDNA reference dataset for forensic investigations and phylogenetic research.

Current pool of ultimate collection of mitochondrial DNA from remnants of Kalash

The mitochondrial DNA (mtDNA) complete control region coverage of 111 individuals from Kalash population of Pakistan has been presented for forensic applications and to infer their genetic parameters. We detected in total 14 different haplotypes with only five unique and nine shared by more than one individual. This population has come up with quite lower haplotype diversity (0.8393) and very higher random match probability (0.1682), and ultimately lower power of discrimination (0.832). Additionally, haplogroup distribution reveals the genetic ancestry of Kalash, mainly from West Eurasia (98.8%) and very little from South Asia (0.9%). Neither African lineages nor East Asian genetic segments were detected among these Kalash. This study will contribute to the database development for forensic applications as well as to track the evolutionary highlights of this ethnic group.

An in-depth analysis of the mitochondrial phylogenetic landscape of Cambodia

Cambodia harbours a variety of human aboriginal populations that have scarcely been studied in terms of genetic diversity of entire mitochondrial genomes. Here we present the matrilineal gene pool of 299 Cambodian refugees from three different ethnic groups (Cham, Khmer, and Khmer Loeu) deriving from 16 Cambodian districts. After establishing a DNA-saving high-throughput strategy for mitochondrial whole-genome Sanger sequencing, a HaploGrep based workflow was used for quality control, haplogroup classification and phylogenetic reconstruction. The application of diverse phylogenetic algorithms revealed an exciting picture of the genetic diversity of Cambodia, especially in relation to populations from Southeast Asia and from the whole world. A total of 224 unique haplotypes were identified, which were mostly classified under haplogroups B5a1, F1a1, or categorized as newly defined basal haplogroups or basal sub-branches of R, N and M clades. The presence of autochthonous maternal lineages could be confirmed as reported in previous studies. The exceptional homogeneity observed between and within the three investigated Cambodian ethnic groups indicates genetic isolation of the whole population. Between ethnicities, genetic barriers were not detected. The mtDNA data presented here increases the phylogenetic resolution in Cambodia significantly, thereby highlighting the need for an update of the current human mtDNA phylogeny.

Mitochondrial DNA control region variation in a population sample from Thailand

Mitochondrial DNA (mtDNA) control region sequences from hair samples of 213 individuals from Thailand were analyzed using Sanger sequencing. A total of 170 different haplotypes were identified, of which 146 occurred only once (unique haplotypes). The dataset showed a random match probability of 0.87% and a haplotype diversity of 0.9960. The samples were assigned to 85 different haplogroups with B5a, F1a1a, and M being the most frequent ones. Pairwise F_ST-values between this and other Southeast and East Asian populations revealed significant but relatively low differences, indicating a close relation. Heteroplasmic positions were observed in 12.2% of hair samples confirming the frequent appearance of heteroplasmic positions in hairs. This dataset will complement existing data as an mtDNA reference for forensic investigations.

Phylogeographic and genome-wide investigations of Vietnam ethnic groups reveal signatures of complex historical demographic movements

The territory of present-day Vietnam was the cradle of one of the world’s earliest civilizations, and one of the first world regions to develop agriculture. We analyzed the mitochondrial DNA (mtDNA) complete control region of six ethnic groups and the mitogenomes from Vietnamese in The 1000 Genomes Project (1000G). Genome-wide data from 1000G (~55k SNPs) were also investigated to explore different demographic scenarios. All Vietnamese carry South East Asian (SEA) haplotypes, which show a moderate geographic and ethnic stratification, with the Mong constituting the most distinctive group. Two new mtDNA clades (M7b1a1f1 and F1f1) point to historical gene flow between the Vietnamese and other neighboring countries. Bayesian-based inferences indicate a time-deep and continuous population growth of Vietnamese, although with some exceptions. The dramatic population decrease experienced by the Cham 700 years ago (ya) fits well with the Nam tiến (“southern expansion”) southwards from their original heartland in the Red River Delta. Autosomal SNPs consistently point to important historical gene flow within mainland SEA, and add support to a main admixture event occurring between Chinese and a southern Asian ancestral composite (mainly represented by the Malay). This admixture event occurred ~800 ya, again coinciding with the Nam tiến.

Understanding influences of culture and history on mtDNA variation and population structure in three populations from Assam, Northeast India

OBJECTIVES

Positioned at the nexus of India, China, and Southeast Asia, Northeast India is presumed to have served as a channel for land-based human migration since the Upper Pleistocene. Assam is the largest state in the Northeast. We characterized the genetic background of three populations and examined the ways in which their population histories and cultural practices have influenced levels of intrasample and intersample variation.

METHODS

We examined sequence data from the mtDNA hypervariable control region and selected diagnostic mutations from the coding region in 128 individuals from three ethnic groups currently living in Assam: two Scheduled tribes (Sonowal Kachari and Rabha), and the non-Scheduled Tai Ahom.

RESULTS

The populations of Assam sampled here express mtDNA lineages indicative of South Asian, Southeast Asian, and East Asian ancestry. We discovered two completely novel haplogroups in Assam that accounted for 6.2% of the lineages in our sample. We also identified a new subhaplogroup of M9a that is prevalent in the Sonowal Kachari of Assam (19.1%), but not present in neighboring Arunachal Pradesh, indicating substantial regional population structuring. Employing a large comparative dataset into a series of multidimensional scaling (MDS) analyses, we saw the Rabha cluster with populations sampled from Yunnan Province, indicating that the historical matrilineality of the Rabha has maintained lineages from Southern China.

CONCLUSION

Assam has undergone multiple colonization events in the time since the initial peopling event, with populations from Southern China and Southeast Asia having the greatest influence on maternal lineages in the region.

EMPOP-quality mtDNA control region sequences from Kashmiri of Azad Jammu & Kashmir, Pakistan

The mitochondrial DNA (mtDNA) control region (nucleotide position 16024-576) sequences were generated through Sanger sequencing method for 317 self-identified Kashmiris from all districts of Azad Jammu & Kashmir Pakistan. The population sample set showed a total of 251 haplotypes, with a relatively high haplotype diversity (0.9977) and a low random match probability (0.54%). The containing matrilineal lineages belonging to three different phylogeographic origins of Western Eurasian (48.9%), South Asian (47.0%) and East Asian (4.1%). The present study was compared to previous data from Pakistan and other worldwide populations (Central Asia, Western Asia, and East & Southeast Asia). The dataset is made available through EMPOP under accession number EMP00679 and will serve as an mtDNA reference database in forensic casework in Pakistan.

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.

Single nucleotide polymorphisms of mitochondrial DNA HVS-I and HVS-II in Chinese Bai ethnic group

For forensic and population genetic purposes, a total of 125 unrelated volunteers' blood samples were collected from Chinese Bai ethnic minority group to analyze sequence variation of two hypervariable segments (HVS-I and HVS-II) in the mitochondrial DNA control region. Comparing the HVS-I and HVS-II sequences of the 125 Chinese Bais to the Anderson reference sequence, we found 86 polymorphic loci in HVS-I and 40 in HVS-II in mitochondrial DNA sequences of the Chinese Bai ethnic minority group, which defined 93 and 53 different haplotypes, respectively. Haplotype diversity and the mean pairwise differences were 0.992 ± 0.003 and 6.553 in HVS-I, and 0.877 ± 0.027 and 2.407 in HVS-II, respectively. We defined four macrohaplogroups R, M, N and D with the proportions ranging from 9.6% to 40.0%. With the analysis of the hypervariable domain from nucleotide 16 180-16 193 in HVS-I, our study revealed new haplotypes of sequence variations. In addition, the Fst metric, phylogenetic tree, and principal component analysis demonstrated a close genetic relationship between the Bai group and Chinese Han populations from South China, Changsha, and Guangdong. The results support that the Bai group is a multiorigin ethnic minority that has merged with the Chinese Han population.

Human settlement history between Sunda and Sahul: a focus on East Timor (Timor-Leste) and the Pleistocenic mtDNA diversity

BACKGROUND

Distinct, partly competing, "waves" have been proposed to explain human migration in(to) today's Island Southeast Asia and Australia based on genetic (and other) evidence. The paucity of high quality and high resolution data has impeded insights so far. In this study, one of the first in a forensic environment, we used the Ion Torrent Personal Genome Machine (PGM) for generating complete mitogenome sequences via stand-alone massively parallel sequencing and describe a standard data validation practice.

RESULTS

In this first representative investigation on the mitochondrial DNA (mtDNA) variation of East Timor (Timor-Leste) population including >300 individuals, we put special emphasis on the reconstruction of the initial settlement, in particular on the previously poorly resolved haplogroup P1, an indigenous lineage of the Southwest Pacific region. Our results suggest a colonization of southern Sahul (Australia) >37 kya, limited subsequent exchange, and a parallel incubation of initial settlers in northern Sahul (New Guinea) followed by westward migrations <28 kya.

CONCLUSIONS

The temporal proximity and possible coincidence of these latter dispersals, which encompassed autochthonous haplogroups, with the postulated "later" events of (South) East Asian origin pinpoints a highly dynamic migratory phase.

Helena, the hidden beauty: Resolving the most common West Eurasian mtDNA control region haplotype by massively parallel sequencing an Italian population sample

The analysis of mitochondrial (mt)DNA is a powerful tool in forensic genetics when nuclear markers fail to give results or maternal relatedness is investigated. The mtDNA control region (CR) contains highly condensed variation and is therefore routinely typed. Some samples exhibit an identical haplotype in this restricted range. Thus, they convey only weak evidence in forensic queries and limited phylogenetic information. However, a CR match does not imply that also the mtDNA coding regions are identical or samples belong to the same phylogenetic lineage. This is especially the case for the most frequent West Eurasian CR haplotype 263G 315.1C 16519C, which is observed in various clades within haplogroup H and occurs at a frequency of 3-4% in many European populations. In this study, we investigated the power of massively parallel complete mtGenome sequencing in 29 Italian samples displaying the most common West Eurasian CR haplotype - and found an unexpected high diversity. Twenty-eight different haplotypes falling into 19 described sub-clades of haplogroup H were revealed in the samples with identical CR sequences. This study demonstrates the benefit of complete mtGenome sequencing for forensic applications to enforce maximum discrimination, more comprehensive heteroplasmy detection, as well as highest phylogenetic resolution.

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.

Epidermal nerve fiber density, oxidative stress, and mitochondrial haplogroups in HIV-infected Thais initiating therapy

OBJECTIVE

We explored associations between mitochondrial DNA (mtDNA) haplogroups, epidermal nerve fiber density (ENFD), and HIV-associated sensory neuropathy (HIV-SN) in a randomized trial of Thai patients initiating antiretroviral therapy (ART).

DESIGN

The South East Asia Research Collaboration with Hawaii 003 study evaluated toxicity of nucleoside reverse transcriptase inhibitors (stavudine vs. zidovudine vs. tenofovir). We present secondary analyses of mtDNA haplogroups and ENFD changes.

METHODS

ENFD, peripheral blood mononuclear cell mitochondrial complex I and IV, and 8-oxo-deoxyguanine (8-oxo-dG) were quantified. Peripheral blood mononuclear cell mtDNA sequences were obtained for haplogroup determination. Multivariate regression of ENFD change was performed.

RESULTS

Paired ENFD was available from 118 patients. Median age, CD4 cell count, and height at entry were 34 years, 172 cells/μl, and 162 cm, respectively. Major haplogroups included M (42%), F (21%), and B (16%). Baseline ENFD, CD4 cell count, randomized ART, and biomarkers did not differ by haplogroup. Haplogroup B patients were older (P=0.02) at baseline, and had an increase in median ENFD (+1.5 vs. -2.9 fibers/mm; P=0.03) and 8-oxo-dG break frequency (+0.05 vs. 0.00; P=0.05) compared to other haplogroups. In a multivariate model, haplogroup B was associated with increased ENFD (β=3.5, P=0.009) at week 24, whereas older age (P=0.02), higher baseline CD4 cell count, (P=0.03), higher complex I level (P=0.03), and higher ENFD (P<0.001) at baseline were all associated with decreased ENFD. Three of the six HIV-SN cases were haplogroup B (P=0.05).

CONCLUSIONS

Thai persons belonging to mtDNA haplogroup B had increased ENFD and 8-oxo-dG on ART, and were more likely to develop HIV-SN. These results suggest that mtDNA variation influences early oxidative damage and ENFD changes.

Large-scale mitochondrial DNA analysis in Southeast Asia reveals evolutionary effects of cultural isolation in the multi-ethnic population of Myanmar

Background

Myanmar is the largest country in mainland Southeast Asia with a population of 55 million people subdivided into more than 100 ethnic groups. Ruled by changing kingdoms and dynasties and lying on the trade route between India and China, Myanmar was influenced by numerous cultures. Since its independence from British occupation, tensions between the ruling Bamar and ethnic minorities increased.

Results

Our aim was to search for genetic footprints of Myanmar’s geographic, historic and sociocultural characteristics and to contribute to the picture of human colonization by describing and dating of new mitochondrial DNA (mtDNA) haplogroups. Therefore, we sequenced the mtDNA control region of 327 unrelated donors and the complete mitochondrial genome of 44 selected individuals according to highest quality standards.

Conclusion

Phylogenetic analyses of the entire mtDNA genomes uncovered eight new haplogroups and three unclassified basal M-lineages. The multi-ethnic population and the complex history of Myanmar were reflected in its mtDNA heterogeneity. Population genetic analyses of Burmese control region sequences combined with population data from neighboring countries revealed that the Myanmar haplogroup distribution showed a typical Southeast Asian pattern, but also Northeast Asian and Indian influences. The population structure of the extraordinarily diverse Bamar differed from that of the Karen people who displayed signs of genetic isolation. Migration analyses indicated a considerable genetic exchange with an overall positive migration balance from Myanmar to neighboring countries. Age estimates of the newly described haplogroups point to the existence of evolutionary windows where climatic and cultural changes gave rise to mitochondrial haplogroup diversification in Asia.

Insight into the peopling of Mainland Southeast Asia from Thai population genetic structure

There is considerable ethno-linguistic and genetic variation among human populations in Asia, although tracing the origins of this diversity is complicated by migration events. Thailand is at the center of Mainland Southeast Asia (MSEA), a region within Asia that has not been extensively studied. Genetic substructure may exist in the Thai population, since waves of migration from southern China throughout its recent history may have contributed to substantial gene flow. Autosomal SNP data were collated for 438,503 markers from 992 Thai individuals. Using the available self-reported regional origin, four Thai subpopulations genetically distinct from each other and from other Asian populations were resolved by Neighbor-Joining analysis using a 41,569 marker subset. Using an independent Principal Components-based unsupervised clustering approach, four major MSEA subpopulations were resolved in which regional bias was apparent. A major ancestry component was common to these MSEA subpopulations and distinguishes them from other Asian subpopulations. On the other hand, these MSEA subpopulations were admixed with other ancestries, in particular one shared with Chinese. Subpopulation clustering using only Thai individuals and the complete marker set resolved four subpopulations, which are distributed differently across Thailand. A Sino-Thai subpopulation was concentrated in the Central region of Thailand, although this constituted a minority in an otherwise diverse region. Among the most highly differentiated markers which distinguish the Thai subpopulations, several map to regions known to affect phenotypic traits such as skin pigmentation and susceptibility to common diseases. The subpopulation patterns elucidated have important implications for evolutionary and medical genetics. The subpopulation structure within Thailand may reflect the contributions of different migrants throughout the history of MSEA. The information will also be important for genetic association studies to account for population-structure confounding effects.

Haplogrouping mitochondrial DNA sequences in Legal Medicine/Forensic Genetics

Haplogrouping refers to the classification of (partial) mitochondrial DNA (mtDNA) sequences into haplogroups using the current knowledge of the worldwide mtDNA phylogeny. Haplogroup assignment of mtDNA control-region sequences assists in the focused comparison with closely related complete mtDNA sequences and thus serves two main goals in forensic genetics: first is the a posteriori quality analysis of sequencing results and second is the prediction of relevant coding-region sites for confirmation or further refinement of haplogroup status. The latter may be important in forensic casework where discrimination power needs to be as high as possible. However, most articles published in forensic genetics perform haplogrouping only in a rudimentary or incorrect way. The present study features PhyloTree as the key tool for assigning control-region sequences to haplogroups and elaborates on additional Web-based searches for finding near-matches with complete mtDNA genomes in the databases. In contrast, none of the automated haplogrouping tools available can yet compete with manual haplogrouping using PhyloTree plus additional Web-based searches, especially when confronted with artificial recombinants still present in forensic mtDNA datasets. We review and classify the various attempts at haplogrouping by using a multiplex approach or relying on automated haplogrouping. Furthermore, we re-examine a few articles in forensic journals providing mtDNA population data where appropriate haplogrouping following PhyloTree immediately highlights several kinds of sequence errors.

Mitochondrial DNA evidence supports northeast Indian origin of the aboriginal Andamanese in the Late Paleolithic

In view of the geographically closest location to Andaman archipelago, Myanmar was suggested to be the origin place of aboriginal Andamanese. However, for lacking any genetic information from this region, which has prevented to resolve the dispute on whether the aboriginal Andamanese were originated from mainland India or Myanmar. To solve this question and better understand the origin of the aboriginal Andamanese, we screened for haplogroups M31 (from which Andaman-specific lineage M31a1 branched off) and M32 among 846 mitochondrial DNAs (mtDNAs) sampled across Myanmar. As a result, two Myanmar individuals belonging to haplogroup M31 were identified, and completely sequencing the entire mtDNA genomes of both samples testified that the two M31 individuals observed in Myanmar were probably attributed to the recent gene flow from northeast India populations. Since no root lineages of haplogroup M31 or M32 were observed in Myanmar, it is unlikely that Myanmar may serve as the source place of the aboriginal Andamanese. To get further insight into the origin of this unique population, the detailed phylogenetic and phylogeographic analyses were performed by including additional 7 new entire mtDNA genomes and 113 M31 mtDNAs pinpointed from South Asian populations, and the results suggested that Andaman-specific M31a1 could in fact trace its origin to northeast India. Time estimation results further indicated that the Andaman archipelago was likely settled by modern humans from northeast India via the land-bridge which connected the Andaman archipelago and Myanmar around the Last Glacial Maximum (LGM), a scenario in well agreement with the evidence from linguistic and palaeoclimate studies.

Southeast Asian diversity: first insights into the complex mtDNA structure of Laos

Background

Vast migrations and subsequent assimilation processes have shaped the genetic composition of Southeast Asia, an area of close contact between several major ethnic groups. To better characterize the genetic variation of this region, we analyzed the entire mtDNA control region of 214 unrelated donors from Laos according to highest forensic quality standards. To detail the phylogeny, we inspected selected SNPs from the mtDNA coding region. For a posteriori data quality control, quasi-median network constructions and autosomal STR typing were performed. In order to describe the mtDNA setup of Laos more thoroughly, the data were subjected to population genetic comparisons with 16 East Asian groups.

Results

The Laos sample exhibited ample mtDNA diversity, reflecting the huge number of ethnic groups listed. We found several new, so far undescribed mtDNA lineages in this dataset and surrounding populations. The Laos population was characteristic in terms of haplotype composition and genetic structure, however, genetic comparisons with other Southeast Asian populations revealed limited, but significant genetic differentiation. Notable differences in the maternal relationship to the major indigenous Southeast Asian ethnolinguistic groups were detected.

Conclusions

In this study, we portray the great mtDNA variety of Laos for the first time. Our findings will contribute to clarify the migration history of the region. They encourage setting up regional and subpopulation databases, especially for forensic applications. The Laotian sequences will be incorporated into the collaborative EMPOP mtDNA database http://www.empop.org upon publication and will be available as the first mtDNA reference data for this country.

HaploGrep: a fast and reliable algorithm for automatic classification of mitochondrial DNA haplogroups

An ongoing source of controversy in mitochondrial DNA (mtDNA) research is based on the detection of numerous errors in mtDNA profiles that led to erroneous conclusions and false disease associations. Most of these controversies could be avoided if the samples' haplogroup status would be taken into consideration. Knowing the mtDNA haplogroup affiliation is a critical prerequisite for studying mechanisms of human evolution and discovering genes involved in complex diseases, and validating phylogenetic consistency using haplogroup classification is an important step in quality control. However, despite the availability of Phylotree, a regularly updated classification tree of global mtDNA variation, the process of haplogroup classification is still time-consuming and error-prone, as researchers have to manually compare the polymorphisms found in a population sample to those summarized in Phylotree, polymorphism by polymorphism, sample by sample. We present HaploGrep, a fast, reliable and straight-forward algorithm implemented in a Web application to determine the haplogroup affiliation of thousands of mtDNA profiles genotyped for the entire mtDNA or any part of it. HaploGrep uses the latest version of Phylotree and offers an all-in-one solution for quality assessment of mtDNA profiles in clinical genetics, population genetics and forensics. HaploGrep can be accessed freely at http://haplogrep.uibk.ac.at.

Application of a west Eurasian-specific filter for quasi-median network analysis: Sharpening the blade for mtDNA error detection

The application of quasi-median networks provides an effective tool to check the quality of mtDNA data. Filtering of highly recurrent mutations prior to network analysis is required to simplify the data set and reduce the complexity of the network. The phylogenetic background determines those mutations that need to be filtered. While the traditional EMPOPspeedy filter was based on the worldwide mtDNA phylogeny, haplogroup-specific filters can more effectively highlight potential errors in data of the respective (sub)-continental region. In this study we demonstrate the performance of a new, west Eurasian filter EMPOPspeedyWE for the fine-tuned examination of data sets belonging to macrohaplogroup N that constitutes the main portion of mtDNA lineages in Europe. The effects on the resulting network of different database sizes, high-quality and flawed data, as well as the examination of a phylogenetically distant data set, are presented by examples. The analyses are based on a west Eurasian etalon data set that was carefully compiled from more than 3500 control region sequences for network purposes. Both, etalon data and the new filter file, are provided through the EMPOP database (www.empop.org).

Evaluating sequence-derived mtDNA length heteroplasmy by amplicon size analysis

Length heteroplasmy (LH) in mitochondrial (mt)DNA is usually observed in homopolymeric tracts and manifest as mixture of various length variants. The generally used difference-coded annotation to report mtDNA haplotypes does not express the degree of LH variation present in a sample, even more so, it is sometimes difficult to establish which length variants are present and clearly distinguishable from background noise. It has therefore become routine practice for some researchers to call the dominant type, the "major molecule", which represents the LH variant that is most abundant in a DNA extract. In the majority of cases a clear single dominant variant can be identified. However, in some samples this interpretation is difficult, i.e. when (almost) equally quantitative LH variants are present or when multiple sequencing primers result in the presentation of different dominant types. To better understand those cases we designed amplicon sizing assays for the five most relevant LH regions in the mtDNA control region (around ntps 16,189, 310, 460, 573, and the AC-repeat between 514 and 524) to determine the ratio of the LH variants by fluorescence based amplicon sizing assays. For difficult LH constellations derived by Sanger sequencing (with Big Dye terminators) these assays mostly gave clear and unambiguous results. In the vast majority of cases we found agreement between the results of the sequence and amplicon analyses and propose this alternative method in difficult cases.

SAM: String-based sequence search algorithm for mitochondrial DNA database queries

The analysis of the haploid mitochondrial (mt) genome has numerous applications in forensic and population genetics, as well as in disease studies. Although mtDNA haplotypes are usually determined by sequencing, they are rarely reported as a nucleotide string. Traditionally they are presented in a difference-coded position-based format relative to the corrected version of the first sequenced mtDNA. This convention requires recommendations for standardized sequence alignment that is known to vary between scientific disciplines, even between laboratories. As a consequence, database searches that are vital for the interpretation of mtDNA data can suffer from biased results when query and database haplotypes are annotated differently. In the forensic context that would usually lead to underestimation of the absolute and relative frequencies. To address this issue we introduce SAM, a string-based search algorithm that converts query and database sequences to position-free nucleotide strings and thus eliminates the possibility that identical sequences will be missed in a database query. The mere application of a BLAST algorithm would not be a sufficient remedy as it uses a heuristic approach and does not address properties specific to mtDNA, such as phylogenetically stable but also rapidly evolving insertion and deletion events. The software presented here provides additional flexibility to incorporate phylogenetic data, site-specific mutation rates, and other biologically relevant information that would refine the interpretation of mitochondrial DNA data. The manuscript is accompanied by freeware and example data sets that can be used to evaluate the new software (http://stringvalidation.org).

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.

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.