Lack of gene-language correlation due to reciprocal female but directional male admixture in Austronesians and non-Austronesians of East Timor.

Nusa Tenggara, including East Timor, located at the crossroad between Island Southeast Asia, Near Oceania, and Australia, are characterized by a complex cultural structure harbouring speakers from two different major linguistic groups of different geographic origins (Austronesian (AN) and non-Austronesian (NAN)). This provides suitable possibilities to study gene-language relationship; however, previous studies from other parts of Nusa Tenggara reported conflicting evidence about gene-language correlation in this region. Aiming to investigate gene-language relationships including sex-mediated aspects in East Timor, we analysed the paternally inherited non-recombining part of the Y chromosome (NRY) and the maternally inherited mitochondrial (mt) DNA in a representative collection of AN- and NAN-speaking groups. Y-SNP (single-nucleotide polymorphism) data were newly generated for 273 samples and combined with previously established Y-STR (short tandem repeat) data of the same samples, and with previously established mtDNA data of 290 different samples with, however, very similar representation of geographic and linguistic coverage of the country. We found NRY and mtDNA haplogroups of previously described putative East/Southeast Asian (E/SEA) and Near Oceanian (NO) origins in both AN and NAN speakers of East Timor, albeit in different proportions, suggesting reciprocal genetic admixture between both linguistic groups for females, but directional admixture for males. Our data underline the dual genetic origin of East Timorese in E/SEA and NO, and highlight that substantial genetic admixture between the two major linguistic groups had occurred, more so via women than men. Our study therefore provides another example where languages and genes do not conform due to sex-biased genetic admixture across major linguistic groups.

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.

Full mtGenome reference data: development and characterization of 588 forensic-quality haplotypes representing three U.S. populations.

Though investigations into the use of massively parallel sequencing technologies for the generation of complete mitochondrial genome (mtGenome) profiles from difficult forensic specimens are well underway in multiple laboratories, the high quality population reference data necessary to support full mtGenome typing in the forensic context are lacking. To address this deficiency, we have developed 588 complete mtGenome haplotypes, spanning three U.S. population groups (African American, Caucasian and Hispanic) from anonymized, randomly-sampled specimens. Data production utilized an 8-amplicon, 135 sequencing reaction Sanger-based protocol, performed in semi-automated fashion on robotic instrumentation. Data review followed an intensive multi-step strategy that included a minimum of three independent reviews of the raw data at two laboratories; repeat screenings of all insertions, deletions, heteroplasmies, transversions and any additional private mutations; and a check for phylogenetic feasibility. For all three populations, nearly complete resolution of the haplotypes was achieved with full mtGenome sequences: 90.3-98.8% of haplotypes were unique per population, an improvement of 7.7-29.2% over control region sequencing alone, and zero haplotypes overlapped between populations. Inferred maternal biogeographic ancestry frequencies for each population and heteroplasmy rates in the control region were generally consistent with published datasets. In the coding region, nearly 90% of individuals exhibited length heteroplasmy in the 12418-12425 adenine homopolymer; and despite a relatively high rate of point heteroplasmy (23.8% of individuals across the entire molecule), coding region point heteroplasmies shared by more than one individual were notably absent, and transversion-type heteroplasmies were extremely rare. The ratio of nonsynonymous to synonymous changes among point heteroplasmies in the protein-coding genes (1:1.3) and average pathogenicity scores in comparison to data reported for complete substitutions in previous studies seem to provide some additional support for the role of purifying selection in the evolution of the human mtGenome. Overall, these thoroughly vetted full mtGenome population reference data can serve as a standard against which the quality and features of future mtGenome datasets (especially those developed via massively parallel sequencing) may be evaluated, and will provide a solid foundation for the generation of complete mtGenome haplotype frequency estimates for forensic applications.

Mitochondrial DNA control region analysis of three ethnic groups in the Republic of Macedonia.

A total of 444 individuals representing three ethnic groups (Albanians, Turks and Romanies) in the Republic of Macedonia were sequenced in the mitochondrial control region. The mtDNA haplogroup composition differed between the three groups. Our results showed relatively high frequencies of haplogroup H12 in Albanians (8.8%) and less in Turks (3.3%), while haplogroups M5a1 and H7a1a were dominant in Romanies (13.7% and 10.3%, respectively) but rare in the former two. This highlights the importance of regional sampling for forensic mtDNA databasing purposes. These population data will be available on EMPOP under accession numbers EMP00644 (Albanians), EMP00645 (Romanies) and EMP00646 (Turks).

Improved visibility of character conflicts in quasi-median networks with the EMPOP NETWORK software.

