Frequency and pattern of heteroplasmy in the control region of human mitochondrial DNA

Nanobiopsy investigation of the subcellular mtDNA heteroplasmy in human tissues

Mitochondrial function is critical to continued cellular vitality and is an important contributor to a growing number of human diseases. Mitochondrial dysfunction is typically heterogeneous, mediated through the clonal expansion of mitochondrial DNA (mtDNA) variants in a subset of cells in a given tissue. To date, our understanding of the dynamics of clonal expansion of mtDNA variants has been technically limited to the single cell-level. Here, we report the use of nanobiopsy for subcellular sampling from human tissues, combined with next-generation sequencing to assess subcellular mtDNA mutation load in human tissue from mitochondrial disease patients. The ability to map mitochondrial mutation loads within individual cells of diseased tissue samples will further our understanding of mitochondrial genetic diseases.

Mitochondrial DNA control region typing from highly degraded skeletal remains by single-multiplex next-generation sequencing

Poor nuclear DNA preservation from highly degraded skeletal remains is the most limiting factor for the genetic identification of individuals. Mitochondrial DNA (mtDNA) typing, and especially of the control region (CR), using next-generation sequencing (NGS), enables retrieval of valuable genetic information in forensic contexts where highly degraded human skeletal remains are the only source of genetic material. Currently, NGS commercial kits can type all mtDNA-CR in fewer steps than the conventional Sanger technique. The PowerSeq CRM Nested System kit (Promega Corporation) employs a nested multiplex-polymerase chain reaction (PCR) strategy to amplify and index all mtDNA-CR in a single reaction. Our study analyzes the success of mtDNA-CR typing of highly degraded human skeletons using the PowerSeq CRM Nested System kit. We used samples from 41 individuals from different time periods to test three protocols (M1, M2, and M3) based on modifications of PCR conditions. To analyze the detected variants, two bioinformatic procedures were compared: an in-house pipeline and the GeneMarker HTS software. The results showed that many samples were not analyzed when the standard protocol (M1) was used. In contrast, the M3 protocol, which includes 35 PCR cycles and longer denaturation and extension steps, successfully recovered the mtDNA-CR from highly degraded skeletal samples. Mixed base profiles and the percentage of damaged reads were both indicators of possible contamination and can provide better results if used together. Furthermore, our freely available in-house pipeline can provide variants concordant with the forensic software.

The complete mitochondrial genome of Trigonisca nataliae (Hymenoptera, Apidae) assemblage reveals heteroplasmy in the control region

The evolution of mitochondrial genomes in the stingless bees is surprisingly dynamic, making them a model system to understand mitogenome structure, function, and evolution. Out of the seven mitogenomes available in this group, five exhibit atypical characteristics, including extreme rearrangements, rapid evolution and complete mitogenome duplication. To further explore the mitogenome diversity in these bees, we utilized isolated mtDNA and Illumina sequencing to assemble the complete mitogenome of Trigonisca nataliae, a species found in Northern Brazil. The mitogenome of T. nataliae was highly conserved in gene content and structure when compared to Melipona species but diverged in the control region (CR). Using PCR amplification, cloning and Sanger sequencing, six different CR haplotypes, varying in size and content, were recovery. These findings indicate that heteroplasmy, where different mitochondrial haplotypes coexist within individuals, occurs in T. nataliae. Consequently, we argue that heteroplasmy might indeed be a common phenomenon in bees that could be associated with variations in mitogenome size and challenges encountered during the assembly process.

DNA and protein analyses of hair in forensic genetics

Hair is one of the most common pieces of biological evidence found at a crime scene and plays an essential role in forensic investigation. Hairs, especially non-follicular hairs, are usually found at various crime scenes, either by natural shedding or by forcible shedding. However, the genetic material in hairs is usually highly degraded, which makes forensic analysis difficult. As a result, the value of hair has not been fully exploited in forensic investigations and trials. In recent years, with advances in molecular biology, forensic analysis of hair has achieved remarkable strides and provided crucial clues in numerous cases. This article reviews recent developments in DNA and protein analysis of hair and attempts to provide a comprehensive solution to improve forensic hair analysis.

Forensic evaluation of mitochondrial DNA heteroplasmy in Gujarat population, India

BACKGROUND

Owing to its high copy number and its small size, mtDNA analysis is the most reliable choice when biological materials from crime scenes are degraded or have mixed STR profiles.

AIM

To examine the occurrence of heteroplasmy along with its frequency and pattern in both HV1 and HV2 regions of the mtDNA among unrelated individuals from India.

SUBJECTS AND METHODS

Mitochondrial DNA control region [hypervariable region one (HV1) and hypervariable region two (HV2)] were analysed in blood and buccal tissues of 104 unrelated individuals from the Indian state of Gujarat.

RESULTS

