mtDNAmanager: a Web-based tool for the management and quality analysis of mitochondrial DNA control-region sequences

Unraveling the Genetic Threads of History: mtDNA HVS-I Analysis Reveals the Ancient Past of the Aburra Valley

This article presents a comprehensive genetic study focused on pre-Hispanic individuals who inhabited the Aburrá Valley in Antioquia, Colombia, between the tenth and seventeenth centuries AD. Employing a genetic approach, the study analyzed maternal lineages using DNA samples obtained from skeletal remains. The results illuminate a remarkable degree of biological diversity within these populations and provide insights into their genetic connections with other ancient and indigenous groups across the American continent. The findings strongly support the widely accepted hypothesis that the migration of the first American settlers occurred through Beringia, a land bridge connecting Siberia to North America during the last Ice Age. Subsequently, these early settlers journeyed southward, crossing the North American ice cap. Of particular note, the study unveils the presence of ancestral lineages from Asian populations, which played a pivotal role in populating the Americas. The implications of these results extend beyond delineating migratory routes and settlement patterns of ancient populations. They also enrich our understanding of the genetic diversity inherent in indigenous populations of the region. By revealing the genetic heritage of pre-Hispanic individuals from the Aburrá Valley, this study offers valuable insights into the history of human migration and settlement in the Americas. Furthermore, it enhances our comprehension of the intricate genetic tapestry that characterizes indigenous communities in the area.

Mitochondrial Haplogroup Classification of Ancient DNA Samples Using Haplotracker

Mitochondrial DNA haplogroup classification is used to study maternal lineage of ancient human populations. The haplogrouping of ancient DNA is not easy because the DNA is usually found in small pieces in limited quantities. We have developed Haplotracker, a straightforward and efficient high-resolution haplogroup classification tool optimized specifically for ancient DNA samples. Haplotracker offers a user-friendly input interface for multiple mitochondrial DNA sequence fragments in a sample. It provides accurate haplogroup classification with full-length mitochondrial genome sequences and provides high-resolution haplogroup predictions for some fragmented control region sequences using a novel algorithm built on Phylotree mtDNA Build 17 (Phylotree) and our haplotype database (n = 118,869). Its performance for accuracy was demonstrated to be high through haplogroup classification using 8,216 Phylotree full-length and control region mitochondrial DNA sequences compared with HaploGrep 2, one of the most accurate current haplogroup classifiers. Haplotracker provides a novel haplogroup tracking solution for fragmented sequences to track subhaplogroups or verify the haplogroups efficiently. Using Haplotracker, we classified mitochondrial haplogroups to the final subhaplogroup level in nine ancient DNA samples extracted from human skeletal remains found in 2,000-year-old elite Xiongnu cemetery in Northeast Mongolia. Haplotracker can be freely accessed at https://haplotracker.cau.ac.kr.

From Forensics to Clinical Research: Expanding the Variant Calling Pipeline for the Precision ID mtDNA Whole Genome Panel

Despite a multitude of methods for the sample preparation, sequencing, and data analysis of mitochondrial DNA (mtDNA), the demand for innovation remains, particularly in comparison with nuclear DNA (nDNA) research. The Applied Biosystems™ Precision ID mtDNA Whole Genome Panel (Thermo Fisher Scientific, USA) is an innovative library preparation kit suitable for degraded samples and low DNA input. However, its bioinformatic processing occurs in the enterprise Ion Torrent Suite™ Software (TSS), yielding BAM files aligned to an unorthodox version of the revised Cambridge Reference Sequence (rCRS), with a heteroplasmy threshold level of 10%. Here, we present an alternative customizable pipeline, the PrecisionCallerPipeline (PCP), for processing samples with the correct rCRS output after Ion Torrent sequencing with the Precision ID library kit. Using 18 samples (3 original samples and 15 mixtures) derived from the 1000 Genomes Project, we achieved overall improved performance metrics in comparison with the proprietary TSS, with optimal performance at a 2.5% heteroplasmy threshold. We further validated our findings with 50 samples from an ongoing independent cohort of stroke patients, with PCP finding 98.31% of TSS's variants (TSS found 57.92% of PCP's variants), with a significant correlation between the variant levels of variants found with both pipelines.

A benchmarking of human mitochondrial DNA haplogroup classifiers from whole-genome and whole-exome sequence data

The mitochondrial genome (mtDNA) is of interest for a range of fields including evolutionary, forensic, and medical genetics. Human mitogenomes can be classified into evolutionary related haplogroups that provide ancestral information and pedigree relationships. Because of this and the advent of high-throughput sequencing (HTS) technology, there is a diversity of bioinformatic tools for haplogroup classification. We present a benchmarking of the 11 most salient tools for human mtDNA classification using empirical whole-genome (WGS) and whole-exome (WES) short-read sequencing data from 36 unrelated donors. We also assessed the best performing tool in third-generation long noisy read WGS data obtained with nanopore technology for a subset of the donors. We found that, for short-read WGS, most of the tools exhibit high accuracy for haplogroup classification irrespective of the input file used for the analysis. However, for short-read WES, Haplocheck and MixEmt were the most accurate tools. Based on the performance shown for WGS and WES, and the accompanying qualitative assessment, Haplocheck stands out as the most complete tool. For third-generation HTS data, we also showed that Haplocheck was able to accurately retrieve mtDNA haplogroups for all samples assessed, although only after following assembly-based approaches (either based on a referenced-based assembly or a hybrid de novo assembly). Taken together, our results provide guidance for researchers to select the most suitable tool to conduct the mtDNA analyses from HTS data.

