Single-Molecule LATE-PCR Analysis of Human Mitochondrial Genomic Sequence Variations

It is thought that changes in mitochondrial DNA are associated with many degenerative diseases, including Alzheimer''s and diabetes. Much of the evidence, however, depends on correlating disease states with changing levels of heteroplasmy within populations of mitochondrial genomes, rather than individual mitochondrial genomes. Thus these measurements are likely to either overestimate the extent of heteroplasmy due to technical artifacts, or underestimate the actual level of heteroplasmy because only the most abundant changes are observable. In contrast, Single Molecule (SM) LATE-PCR analysis achieves efficient amplification of single-stranded amplicons from single target molecules. The product molecules, in turn, can be accurately sequenced using a convenient Dilute-‘N’-Go protocol, as shown here. Using these novel technologies we have rigorously analyzed levels of mitochondrial genome heteroplasmy found in single hair shafts of healthy adult individuals. Two of the single molecule sequences (7% of the samples) were found to contain mutations. Most of the mtDNA sequence changes, however, were due to the presence of laboratory contaminants. Amplification and sequencing errors did not result in mis-identification of mutations. We conclude that SM-LATE-PCR in combination with Dilute-‘N’-Go Sequencing are convenient technologies for detecting infrequent mutations in mitochondrial genomes, provided great care is taken to control and document contamination. We plan to use these technologies in the future to look for age, drug, and disease related mitochondrial genome changes in model systems and clinical samples.     

Molecular characterization of mitochondrial Amerindian haplogroups and the amelogenin gene in human ancient DNA from three archaeological sites in Lambayeque - Peru

Important pre-Inca civilizations, known by their great political and religious structures, inhabited the northern coast of Peru. Archeological and anthropological studies have shown that people from these villages have hierarchical strata, but the genetic structure has been poorly studied. Here, we aimed to perform a molecular characterization of the Amerindian maternal lineages and the amelogenin gene in skeletons collected from three archeological sites in Lambayeque. Ancient DNA (aDNA) samples were analyzed with conventional PCR to assess the nine-base pair (9 bp) deletion corresponding to mitochondrial haplogroup B and the identification of haplogroups A, C, and D were obtained with PCR-RFLP experiments. The sex was characterized via amplification of the AMEL(X/Y) locus. Haplogroup frequencies were compared with available data from other ancient and modern civilizations from the Peruvian coast and highlands using statistical methods. Our results showed that haplogroup C had the highest frequency, while haplogroup B showed variable diversity in the analyzed populations. The meta-analysis revealed a positive correlation among some coastal villages. We concluded that ancient populations analyzed in our study showed the presence of four Amerindian mitochondrial haplogroups, which is consistent with previous studies.     

A Complete Mitochondrial Genome Sequence from a Mesolithic Wild Aurochs (Bos primigenius)

Background

The derivation of domestic cattle from the extinct wild aurochs (Bos primigenius) has been well-documented by archaeological and genetic studies. Genetic studies point towards the Neolithic Near East as the centre of origin for Bos taurus, with some lines of evidence suggesting possible, albeit rare, genetic contributions from locally domesticated wild aurochsen across Eurasia. Inferences from these investigations have been based largely on the analysis of partial mitochondrial DNA sequences generated from modern animals, with limited sequence data from ancient aurochsen samples. Recent developments in DNA sequencing technologies, however, are affording new opportunities for the examination of genetic material retrieved from extinct species, providing new insight into their evolutionary history. Here we present DNA sequence analysis of the first complete mitochondrial genome (16,338 base pairs) from an archaeologically-verified and exceptionally-well preserved aurochs bone sample.

Methodology

DNA extracts were generated from an aurochs humerus bone sample recovered from a cave site located in Derbyshire, England and radiocarbon-dated to 6,738±68 calibrated years before present. These extracts were prepared for both Sanger and next generation DNA sequencing technologies (Illumina Genome Analyzer). In total, 289.9 megabases (22.48%) of the post-filtered DNA sequences generated using the Illumina Genome Analyzer from this sample mapped with confidence to the bovine genome. A consensus B. primigenius mitochondrial genome sequence was constructed and was analysed alongside all available complete bovine mitochondrial genome sequences.

Conclusions