A high frequency of point heteroplasmy (PH) and length heteroplasmy (LH) was revealed. PH was detected in 7.69% of the population, with a higher frequency observed in blood than in buccal samples. However, there were no statistically significant differences in PH between the two tissues (Chi-square = 0.552, p ≥ 0.05). A total of six PH positions were detected: three at HV1, and another three at HV2. The studied population showed 46.15% LH in the HV1 and HV2 regions of both tissues. The LH positions observed in the Gujarat population were the same as those previously reported at HV1 np16184-16193 and HV2 np303-315.

CONCLUSIONS

Our findings suggest that differences in the pattern of heteroplasmy found in different tissues can complicate the forensic analysis, on the other hand, the probability of a match between the questioned and reference samples increases when the heteroplasmy is identical in both tissues. Variability of PH among persons and even within tissues recommends analysing multiple tissue samples before drawing a conclusion in forensic mtDNA analyses.

Pedigree derived mutation rate across the entire mitochondrial genome of the Norfolk Island population

Estimates of mutation rates for various regions of the human mitochondrial genome (mtGenome) vary widely, depending on whether they are inferred using a phylogenetic approach or obtained directly from pedigrees. Traditionally, only the control region, or small portions of the coding region have been targeted for analysis due to the cost and effort required to produce whole mtGenome Sanger profiles. Here, we report one of the first pedigree derived mutation rates for the entire human mtGenome. The entire mtGenome from 225 individuals originating from Norfolk Island was analysed to estimate the pedigree derived mutation rate and compared against published mutation rates. These individuals were from 45 maternal lineages spanning 345 generational events. Mutation rates for various portions of the mtGenome were calculated. Nine mutations (including two transitions and seven cases of heteroplasmy) were observed, resulting in a rate of 0.058 mutations/site/million years (95% CI 0.031-0.108). These mutation rates are approximately 16 times higher than estimates derived from phylogenetic analysis with heteroplasmy detected in 13 samples (n = 225, 5.8% individuals). Providing one of the first pedigree derived estimates for the entire mtGenome, this study provides a better understanding of human mtGenome evolution and has relevance to many research fields, including medicine, anthropology and forensics.

Evolution, systematics, and the unnatural history of mitochondrial DNA

The tenets underlying the use of mtDNA in phylogenetic and systematic analyses are strict maternal inheritance, clonality, homoplasmy, and difference due to mutation: that is, there are species-specific mtDNA sequences and phylogenetic reconstruction is a matter of comparing these sequences and inferring closeness of relatedness from the degree of sequence similarity. Yet, how mtDNA behavior became so defined is mysterious. Even though early studies of fertilization demonstrated for most animals that not only the head, but the sperm's tail and mitochondria-bearing midpiece penetrate the egg, the opposite - only the head enters the egg - became fact, and mtDNA conceived as maternally transmitted. When midpiece/tail penetration was realized as true, the conceptions 'strict maternal inheritance', etc., and their application to evolutionary endeavors, did not change. Yet there is mounting evidence of paternal mtDNA transmission, paternal and maternal combination, intracellular recombination, and intra- and intercellular heteroplasmy. Clearly, these phenomena impact the systematic and phylogenetic analysis of mtDNA sequences.

Genomic Landscape of the Mitochondrial Genome in the United Arab Emirates Native Population

In order to assess the genomic landscape of the United Arab Emirates (UAE) mitogenome, we sequenced and analyzed the complete genomes of 232 Emirate females mitochondrial DNA (mtDNA) within and compared those to Africa. We investigated the prevalence of haplogroups, genetic variation, heteroplasmy, and demography among the UAE native population with diverse ethnicity and relatively high degree of consanguinity. We identified 968 mtDNA variants and high-resolution 15 haplogroups. Our results show that the UAE population received enough gene flow from Africa represented by the haplogroups L, U6, and M1, and that 16.8% of the population has an eastern provenance, depicted by the U haplogroup and the M Indian haplogroup (12%), whereas western Eurasian and Asian haplogroups (R, J, and K) represent 11 to 15%. Interestingly, we found an ancient migration present through the descendant of L (N1 and X) and other sub-haplogroups (L2a1d and L4) and (L3x1b), which is one of the oldest evolutionary histories outside of Africa. Our demographic analysis shows no population structure among populations, with low diversity and no population differentiation. In addition, we show that the transmission of mtDNA in the UAE population is under purifying selection with hints of diversifying selection on ATP8 gene. Last, our results show a population bottleneck, which coincides with the Western European contact (1400 ybp). Our study of the UAE mitogenomes suggest that several maternal lineage migratory episodes liking African–Asian corridors occurred since the first modern human emerges out of Africa.

The role of control region mitochondrial DNA mutations in cardiovascular disease: stroke and myocardial infarction