AIM: To provide a valuable tool for graphical representation of mitochondrial DNA (mtDNA) data that enables visual emphasis on complex substructures within the network to highlight possible ambiguities and errors. METHOD: We applied the new NETWORK graphical user interface, available via EMPOP (European DNA Profiling Group Mitochondrial DNA Population Database; www.empop.org) by means of two mtDNA data sets that were submitted for quality control. RESULTS: The quasi-median network torsi of the two data sets resulted in complex reticulations, suggesting ambiguous data. To check the corresponding raw data, accountable nodes and connecting branches of the network could be identified by highlighting induced subgraphs with concurrent dimming of their complements. This is achieved by accentuating the relevant substructures in the network: mouse clicking on a node displays a list of all mtDNA haplotypes included in that node; the selection of a branch specifies the mutation(s) connecting two nodes. It is indicated to evaluate these mutations by means of the raw data. CONCLUSION: Inspection of the raw data confirmed the presence of phantom mutations due to suboptimal electrophoresis conditions and data misinterpretation. The network software proved to be a powerful tool to highlight problematic data and guide quality control of mtDNA data tables.

Development of forensic-quality full mtGenome haplotypes: success rates with low template specimens.

Forensic mitochondrial DNA (mtDNA) testing requires appropriate, high quality reference population data for estimating the rarity of questioned haplotypes and, in turn, the strength of the mtDNA evidence. Available reference databases (SWGDAM, EMPOP) currently include information from the mtDNA control region; however, novel methods that quickly and easily recover mtDNA coding region data are becoming increasingly available. Though these assays promise to both facilitate the acquisition of mitochondrial genome (mtGenome) data and maximize the general utility of mtDNA testing in forensics, the appropriate reference data and database tools required for their routine application in forensic casework are lacking. To address this deficiency, we have undertaken an effort to: (1) increase the large-scale availability of high-quality entire mtGenome reference population data, and (2) improve the information technology infrastructure required to access/search mtGenome data and employ them in forensic casework. Here, we describe the application of a data generation and analysis workflow to the development of more than 400 complete, forensic-quality mtGenomes from low DNA quantity blood serum specimens as part of a U.S. National Institute of Justice funded reference population databasing initiative. We discuss the minor modifications made to a published mtGenome Sanger sequencing protocol to maintain a high rate of throughput while minimizing manual reprocessing with these low template samples. The successful use of this semi-automated strategy on forensic-like samples provides practical insight into the feasibility of producing complete mtGenome data in a routine casework environment, and demonstrates that large (>2kb) mtDNA fragments can regularly be recovered from high quality but very low DNA quantity specimens. Further, the detailed empirical data we provide on the amplification success rates across a range of DNA input quantities will be useful moving forward as PCR-based strategies for mtDNA enrichment are considered for targeted next-generation sequencing workflows.

Mass spectrometric base composition profiling: Implications for forensic mtDNA databasing.

In forensic genetics mitochondrial DNA (mtDNA) is usually analyzed by direct Sanger-type sequencing (STS). This method is known to be laborious and sometimes prone to human error. Alternative methods have been proposed that lead to faster results. Among these are methods that involve mass-spectrometry resulting in base composition profiles that are, by definition, less informative than the full nucleotide sequence. Here, we applied a highly automated electrospray ionization mass spectrometry (ESI-MS) system (PLEX-ID) to an mtDNA population study to compare its performance with respect to throughput and concordance to STS. We found that the loss of information power was relatively low compared to the gain in speed and analytical standardization. The detection of point and length heteroplasmy turned out to be roughly comparable between the technologies with some individual differences related to the processes. We confirm that ESI-MS provides a valuable platform for analyzing mtDNA variation that can also be applied in the forensic context.

Concept for estimating mitochondrial DNA haplogroups using a maximum likelihood approach (EMMA).

The assignment of haplogroups to mitochondrial DNA haplotypes contributes substantial value for quality control, not only in forensic genetics but also in population and medical genetics. The availability of Phylotree, a widely accepted phylogenetic tree of human mitochondrial DNA lineages, led to the development of several (semi-)automated software solutions for haplogrouping. However, currently existing haplogrouping tools only make use of haplogroup-defining mutations, whereas private mutations (beyond the haplogroup level) can be additionally informative allowing for enhanced haplogroup assignment. This is especially relevant in the case of (partial) control region sequences, which are mainly used in forensics. The present study makes three major contributions toward a more reliable, semi-automated estimation of mitochondrial haplogroups. First, a quality-controlled database consisting of 14,990 full mtGenomes downloaded from GenBank was compiled. Together with Phylotree, these mtGenomes serve as a reference database for haplogroup estimates. Second, the concept of fluctuation rates, i.e. a maximum likelihood estimation of the stability of mutations based on 19,171 full control region haplotypes for which raw lane data is available, is presented. Finally, an algorithm for estimating the haplogroup of an mtDNA sequence based on the combined database of full mtGenomes and Phylotree, which also incorporates the empirically determined fluctuation rates, is brought forward. On the basis of examples from the literature and EMPOP, the algorithm is not only validated, but both the strength of this approach and its utility for quality control of mitochondrial haplotypes is also demonstrated.