For all nucleotide positions where both Sanger and Illumina Genome Analyzer sequencing methods gave high-confidence calls, no discrepancies were observed. Sequence analysis reveals evidence of heteroplasmy in this sample and places this mitochondrial genome sequence securely within a previously identified aurochsen haplogroup (haplogroup P), thus providing novel insights into pre-domestic patterns of variation. The high proportion of authentic, endogenous aurochs DNA preserved in this sample bodes well for future efforts to determine the complete genome sequence of a wild ancestor of domestic cattle.     

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.     

Simultaneous A8344G heteroplasmy and mitochondrial DNA copy number quantification in myoclonus epilepsy and ragged-red fibers (MERRF) syndrome by a multiplex molecular beacon based real-time fluorescence PCR

The association of a particular mitochondrial DNA (mtDNA) mutation with different clinical phenotypes is a well-known feature of mitochondrial diseases. A simple genotype–phenotype correlation has not been found between mutation load and disease expression. Tissue and intercellular mosaicism as well as mtDNA copy number are thought to be responsible for the different clinical phenotypes. As disease expression of mitochondrial tRNA mutations is mostly in postmitotic tissues, studies to elucidate disease mechanisms need to be performed on patient material. Heteroplasmy quantitation and copy number estimation using small patient biopsy samples has not been reported before, mainly due to technical restrictions. In order to resolve this problem, we have developed a robust assay that utilizes Molecular Beacons to accurately quantify heteroplasmy levels and determine mtDNA copy number in small samples carrying the A8344G tRNA^Lys mutation. It provides the methodological basis to investigate the role of heteroplasmy and mtDNA copy number in determining the clinical phenotypes.     

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. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Evidence for multi-copy Mega-NUMTs in the human genome

The maternal mode of mitochondrial DNA (mtDNA) inheritance is central to human genetics. Recently, evidence for bi-parental inheritance of mtDNA was claimed for individuals of three pedigrees that suffered mitochondrial disorders. We sequenced mtDNA using both direct Sanger and Massively Parallel Sequencing in several tissues of eleven maternally related and other affiliated healthy individuals of a family pedigree and observed mixed mitotypes in eight individuals. Cells without nuclear DNA, i.e. thrombocytes and hair shafts, only showed the mitotype of haplogroup (hg) V. Skin biopsies were prepared to generate ρ° cells void of mtDNA, sequencing of which resulted in a hg U4c1 mitotype. The position of the Mega-NUMT sequence was determined by fluorescence in situ hybridization and two different quantitative PCR assays were used to determine the number of contributing mtDNA copies. Thus, evidence for the presence of repetitive, full mitogenome Mega-NUMTs matching haplogroup U4c1 in various tissues of eight maternally related individuals was provided. Multi-copy Mega-NUMTs mimic mixtures of mtDNA that cannot be experimentally avoided and thus may appear in diverse fields of mtDNA research and diagnostics. We demonstrate that hair shaft mtDNA sequencing provides a simple but reliable approach to exclude NUMTs as source of misleading results.     

Mitochondrial Haplogroup T Is Associated with Obesity in Austrian Juveniles and Adults

Background

Recent publications have reported contradictory data regarding mitochondrial DNA (mtDNA) variation and its association with body mass index. The aim of the present study was to compare the frequencies of mtDNA haplogroups as well as control region (CR) polymorphisms of obese juveniles (n = 248) and obese adults (n = 1003) versus normal weight controls (n_juvenile = 266, n_adults = 595) in a well-defined, ethnically homogenous, age-matched comparative cohort of Austrian Caucasians.

Methodology and Principal Findings

Using SNP analysis and DNA sequencing, we identified the nine major European mitochondrial haplogroups and CR polymorphisms. Of these, only the T haplogroup frequency was increased in the juvenile obese cohort versus the control subjects [11.7% in obese vs. 6.4% in controls], although statistical significance was lost after adjustment for sex and age. Similar data were observed in a local adult cohort, in which haplogroup T was found at a significantly higher frequency in the overweight and obese subjects than in the normal weight group [9.7% vs. 6.2%, p = 0.012, adjusted for sex and age]. When all obese subjects were considered together, the difference in the frequency of haplogroup T was even more clearly seen [10.1% vs. 6.3%, p = 0.002, OR (95% CI) 1.71 (1.2–2.4), adjusted for sex and age]. The frequencies of the T haplogroup-linked CR polymorphisms C16294T and the C16296T were found to be elevated in both the juvenile and the adult obese cohort compared to the controls. Nevertheless, no mtDNA haplogroup or CR polymorphism was robustly associated with any of several investigated metabolic and cardiovascular parameters (e.g., blood pressure, blood glucose concentration, triglycerides, cholesterol) in all obese subjects.