Massively parallel sequencing of human skeletal remains in Vietnam using the precision ID mtDNA control region panel on the Ion S5™ system

Mitochondrial DNA (mtDNA) analysis using Sanger sequencing has been a routine practice for the identification of human skeletal remains. However, this process is usually challenging since DNA from the remains is highly degraded and at low concentration. Recently, the advent and implementation of massively parallel sequencing (MPS) have been offered the ability to improve mtDNA sequence data for forensic analysis. To assess the utility of the Ion S5™ system - an MPS platform for mtDNA analysis in challenging samples, we sequenced the mitochondrial control region of 52 age-old skeletal remains. Using the Precision ID mtDNA Control Region Panel, 50 full and two partial control region haplotypes at relatively high mean coverage of 2494 × were achieved for variant calling. Further variant analysis at 10% threshold for point heteroplasmy showed high degradation degree in terms of DNA damage in our bone samples. A higher point heteroplasmy threshold of 20% was required to diminish most of background noise caused by the damage. The results from this study indicated the potential application of the Ion S5™ system in sequencing degraded samples in Vietnam and provided valuable data sources for forensic analyses in the future.

Genetic characterisation of the North-West Indian populations: analysis of mitochondrial DNA control region variations

BACKGROUND

Human mitochondrial DNA presents several interesting characteristics, making it a favourable tool in the field of molecular anthropology, medical genetics, population history, and forensic science.

AIM

The present study investigated the mitochondrial DNA (mtDNA) control region variations in diverse ethnic groups of North-West India for which population data is insufficient.

SUBJECTS AND METHODS

The complete mtDNA control regions of 197 unrelated (for up to three generations) healthy individuals belonging to different ethnic groups of North-West India were sequenced. The haplotype frequencies, haplogroup distribution, and pairwise F_ST values between the studied and other worldwide populations were generated to study patterns of variation in human mtDNA.

RESULTS

The results ascertained high gene diversity (0.998) in the studied maternal lineages, identifying 166 distinct haplotypes, of which 158 were unique and characterised by 117 variable sites. Three haplogroups: M3, M30, and U7 were observed to be the most prevalent, and phylogeographically a total of 55.86% of sequences were characterised into South Asian, followed by West Eurasian (40.18%) and East Asian (3.96%), ancestry haplogroups. Pairwise genetic differentiation comparisons revealed maternal homogeneity in the studied groups. No population substructure was detected within the North-West Indian populations.

CONCLUSION

The results of this preliminary study will contribute to an existing database of mtDNA variations of the Indian population and facilitate prospective studies investigating population genetics and human diseases.

Contaminations in (meta)genome data: An open issue for the scientific community

In recent years, the high throughput and the low cost of next-generation sequencing (NGS) technologies have led to an increase of the amount of (meta)genomic data, revolutionizing genomic research studies. However, the quality of sequencing data could be affected by experimental errors derived from defective methods and protocols. This represents a serious problem for the scientific community with a negative impact on the correctness of studies that involve genomic sequence analysis. As a countermeasure, several alignment and taxonomic classification tools have been developed to uncover and correct errors. In this critical review some of these integrated software tools and pipelines used to detect contaminations in reference genome databases and sequenced samples are reported. In particular, case studies of bacterial contaminations, contaminations of human origin, mitochondrial contaminations of ancient DNA, and cross contaminations are examined.

Development of a Novel ex vivo Nasal Epithelial Cell Model Supporting Colonization With Human Nasal Microbiota

The nasal mucosa provides first line defense against inhaled pathogens while creating a unique microenvironment for bacterial communities. Studying the impact of microbiota in the nasal cavity has been difficult due to limitations with current models including explant cultures, primary cells, or neoplastic cell lines. Most notably, none have been shown to support reproducible colonization by bacterial communities from human donors. Therefore, to conduct controlled studies of the human nasal ecosystem, we have developed a novel ex vivo mucosal model that supports bacterial colonization of a cultured host mucosa created by immortalized human nasal epithelial cells (NEC). For this model, immortalized NEC established from 5 male and 5 female donors were cultured with an air-interfaced, apical surface on a porous transwell membrane. NEC were grown from nasal turbinate tissues harvested from willed bodies or from discarded tissue collected during sinonasal procedures. Immortalized cells were evaluated through molecular verification of cell type, histological confirmation of tissue differentiation including formation of tight junctions, NEC multilayer viability, metabolism, physiology and imaging of the luminal surface by scanning electron microscopy. Results showed proper differentiation and multilayer formation at seven to 10 days after air interface that was maintained for up to 3 weeks. The optimized mucosal cultures created an environment necessary to sustain colonization by nasal microbiomes (NMBs) that were collected from healthy volunteers, cryogenically preserved and characterized with customized quantitative polymerase chain reaction (qPCR) arrays. Polymicrobial communities of nasal bacteria associated with healthy and inflamed states were consistently reproduced in matured NEC co-cultures by transplant of NMBs from multiple community types. The cultured NMBs were stable after an initial period of bacterial replication and equilibration. This novel ex vivo culture system is the first model that supports controlled cultivation of NMBs, allowing for lab-based causation studies and further experimentation to explore the complexities of host-microbe and microbe-microbe interactions.

