Gabriela Huber

The initial peopling of the Americas: a growing number of founding mitochondrial genomes from Beringia.

Pan-American mitochondrial DNA (mtDNA) haplogroup C1 has been recently subdivided into three branches, two of which (C1b and C1c) are characterized by ages and geographical distributions that are indicative of an early arrival from Beringia with Paleo-Indians. In contrast, the estimated ages of C1d--the third subset of C1--looked too young to fit the above scenario. To define the origin of this enigmatic C1 branch, we completely sequenced 63 C1d mitochondrial genomes from a wide range of geographically diverse, mixed, and indigenous American populations. The revised phylogeny not only brings the age of C1d within the range of that of its two sister clades, but reveals that there were two C1d founder genomes for Paleo-Indians. Thus, the recognized maternal founding lineages of Native Americans are at least 15, indicating that the overall number of Beringian or Asian founder mitochondrial genomes will probably increase extensively when all Native American haplogroups reach the same level of phylogenetic and genomic resolution as obtained here for C1d.

Evaluation of next generation mtGenome sequencing using the Ion Torrent Personal Genome Machine (PGM).

Insights into the human mitochondrial phylogeny have been primarily achieved by sequencing full mitochondrial genomes (mtGenomes). In forensic genetics (partial) mtGenome information can be used to assign haplotypes to their phylogenetic backgrounds, which may, in turn, have characteristic geographic distributions that would offer useful information in a forensic case. In addition and perhaps even more relevant in the forensic context, haplogroup-specific patterns of mutations form the basis for quality control of mtDNA sequences. The current method for establishing (partial) mtDNA haplotypes is Sanger-type sequencing (STS), which is laborious, time-consuming, and expensive. With the emergence of Next Generation Sequencing (NGS) technologies, the body of available mtDNA data can potentially be extended much more quickly and cost-efficiently. Customized chemistries, laboratory workflows and data analysis packages could support the community and increase the utility of mtDNA analysis in forensics. We have evaluated the performance of mtGenome sequencing using the Personal Genome Machine (PGM) and compared the resulting haplotypes directly with conventional Sanger-type sequencing. A total of 64 mtGenomes (>1 million bases) were established that yielded high concordance with the corresponding STS haplotypes (<0.02% differences). About two-thirds of the differences were observed in or around homopolymeric sequence stretches. In addition, the sequence alignment algorithm employed to align NGS reads played a significant role in the analysis of the data and the resulting mtDNA haplotypes. Further development of alignment software would be desirable to facilitate the application of NGS in mtDNA forensic genetics.

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.

Mitochondrial Haplogroup U5b3: A Distant Echo of the Epipaleolithic in Italy and the Legacy of the Early Sardinians

There are extensive data indicating that some glacial refuge zones of southern Europe (Franco-Cantabria, Balkans, and Ukraine) were major genetic sources for the human recolonization of the continent at the beginning of the Holocene. Intriguingly, there is no genetic evidence that the refuge area located in the Italian Peninsula contributed to this process. Here we show, through phylogeographic analyses of mitochondrial DNA (mtDNA) variation performed at the highest level of molecular resolution (52 entire mitochondrial genomes), that the most likely homeland for U5b3-a haplogroup present at a very low frequency across Europe-was the Italian Peninsula. In contrast to mtDNA haplogroups that expanded from other refugia, the Holocene expansion of haplogroup U5b3 toward the North was restricted by the Alps and occurred only along the Mediterranean coasts, mainly toward nearby Provence (southern France). From there, approximately 7,000-9,000 years ago, a subclade of this haplogroup moved to Sardinia, possibly as a result of the obsidian trade that linked the two regions, leaving a distinctive signature in the modern people of the island. This scenario strikingly matches the age, distribution, and postulated geographic source of a Sardinian Y chromosome haplogroup (I2a2-M26), a paradigmatic case in the European context of a founder event marking both female and male lineages.

Massively parallel sequencing of complete mitochondrial genomes from hair shaft samples

Though shed hairs are one of the most commonly encountered evidence types, they are among the most limited in terms of DNA quantity and quality. As a result, DNA testing has historically focused on the recovery of just about 600 base pairs of the mitochondrial DNA control region. Here, we describe our success in recovering complete mitochondrial genome (mtGenome) data (∼16,569bp) from single shed hairs. By employing massively parallel sequencing (MPS), we demonstrate that particular hair samples yield DNA sufficient in quantity and quality to produce 2-3kb mtGenome amplicons and that entire mtGenome data can be recovered from hair extracts even without PCR enrichment. Most importantly, we describe a small amplicon multiplex assay comprised of sixty-two primer sets that can be routinely applied to the compromised hair samples typically encountered in forensic casework. In all samples tested here, the MPS data recovered using any one of the three methods were consistent with the control Sanger sequence data developed from high quality known specimens. Given the recently demonstrated value of complete mtGenome data in terms of discrimination power among randomly sampled individuals, the possibility of recovering mtGenome data from the most compromised and limited evidentiary material is likely to vastly increase the utility of mtDNA testing for hair evidence.

Amerindian mitochondrial DNA haplogroups predominate in the population of Argentina: towards a first nationwide forensic mitochondrial DNA sequence database