Recent studies associated certain type of cardiovascular disease (CVD) with specific mitochondrial DNA (mtDNA) defects, mainly driven by the central role of mitochondria in cellular metabolism. Considering the importance of the control region (CR) on the regulation of the mtDNA gene expression, the aim of the present study was to investigate the role of mtDNA CR mutations in two CVDs: stroke and myocardial infarction (MI). MtDNA CR mutations (both fixed and in heteroplasmy) were analysed in two demographically-matched case-control samples, using 154 stroke cases, 211 MI cases and their corresponding control individuals. Significant differences were found, reporting mutations m.16145 G > A and m.16311 T > C as potential genetic risk factors for stroke (conditional logistic regression: p = 0.038 and p = 0.018, respectively), whereas the m.72 T > C, m.73 A > G and m.16356 T > C mutations could act as possible beneficial genetic factors for MI (conditional logistic regression: p = 0.001, p = 0.009 and p = 0.016, respectively). Furthermore, our findings also showed a high percentage of point heteroplasmy in MI controls (logistic regression: p = 0.046; OR = 0.209, 95% CI [0.045-0.972]). These results demonstrate the possible role of mtDNA mutations in the CR on the pathogenesis of stroke and MI, and show the importance of including this regulatory region in genetic association studies.

Epigenetic Modifications Compromise Mitochondrial DNA Quality Control in the Development of Diabetic Retinopathy

Purpose

Diabetes causes dysfunction in the retinal mitochondria and increases base mismatches in their DNA (mtDNA). The enzyme responsible for repairing the base mismatches, MutL homolog 1 (Mlh1), is compromised. Diabetes also favors many epigenetic modifications and activates DNA methylation machinery, and Mlh1 has a CpG-rich promoter. Our aim is to identify the molecular mechanism responsible for impaired mtDNA mismatch repair in the pathogenesis of diabetic retinopathy.

Methods

Human retinal endothelial cells, incubated in 20 mM glucose, were analyzed for mitochondrial localization of Mlh1 by an immunofluorescence technique, Mlh1 promoter DNA methylation by the methylated DNA capture method, and the binding of Dnmt1 and transcriptional factor Sp1 by chromatin immunoprecipitation. The results were confirmed in retinal microvessels from streptozotocin-induced diabetic mice, with or without Dnmt inhibitors, and from human donors with diabetic retinopathy.

Results

Compared with cells in 5 mM glucose, high glucose decreased Mlh1 mitochondrial localization, and its promoter DNA was hypermethylated with increased Dnmt-1 binding and decreased Sp1 binding. Dnmt inhibitors attenuated Mlh1 promoter hypermethylation and prevented a decrease in its gene transcripts and an increase in mtDNA mismatches. The administration of Dnmt inhibitors in mice ameliorated a diabetes-induced increase in Mlh1 promoter hypermethylation and a decrease in its gene transcripts. Similar decreases in Mlh1 gene transcripts and its promoter DNA hypermethylation were observed in human donors.

Conclusions

Thus, as a result of the epigenetic modifications of the Mlh1 promoter, its transcription is decreased, and decreased mitochondrial accumulation fails to repair mtDNA mismatches. Therapies targeted to halt DNA methylation have the potential to prevent/halt mtDNA damage and the development of diabetic retinopathy.

Estimation of the number of contributors to mixed samples of DNA by mitochondrial DNA analyses using massively parallel sequencing

We evaluated whether the number of contributors to mixed DNA samples can be estimated by analyzing the D-loop of mitochondrial DNA using massively parallel sequencing. The A- (positions 16,209-16,400) and B- (positions 30-284) amplicons in hypervariable regions 1 and 2, respectively, were sequenced using MiSeq with 2 × 251 cycles. Sequence extraction and trimming were performed using CLC Genomics Workbench 11 and the number of observed haplotypes was counted for each amplicon type using Microsoft Excel. The haplotype ratios were calculated by dividing the number of counted reads of the corresponding haplotype by the total number of sequence reads. Haplotypes that were over the threshold (5%) were defined as positive haplotypes. The number of larger positive haplotypes in either of the two amplicon types was defined as the number of contributors. Samples were collected from seven individuals. Seventeen mixed samples were prepared by mixing DNA from two to five contributors at various ratios. The number of contributors was correctly estimated from almost all of the mixed samples containing equal amounts of DNA from two to five people. In mixed samples of two or three people, the minor components were detected down to a ratio of 20:1 or 8:2:1. However, heteroplasmy, base deletions, and sharing of the same haplotypes caused incorrect estimations of the number of contributors. Although this method still has room for improvement, it may be useful for estimating the number of contributors in a mixed sample, as it does not rely on forensic mathematics.

Exploration of mtDNA control region sequences in Chinese Tibetan Mastiffs

The control region of mitochondrial DNA (mtDNA) was obtained from 40 purebred Chinese Tibetan Mastiffs (TMs). Sequence structure and genetic diversity were analyzed, and a phylogenetic tree was constructed. The TM mtDNA control region was composed of ETAS (extended termination associated sequences), CD (a central domain) and CSBs (conserved sequenced blocks) and sequence length showed some diversity, which was mainly caused by the number of 10 nucleotide repeat units [5'-GTA CAC GT (G/A) C-3'] between CSB I and CSB II, which ranged from 27 to 35 among individuals. Seventy-five polymorphic sites were identified, which defined 37 haplotypes; the haplotype diversity was 0.990, and the nucleotide diversity was 1.201. Based on the control region sequences, Chinese TMs were divided into three categories, which were consistent with the origin and geographical classification of TMs. Phylogenetic analysis of 538-bp HVR-I sequences revealed that TMs were most closely related to Labrador Retrievers.