A review of bioinformatic methods for forensic DNA analyses

Short tandem repeats, single nucleotide polymorphisms, and whole mitochondrial analyses are three classes of markers which will play an important role in the future of forensic DNA typing. The arrival of massively parallel sequencing platforms in forensic science reveals new information such as insights into the complexity and variability of the markers that were previously unseen, along with amounts of data too immense for analyses by manual means. Along with the sequencing chemistries employed, bioinformatic methods are required to process and interpret this new and extensive data. As more is learnt about the use of these new technologies for forensic applications, development and standardization of efficient, favourable tools for each stage of data processing is being carried out, and faster, more accurate methods that improve on the original approaches have been developed. As forensic laboratories search for the optimal pipeline of tools, sequencer manufacturers have incorporated pipelines into sequencer software to make analyses convenient. This review explores the current state of bioinformatic methods and tools used for the analyses of forensic markers sequenced on the massively parallel sequencing (MPS) platforms currently most widely used.

The Sommersdorf mummies-An interdisciplinary investigation on human remains from a 17th-19th century aristocratic crypt in southern Germany

Sommersdorf Castle (Bavaria, Germany) is a medieval castle complex which has been inhabited by the aristocratic family von Crailsheim. The deceased were entombed in a crypt located in the parapets underneath the castle’s church, resulting in mummification of the bodies. Based on the family chronicle and oral history, identities have been ascribed to the mummies. The aim of the study is therefore to test the accuracy of the historical records in comparison to archaeological, anthropological and genetic data. Today, the crypt houses eleven wooden coffins from the 17^th to 19^th century AD. In ten of these, mummified and scattered human remains were found. Archive records were studied in order to identify names, ancestry, titles, occupation, date of birth and death, and place of interment of the individuals. The coffins were visually inspected and dated by typo-chronology, and the mummified and scattered skeletal remains were subjected to a physical anthropological examination. In total, the crypt contains the remains of a minimum number of nine individuals, among them three adult males, five adult females and one infant. A detailed scientific examination, including prior conservation, ancient DNA analyses, and computed tomography (CT), was performed on five mummies. By means of the CT data age at death, sex, body height, pathologies, and anatomical variants were investigated. CT analysis further showed that the bodies were naturally mummified. Mitochondrial DNA analyses revealed that the tested individuals are not maternally related. In addition, health, living conditions and circumstances of death of the entombed individuals could be highlighted. Being confronted with the strengths, weaknesses and limitations of each methodological approach, probable identification was achieved in two cases.

MitoSuite: a graphical tool for human mitochondrial genome profiling in massive parallel sequencing

Recent rapid advances in high-throughput, next-generation sequencing (NGS) technologies have promoted mitochondrial genome studies in the fields of human evolution, medical genetics, and forensic casework. However, scientists unfamiliar with computer programming often find it difficult to handle the massive volumes of data that are generated by NGS. To address this limitation, we developed MitoSuite, a user-friendly graphical tool for analysis of data from high-throughput sequencing of the human mitochondrial genome. MitoSuite generates a visual report on NGS data with simple mouse operations. Moreover, it analyzes high-coverage sequencing data but runs on a stand-alone computer, without the need for file upload. Therefore, MitoSuite offers outstanding usability for handling massive NGS data, and is ideal for evolutionary, clinical, and forensic studies on the human mitochondrial genome variations. It is freely available for download from the website https://mitosuite.com.

Mitochondrial super-haplogroup U diversity in Serbians

BACKGROUND

Available mitochondrial (mtDNA) data demonstrate genetic differentiation among South Slavs inhabiting the Balkan Peninsula. However, their resolution is insufficient to elucidate the female-specific aspects of the genetic history of South Slavs, including the genetic impact of various migrations which were rather common within the Balkans, a region having a turbulent demographic history.

AIM

The aim was to thoroughly study complete mitogenomes of Serbians, a population linking westward and eastward South Slavs.

SUBJECTS AND METHODS

Forty-six predominantly Serbian super-haplogroup U complete mitogenomes were analysed phylogenetically against ∼4000 available complete mtDNAs of modern and ancient Western Eurasians.

RESULTS

Serbians share a number of U mtDNA lineages with Southern, Eastern-Central and North-Western Europeans. Putative Balkan-specific lineages (e.g. U1a1c2, U4c1b1, U5b3j, K1a4l and K1a13a1) and lineages shared among Serbians (South Slavs) and West and East Slavs were detected (e.g. U2e1b1, U2e2a1d, U4a2a, U4a2c, U4a2g1, U4d2b and U5b1a1).

CONCLUSION

The exceptional diversity of maternal lineages found in Serbians may be associated with the genetic impact of both autochthonous pre-Slavic Balkan populations whose mtDNA gene pool was affected by migrations of various populations over time (e.g. Bronze Age pastoralists) and Slavic and Germanic newcomers in the early Middle Ages.

Molecular Genealogy of a Mongol Queen's Family and Her Possible Kinship with Genghis Khan

