Liane Fendt

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.

Sequencing strategy for the whole mitochondrial genome resulting in high quality sequences.

BackgroundIt has been demonstrated that a reliable and fail-safe sequencing strategy is mandatory for high-quality analysis of mitochondrial (mt) DNA, as the sequencing and base-calling process is prone to error. Here, we present a high quality, reliable and easy handling manual procedure for the sequencing of full mt genomes that is also appropriate for laboratories where fully automated processes are not available.ResultsWe amplified whole mitochondrial genomes as two overlapping PCR-fragments comprising each about 8500 bases in length. We developed a set of 96 primers that can be applied to a (manual) 96 well-based technology, which resulted in at least double strand sequence coverage of the entire coding region (codR).ConclusionThis elaborated sequencing strategy is straightforward and allows for an unambiguous sequence analysis and interpretation including sometimes challenging phenomena such as point and length heteroplasmy that are relevant for the investigation of forensic and clinical samples.

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.

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.

Polymorphisms of mtDNA control region in Tunisian and Moroccan populations: An enrichment of forensic mtDNA databases with Northern Africa data

Current forensic mitochondrial (mt)DNA databases are limited in representative population data of African origin. We investigated HVS-I/HVS-II sequences of 120 Tunisian and Moroccan healthy male donors applying stringent quality criteria to assure high quality of the data and phylogenetic alignment and notation of the sequences. Among 64 Tunisians, 56 different haplotypes were observed and the most common haplotype (16187T 16189C 16223T 16264T 16270T 16278T 16293G 16311C 73G 152C 182T 185T 195C 247A 263G 309.1C 315.1C; haplogroup (hg) L1b) was shared by four individuals. 56 Moroccans could be assigned to 52 different haplotypes where the most common haplotype was of West Eurasian origin with the hg H sequence motif 263G 315.1C and variations in the HVS-II polyC-stretch (309.1C 309.2C) shared by six samples. The majority of the observed haplotypes belong to the west Eurasian phylogeny (50% in Tunisians and 62.5% in Moroccans). Our data are consistent with the current phylogeographic knowledge displaying the occurrence of sub-Saharan haplogroup L sequences, found in 48.4% of Tunisians and 25% of Moroccans as well as the presence of the two re-migrated haplogroups U6 (7.8% and 1.8% in Tunisians and Moroccans, respectively) and M1 (1.6% in Tunisians and 8.9% in Moroccans).

MtDNA diversity of Ghana: a forensic and phylogeographic view

West Africa is characterized by a migration history spanning more than 150,000 years. Climate changes but also political circumstances were responsible for several early but also recent population movements that shaped the West African mitochondrial landscape. The aim of the study was to establish a Ghanaian mtDNA dataset for forensic purposes and to investigate the diversity of the Ghanaian population sample with respect to surrounding populations. We sequenced full mitochondrial control regions of 193 Akan people from Ghana and excluded two apparently close maternally related individuals due to preceding kinship testing. The remaining dataset comprising 191 sequences was applied as etalon for quasi-median network analysis and was subsequently combined with 99 additional control region sequences from surrounding West African countries. All sequences were incorporated into the EMPOP database enriching the severely underrepresented African mtDNA pool. For phylogeographic considerations, the Ghanaian haplotypes were compared to those of 19 neighboring populations comprising a total number of 6198 HVS1 haplotypes. We found extensive genetic admixture between the Ghanaian lineages and those from adjacent populations diminishing with geographical distance. The extent of genetic admixture reflects the long but also recent history of migration waves within West Africa mainly caused by changing environmental conditions. Also, evidence for potential socio-economical influences such as trade routes is provided by the occurrence of U6b and U6d sequences found in Dubai but also in Tunisia leading to the African West Coast via Mauritania and Senegal but also via Niger, Nigeria to Cameroon.

Validation of Next-Generation Sequencing of Entire Mitochondrial Genomes and the Diversity of Mitochondrial DNA Mutations in Oral Squamous Cell Carcinoma

Background Oral squamous cell carcinoma (OSCC) is mainly caused by smoking and alcohol abuse and shows a five-year survival rate of ~50%. We aimed to explore the variation of somatic mitochondrial DNA (mtDNA) mutations in primary oral tumors, recurrences and metastases. Methods We performed an in-depth validation of mtDNA next-generation sequencing (NGS) on an Illumina HiSeq 2500 platform for its application to cancer tissues, with the goal to detect low-level heteroplasmies and to avoid artifacts. Therefore we genotyped the mitochondrial genome (16.6 kb) from 85 tissue samples (tumors, recurrences, resection edges, metastases and blood) collected from 28 prospectively recruited OSCC patients applying both Sanger sequencing and high-coverage NGS (~35,000 reads per base). Results We observed a strong correlation between Sanger sequencing and NGS in estimating the mixture ratio of heteroplasmies (r = 0.99; p<0.001). Non-synonymous heteroplasmic variants were enriched among cancerous tissues. The proportions of somatic and inherited variants in a given gene region were strongly correlated (r = 0.85; p<0.001). Half of the patients shared mutations between benign and cancerous tissue samples. Low level heteroplasmies (<10%) were more frequent in benign samples compared to tumor samples, where heteroplasmies >10% were predominant. Four out of six patients who developed a local tumor recurrence showed mutations in the recurrence that had also been observed in the primary tumor. Three out of five patients, who had tumor metastases in the lymph nodes of their necks, shared mtDNA mutations between primary tumors and lymph node metastases. The percentage of mutation heteroplasmy increased from the primary tumor to lymph node metastases. Conclusions We conclude that Sanger sequencing is valid for heteroplasmy quantification for heteroplasmies ≥10% and that NGS is capable of reliably detecting and quantifying heteroplasmies down to the 1%-level. The finding of shared mutations between primary tumors, recurrences and metastasis indicates a clonal origin of malignant cells in oral cancer.

Identification of West Eurasian mitochondrial haplogroups by mtDNA SNP screening: Results of the 2006-2007 EDNAP collaborative exercise

The European DNA Profiling (EDNAP) Group performed a collaborative exercise on a mitochondrial (mt) DNA screening assay that targeted 16 nucleotide positions in the coding region and allowed for the discrimination of major west Eurasian mtDNA haplogroups. The purpose of the exercise was to evaluate the stability and reproducibility of the self-developed multiplex-PCR and multiplex-single base extension kit by blind-testing saliva and hair shaft samples provided by the organizing laboratory. The overall success rate in obtaining useful results was high given that some of the participating laboratories had no previous experience with the technology and/or mtDNA analysis. The results of this collaborative exercise stimulate the expansion of screening methods in forensic laboratories to increase efficiency and performance of mtDNA typing, and thus demonstrates that mtDNA SNP typing is a powerful tool for forensic casework analysis.

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