Phylogenetic analysis of Tibetan mastiffs based on mitochondrial hypervariable region I

Recently, the number of Tibetan mastiffs, which is a precious germplasm resource and cultural heritage, is decreasing sharply. Therefore, the genetic diversity of Tibetan mastiffs needs to be studied to clarify its phylogenetics relationships and lay the foundation for resource protection, rational development and utilization of Tibetan mastiffs. We sequenced hypervariable region I of mitochondrial DNA (mtDNA) of 110 individuals from Tibet region and Gansu province. A total of 12 polymorphic sites were identified which defined eight haplotypes of which H4 and H8 were unique to Tibetan population with H8 being identified first. The haplotype diversity (Hd: 0.808), nucleotide diversity (Pi: 0.603%), the average number of nucleotide difference (K: 3.917) of Tibetan mastiffs from Gansu were higher than those from Tibet region (Hd: 0.794; Pi: 0.589%; K: 3.831), which revealed higher genetic diversity in Gansu. In terms of total population, the genetic variation was low. The median-joining network and phylogenetic tree based on the mtDNA hypervariable region I showed that Tibetan mastiffs originated from grey wolves, as the other domestic dogs and had different history of maternal origin. The mismatch distribution analysis and neutrality tests indicated that Tibetan mastiffs were in genetic equilibrium or in a population decline.

Length heteroplasmy of the polyC-polyT-polyC stretch in the dog mtDNA control region

Previously, the mitochondrial control region of 214 Belgian dogs was sequenced. Analysis of this data indicated length heteroplasmy of the polyT stretch in the polyC-polyT-polyC stretch from positions 16661 to 16674. Nine polyC-polyT-polyC haplotype combinations were observed, consisting of seven major haplotypes (highest signal intensity) combined with minor haplotypes (lower signal intensity) one T shorter than the major haplotype in all but three dogs. The longer the polyT stretch, the smaller was the difference in signal intensity between the major and minor haplotype peaks. Additional sequencing, cloning, and PCR trap experiments were performed to further study the intra-individual variation of this mitochondrial DNA (mtDNA) region. Cloning experiments demonstrated that the proportion of clones displaying the minor haplotypes also increased with the length of the polyT stretch. Clone amplification showed that in vitro polymerase errors might contribute to the length heteroplasmy of polyT stretches with at least 10 Ts. Although major and minor polyC-polyT-polyC haplotypes did not differ intra-individually within and between tissues in this study, interpretation of polyT stretch variation should be handled with care in forensic casework.

Primer effect in the detection of mtDNA heteroplasmy: insights from horse Cytochrome b gene

Heteroplasmy, the presence of bi-allelic mtDNA types within an individual, has been previously detected in the D-loop region and Cytochrome b gene (Cytb) of mitochondrial DNA using PCR-RFLP. However, heteroplasmy was absent in thousands of equine mtDNA sequences deposited in GenBank. To address whether heteroplasmy widely exists in mitochondria of Chinese indigenous horses, we generated the data set of the target sites in Cytb region with Sanger sequencing and PCR-RFLP method as well. In this study, 23 heteroplasmic individuals were detected in 430 Chinese local horses. Both site and length heteroplasmy were identified in horse Cytb, especially in Xinihe and Ujumqin breeds. Our data provide evidence that the forward and reverse primers seem to produce a similar approximation to the proportion of mutation base call. However, locations of primers affected the proper detection of mtDNA heteroplasmy. The data obtained in this study highlight the importance of the primers in the accurate detection of heteroplasmy.

Frequency and pattern of heteroplasmy in the complete human mitochondrial genome

Determining the levels of human mitochondrial heteroplasmy is of utmost importance in several fields. In spite of this, there are currently few published works that have focused on this issue. In order to increase the knowledge of mitochondrial DNA (mtDNA) heteroplasmy, the main goal of this work is to investigate the frequency and the mutational spectrum of heteroplasmy in the human mtDNA genome. To address this, a set of nine primer pairs designed to avoid co-amplification of nuclear DNA (nDNA) sequences of mitochondrial origin (NUMTs) was used to amplify the mitochondrial genome in 101 individuals. The analysed individuals represent a collection with a balanced representation of genders and mtDNA haplogroup distribution, similar to that of a Western European population. The results show that the frequency of heteroplasmic individuals exceeds 61%. The frequency of point heteroplasmy is 28.7%, with a widespread distribution across the entire mtDNA. In addition, an excess of transitions in heteroplasmy were detected, suggesting that genetic drift and/or selection may be acting to reduce its frequency at population level. In fact, heteroplasmy at highly stable positions might have a greater impact on the viability of mitochondria, suggesting that purifying selection must be operating to prevent their fixation within individuals. This study analyses the frequency of heteroplasmy in a healthy population, carrying out an evolutionary analysis of the detected changes and providing a new perspective with important consequences in medical, evolutionary and forensic fields.