Members of the Mongol imperial family (designated the Golden family) are buried in a secret necropolis; therefore, none of their burial grounds have been found. In 2004, we first discovered 5 graves belonging to the Golden family in Tavan Tolgoi, Eastern Mongolia. To define the genealogy of the 5 bodies and the kinship among them, SNP and/or STR profiles of mitochondria, autosomes, and Y chromosomes were analyzed. Four of the 5 bodies were determined to carry the mitochondrial DNA haplogroup D4, while the fifth carried haplogroup CZ, indicating that this individual had no kinship with the others. Meanwhile, Y-SNP and Y-STR profiles indicate that the males examined belonged to the R1b-M343 haplogroup. Thus, their East Asian D4 or CZ matrilineal and West Eurasian R1b-M343 patrilineal origins reveal genealogical admixture between Caucasoid and Mongoloid ethnic groups, despite a Mongoloid physical appearance. In addition, Y chromosomal and autosomal STR profiles revealed that the four D4-carrying bodies bore the relationship of either mother and three sons or four full siblings with almost the same probability. Moreover, the geographical distribution of R1b-M343-carrying modern-day individuals demonstrates that descendants of Tavan Tolgoi bodies today live mainly in Western Eurasia, with a high frequency in the territories of the past Mongol khanates. Here, we propose that Genghis Khan and his family carried Y-haplogroup R1b-M343, which is prevalent in West Eurasia, rather than the Y-haplogroup C3c-M48, which is prevalent in Asia and which is widely accepted to be present in the family members of Genghis Khan. Additionally, Tavan Tolgoi bodies may have been the product of marriages between the lineage of Genghis Khan’s Borjigin clan and the lineage of either the Ongud or Hongirad clans, indicating that these individuals were members of Genghis Khan’s immediate family or his close relatives.

The Use and Effectiveness of Triple Multiplex System for Coding Region Single Nucleotide Polymorphism in Mitochondrial DNA Typing of Archaeologically Obtained Human Skeletons from Premodern Joseon Tombs of Korea

Previous study showed that East Asian mtDNA haplogroups, especially those of Koreans, could be successfully assigned by the coupled use of analyses on coding region SNP markers and control region mutation motifs. In this study, we tried to see if the same triple multiplex analysis for coding regions SNPs could be also applicable to ancient samples from East Asia as the complementation for sequence analysis of mtDNA control region. By the study on Joseon skeleton samples, we know that mtDNA haplogroup determined by coding region SNP markers successfully falls within the same haplogroup that sequence analysis on control region can assign. Considering that ancient samples in previous studies make no small number of errors in control region mtDNA sequencing, coding region SNP analysis can be used as good complimentary to the conventional haplogroup determination, especially of archaeological human bone samples buried underground over long periods.

Mitochondrial DNA variation in the Viking age population of Norway

The medieval Norsemen or Vikings had an important biological and cultural impact on many parts of Europe through raids, colonization and trade, from about AD 793 to 1066. To help understand the genetic affinities of the ancient Norsemen, and their genetic contribution to the gene pool of other Europeans, we analysed DNA markers in Late Iron Age skeletal remains from Norway. DNA was extracted from 80 individuals, and mitochondrial DNA polymorphisms were detected by next-generation sequencing. The sequences of 45 ancient Norwegians were verified as genuine through the identification of damage patterns characteristic of ancient DNA. The ancient Norwegians were genetically similar to previously analysed ancient Icelanders, and to present-day Shetland and Orkney Islanders, Norwegians, Swedes, Scots, English, German and French. The Viking Age population had higher frequencies of K*, U*, V* and I* haplogroups than their modern counterparts, but a lower proportion of T* and H* haplogroups. Three individuals carried haplotypes that are rare in Norway today (U5b1b1, Hg A* and an uncommon variant of H*). Our combined analyses indicate that Norse women were important agents in the overseas expansion and settlement of the Vikings, and that women from the Orkneys and Western Isles contributed to the colonization of Iceland.

Joseon funerary texts tested using ancient DNA analysis of a Korean mummy

In Korea, ancient DNA (aDNA) analysis has been applied to investigations into the genetic affiliations of mummies found in Joseon Dynasty tombs (1392-1910 CE), becoming now indispensable tool for researches studying human remains from archaeological sites. In the course of our recent examinations on a Korean mummy of Joseon Dynasty, we discovered many teeth contained in a pouch. And in fact, the historical literature on the topic of Joseon funerals contain general accounts of pouches in which an individual's lost teeth were collected over the course of a lifetime and, after death, placed in the coffin with the body. To test the veracity of the historical texts, the present study undertook aDNA analyses and compared the results between specifically questioned (Q) samples (teeth) and known (K) samples (brain and bone) from the mummy to ensure that they came from the same individual. Although the Q-K comparison of autosomal short tandem repeat results did not show full concordance due to allelic drop-outs in some loci, our statistical calculation indicated that the teeth in the pouch are highly likely those of the mummy. Additionally, Q-K comparison of mitochondrial DNA sequence results showed 100% matches between samples. There results, in short, could not gainsay the conjecture that the teeth samples originated from the person buried in the tomb; and if so, he must have kept his teeth for a long time after their loss. As the application of aDNA analysis to Korean mummy studies develops, there will be other opportunities to test historical documents, particularly those referring to funerary rites.