Forensic and phylogeographic characterisation of mtDNA lineages from Somalia.

The African mitochondrial (mt) phylogeny is coarsely resolved but the majority of population data generated so far is limited to the analysis of the first hypervariable segment (HVS-1) of the control region (CR). Therefore, this study aimed on the investigation of the entire CR of 190 unrelated Somali individuals to enrich the severely underrepresented African mtDNA pool. The majority (60.5 %) of the haplotypes were of sub-Saharan origin with L0a1d, L2a1h and L3f being the most frequently observed haplogroups. This is in sharp contrast to previous data reported from the Y-chromosome, where only about 5 % of the observed haplogroups were of sub-Saharan provenance. We compared the genetic distances based on population pairwise F (st) values between 11 published East, Central and North African as well as western Asian populations and the Somali sequences and displayed them in a multi-dimensional scaling plot. Genetic proximity evidenced by clustering roughly reflected the relative geographic location of the populations. The sequences will be included in the EMPOP database ( www.empop.org ) under accession number EMP00397 upon publication (Parson and Dür Forensic Sci Int Genet 1:88-92, 2007).

Mitochondrial DNA control region data from indigenous Angolan Khoe-San lineages.

Here we provide 129 complete mitochondrial control region sequences of indigenous Khoe-San individuals from Angola to contribute to the still underrepresented pool of data from Africa. The dataset consists of exclusively African lineages with a majority of Sub-Saharan haplogroups. The probability of a random match was calculated as 0.09. The data set comprises 21 haplotypes occurring more than once and 17 unique haplotypes. Upon publication, haplotypes were incorporated in the EMPOP database (www.empop.org; EMP00069) [1].

Rapid coastal spread of First Americans: novel insights from South America's Southern Cone mitochondrial genomes.

It is now widely agreed that the Native American founders originated from a Beringian source population ~15-18 thousand years ago (kya) and rapidly populated all of the New World, probably mainly following the Pacific coastal route. However, details about the migration into the Americas and the routes pursued on the continent still remain unresolved, despite numerous genetic, archaeological, and linguistic investigations. To examine the pioneering peopling phase of the South American continent, we screened literature and mtDNA databases and identified two novel mitochondrial DNA (mtDNA) clades, here named D1g and D1j, within the pan-American haplogroup D1. They both show overall rare occurrences but local high frequencies, and are essentially restricted to populations from the Southern Cone of South America (Chile and Argentina). We selected and completely sequenced 43 D1g and D1j mtDNA genomes applying highest quality standards. Molecular and phylogeographic analyses revealed extensive variation within each of the two clades and possibly distinct dispersal patterns. Their age estimates agree with the dating of the earliest archaeological sites in South America and indicate that the Paleo-Indian spread along the entire longitude of the American double continent might have taken even <2000 yr. This study confirms that major sampling and sequencing efforts are mandatory for uncovering all of the most basal variation in the Native American mtDNA haplogroups and for clarification of Paleo-Indian migrations, by targeting, if possible, both the general mixed population of national states and autochthonous Native American groups, especially in South America.

Timing and deciphering mitochondrial DNA macro-haplogroup R0 variability in Central Europe and Middle East.

BACKGROUND: Nearly half of the West Eurasian assemblage of human mitochondrial DNA (mtDNA) is fractioned into numerous sub-lineages of the predominant haplogroup (hg) R0. Several hypotheses have been proposed on the origin and the expansion times of some R0 sub-lineages, which were partially inconsistent with each other. Here we describe the phylogenetic structure and genetic variety of hg R0 in five European populations and one population from the Middle East. RESULTS: Our analysis of 1,350 mtDNA haplotypes belonging to R0, including entire control region sequences and 45 single nucleotide polymorphisms from the coding region, revealed significant differences in the distribution of different sub-hgs even between geographically closely located regions. Estimates of coalescence times that were derived using diverse algorithmic approaches consistently affirmed that the major expansions of the different R0 hgs occurred in the terminal Pleistocene and early Holocene. CONCLUSION: Given an estimated coalescence time of the distinct lineages of 10-18 kya, the differences in the distributions could hint to either limited maternal gene flow after the Last Glacial Maximum due to the alpine nature of the regions involved or to a stochastic loss of diversity due to environmental events and/or disease episodes occurred at different times and in distinctive regions. Our comparison of two different ways of obtaining the timing of the most recent common ancestor confirms that the time of a sudden expansion can be adequately recovered from control region data with valid confidence intervals. For reliable estimates, both procedures should be applied in order to cross-check the results for validity and soundness.