Heteroplasmy in the mitochondrial genomes of human lice and ticks revealed by high throughput sequencing

The typical mitochondrial (mt) genomes of bilateral animals consist of 37 genes on a single circular chromosome. The mt genomes of the human body louse, Pediculus humanus, and the human head louse, Pediculus capitis, however, are extensively fragmented and contain 20 minichromosomes, with one to three genes on each minichromosome. Heteroplasmy, i.e. nucleotide polymorphisms in the mt genome within individuals, has been shown to be significantly higher in the mt cox1 gene of human lice than in humans and other animals that have the typical mt genomes. To understand whether the extent of heteroplasmy in human lice is associated with mt genome fragmentation, we sequenced the entire coding regions of all of the mt minichromosomes of six human body lice and six human head lice from Ethiopia, China and France with an Illumina HiSeq platform. For comparison, we also sequenced the entire coding regions of the mt genomes of seven species of ticks, which have the typical mitochondrial genome organization of bilateral animals. We found that the level of heteroplasmy varies significantly both among the human lice and among the ticks. The human lice from Ethiopia have significantly higher level of heteroplasmy than those from China and France (Pt<0.05). The tick, Amblyomma cajennense, has significantly higher level of heteroplasmy than other ticks (Pt<0.05). Our results indicate that heteroplasmy level can be substantially variable within a species and among closely related species, and does not appear to be determined by single factors such as genome fragmentation.

Primer effect in the detection of mitochondrial DNA point heteroplasmy by automated sequencing

The correct detection of mitochondrial DNA (mtDNA) heteroplasmy by automated sequencing presents methodological constraints. The main goals of this study are to investigate the effect of sense and distance of primers in heteroplasmy detection and to test if there are differences in the accurate determination of heteroplasmy involving transitions or transversions. A gradient of the heteroplasmy levels was generated for mtDNA positions 9477 (transition G/A) and 15,452 (transversion C/A). Amplification and subsequent sequencing with forward and reverse primers, situated at 550 and 150 bp from the heteroplasmic positions, were performed. Our data provide evidence that there is a significant difference between the use of forward and reverse primers. The forward primer is the primer that seems to give a better approximation to the real proportion of the variants. No significant differences were found concerning the distance at which the sequencing primers were placed neither between the analysis of transitions and transversions. The data collected in this study are a starting point that allows to glimpse the importance of the sequencing primers in the accurate detection of point heteroplasmy, providing additional insight into the overall automated sequencing strategy.

Site heteroplasmy in the mitochondrial cytochrome b gene of the sterlet sturgeon Acipenser ruthenus

Sturgeons are fish species with a complex biology. They are also characterized by complex aspects including polyploidization and easiness of hybridization. As with most of the Ponto-Caspian sturgeons, the populations of Acipenser ruthenus from the Danube have declined drastically during the last decades. This is the first report on mitochondrial point heteroplasmy in the cytochrome b gene of this species. The 1141 bp sequence of the cytb gene in wild sterlet sturgeon individuals from the Lower Danube was determined, and site heteroplasmy evidenced in three of the 30 specimens collected. Two nucleotide sequences were identified in these heteroplasmic individuals. The majority of the heteroplasmic sites are synonymous and do not modify the sequence of amino acids in cytochrome B protein. To date, several cases of point heteroplasmy have been reported in animals, mostly due to paternal leakage of mtDNA. The presence of specific point heteroplasmic sites might be interesting for a possible correlation with genetically distinct groups in the Danube River.

A twin study of mitochondrial DNA polymorphisms shows that heteroplasmy at multiple sites is associated with mtDNA variant 16093 but not with zygosity

The mitochondrial theory of ageing proposes that damage to mitochondria and diminished mitochondrial DNA (mtDNA) repair are major contributors to cellular dysfunction and age-related diseases. We investigate the prevalence of heteroplasmy in the mtDNA control region in buccal swab and blood derived samples for 178 women from the TwinsUK cohort (41 DZ pair 39 MZ pairs, 18 singletons, mean age 57.5 range 28–82) and its relationship to age, BMI and fasting insulin and glucose serum levels. The overall estimated prevalence of heteroplasmy for both tissues in the control region measured for 37 sites was 17%. The prevalence of heteroplasmy was higher among the older half of the study subjects than in the younger half (23% vs 10% p<0.03), primarily reflecting the increase in the prevalence of a heteroplasmic dinucleotide CA repeat in variable region II (VRII) with age. The VRII 523–524 heteroplasmic site (heteroplasmic in 25 subjects) was also associated with a decrease in BMI. In addition, concordance rates for common heteroplasmy were observed to be near complete for both dizygotic (DZ = 94%) and monozygotic twin pairs (MZ = 100%), consistent with previous reports that suggest variation in heteroplasmy rates between generations are determined by bottlenecks in maternal transmission of mitochondria. Differences in the prevalence of heteroplasmy were observed overall between samples derived from buccal swabs (19%) and blood (15%, p<0.04). These were particularly marked at position 16093 of hypervariable region I (HVI, 7% vs 0%, respectively, p<4×10⁻¹¹). The presence of the C allele at position 16093 in blood was associated with the presence of heteroplasmy in buccal swabs at this position (p = 3.5×10⁻¹⁴) and also at VRII (p = 2×10⁻⁴) suggesting a possible predisposing role for this site in the accumulation of heteroplasmy. Our data indicate that BMI is potentially associated with control region heteroplasmy.

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.