Breast cancer risk and genetic ancestry: a case-control study in Uruguay

Background

Uruguay exhibits one of the highest rates of breast cancer in Latin America, similar to those of developed nations, the reasons for which are not completely understood. In this study we investigated the effect that ancestral background has on breast cancer susceptibility among Uruguayan women.

Methods

We carried out a case–control study of 328 (164 cases, 164 controls) women enrolled in public hospitals and private clinics across the country. We estimated ancestral proportions using a panel of nuclear and mitochondrial ancestry informative markers (AIMs) and tested their association with breast cancer risk.

Results

Nuclear individual ancestry in cases was (mean ± SD) 9.8 ± 7.6% African, 13.2 ± 10.2% Native American and 77.1 ± 13.1% European, and in controls 9.1 ± 7.5% African, 14.7 ± 11.2% Native American and 76.2 ± 14.2% European. There was no evidence of a difference in nuclear or mitochondrial ancestry between cases and controls. However, European mitochondrial haplogroup H was associated with breast cancer (OR = 2.0; 95% CI 1.1, 3.5).

Conclusions

We have not found evidence that overall genetic ancestry differs between breast cancer patients and controls in Uruguay but we detected an association of the disease with a European mitochondrial lineage, which warrants further investigation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12905-015-0171-8) contains supplementary material, which is available to authorized users.

Mitochondrial DNA perspective of Serbian genetic diversity

Although south-Slavic populations have been studied to date from various aspects, the population of Serbia, occupying the central part of the Balkan Peninsula, is still genetically understudied at least at the level of mitochondrial DNA (mtDNA) variation. We analyzed polymorphisms of the first and the second mtDNA hypervariable segments (HVS-I and HVS-II) and informative coding-region markers in 139 Serbians to shed more light on their mtDNA variability, and used available data on other Slavic and neighboring non-Slavic populations to assess their interrelations in a broader European context. The contemporary Serbian mtDNA profile is consistent with the general European maternal landscape having a substantial proportion of shared haplotypes with eastern, central, and southern European populations. Serbian population was characterized as an important link between easternmost and westernmost south-Slavic populations due to the observed lack of genetic differentiation with all other south-Slavic populations and its geographical positioning within the Balkan Peninsula. An increased heterogeneity of south Slavs, most likely mirroring turbulent demographic events within the Balkan Peninsula over time (i.e., frequent admixture and differential introgression of various gene pools), and a marked geographical stratification of Slavs to south-, east-, and west-Slavic groups, were also found. A phylogeographic analyses of 20 completely sequenced Serbian mitochondrial genomes revealed not only the presence of mtDNA lineages predominantly found within the Slavic gene pool (U4a2a*, U4a2a1, U4a2c, U4a2g, HV10), supporting a common Slavic origin, but also lineages that may have originated within the southern Europe (H5*, H5e1, H5a1v) and the Balkan Peninsula in particular (H6a2b and L2a1k).

Mitochondrial DNA diversity in the African American population

Genetic polymorphism along mitochondrial DNA (mtDNA) defines population-specific signatures called mtDNA haplogroups. Estimation of mtDNA haplogroup distribution may be prone to errors, notably if the study sample is not drawn from a multicenter cohort. Here, we report on mtDNA diversity in a sample of African American individuals (n = 343) enrolled in a multicenter cohort. Sequencing of the hypervariable regions I and II of the D-loop control region showed that the most common mitochondrial variants are 73G, 146C, 150T, 152C, 189G, 16278T, and 16311C. In agreement with the published data, we observed 17 common mtDNA haplogroups: L0, L1, L1b, L1c, L2, L2a, L2b, L2c, L2e, L3, L3b, L3d, L3e, L3f, L3h, L3x, and L4. The most commonly observed haplogroup is L2a (19.8%), followed by L1b (10.2%). Overall, the observed mtDNA haplogroup distribution in our study is similar to those published for the African American and the African populations.

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.

HAPLOFIND: a new method for high-throughput mtDNA haplogroup assignment

Deep sequencing technologies are completely revolutionizing the approach to DNA analysis. Mitochondrial DNA (mtDNA) studies entered in the "postgenomic era": the burst in sequenced samples observed in nuclear genomics is expected also in mitochondria, a trend that can already be detected checking complete mtDNA sequences database submission rate. Tools for the analysis of these data are available, but they fail in throughput or in easiness of use. We present here a new pipeline based on previous algorithms, inherited from the "nuclear genomic toolbox," combined with a newly developed algorithm capable of efficiently and easily classify new mtDNA sequences according to PhyloTree nomenclature. Detected mutations are also annotated using data collected from publicly available databases. Thanks to the analysis of all freely available sequences with known haplogroup obtained from GenBank, we were able to produce a PhyloTree-based weighted tree, taking into account each haplogroup pattern conservation. The combination of a highly efficient aligner, coupled with our algorithm and massive usage of asynchronous parallel processing, allowed us to build a high-throughput pipeline for the analysis of mtDNA sequences that can be quickly updated to follow the ever-changing nomenclature. HaploFind is freely accessible at the following Web address: https://haplofind.unibo.it.

