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

Blind study evaluation illustrates utility of the Ion PGM™ system for use in human identity DNA typing

Aim

To perform a blind study to assess the capability of the Ion Personal Genome Machine™ (PGM) system to sequence forensically relevant genetic marker panels and to characterize unknown individuals for ancestry and possible relatedness.

Methods

Twelve genomic samples were provided by a third party for blinded genetic analysis. For these 12 samples, the mitochondrial genome and three PGM™ panels containing human identity single nucleotide polymorphisms (SNPs), ancestry informative SNPs, and short tandem repeats (STRs) were sequenced on the PGM™ system and analyzed.

Results

All four genetic systems were run and analyzed on the PGM™ system in a reasonably quick time frame. Completeness of genetic profiles, depth of coverage, strand balance, and allele balance were informative metrics that illustrated the quality and reliability of the data produced. SNP genotypes allowed for identification of sex, paternal lineage, and population ancestry. STR genotypes were shown to be in complete concordance with genotypes generated by standard capillary electrophoresis-based technologies. Variants in the mitochondrial genome data provided information on population background and maternal relationships.

Conclusion

All results from analysis of the 12 genomic samples were consistent with sample information provided by the sample providers at the end of the blinded study. The relatively easy identification of intra-STR allele SNPs offered the potential for increased discrimination power. The promising nature of these results warrants full validation studies of this massively parallel sequencing technology and its further development for forensic data analysis.

The mitochondrial landscape of African Americans: An examination of more than 2500 control region haplotypes from 22 U.S. locations

The mitochondrial DNA (mtDNA) control region (16024-576) was Sanger-sequenced for a total of 2563 self-identified African Americans, using automated processing techniques and data review standards exceeding guidelines for forensic applications. Genetic diversity ranged from 0.9952 to 0.9998 in 22 population samples from 20 different states. Haplogroups of African ancestry, found in 82.48% of individuals overall, were most concentrated in the Southeast U.S. and decreased to the north and west. West African and West Central African haplotypes were well-represented in the population samples, especially in the southern U.S. states, while East African haplogroups were observed in low-frequency clusters in a handful of locations across the country. East Asian, Native American, and West Eurasian admixture was present in 3.16%, 2.93%, and 11.43% of samples, respectively. While some geographic substructure was detected across the population samples as clines in admixture frequencies, 20 of the 22 population samples were found to be statistically indistinguishable by pairwise comparisons and AMOVA calculations. Datasets from Hawaii and Idaho, however, were clear outliers. Overall, these more than 2500 control region sequences represent the most comprehensive regional sampling of African American mtDNA diversity to date, and are suitable for use in a forensic mtDNA database. The population data are made available via EMPOP (www.empop.org) and GenBank.

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.

mtDNA analysis of 174 Eurasian populations using a new iterative rank correlation method

In this study, we analyse 27-dimensional mtDNA haplogroup distributions of 174 Eurasian, North-African and American populations, including numerous ancient data as well. The main contribution of this work was the description of the haplogroup distribution of recent and ancient populations as compounds of certain hypothetic ancient core populations immediately or indirectly determining the migration processes in Eurasia for a long time. To identify these core populations, we developed a new iterative algorithm determining clusters of the 27 mtDNA haplogroups studied having strong rank correlation among each other within a definite subset of the populations. Based on this study, the current Eurasian populations can be considered as compounds of three early core populations regarding to maternal lineages. We wanted to show that a simultaneous analysis of ancient and recent data using a new iterative rank correlation algorithm and the weighted SOC learning technique may reveal the most important and deterministic migration processes in the past. This technique allowed us to determine geographically, historically and linguistically well-interpretable clusters of our dataset having a very specific, hardly classifiable structure. The method was validated using a 2-dimensional stepping stone model.

Admixture and genetic diversity distribution patterns of non-recombining lineages of Native American ancestry in Colombian populations

Genetic diversity of present American populations results from very complex demographic events involving different types and degrees of admixture. Through the analysis of lineage markers such as mtDNA and Y chromosome it is possible to recover the original Native American haplotypes, which remained identical since the admixture events due to the absence of recombination. However, the decrease in the effective population sizes and the consequent genetic drift effects suffered by these populations during the European colonization resulted in the loss or under-representation of a substantial fraction of the Native American lineages. In this study, we aim to clarify how the diversity and distribution of uniparental lineages vary with the different demographic characteristics (size, degree of isolation) and the different levels of admixture of extant Native groups in Colombia. We present new data resulting from the analyses of mtDNA whole control region, Y chromosome SNP haplogroups and STR haplotypes, and autosomal ancestry informative insertion-deletion polymorphisms in Colombian individuals from different ethnic and linguistic groups. The results demonstrate that populations presenting a high proportion of non-Native American ancestry have preserved nevertheless a substantial diversity of Native American lineages, for both mtDNA and Y chromosome. We suggest that, by maintaining the effective population sizes high, admixture allowed for a decrease in the effects of genetic drift due to Native population size reduction and thus resulting in an effective preservation of the Native American non-recombining lineages.

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.

Evidence for frequent and tissue-specific sequence heteroplasmy in human mitochondrial DNA