The link between mitochondrial DNA hypervariable segment I heteroplasmy and ageing among genetically unrelated Latvians

Various studies have demonstrated that mitochondrial DNA (mtDNA) heteroplasmy tends to increase with age and that the observed frequency of heteroplasmy among populations mostly depends on the way it is measured. Therefore, we investigated age-related association on the presence of mtDNA heteroplasmy within the hypervariable segment 1 (HVS-I) in a selected study group. The study group consisted of 300 maternally unrelated Latvians ranging in age from 18 to over 90 years. To determine the optimal method for mtDNA heteroplasmy detection, three approaches were used: (i) SURVEYOR Mutation Detection Kit, (ii) sequencing and (iii) denaturing gradient-gel electrophoresis (DGGE). Among the studied individuals, 30.3% were found to be heteroplasmic. The distribution of heteroplasmy statistically significantly increased with individuals' age (17%; 95% confidence interval [CI] 0.095-0.244 in the 18-40 year age group vs. 39%; [CI] 0.294-0.487 in the >90 year age group). Heteroplasmy occurred in a total of 21 different positions within HVS-I, and was the most frequent at fast-mutated positions 16189, 16304 and 16311. The results indicate that heteroplasmy in HVS-I is relatively common and occurs in a broad spectrum of sites. The above is supported by evidence to eventual increase of the probability of heteroplasmy with age due to specific mitochondrial haplogroup background.

Evaluating sequence-derived mtDNA length heteroplasmy by amplicon size analysis

Length heteroplasmy (LH) in mitochondrial (mt)DNA is usually observed in homopolymeric tracts and manifest as mixture of various length variants. The generally used difference-coded annotation to report mtDNA haplotypes does not express the degree of LH variation present in a sample, even more so, it is sometimes difficult to establish which length variants are present and clearly distinguishable from background noise. It has therefore become routine practice for some researchers to call the dominant type, the "major molecule", which represents the LH variant that is most abundant in a DNA extract. In the majority of cases a clear single dominant variant can be identified. However, in some samples this interpretation is difficult, i.e. when (almost) equally quantitative LH variants are present or when multiple sequencing primers result in the presentation of different dominant types. To better understand those cases we designed amplicon sizing assays for the five most relevant LH regions in the mtDNA control region (around ntps 16,189, 310, 460, 573, and the AC-repeat between 514 and 524) to determine the ratio of the LH variants by fluorescence based amplicon sizing assays. For difficult LH constellations derived by Sanger sequencing (with Big Dye terminators) these assays mostly gave clear and unambiguous results. In the vast majority of cases we found agreement between the results of the sequence and amplicon analyses and propose this alternative method in difficult cases.

Segregation of point mutation heteroplasmy in the control region of dog mtDNA studied systematically in deep generation pedigrees

Heteroplasmy, the presence of two or more variants in an organism, may render mitochondrial DNA (mtDNA)-based individual identification challenging in forensic analysis. However, the variation of heteroplasmic proportions and the segregation of heteroplasmic variants through generations and within families have not been systematically described at a large scale in animals such as the domestic dog. Therefore, we performed the largest study to date in domestic dogs and screened a 582-bp-long fragment of the mtDNA control region in 180 individuals in 58 pedigrees for signs of heteroplasmy. We identified three pedigrees (5.17%) with heteroplasmic point mutations. To follow the segregation of the point mutations, we then analyzed 131 samples from these three independent pedigrees and found significant differences in heteroplasmy between generations and among siblings. Frequently (10% of cases), the proportion of one base changed from 0–10% to 80–90% (as judged from Sanger electropherograms) between generations and varied to a similar extent among siblings. We included also a literature review of heteroplasmic and potential mutational hot spot positions in the studied region which showed that all heteroplasmic positions appear to be mutational hot spots. Thus, although heteroplasmy may be used to increase the significance of a match in forensic case work, it may also cause erroneous exclusion of related individuals because of sharp switches from one state to the other within a single generation or among siblings especially in the presented mutational hot spots.

Mitochondrial DNA patterns in the Macaronesia islands: Variation within and among archipelagos