An update to MitoTool: using a new scoring system for faster mtDNA haplogroup determination

The determination of human mitochondrial DNA (mtDNA) haplogroups is not only crucial in anthropological and forensic studies, but is also helpful in the medical field to prevent establishment of wrong disease associations. In recent years, high-throughput technologies and the huge amounts of data they create, as well as the regular updates to the mtDNA phylogenetic tree, mean that there is a need for an automated approach which can make a speedier determination of haplogroups than can be made by using the traditional manual method. Here, we update the MitoTool (www.mitotool.org) by incorporating a novel scoring system for the determination of mtDNA into haplogroups, which has advantages on speed, accuracy and ease of implementation. In order to make the access to MitoTool easier, we also provide a stand-alone version of the program that will run on a local computer and this version is freely available at the MitoTool website.

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.

mtDNAoffice: a software to assign human mtDNA macro haplogroups through automated analysis of the protein coding region

We describe a fast, automated process to determine distances between mtDNA sequences allowing their subsequent clustering and haplogroup assignment that may increase the speed of data analysis and avoid human errors. In order to avoid complexities/ambiguities resulting from recurrence and insertion/deletion phenomena and thus improving evolutionary signal-to-noise ratio, protein coding sequences were compared using a vectorial representation method, and the corresponding genetic distance matrix was used for the construction of a neighbor-joining/UPGMA tree or an MDS graphic, which generally agrees with the consensus mtDNA phylogeny. mtDNAoffice software, detailed instructions and example files are freely available on the web at http://www.portugene.com/SupMat/setupmtDNAoffice.rar.

Mitochondrial haplogroups and polymorphisms reveal no association with sporadic prostate cancer in a southern European population

BACKGROUND

It is known that mitochondria play an important role in certain cancers (prostate, renal, breast, or colorectal) and coronary disease. These organelles play an essential role in apoptosis and the production of reactive oxygen species; in addition, mtDNA also reveals the history of populations and ancient human migration. All these events and variations in the mitochondrial genome are thought to cause some cancers, including prostate cancer, and also help us to group individuals into common origin groups. The aim of the present study is to analyze the different haplogroups and variations in the sequence in the mitochondrial genome of a southern European population consisting of subjects affected (n = 239) and non-affected (n = 150) by sporadic prostate cancer.

METHODOLOGY AND PRINCIPAL FINDINGS

Using primer extension analysis and DNA sequencing, we identified the nine major European haplogroups and CR polymorphisms. The frequencies of the haplogroups did not differ between patients and control cohorts, whereas the CR polymorphism T16356C was significantly higher in patients with PC compared to the controls (p = 0.029). PSA, staging, and Gleason score were associated with none of the nine major European haplogroups. The CR polymorphisms G16129A (p = 0.007) and T16224C (p = 0.022) were significantly associated with Gleason score, whereas T16311C (p = 0.046) was linked with T-stage.

CONCLUSIONS AND SIGNIFICANCE

Our results do not suggest that mtDNA haplogroups could be involved in sporadic prostate cancer etiology and pathogenesis as previous studies performed in middle Europe population. Although some significant associations have been obtained in studying CR polymorphisms, further studies should be performed to validate these results.

Ancient DNA from an Early Neolithic Iberian population supports a pioneer colonization by first farmers

The Neolithic transition has been widely debated particularly regarding the extent to which this revolution implied a demographic expansion from the Near East. We attempted to shed some light on this process in northeastern Iberia by combining ancient DNA (aDNA) data from Early Neolithic settlers and published DNA data from Middle Neolithic and modern samples from the same region. We successfully extracted and amplified mitochondrial DNA from 13 human specimens, found at three archaeological sites dated back to the Cardial culture in the Early Neolithic (Can Sadurní and Chaves) and to the Late Early Neolithic (Sant Pau del Camp). We found that haplogroups with a low frequency in modern populations-N* and X1-are found at higher frequencies in our Early Neolithic population (∼31%). Genetic differentiation between Early and Middle Neolithic populations was significant (F(ST) ∼0.13, P<10(-5)), suggesting that genetic drift played an important role at this time. To improve our understanding of the Neolithic demographic processes, we used a Bayesian coalescence-based simulation approach to identify the most likely of three demographic scenarios that might explain the genetic data. The three scenarios were chosen to reflect archaeological knowledge and previous genetic studies using similar inferential approaches. We found that models that ignore population structure, as previously used in aDNA studies, are unlikely to explain the data. Our results are compatible with a pioneer colonization of northeastern Iberia at the Early Neolithic characterized by the arrival of small genetically distinctive groups, showing cultural and genetic connections with the Near East.

DNA typing for the identification of old skeletal remains from Korean War victims

The identification of missing casualties of the Korean War (1950-1953) has been performed using mitochondrial DNA (mtDNA) profiles, but recent advances in DNA extraction techniques and approaches using smaller amplicons have significantly increased the possibility of obtaining DNA profiles from highly degraded skeletal remains. Therefore, 21 skeletal remains of Korean War victims and 24 samples from biological relatives of the supposed victims were selected based on circumstantial evidence and/or mtDNA-matching results and were analyzed to confirm the alleged relationship. Cumulative likelihood ratios were obtained from autosomal short tandem repeat, Y-chromosomal STR, and mtDNA-genotyping results, and mainly confirmed the alleged relationship with values over 10⁵. The present analysis emphasizes the value of mini- and Y-STR systems as well as an efficient DNA extraction method in DNA testing for the identification of old skeletal remains.