Conclusions and Significance

By investigation of this large ethnically and geographically homogenous cohort of Middle European Caucasians, only mtDNA haplogroup T was identified as an obesity risk factor.     

Influence of cell distribution and diabetes status on the association between mitochondrial DNA copy number and aging phenotypes in the InCHIANTI study

Mitochondrial DNA copy number (mtDNA‐CN) estimated in whole blood is a novel marker of mitochondrial mass and function that can be used in large population‐based studies. Analyses that attempt to relate mtDNA‐CN to specific aging phenotypes may be confounded by differences in the distribution of blood cell types across samples. Also, low or high mtDNA‐CN may have a different meaning given the presence of diseases associated with mitochondrial damage. We evaluated the impact of blood cell type distribution and diabetes status on the association between mtDNA‐CN and aging phenotypes, namely chronologic age, interleukin‐6, hemoglobin, and all‐cause mortality, among 672 participants of the InCHIANTI study. After accounting for white blood cell count, platelet count, and white blood cell proportions in multivariate models, associations of mtDNA‐CN with age and interleukin‐6 were no longer statistically significant. Evaluation of a statistical interaction by diabetes status suggested heterogeneity of effects in the analysis of mortality (P < 0.01). The magnitude and direction of associations between mtDNA‐CN estimated from blood samples and aging phenotypes are influenced by the sample cell type distribution and disease status. Therefore, accounting for these factors may aid understanding of the relevance of mitochondrial DNA copy number to health and aging.     

Heteroplasmy analysis in the Polish patients with 11778A mutation responsible for Leber hereditary optic neuropathy

We have analysed the heteroplasmy level in 11 individuals from 3 families harbouring the mitochondrial 11778A mutation responsible for Leber hereditary optic neuropathy using last cycle hot PCR. The mutation level exceeded 90% both in affected and in unaffected individuals. We also checked whether any of the families belonged to the J haplogroup of mitochondrial DNA and obtained a negative result.     

Demographic history of Canary Islands male gene-pool: replacement of native lineages by European

Background  

The origin and prevalence of the prehispanic settlers of the Canary Islands has attracted great multidisciplinary interest. However, direct ancient DNA genetic studies on indigenous and historical 17^th–18^th century remains, using mitochondrial DNA as a female marker, have only recently been possible. In the present work, the analysis of Y-chromosome polymorphisms in the same samples, has shed light on the way the European colonization affected male and female Canary Island indigenous genetic pools, from the conquest to present-day times.     

An economical mtDNA SNP assay detecting different mitochondrial haplogroups in identical HVR 1 samples of Caucasian ancestry

BackgroundWe had sequenced 329 Caucasian samples in Hypervariable Region 1 (HVR 1) and found that they belong to eleven different mitochondrial DNA (mtDNA) haplotypes. The sample set was further analysed by an mtDNA assay examining 32 single nucleotide polymorphisms (SNPs) for haplogroup discrimination.In a validation study on 160 samples of different origin it was shown that these SNPs were able to discriminate between the evolved superhaplogroups worldwide (L, M and N) and between the nine most common Caucasian haplogroups (H, I, J, K, T, U, V, W and X).ResultsThe 32 mtDNA SNPs comprised 42 different SNP haplotypes instead of only eleven haplotypes after HVR 1 sequencing. The assay provided stable results in a range of 5 ng genomic DNA down to virtually no genomic DNA per reaction. It was possible to detect samples of African, Asian and Eurasian ancestry, respectively.DiscussionThe 32 mtDNA SNP assay is a helpful adjunct to further distinguish between identical HVR 1 sequences of Caucasian origin. Our results suggest that haplogroup prediction using HVR 1 sequencing provides instable results. The use of coding region SNPs for haplogroup assignment is more suited than using HVR 1 haplotypes.     

Assessing Mitochondrial DNA Variation and Copy Number in Lymphocytes of ~2,000 Sardinians Using Tailored Sequencing Analysis Tools