The study presents South American mitochondrial DNA (mtDNA) data from selected north (N = 98), central (N = 193) and south (N = 47) Argentinean populations. Sequence analysis of the complete mtDNA control region (CR, 16024–576) resulted in 288 unique haplotypes ignoring C-insertions around positions 16193, 309, and 573; the additional analysis of coding region single nucleotide polymorphisms enabled a fine classification of the described lineages. The Amerindian haplogroups were most frequent in the north and south representing more than 60% of the sequences. A slightly different situation was observed in central Argentina where the Amerindian haplogroups represented less than 50%, and the European contribution was more relevant. Particular clades of the Amerindian subhaplogroups turned out to be nearly region-specific. A minor contribution of African lineages was observed throughout the country. This comprehensive admixture of worldwide mtDNA lineages and the regional specificity of certain clades in the Argentinean population underscore the necessity of carefully selecting regional samples in order to develop a nationwide mtDNA database for forensic and anthropological purposes. The mtDNA sequencing and analysis were performed under EMPOP guidelines in order to attain high quality for the mtDNA database.

Accumulation of mutations over the entire mitochondrial genome of breast cancer cells obtained by tissue microdissection.

The occurrence of heteroplasmy and mixtures is technically challenging for the analysis of mitochondrial DNA. More than that, observed mutations need to be carefully interpreted in the light of the phylogeny as mitochondrial DNA is a uniparental marker reflecting human evolution. Earlier attempts to explain the role of mtDNA in cancerous tissues led to substantial confusion in medical genetics mainly due to the presentation of low sequence data quality and misinterpretation of mutations representing a particular haplogroup background rather than being cancer-specific. The focus of this study is to characterize the extent and level of mutations in breast cancer samples obtained by tissue microdissection by application of an evaluated full mtDNA genome sequencing protocol. We amplified and sequenced the complete mitochondrial genomes of microdissected breast cancer cells of 15 patients and compared the results to those obtained from paired non-cancerous breast tissue derived from the same patients. We observed differences in the heteroplasmic states of substitutions between cancerous and normal cells, one of which was affecting a position that has been previously reported in lung cancer and another one that has been identified in 16 epithelial ovarian tumors, possibly indicating functional relevance. In the coding region, we found full transitions in two cancerous mitochondrial genomes and 12 heteroplasmic substitutions as compared to the non-cancerous breast cells. We identified somatic mutations over the entire mtDNA of human breast cancer cells potentially impairing the mitochondrial OXPHOS system.

Increasing the discrimination power of forensic STR testing by employing high-performance mass spectrometry, as illustrated in indigenous South African and Central Asian populations

Short tandem repeat (STR) typing has become the standard technique in forensic methodology for the identification of unknown samples. National DNA databases have been established that contain STR genotypes for intelligence purposes. Due to their success, national DNA databases have been growing so fast that the number of advantageous matches may become a logistic problem for the analysts. This is especially true for partial STR profiles as they display reduced discrimination power. To overcome this drawback, modified versions (so-called mini-STRs) of existing loci were introduced as well as new loci to improve the information content of (partial) STR profiles. We pursue an alternative approach that makes use of nucleotide variation within the amplified STR fragments, which can be discerned by mass spectrometry. We have developed an assay that determines molecular masses from crude STR amplicons which were purified and separated by a liquid chromatographic system directly hyphenated to an electrospray ionization mass spectrometer. We present here new population data of forensically relevant STRs in Khoisan and Yakut populations. These autochthonous groups were selected as they may harbor additional STR alleles that are rare or unobserved in modern humans from cosmopolitan areas, especially for the Khoisan, which are known to represent a very ancient human population. The analysis of the molecular mass of STRs offered a widened spectrum of allele variability escorted by enhanced forensic use. Thus, established STR data derived from fragment size analysis can still be used in casework or in the context of intelligence databasing.

The next generation of DNA profiling--STR typing by multiplexed PCR--ion-pair RP LC-ESI time-of-flight MS.

For the first time a multiplexed PCR approach suitable for mass spectrometric STR allele identification is presented. Thirteen forensically important STR markers (vWA, D21S11, D3S1358, D16S539, D8S1179, D7S820, D13S317, D5S818, TPOX, CSF1PO, D2S441, D10S1248, and D22S1045) and the gender typing locus amelogenin were simultaneously amplified. Ion‐pair reversed‐phase high‐performance liquid chromatography electrospray‐ionization time‐of‐flight mass spectrometry (ICEMS) was applied for genotyping, and allowed for highly efficient characterization of multiple PCR amplicons. Compared with electrophoretic sizing ICEMS enabled for the simultaneous detection of length and nucleotide variations. Thus, the obtained amount of biological information present within STR profiles was significantly increased even though the compatibility of typing results with electrophoretically generated data(bases) was maintained. Other advantages of the ICEMS platform included the abandonment of internal size standards, allelic ladders, and any kind of spectral calibration. The 14‐plex PCR was tailor‐made for ICEMS analysis by designing primer pairs that bind close to the repeat region, by using a proof reading polymerase for amplification, and by implementing molecular mass modifiers for prevention of molecular mass overlaps. In a series of experiments, the performance of the multiplexed PCR‐ICEMS assay was evaluated. The ICEMS‐based DNA profiling assay was found to be competitive regarding detection sensitivity and analyzability of degraded and casework samples with commercially available electrophoretic typing approaches, which suggests that multiplexed PCR‐ICEMS assays could represent a valuable tool for (forensic) genetics.

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).