Forensic and genetic characterization of mtDNA from Pathans of Pakistan

Complete mitochondrial control region data were generated for 230 unrelated Pathans from North West Frontier Province and Federally Administered Tribal Areas of Pakistan. To confirm data quality and to explore the genetic structure of Pathans, mitochondrial DNA haplogroup affiliation was determined by shared haplogroup-specific polymorphisms in the control region and by the analysis of diagnostic coding region single-nucleotide polymorphisms using a multiplex system for the assignment of eight haplogroups: M, N1'5, W, R, R0, T, J, and U. Sequence comparison revealed that 193 haplotypes were defined by 215 variable sites when major insertions were ignored at nucleotide positions 16193, 309, and 573. From a phylogenetic perspective, Pathans have a heterogeneous origin, displaying a high percentage of West Eurasian haplogroups followed by haplogroups native to South Asia and a small fraction from East Asian lineages. In population comparisons, this ethnic group differed significantly from several other ethnic groups from Pakistan and surrounding countries. These results suggest that frequency estimates for mtDNA haplotypes should be determined for endogamous ethnic groups individually instead of pooling data for these subpopulations into a single dataset for the Pakistani population. Data presented here may contribute to the accuracy of forensic mtDNA comparisons in the Pathans of Pakistan.

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.

HaploSearch: a tool for haplotype-sequence two-way transformation

Comparison of available mitochondrial DNA data is some times hindered by the data presentation format. HaploSearch is a simple tool for transforming DNA sequences into haplotype data and vice versa, speeding up the manipulation of large datasets. Although designed for mitochondrial DNA, HaploSearch could be used with any kind of DNA type. HaploSearch program, detailed software instructions and example files are freely available on the web http://www.haplosite.com/haplosearch/.

MitoTool: a web server for the analysis and retrieval of human mitochondrial DNA sequence variations

MitoTool, a web-based bioinformatics platform, is designed for deciphering human mitochondrial DNA (mtDNA) data in batch mode. The platform has advantages in (i) parsing diverse types of mtDNA data; (ii) automatically classifying haplogroup according to mtDNA sequences or variants; (iii) discovering possibly missing variants of the samples with claimed haplogroups status; (iv) estimating the evolutionary conservation index, protein coding effect and potential pathogenicity of certain substitutions; (v) performing statistical analysis for haplogroup distribution frequency between case and control groups. Furthermore, it offers an integrated database for retrieving five types of mitochondrion-related information. The MitoTool is freely accessed at http://www.mitotool.org.

Mitochondrial DNA variation in an isolated area of Central Italy

BACKGROUND

The genetic variation in Italy is the result of ancient population movement, demographic change, and geography. The increasing possibility of studying the maternal genetic structure of selected Italian population samples at a high level of phylogenetic resolution provides a particularly useful model to assess the presence of genetic traces of the ancient people who lived in Italy in pre-Roman times in present populations

AIM

In this study we reconstructed the genetic maternal history of Jenne and Vallepietra, two mountain communities in the Aniene Valley in the Simbruini Mountains near Rome. Both communities have been spared external invasion due to their geographic location, which very likely preserved the genetic pool of these autochthonous populations.

SUBJECTS AND METHODS

The study population (124 individuals from Jenne and Vallepietra) were investigated for D-loop mtDNA hypervariable segments I (HVS-I) and II (HVS-II) and for informative single nucleotide polymorphisms (SNPs) within the coding region. The detected haplotypes were then compared with those of other Italian, European and Mediterranean populations.

RESULTS

The distribution of mtDNA diversity in Jenne and Vallepietra, although similar to that found in other European populations, shows a basic variability and the typical signs of a certain degree of isolation between them and other populations analysed; in particular, the Vallepietra sample showed an unusually high frequency (71.3%) of mtDNA haplogroups which are typical of Near Eastern and South-Western Asian populations.

CONCLUSION

The high degree of differentiation between the two villages is intriguing, since it suggests a low level of gene flow between them, despite their close geographic proximity and shared linguistic features. The degree of their genetic isolation, also in comparison to other Italian, European and Mediterranean populations, is consistent with isolation among geographically separated populations.

Reconstructing the population history of Nicaragua by means of mtDNA, Y-chromosome STRs, and autosomal STR markers

Before the arrival of the Spaniards in Nicaragua, diverse Native American groups inhabited the territory. In colonial times, Native Nicaraguan populations interacted with Europeans and slaves from Africa. To ascertain the extent of this genetic admixture and provide genetic evidence about the origin of the Nicaraguan ancestors, we analyzed the mitochondrial control region (HVSI and HVSII), 17 Y chromosome STRs, and 15 autosomal STRs in 165 Mestizo individuals from Nicaragua. To carry out interpopulation comparisons, HVSI sequences from 29 American populations were compiled from the literature. The results reveal a close relationship between Oto-manguean, Uto-Aztecan, Mayan groups from Mexico, and a Chibchan group to Nicaraguan lineages. The Native American contribution to present-day Nicaraguan Mestizos accounts for most of the maternal lineages, whereas the majority of Nicaraguan Y chromosome haplogroups can be traced back to a West Eurasian origin. Pairwise Fst distances based on Y-STRs between Nicaragua and European, African and Native American populations show that Nicaragua is much closer to Europeans than the other populations. Additionally, admixture proportions based on autosomal STRs indicate a predominantly Spanish contribution. Our study reveals that the Nicaraguan Mestizo population harbors a high proportion of European male and Native American female substrate. Finally, the amount of African ancestry is also interesting, probably because of the contribution of Spanish conquerors with North African genetic traces or that of West African slaves.