Macaronesia covers four Atlantic archipelagos: the Azores, Madeira, the Canary Islands, and the Cape Verde islands. When discovered by Europeans in the 15th century, only the Canaries were inhabited. Historical reports highlight the impact of Iberians on settlement in Macaronesia. Although important differences in their settlement are documented, its influence on their genetic structures and relationships has yet to be ascertained. In this study, the hypervariable region I (HVRI) sequence and coding region polymorphisms of mitochondrial DNA (mtDNA) in 623 individuals from the Azores (120) and Canary Islands (503) were analyzed. Combined with published data, these give a total of 1,542 haplotypes from Macaronesia and 1,067 from the Iberian Peninsula. The results obtained indicate that Cape Verde is the most distinctive archipelago, with an mtDNA pool composed almost exclusively of African lineages. However, the other archipelagos present an mtDNA profile dominated by the presence of West-Eurasian mtDNA haplogroups with African lineages present in varying proportions. Moreover, no signs of integration of typical Canarian U6 lineages in the other archipelagos were detected. The four Macaronesia archipelagos currently have differentiated genetic profiles, and the Azores present the highest intra-archipelago differentiation and the lowest values of diversity. The analyses performed show that the present-day genetic profile of the Macaronesian archipelagos was mainly determined by the initial process of settlement and further microdifferentiation probably as a consequence of the small population size of some islands. Moreover, contacts between archipelagos seem to have had a low impact on the mtDNA genetic pool of each archipelago.

A western Eurasian male is found in 2000-year-old elite Xiongnu cemetery in Northeast Mongolia

We analyzed mitochondrial DNA (mtDNA), Y-chromosome single nucleotide polymorphisms (Y-SNP), and autosomal short tandem repeats (STR) of three skeletons found in a 2,000-year-old Xiongnu elite cemetery in Duurlig Nars of Northeast Mongolia. This study is one of the first reports of the detailed genetic analysis of ancient human remains using the three types of genetic markers. The DNA analyses revealed that one subject was an ancient male skeleton with maternal U2e1 and paternal R1a1 haplogroups. This is the first genetic evidence that a male of distinctive Indo-European lineages (R1a1) was present in the Xiongnu of Mongolia. This might indicate an Indo-European migration into Northeast Asia 2,000 years ago. Other specimens are a female with mtDNA haplogroup D4 and a male with Y-SNP haplogroup C3 and mtDNA haplogroup D4. Those haplogroups are common in Northeast Asia. There was no close kinship among them. The genetic evidence of U2e1 and R1a1 may help to clarify the migration patterns of Indo-Europeans and ancient East-West contacts of the Xiongnu Empire. Artifacts in the tombs suggested that the Xiongnu had a system of the social stratification. The West Eurasian male might show the racial tolerance of the Xiongnu Empire and some insight into the Xiongnu society.

Detecting heteroplasmy from high-throughput sequencing of complete human mitochondrial DNA genomes

Heteroplasmy, the existence of multiple mtDNA types within an individual, has been previously detected by using mostly indirect methods and focusing largely on just the hypervariable segments of the control region. Next-generation sequencing technologies should enable studies of heteroplasmy across the entire mtDNA genome at much higher resolution, because many independent reads are generated for each position. However, the higher error rate associated with these technologies must be taken into consideration to avoid false detection of heteroplasmy. We used simulations and phiX174 sequence data to design criteria for accurate detection of heteroplasmy with the Illumina Genome Analyzer platform, and we used artificial mixtures and replicate data to test and refine the criteria. We then applied these criteria to mtDNA sequence reads for 131 individuals from five Eurasian populations that had been generated via a parallel tagged approach. We identified 37 heteroplasmies at 10% frequency or higher at 34 sites in 32 individuals. The mutational spectrum does not differ between heteroplasmic mutations and polymorphisms in the same individuals, but the relative mutation rate at heteroplasmic mutations is significantly higher than that estimated for all mutable sites in the human mtDNA genome. Moreover, there is also a significant excess of nonsynonymous mutations observed among heteroplasmies, compared to polymorphism data from the same individuals. Both mutation-drift and negative selection influence the fate of heteroplasmies to determine the polymorphism spectrum in humans. With appropriate criteria for avoiding false positives due to sequencing errors, next-generation technologies can provide novel insights into genome-wide aspects of mtDNA heteroplasmy.

Evaluating length heteroplasmy in the human mitochondrial DNA control region

We present allelic data for three known and one new C-tract in the human mitochondrial DNA (mtDNA) control region, and we measure intergenerational mutation rates at such C-tracts. In detail, in a sample of 1,172 mtDNA sequences, we demonstrate the existence of an instability threshold of eight consecutive cytosines, at and above which the phenomenon of length heteroplasmy arises. To determine mutation rates, we draw on mtDNA sequences in up to four generations of 248 pedigrees for families living in high or low-radiation environmental conditions. The high-radiation sample gives the most conservative (fastest) mutation rate likely to be encountered in any forensic context. We find that the C-tract mutation rate is up to 6% per generation, and we observe an excess of cytosine gains over losses. Case studies and guidelines for evaluating mtDNA heteroplasmy are provided.

Heteroplasmic mitochondrial DNA mutations in normal and tumour cells