DNA sequencing identifies common and rare genetic variants for association studies, but studies typically focus on variants in nuclear DNA and ignore the mitochondrial genome. In fact, analyzing variants in mitochondrial DNA (mtDNA) sequences presents special problems, which we resolve here with a general solution for the analysis of mtDNA in next-generation sequencing studies. The new program package comprises 1) an algorithm designed to identify mtDNA variants (i.e., homoplasmies and heteroplasmies), incorporating sequencing error rates at each base in a likelihood calculation and allowing allele fractions at a variant site to differ across individuals; and 2) an estimation of mtDNA copy number in a cell directly from whole-genome sequencing data. We also apply the methods to DNA sequence from lymphocytes of ~2,000 SardiNIA Project participants. As expected, mothers and offspring share all homoplasmies but a lesser proportion of heteroplasmies. Both homoplasmies and heteroplasmies show 5-fold higher transition/transversion ratios than variants in nuclear DNA. Also, heteroplasmy increases with age, though on average only ~1 heteroplasmy reaches the 4% level between ages 20 and 90. In addition, we find that mtDNA copy number averages ~110 copies/lymphocyte and is ~54% heritable, implying substantial genetic regulation of the level of mtDNA. Copy numbers also decrease modestly but significantly with age, and females on average have significantly more copies than males. The mtDNA copy numbers are significantly associated with waist circumference (p-value = 0.0031) and waist-hip ratio (p-value = 2.4×10^-5), but not with body mass index, indicating an association with central fat distribution. To our knowledge, this is the largest population analysis to date of mtDNA dynamics, revealing the age-imposed increase in heteroplasmy, the relatively high heritability of copy number, and the association of copy number with metabolic traits.     

Sensory Ataxic Neuropathy in Golden Retriever Dogs Is Caused by a Deletion in the Mitochondrial tRNATyr Gene

Sensory ataxic neuropathy (SAN) is a recently identified neurological disorder in golden retrievers. Pedigree analysis revealed that all affected dogs belong to one maternal lineage, and a statistical analysis showed that the disorder has a mitochondrial origin. A one base pair deletion in the mitochondrial tRNA^Tyr gene was identified at position 5304 in affected dogs after re-sequencing the complete mitochondrial genome of seven individuals. The deletion was not found among dogs representing 18 different breeds or in six wolves, ruling out this as a common polymorphism. The mutation could be traced back to a common ancestor of all affected dogs that lived in the 1970s. We used a quantitative oligonucleotide ligation assay to establish the degree of heteroplasmy in blood and tissue samples from affected dogs and controls. Affected dogs and their first to fourth degree relatives had 0–11% wild-type (wt) sequence, while more distant relatives ranged between 5% and 60% wt sequence and all unrelated golden retrievers had 100% wt sequence. Northern blot analysis showed that tRNA^Tyr had a 10-fold lower steady-state level in affected dogs compared with controls. Four out of five affected dogs showed decreases in mitochondrial ATP production rates and respiratory chain enzyme activities together with morphological alterations in muscle tissue, resembling the changes reported in human mitochondrial pathology. Altogether, these results provide conclusive evidence that the deletion in the mitochondrial tRNA^Tyr gene is the causative mutation for SAN.     

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.     

Haplotype-matched controls as a tool to discriminate polymorphisms from pathogenic mutations in mtDNA

We have studied the pathogenic role of 10044A-->G, a heteroplasmic mitochondrial DNA (mtDNA) mutation that has been suggested to be pathogenic in one family with severe pediatric morbidity. We found the mutation at an average frequency of 1.9% among 259 individuals including healthy controls. The mutation appeared to be heteroplasmic by restriction fragment analysis but analysis of subcloned polymerase chain reaction fragments confirmed homoplasmy. The polymorphic nature of 10044A-->G was verified by demonstrating exclusive association with a rare mtDNA haplotype within haplogroup H. We suggest that the evaluation of putatively pathogenic mutations in mtDNA should include the analysis of a sufficient number of haplotype-matched control samples and that the heteroplasmy should be verified by cloning.     

A quantitative evaluation of two methods for preserving hair samples

Hair samples are an increasingly important DNA source for wildlife studies, yet optimal storage methods and DNA degradation rates have not been rigorously evaluated. We tested amplification success rates over a one‐year storage period for DNA extracted from brown bear (Ursus arctos) hair samples preserved using silica desiccation and ?20 °C freezing. For three nuclear DNA microsatellites, success rates decreased significantly after a six‐month time point, regardless of storage method. For a 1000 bp mitochondrial fragment, a similar decrease occurred after a two‐week time point. Minimizing delays between collection and DNA extraction will maximize success rates for hair‐based noninvasive genetic sampling projects.     