Mitochondrial point heteroplasmy is a common event observed not only in patients with mitochondrial diseases but also in healthy individuals. We here report a comprehensive investigation of heteroplasmy occurrence in human including the whole mitochondrial control region from nine different tissue types of 100 individuals. Sanger sequencing was used as a standard method and results were supported by cloning, minisequencing, and massively parallel sequencing. Only 12% of all individuals showed no heteroplasmy, whereas 88% showed at least one heteroplasmic position within the investigated tissues. In 66% of individuals up to 8 positions were affected. The highest relative number of heteroplasmies was detected in muscle and liver (79%, 69%), followed by brain, hair, and heart (36.7%-30.2%). Lower percentages were observed in bone, blood, lung, and buccal cells (19.8%-16.2%). Accumulation of position-specific heteroplasmies was found in muscle (positions 64, 72, 73, 189, and 408), liver (position 72) and brain (partial deletion at position 71). Deeper analysis of these specific positions in muscle revealed a non-random appearance and position-specific dependency on age. MtDNA heteroplasmy frequency and its potential functional importance have been underestimated in the past and its occurrence is ubiquitous and dependent at least on age, tissue, and position-specific mutation rates.

Phy-Mer: a novel alignment-free and reference-independent mitochondrial haplogroup classifier

MOTIVATION

All current mitochondrial haplogroup classification tools require variants to be detected from an alignment with the reference sequence and to be properly named according to the canonical nomenclature standards for describing mitochondrial variants, before they can be compared with the haplogroup determining polymorphisms. With the emergence of high-throughput sequencing technologies and hence greater availability of mitochondrial genome sequences, there is a strong need for an automated haplogroup classification tool that is alignment-free and agnostic to reference sequence.

RESULTS

We have developed a novel mitochondrial genome haplogroup-defining algorithm using a k-mer approach namely Phy-Mer. Phy-Mer performs equally well as the leading haplogroup classifier, HaploGrep, while avoiding the errors that may occur when preparing variants to required formats and notations. We have further expanded Phy-Mer functionality such that next-generation sequencing data can be used directly as input.

AVAILABILITY AND IMPLEMENTATION

Phy-Mer is publicly available under the GNU Affero General Public License v3.0 on GitHub (https://github.com/danielnavarrogomez/phy-mer).

CONTACT

Xiaowu_Gai@meei.harvard.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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

Portuguese mitochondrial DNA genetic diversity-An update and a phylogenetic revision

In recent years a large amount of mitochondrial population data for forensic purposes has been produced. Current efforts are focused at increasing the number of studied populations while generating updated genetic information of forensic quality. However, complete mitochondrial control region sequences are still scarce for most populations and even more so for complete mitochondrial genomes. In the case of Portugal, previous population genetics studies have already revealed the general portrait of HVS-I and HVS-II mitochondrial diversity, becoming now important to update and expand the mitochondrial region analysed. Accordingly, a total of 292 complete control region sequences from continental Portugal were obtained, under a stringent experimental design to ensure the quality of data through double sequencing of each target region. Furthermore, H-specific coding region SNPs were examined to detail haplogroup classification and complete mitogenomes were obtained for all sequences belonging to haplogroups U4 and U5. In general, a typical Western European haplogroup composition was found in mainland Portugal, associated to high level of mitochondrial genetic diversity. Within the country, no signs of substructure were detected. The typing of extra coding region SNPs has provided the refinement or confirmation of the previous classification obtained with EMMA tool in 96% of the cases. Finally, it was also possible to enlarge haplogroup U phylogeny with 28 new U4 and U5 mitogenomes.

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.

mitoSAVE: mitochondrial sequence analysis of variants in Excel

The mitochondrial genome (mtGenome) contains genetic information amenable to numerous applications such as medical research, population and evolutionary studies, and human identity testing. However, inconsistent nomenclature assignment makes haplotype comparison difficult and can lead to false exclusion of potentially useful profiles. Massively Parallel Sequencing (MPS) is a platform for sequencing large datasets and potentially whole populations with relative ease. However, the data generated are not easily parsed and interpreted. With this in mind, mitoSAVE has been developed to enable fast conversion of Variant Call Format (VCF) files. mitoSAVE is an Excel-based workbook that converts data within the VCF into mtDNA haplotypes using phylogenetically-established nomenclature as well as rule-based alignments consistent with current forensic standards. mitoSAVE is formatted for human mitochondrial genome; however, it can easily be adapted to support other reasonably small genomes.

Strong genetic admixture in the Altai at the Middle Bronze Age revealed by uniparental and ancestry informative markers

The Altai Mountains have been a long-term boundary zone between the Eurasian Steppe populations and South and East Asian populations. To disentangle some of the historical population movements in this area, 14 ancient human specimens excavated in the westernmost part of the Mongolian Altai were studied. Thirteen of them were dated from the Middle to the End of the Bronze Age and one of them to the Eneolithic period. The environmental conditions encountered in this region led to the good preservation of DNA in the human remains. Therefore, a multi-markers approach was adopted for the genetic analysis of identity, ancestry and phenotype markers. Mitochondrial DNA analyses revealed that the ancient Altaians studied carried both Western (H, U, T) and Eastern (A, C, D) Eurasian lineages. In the same way, the patrilineal gene pool revealed the presence of different haplogroups (Q1a2a1-L54, R1a1a1b2-Z93 and C), probably marking different origins for the male paternal lineages. To go further in the search of the origin of these ancient specimens, phenotypical characters (i.e. hair and eye color) were determined. For this purpose, we adapted the HIrisPlex assay recently described to MALDI-TOF mass spectrometry. In addition, some ancestry informative markers were analyzed with this assay. The results revealed mixed phenotypes among this group confirming the probable admixed ancestry of the studied Altaian population at the Middle Bronze Age. The good results obtained from ancient DNA samples suggest that this approach might be relevant for forensic casework too.

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.

Reprint of: 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 64mtGenomes (>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.

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.