The presence of hundreds of copies of mitochondrial (mt) DNA in each human cell poses a challenge for complete characterization of mtDNA genomes by conventional sequencing technologies¹. Here, we describe digital sequencing of mtDNA genomes using massively parallel sequencing-by-synthesis. Though the mtDNA of human cells is considered to be homogeneous, we found widespread heterogeneity (heteroplasmy) in the mtDNA of normal human cells. Moreover, the frequency of heteroplasmic variants among different tissues of the same individual varied considerably. In addition to the variants identified in normal tissues, cancer cells harbored additional homoplasmic and heteroplasmic mutations that could also be detected in patient plasma. These studies provide new insights into the nature and variability of mtDNA sequences and have intriguing implications for mitochondrial processes during embryogenesis, cancer biomarker development, and forensic analysis. In particular, they demonstrate that individual humans are characterized by a complex mixture of related mitochondrial genotypes rather than a single genotype.

Fluorescent duplex allele-specific PCR and amplicon melting for rapid homogeneous mtDNA haplogroup H screening and sensitive mixture detection

BACKGROUND

For large scale studies aiming at a better understanding of mitochondrial DNA (mtDNA), sequence variation in particular mt haplogroups (hgs) and population structure, reliable low-cost high-throughput genotyping assays are needed. Furthermore, methods facilitating sensitive mixture detection and relative quantification of allele proportions are indispensable for the study of heteroplasmy, mitochondrial sequence evolution, and mitochondrial disorders. Here the properties of a homogeneous competitive duplex allele specific PCR (ARMS) assay were scrutinized in the light of these requirements.

METHODOLOGY/PRINCIPAL FINDINGS

A duplex ARMS assay amplifying either the ancestral mtDNA 2706G allele (non-hg H samples) or the derived 7028C allele (hg H samples) in the presence of SYBR Green fluorescent reporter dye was developed and characterized. Product detection, allele calling, and hg inference were based on the amplicon-characteristic melting-point temperatures obtained with on-line post-PCR fluorescent dissociation curve analysis (DCA). The analytical window of the assay covered at least 5 orders of magnitude of template DNA input with a detection limit in the low picogram range of genomic DNA. A set of forensically relevant test specimens was analyzed successfully. The presence of mtDNA mixtures was detected over a broad range of input DNA amounts and mixture ratios, and the estimation of allele proportions in samples with known total mtDNA content was feasible with limitations. A qualified DNA analyst successfully analyzed approximately 2,200 DNA extracts within three regular working days, without using robotic lab-equipment. By performing the amplification on-line, the assay also facilitated absolute mtDNA quantification.

CONCLUSIONS

Although this assay was developed just for a particular purpose, the approach is general in that it is potentially suitable in a broad variety of assay-layouts for many other applications, including the analysis of mixtures. Homogeneous ARMS-DCA is a valuable tool for large-volume studies targeting small numbers of single nucleotide polymorphisms (SNPs).

Human mitochondrial DNA complete amplification and sequencing: a new validated primer set that prevents nuclear DNA sequences of mitochondrial origin co-amplification

To date, there are no published primers to amplify the entire mitochondrial DNA (mtDNA) that completely prevent the amplification of nuclear DNA (nDNA) sequences of mitochondrial origin. The main goal of this work was to design, validate and describe a set of primers, to specifically amplify and sequence the complete human mtDNA, allowing the correct interpretation of mtDNA heteroplasmy in healthy and pathological samples. Validation was performed using two different approaches: (i) Basic Local Alignment Search Tool and (ii) amplification using isolated nDNA obtained from sperm cells by differential lyses. During the validation process, two mtDNA regions, with high similarity with nDNA, represent the major problematic areas for primer design. One of these could represent a non-published nuclear DNA sequence of mitochondrial origin. For two of the initially designed fragments, the amplification results reveal PCR artifacts that can be attributed to the poor quality of the DNA. After the validation, nine overlapping primer pairs to perform mtDNA amplification and 22 additional internal primers for mtDNA sequencing were obtained. These primers could be a useful tool in future projects that deal with mtDNA complete sequencing and heteroplasmy detection, since they represent a set of primers that have been tested for the non-amplification of nDNA.

Investigation of heteroplasmy in the human mitochondrial DNA control region: a synthesis of observations from more than 5000 global population samples

Instances of point and length heteroplasmy in the mitochondrial DNA control region were compiled and analyzed from over 5,000 global human population samples. These data represent observations from a large and broad population sample, representing nearly 20 global populations. As expected, length heteroplasmy was frequently observed in the HVI, HVII and HVIII C-stretches. Length heteroplasmy was also observed in the AC dinucleotide repeat region, as well as other locations. Point heteroplasmy was detected in approximately 6% of all samples, and while the vast majority of heteroplasmic samples comprised two molecules differing at a single position, samples exhibiting two and three mixed positions were also observed in this data set. In general, the sites at which heteroplasmy was most commonly observed correlated with reported control region mutational hotspots. However, for some sites, observations of heteroplasmy did not mirror established mutation rate data, suggesting the action of other mechanisms, both selective and neutral. Interestingly, these data indicate that the frequency of heteroplasmy differs between particular populations, perhaps reflecting variable mutation rates among different mtDNA lineages and/or artifacts of particular population groups. The results presented here contribute to our general understanding of mitochondrial DNA control region heteroplasmy and provide additional empirical information on the mechanisms contributing to mtDNA control region mutation and evolution.