No relationship found between point heteroplasmy in mitochondrial DNA control region and age range, sex and haplogroup in human hairs

The analysis of heteroplasmy (presence of more than one type of mitochondrial DNA in an individual) is used as a tool in human identification studies, anthropology, and most currently in studies that relate heteroplasmy with longevity. The frequency of heteroplasmy and its correlation with age has been analyzed using different tissues such as blood, muscle, heart, bone and brain and in different regions of mitochondrial DNA, but this analysis had never been performed using hair samples. In this study, samples of hair were sequenced in order to ascertain whether the presence or not of heteroplasmy varied according to age, sex and origin of haplogroup individuals. The samples were grouped by age (3 groups), gender (male and female) and haplogroup of origin (European, African and Native American), and analyzed using the chi-square statistical test (χ(2)). Based in statistical results obtained, we conclude that there is no relationship between heteroplasmy and sex, age and haplogroup origin using hair samples.

eCOMPAGT integrates mtDNA: import, validation and export of mitochondrial DNA profiles for population genetics, tumour dynamics and genotype-phenotype association studies

Background

Mitochondrial DNA (mtDNA) is widely being used for population genetics, forensic DNA fingerprinting and clinical disease association studies. The recent past has uncovered severe problems with mtDNA genotyping, not only due to the genotyping method itself, but mainly to the post-lab transcription, storage and report of mtDNA genotypes.

Description

eCOMPAGT, a system to store, administer and connect phenotype data to all kinds of genotype data is now enhanced by the possibility of storing mtDNA profiles and allowing their validation, linking to phenotypes and export as numerous formats. mtDNA profiles can be imported from different sequence evaluation programs, compared between evaluations and their haplogroup affiliations stored. Furthermore, eCOMPAGT has been improved in its sophisticated transparency (support of MySQL and Oracle), security aspects (by using database technology) and the option to import, manage and store genotypes derived from various genotyping methods (SNPlex, TaqMan, and STRs). It is a software solution designed for project management, laboratory work and the evaluation process all-in-one.

Conclusions

The extended mtDNA version of eCOMPAGT was designed to enable error-free post-laboratory data handling of human mtDNA profiles. This software is suited for small to medium-sized human genetic, forensic and clinical genetic laboratories. The direct support of MySQL and the improved database security options render eCOMPAGT a powerful tool to build an automated workflow architecture for several genotyping methods. eCOMPAGT is freely available at http://dbis-informatik.uibk.ac.at/ecompagt.

Human mitochondrial DNA variation in Southern Italy

BACKGROUND

Since prehistoric times Southern Italy has been a cultural crossroads of the Mediterranean basin. Genetic data on the peoples of Basilicata and Calabria are scarce and, particularly, no records on mtDNA variability have been published.

AIM

In this study mtDNA haplotypes of populations from Basilicata, Calabria and Sicily are compared with those of other Italian and Mediterranean populations, so as to investigate their genetic relationships.

SUBJECTS AND METHODS

A total of 341 individuals was analysed for mtDNA in order to provide their classification into haplogroups. Multivariate analysis was used to compare the studied populations with other Mediterranean samples; median-joining network analysis was applied to observe the relationship between the major lineages of the Southern Italians.

RESULTS

The haplogroup distribution in the Southern Italian samples falls within the typical pattern of mtDNA variability of Western Eurasia. The comparison with other Mediterranean countries showed a substantial homogeneity of the area, which is probably related to the historic contact through the Mediterranean Sea.

CONCLUSION

The mtDNA analysis demonstrated that Southern Italy displays a typical pattern of Mediterranean basin variability, even though it appears plausible that Southern Italy was less affected by the effects of the Late Glacial Maximum, which reduced genetic diversity in Europe.

A new deep branch of eurasian mtDNA macrohaplogroup M reveals additional complexity regarding the settlement of Madagascar

Background

Current models propose that mitochondrial DNA macrohaplogroups M and N evolved from haplogroup L3 soon after modern humans left Africa. Increasingly, however, analysis of isolated populations is filling in the details of, and in some cases challenging, aspects of this general model.

Results

Here, we present the first comprehensive study of three such isolated populations from Madagascar: the Mikea hunter-gatherers, the neighbouring Vezo fishermen, and the Merina central highlanders (n = 266). Complete mitochondrial DNA genome sequences reveal several unresolved lineages, and a new, deep branch of the out-of-Africa founder clade M has been identified. This new haplogroup, M23, has a limited global distribution, and is restricted to Madagascar and a limited range of African and Southwest Asian groups.

Conclusions

The geographic distribution, phylogenetic placement and molecular age of M23 suggest that the colonization of Madagascar was more complex than previously thought.