Lack of Structural Variation but Extensive Length Polymorphisms and Heteroplasmic Length Variations in the Mitochondrial DNA Control Region of Highly Inbred Crested Ibis,Nipponia nippon

The animal mitochondrial DNA (mtDNA) length polymorphism and heteroplasmy are accepted to be universal. Here we report the lack of structural variation but the presence of length polymorphism as well as heteroplasmy in mtDNA control region of an endangered avian species – the Crested Ibis (Nipponia nippon). The complete control region was directly sequenced while the distribution pattern and inheritance of the length variations were examined using both direct sequencing and genotyping of the PCR fragments from captive birds with pedigrees, wild birds and a historical specimen. Our results demonstrated that there was no structural variation in the control region, however, different numbers of short tandem repeats with an identical motif of CA₃CA₂CA₃ at the 3′-end of the control region determined the length polymorphisms among and heteroplasmy within individual birds. There were one to three predominant fragments in every bird; nevertheless multiple minor fragments coexist in all birds. These extremely high polymorphisms were suggested to have derived from the ‘replication slippage’ of a perfect microsatellite evolution following the step-wise mutational model. The patterns of heteroplasmy were found to be shifted between generations and among siblings but rather stable between blood and feather samples. This study provides the first evidence of a very extensive mtDNA length polymorphism and heteroplasmy in the highly inbred Crested Ibis which carries an mtDNA genome lack of structural genetic diversity. The analysis of pedigreed samples also sheds light on the transmission of mtDNA length heteroplasmy in birds following the genetic bottleneck theory. Further research focusing on the generation and transmission of particular mtDNA heteroplasmy patterns in single germ line of Crested Ibis is encouraged by this study.     

Heteroplasmy and Ancient Translocation of Mitochondrial DNA to the Nucleus in the Chinese Horseshoe Bat (Rhinolophus sinicus) Complex

The utility and reliability of mitochondrial DNA sequences in phylogenetic and phylogeographic studies may be compromised by widespread and undetected nuclear mitochondrial copies (numts) as well as heteroplasmy within individuals. Both numts and heteroplasmy are likely to be common across diverse taxa yet few studies have characterised their frequencies and variation at the intra-specific level. Here we report the presence of both numts and heteroplasmy in the mitochondrial control region of the Chinese horseshoe bat Rhinolophus sinicus. In total we generated 123 sequences from 18 bats, which contained two different numt clades (i.e. Numt-1 and Numt-2) and one mtDNA clade. The sequence divergence between Numt-1 and Numt-2 was 16.8% and each numt type was found in all four R. sinicus taxa, suggesting either two ancient translocations of mitochondrial DNA into the nucleus from the same source taxon, or a single translocation from different source taxa that occurred before the split of R. sinicus into different lineages. Within the mtDNA clade, phylogenetic relationships among the four taxa of R. sinicus were similar to those seen in previous results. Based on PCR comparisons, heteroplasmy was inferred between almost all individuals of R. sinicus with respect to sequence variation. Consistent with introgression of mtDNA between Central sinicus and septentrionalis, individuals from these two taxa exhibited similar signatures of repeated sequences in the control region. Our study highlights the importance of testing for the presence of numts and heteroplasmy when applying mtDNA markers to phylogenetic studies.     

A family exhibiting heteroplasmy in the human mitochondrial DNA control region reveals both somatic mosaicism and pronounced segregation of mitotypes

A family exhibiting heteroplasmy at position 16 355 in hypervariable region I of the human mtDNA control region has been identified. This family consists of a mother, daughter, and son. DNA samples extracted from blood stains, buccal swabs, and hairs from these individuals were amplified by PCR and sequenced utilizing fluoresence-labeled dye terminator chemistry in an automated DNA sequencer. In both the daughter and mother, heteroplasmy was observed in DNA extracted from blood stains, buccal swabs, and hairs. In the blood stains, the proportion of cytosine was greater than thymine in both individuals. Buccal swab extracts showed a more balanced contribution from the two nucleotides. Telogenic hair root and hair shaft samples exhibited a wide range of nucleotide contributions at this position, from predominately cytosine in some samples to predominately thymine in others. The apparent stochastic segregation of mitotypes in hair samples is discussed from a forensic viewpoint, and the mechanism of mtDNA heteroplasmy is considered. Received: 6 November 1996 / Accepted: 13 February 1997