Page 1
Results of a collaborative study of the EDNAP group regardingmitochondrial DNA heteroplasmy and segregation in hair shafts
G. Tullya,*, S.M. Barritti, K. Bendere, E. Brignonb, C. Capellij, N. Dimo-Simoninc,C. Eichmannf, C.M. Ernsti, C. Lamberta, M.V. Lareud, B. Ludesb, B. Mevagh,
W. Parsonf, H. Pfeifferg, A. Salasd, P.M. Schneidere, E. Staalstromh
aThe Forensic Science ServiceTM, Trident Court, Solihull Parkway, Birmingham Business Park, Solihull B37 7YN, UKbCODGENE Laboratory 11, Rue Humann 67085, Strasbourg, France
cInstitut Universitaire de Medecine Legale, Rue de Bugnon 21, CH-1005 Lausanne, SwitzerlanddFaculty of Medicine, Institute of Legal Medicine, c/San Fransisco s/n, 15705 Santiago de Compestela, Spain
eInstitute of Legal Medicine, University of Mainz, Am Pulverturm 3, D 55131 Mainz, GermanyfInstitut fuer Gerichtliche Medizin, DNA-Labor, Universitaet Innsbruck, Muellerstrasse 44, A-6020 Innsbruck, Austria
gInstitut fuer Rechtsmedizjn der Universitaet Muenster, Von-Esmarch Strasse 62, D-48149 Muenster, GermanyhRettsmedisinsk Institutt, Rikshospitalet, 0027 Oslo, Norway
iDepartment of the Army, Armed Forces Institute of Pathology, DNA Identification Laboratory,
1413 Research Boulevard, Rockville, MD 209850, USAjIstituto di Midicina Legale, Universita Cattoilca del sacro Cuore, Largo F. Vito 1, 00168 Rome, Italy
Received 26 November 2002; received in revised form 5 March 2003; accepted 4 April 2003
Abstract
A collaborative exercise was carried out by the European DNA Profiling Group (EDNAP) in order to evaluate the distribution
of mitochondrial DNA (mtDNA) heteroplasmy amongst the hairs of an individual who displays point heteroplasmy in blood and
buccal cells. A second aim of the exercise was to study reproducibility of mtDNA sequencing of hairs between laboratories using
differing chemistries, further to the first mtDNA reproducibility study carried out by the EDNAP group [Forensic Sci. Int. 97
(1998) 165]. Laboratories were asked to type 2 sections from each of 10 hairs, such that each hair was typed by at least two
laboratories. Ten laboratories participated in the study, and a total of 55 hairs were typed.
The results showed that the C/T point heteroplasmy observed in blood and buccal cells at position 16234 segregated
differentially between hairs, such that some hairs showed only C, others only T and the remainder, C/T heteroplasmy at varying
ratios. Additionally, differential segregation of heteroplasmic variants was confirmed in independent extracts at positions 16093
and the poly(C) tract at 302–309, whilst a complete A–G transition was confirmed at position 16129 in one hair.
Heteroplasmy was observed at position 16195 on both strands of a single extract from one hair segment, but was not observed
in the extracts from any other segment of the same hair. Similarly, heteroplasmy at position 16304 was observed on both strands
of a single extract from one hair. Additional variants at positions 73, 249 and the HVII poly(C) region were reported by one
laboratory; as these were not confirmed in independent extracts, the possibility of contamination cannot be excluded.
Additionally, the electrophoresis and detection equipment used by this laboratory was different to those of the other
laboratories, and the discrepancies at position 249 and the HVII poly(C) region appear to be due to reading errors that may be
associated with this technology.
The results, and their implications for forensic mtDNA typing, are discussed in the light of the biology of hair formation.
# 2003 Published by Elsevier Ireland Ltd.
Keywords: Mitochondrial DNA; Heteroplasmy; Mutation; Segregation; Hair
Forensic Science International 140 (2004) 1–11
* Corresponding author. Tel.: þ44-121-399-5041; fax: þ44-121-622-2051.
E-mail address: [email protected] (G. Tully).
0379-0738/$ – see front matter # 2003 Published by Elsevier Ireland Ltd.
doi:10.1016/S0379-0738(03)00181-6
Page 2
1. Introduction
Mitochondrial DNA (mtDNA) analysis is widely used in
forensic examination of degraded remainsandsamples such as
hair shaftsdue to itshighcopy number [2]and hence its relative
abundance. The control region of human mtDNA, which is the
region most commonly analysed for forensic purposes, dis-
plays considerable sequence variation in human populations
(e.g. [3–9]). Comparison of mtDNA sequences from a ques-
tioned and a known sample can thus provide significant evi-
dence to exclude an individual as the source of an evidential
sample, or to support an association between the two. It is well
established that mtDNA mutates at a substantially higher rate
than that does nuclear DNA (e.g. [10–15]). Consequently, it is
not uncommon for differences to be observed in the mtDNA
sequence when comparing close maternal relatives (such as
mother and child) [15]. Substitution has also been observed in
somatic tissues, presumably due to segregation of existing
heteroplasmy within the individual. This means that differ-
ences may be observed between different hairs and/or tissues
within an individual (e.g. [16–19]). The clonal nature of hair
follicles and the high energy requirements of keratinizing hair
shaft cells are two features of hair histogenesis that could
contribute to the high observation of segregation of hetero-
plasmic variants in mtDNA from hair shafts [20].
Interpretation of mtDNA results in forensic casework
must take account of the possibility of sequence differences
within an individual [21–23]. This interpretation is aided by
accumulation of data concerning the occurrence and segre-
gation of mutation. The aim of this exercise was therefore
to add to the body of data concerning segregation of hetero-
plasmic mtDNA mutations in hairs. Additionally, the
exercise provided an opportunity to add to the earlier
reproducibility study undertaken by this group [1].
2. Materials and methods
2.1. Samples
Reference blood and buccal scrape samples from a donor
were extracted using Chelex100 resin (BioRad, Hemel
Hempstead, UK) following published protocols [24].
Multiple hairs with roots were obtained from the donor
and numbered sequentially. Each was cut into at least
5 cm � 2 cm sections, with the root section being given
the suffix ‘‘.1’’, the next section ‘‘.2’’ and so on.
Two sections from each of 10 hairs were sent to each
participating laboratory, such that each hair was analysed by
at least two laboratories; the root sections were retained for
subsequent analysis where required. Each laboratory used
their own protocols for the analysis.
2.2. Methods
The methods employed for washing the hairs, extracting,
amplifying and sequencing the DNA and for electrophoresis
and analysis of the sequenced products are summarised in
Table 1.
As shown in the table, a wide variety of methodologies
were employed. The wash procedures applied to the hairs
before digestion ranged from mild washing using only
ethanol and/or water to more vigorous washing using deter-
gents and/or proteinase K (pK) digestion from the external
surface of the hair.
For extraction, chemical digestion was employed, fol-
lowed by organic or resin-based purification and concentra-
tion by ethanol precipitation or spin columns. Where the root
section was analysed, a preferential extraction method was
used to first remove sheath material by a mild pK digestion
and washing, before the hair material itself was chemically
digested (pK, DTT, SDS) and the DNA organically purified
and ethanol precipitated.
PCR strategies included direct amplification of either
single or overlapping fragments, semi-nested PCR or nested
PCR. The majority of laboratories used an electrophoretic
assessment stage prior to sequencing, with some laboratories
also employing a spin column purification step.
Both dye primer and dye terminator sequencing chemis-
tries were employed, with a variety of dye chemistries.
Electrophoresis was performed using both slab-gels and
capillaries; all laboratories used real-time fluorescent
detection.
3. Results and discussion
The results are summarised in Table 2.
In instances where sequence was not obtained by two
different laboratories, or where the hair seemed to show a
graduation of heteroplasmy along its length, the root section
of the hair was extracted and analysed. For this analysis, the
sheath material was digested using proteinase K and dis-
carded prior to digestion of the hair itself using additional pK
plus dithiothreitol (DTT). The rationale of removing the soft
tissue adhering to the root in this manner was to remove the
possibility of mesoderm and endoderm cells from such
tissue, which may contain a different mtDNA population
[20], contributing to the mtDNA analysed.
The existence and differential segregation of sequence
heteroplasmy at positions 16093 and 16234 and of length
heteroplasmy at the polycytosine region between 302 and
310 was confirmed by the present study. Within some hairs,
graduation of heteroplasmy along the length of the hair was
apparent; for example, in hair 17, the section closest to the
root showed a T at position 16093, the next section showed a
T with a small C underneath, the next showed C at greater
peak height than T and the section nearest the tip showed C
with a small T beneath (Fig. 1). In other hairs, the change
along the length of the hair did not show a gradual variation;
an example is seen at position 16234 in hair 43, where
apparently homoplasmic T is observed in the two sections
nearest the root and the two sections at the tip, whilst the
2 G. Tully et al. / Forensic Science International 140 (2004) 1–11
Page 3
Table 1
Methods employed by participating laboratories
Laboratory Hair washing method Hair extraction method PCR primers PCR strategy Post-PCR assessment
and/or purification
Sequencing
strategy
Electrophoresis
instrument
Analysis software
1 2� 10% SDS washes, 2� extraction
buffer washes, all with vortexing
Chemical digestion
(pK, SDS and DTT) then
phenol chloroform extraction
and EtOH precipitation
First round: L15933/H00575;
second round: L15997/H16401;
L00029/H00408
Nested PCR
(2 � 25 cycles)
and direct PCR
(35 cycles)
Agarose gel electrophoresis
to quantify
BigDye primer
cycle sequencing
377 Sequence analysis v. 3.0,
sequence navigator
2 1% SDS washes, water washes,
ethanol washes, all with vortexing
Chemical digestion (DTT, pK)
then phenol chloroform
extraction and purification
with microcon
L15997; H16401; L29; H408 Direct PCR Agarose gel electrophoresis
followed by microcon
(Millipore) purification
dRhodamine dye
primer cycle
sequencing
373A Sequence analysis v. 3.0,
sequence navigator v. 1.0.1
3 1� ethanol wash, 2� water washes,
all with vortexing
Chelex, pK and DTT with
Centricon 100 purification
and concentration
L15926/H00580;
M13(-21)L15997/M13revH16401;
M13revL00029/M13(-21)H00408
Nested PCR
(2 � 30 cycles)
Agarose gel electrophoresis
to quantify and spin
column purification
BigDye terminator
cycle sequencing
(nested sequencing)
310 Sequence analysis v. 3.0,
sequence navigator v. 1.02b3
4 1� 80% ethanol wash, 1� water,
all with vortexing
Chelex, pK and DTT L15997/H16236; L16159/H16401;
L047/H285; L172/H408
Overlapping
fragments
PAGE assessment and spin
column purification
dRhodamine
terminator cycle
sequencing
377 Sequence analysis,
sequence navigator
5 1� ethanol wash, 2� water washes,
all with vortexing
Chemical digestion (pK, SDS
and DTT) then phenol
chloroform extraction and
amicon30 concentration
L15997/H16401; L291/H408 Direct PCR
(35 cycles)
Agarose gel electrophoresis
to quantify and spin
column purification
BigDye terminator
cycle sequencing
310 Sequence analysis v. 2.1.2,
sequence navigator v. 1.0.1
6 Mild pK digestion Chemical digestion (pK,
SDS and DTT) then phenol
chloroform extraction and EtOH
precipitation
L15971/H16410; L15/H484 Direct PCR
(35 cycles)
PAGE gel (sliver stained)
for assessment
BigDye terminator
cycle sequencing
310/3700 Sequence analysis v. 3.0
and v. 3.6
7 30 min 70% ethanol followed by
30 min HPLC-water
Chemical digestion proK, DTT,
PCR buffer, HPLC-water
phenol chloroform extraction,
microcon
L15997/H16401; L0029/H00389 Direct PCR
(30 cycles)
Qiagen BigDye terminator
cycle sequencing
310 Sequence analysis v. 3.0,
sequence navigator
8 1� water wash, 5% Tergazym 60 8Cfor 20 min, 1� water 37 8C for
10 min, 1� ethanol left to air dry
Chemical digestion (DTT, pK,
Laureth 10), phenol chloroform
extraction and EtOH precipitation
L15997/H16405; H16169,
H16209; L00029; H00408;
H00266
Semi-nested Agarose gel electrophoresis
and spin column purification
TaqFS dye primer
cycle sequencing
377 Sequence analysis v. 3.0,
sequence navigator
9 4� wash with 5% detergent solution,
wash with ethanol and water
Micro tissue grinder and chemical
digestion (proK, DTT) organic
extraction, Centricon 30 concentration
L15989/H16251; L16190/H16410;
L15/H285; L155/H389
Overlapping
fragments
(38 cycles)
Agarose gel electrophoresis
and spin column purification
BigDye terminator
cycle sequencing
377 Sequence analysis v. 3.0,
sequence navigator
10 1� ethanol wash Chemical digestion (pK,
SDS and DTT) and silica purification
L15990/H16401; L34/H370 Direct PCR
(38 cycles)
Agarose gel electrophoresis TaqFS dye primer
cycle sequencing
ALF express ALF win 1.0
G.
Tu
llyet
al./F
oren
sicS
cience
Intern
atio
na
l1
40
(20
04
)1
–1
13
Page 4
Table 2
Results, expressed as differences from the Anderson reference sequence (‘‘AND’’) [33]
4 G. Tully et al. / Forensic Science International 140 (2004) 1–11
Page 5
Table 2 (Continued )
G. Tully et al. / Forensic Science International 140 (2004) 1–11 5
Page 6
Table 2 (Continued )
6 G. Tully et al. / Forensic Science International 140 (2004) 1–11
Page 7
intermediate section has C as the major component with only
a small T peak underneath (Fig. 2). The same laboratory,
using the same methodology, sequenced two of the mid-
sections that demonstrated the dramatic shift in proportions.
This switch is therefore unlikely to be attributable to differ-
ences in sequencing chemistry. Linch et al. in their review of
hair histogenesis for the mtDNA forensic scientist [20],
highlight that a significant proportion of the mitochondria
in hair shafts are likely to originate from a relatively small
number of melanocytes, the remainder originating from
keratinocytes. Melanocytes, whilst not undergoing mitosis,
continue to produce mitochondria until later in life, when
grey hairs are produced. They therefore represent a bottle-
neck, through which heteroplasmic ratios could shift mark-
edly, both within and between individual hair shafts [20].
Intuitively, gradual changes in heteroplasmic ratios along the
length of a hair shaft would be expected rather than appar-
ently homoplasmic switches between one base and another
in adjacent sections of hair, arbitrarily cut into 2 cm lengths.
However, present sequencing technologies do not detect all
of the variability present, as has been demonstrated by using
more sensitive techniques such as DGGE [25]. Thus, if the
detection threshold of sequencing technologies is in the
order of 20%, apparently homoplasmic changes between
Table 2 (Continued )
The sequence referred to as ‘‘REF’’ was determined from blood and saliva, both of which gave the same result. Lower case indicates that
sequence was determined from only one strand. In such instances, as per the ISFG recommendations [22], if heteroplasmy was suspected, the
position was marked as not determined. In order to clarify the heteroplasmic proportions observed, the following nomenclature was adopted.
C�: C peak with small T peak underneath; C > T: C peak stronger than T peak; C � T: C and T peaks of approximately equal intensity; C ¼ T:
C and T peaks of equal intensity; T > C: T peak stronger than C peak, etc. NI was used to denote the absence of an inserted base; (�) was used
to denote an undetermined base; DEL was used to denote a deletion.
G. Tully et al. / Forensic Science International 140 (2004) 1–11 7
Page 8
adjacent sections of hair can be explained by an approxi-
mately 60% shift in the proportions of the heteroplasmic
variants (Fig. 3). Comparing fragments within each hair, a
shift in the proportions of heteroplasmic variants was
observed within 37 hairs, and an apparently homoplasmic
change observed within 11. Such changes in heteroplasmic
proportions within individual hairs has been documented
previously [19]. Wherever possible, analysis of two separate
segments of an evidential hair is considered by the EDNAP
group to be optimal [23]. Whilst this recommendation was
primarily aimed at ensuring that erroneous results due to
contamination are never reported, the results of the present
study show that sequencing more than one fragment of a hair
can give additional information to the forensic scientist to
aid interpretation. In only 13 of the hairs from this study
would a single fragment from each hair have given as much
information about the occurrence of heteroplasmy within the
donor individual as did the multiple hair segments.
An apparently complete base change was observed at
position 16129 in hair 39; no heteroplasmy was detected in
any other hair at this position. One potential explanation for
this result is that heteroplasmy at position 16129 exists
within the donor individual at a low level, and has
segregated, via a bottleneck with a very low number of
segregating units, into the melanocytes contributing mtDNA
to this one hair. An alternative explanation is that a de novo
mutation has arisen during the course of the formation of this
hair; study of maternal family members of the donor indi-
vidual would assist in determining which of these possibi-
lities is the more likely. In concordance with family studies,
phylogenetic analyses have shown that position 16129 is a
rapidly mutating base position [26,27]; similarly, position
16093 is a quickly mutating position, and 16148 has a higher
than average mutation rate. Moreover, positions 16093 and
16129 have been identified as by far and away the most
prominent sites for heteroplasmy in HV1 [25], to an extent
even in excess of their high relative mutation rates that have
been inferred phylogenetically. The demonstration that 1 of
55 hairs shows an apparently complete change at this posi-
tion underlines the need for careful consideration in the
interpretation of mtDNA evidence, especially in cases where
only one or a few bases differ between the reference and
unknown samples. Sequencing of multiple reference hairs
has been suggested in such instances (e.g. [21]); whilst often
providing useful information, the result in hair 39 shows that
even if additional hairs are typed, instances are likely to
Fig. 1. Gradual variation in the distribution of heteroplasmy at position 16093 along the length of hair 17. Section 17.2 is the section 2 cm
proximal to the root; sections 0.3–0.5 are progressively further away from the root. The reference blood and saliva show C and T in
approximately equal proportions, whilst this hair shows an increasing proportion of C along its length. In section 17.3, electrophoretic
‘‘shouldering’’ of each peak is present, but the heteroplasmy is distinguishable from this shouldering as it is raised in height, and is not
adjacent to a blue peak. This shouldering was not present in the sequence from the reverse strand, but the heteroplasmy was clearly
distinguishable.
8 G. Tully et al. / Forensic Science International 140 (2004) 1–11
Page 9
Fig. 2. Variation in the distribution of heteroplasmy at position 16234 along the length of hair 43. Section 43.1 is the root section, with
sections 43.2–43.5 being progressively further away from the root. The reference blood and saliva are apparently homoplasmic for the C
variant, section 43.3 shows C > T heteroplasmy, whilst all other sections are apparently homoplasmic T. Sections 43.3 and 43.4 were
sequenced by the same laboratory, using the same chemistry.
Fig. 3. Theoretical heteroplasmic proportions at position 16234 in the fragments of hair 43. Assuming a detection threshold of 20%, any point
above the upper dashed line would be observed as a C base, whilst any point below the lower dashed line would be observed as a T base; points
between the lines would be observed to exhibit C/T heteroplasmy.
G. Tully et al. / Forensic Science International 140 (2004) 1–11 9
Page 10
occur in which the heteroplasmy or base change is
not sufficiently common within the individual to be detected.
Thus, although any difference between unknown and
reference sample must diminish the strength of the evidence
[23], a single such difference cannot provide conclusive
evidence that the two samples are from different maternal
origins.
An intriguing suggestion of Linch et al. is the potential
link between degree of pigmentation and degree of hetero-
plasmy [20]. This link is postulated to be due to the
melanocyte contribution of mitochondria to hair shaft being
coincident with contribution of melanosomes, although the
authors note that there have been no reports in the literature
relating degree of pigmentation to observed mtDNA het-
eroplasmy. The donor of the hairs examined in the present
study is of Asian descent (belonging to the typically Asian
haplogroup F1a, characterised in the hypervariable regions
by positions 16129, 16172, 16304, 249del [28]) and as such
has extremely darkly pigmented hair. The EDNAP group
is currently undertaking a similar exercise, using an indi-
vidual in whose reference sample no heteroplasmy
was evident. This second donor individual has mid brown
hair; this and similar studies may cast more light on the
hypothesised link between pigmentation and heteroplasmy
in hair shafts. A range of hair colours was studied by
Huhne et al., but no heteroplasmy was recorded in their
study [29].
Heteroplasmy has been observed at multiple positions
within a single individual in several previous studies
[18,19,25,30]. The findings of Grzybowski [30] have been
challenged previously [31,32], largely on the basis of poten-
tial contamination or amplification of nuclear pseudogenes.
With respect to contamination, the EDNAP group has pre-
viously stated that ‘‘in the event of duplication failing, or
where insufficient sample is present, the report should
clearly state the limitations of an unconfirmed result. It is
possible to have apparently clear negative controls and yet to
obtain an erroneous result due to contamination’’ [23].
Following this guideline, the potential heteroplasmy at
position 16095 in hair fragment 34.2 (this position is not
a common variant in human population studies), at position
16148 in hair fragment 44.3, at position 73 in fragment 60.3,
and the potential base changes at positions 249 in fragments
58.3 and 60.3, and position 315.1 in fragments 57.3, 58.4,
59.3 and 60.3 would be treated with caution. The results
from hair fragment 60.3 are worthy of particular considera-
tion: three differences to the results from the other two
fragments of this hair were observed. These results may
suggest that contamination has been encountered in the
analysis of this hair fragment. Overall, the results from
laboratory 10 show several anomalies in comparison to
the results from laboratories 1–9 (positions 73, 249,
315.1) and the possibility of contamination cannot be ruled
out. This contamination may have been introduced during
the analysis, or may have been present on the surface of the
hair and have been inefficiently removed; of all the methods
used for washing the hairs, that of laboratory 10 was the least
stringent. More detailed examination of the electrophero-
grams from laboratory 10 has highlighted that the resolution
was poor in some regions, particularly around position 249
and the poly(C) stretch in HVII. This laboratory was the only
laboratory in this study using sequencing technology requir-
ing the four reactions (A, C, G and T) to be electrophoresed
in separate lanes; such technology increases the potential for
error in sequencing, particularly in polynucleotide regions,
and should not therefore be used to assess the length of
polynucleotide stretches in forensic analyses.
Nevertheless, the results from the remaining laboratories
show a high level of consistency. In particular, it is note-
worthy that despite the relatively high levels of heteroplasmy
observed in this study, the maximum number of homoplas-
mic differences between any single hair and the reference
blood and saliva was found to be one (excluding results from
laboratory 10). While it is clear that in rare cases two or more
homoplasmic differences would be possible, this study
demonstrates that the magnitude of mutational segregation
within hair samples is tractable, and follows predictable
patterns of increasing rarity with increasing heteroplasmic
and homoplasmic differences. If interpreted with care, fol-
lowing published guidelines [21–23], such results can pro-
vide valuable and reliable information to assist in forensic
examinations.
Acknowledgements
The authors wish to thank all of the members of the
EDNAP group for useful comments and discussion. Parti-
cular thanks for helpful comments on the manuscript are due
to Angel Carracedo and Tom Parsons. The EDNAP group
worked in the period 1997–2000 within the framework of the
Standardisation of DNA Profiling in Europe (STADNAP)
consortium, a network project of the European Commis-
sion—DG XII programme ‘‘Standards, Measurement and
Testing’’ (Contract 97-7506).
References
[1] A. Carracedo, E. D’Aloja, B. Dupuy, A. Jangblad, M.
Karjalainen, C. Lambert, W. Parson, H. Pfeiffer, H. Pfitzinger,
M. Sabatier, D. Syndercombe Court, C. Vide, Reproducibility
of mtDNA analysis between laboratories: a report of the
European DNA profiling group (EDNAP), Forensic Sci. Int.
97 (1998) 165–170.
[2] D. Bogenhagen, D. Clayton, The number of mitochondrial
deoxyribonucleic acid genomes in mouse L and human HeLa
cells, J. Biol. Chem. 249 (1974) 7991–7995.
[3] C. Aquadro, B. Greenberg, Human mitochondrial DNA
variation and evolution: analysis of nucleotide sequences
from seven individuals, Genetics 103 (1983) 287–312.
[4] L. Vigilant, R. Pennington, H. Harpending, T. Kocher, A.
Wilson, Mitochondrial DNA sequences in single hairs from a
10 G. Tully et al. / Forensic Science International 140 (2004) 1–11
Page 11
southern African population, Proc. Natl. Acad. Sci. U.S.A. 86
(1989) 9350–9354.
[5] M. Stoneking, D. Hedgecock, R. Higuchi, L. Vigilant, H.
Erlich, Population variation of human mtDNA control region
sequences detected by enzymatic amplification and sequence-
specific oligonucleotide probes, Am. J. Hum. Genet. 48
(1991) 370–382.
[6] C. Herrnstadt, J. Elson, E. Fahy, G. Preston, D. Turnbull, C.
Anderson, S. Ghosh, J. Olefsky, M. Beal, R. Davis, N.
Howell, Reduced-median-network analysis of complete
mitochondrial DNA coding-region sequences for the major
African, Asian, and European haplogroups, Am. J. Hum.
Genet. 70 (2002) 1152–1171.
[7] T. Melton, M. Wilson, M. Batzer, M. Stoneking, Extent of
heterogeneity in mitochondrial DNA of European popula-
tions, J. Forensic Sci. 42 (1997) 437–446.
[8] P. Forster, F. Call, A. Roe Metspalu, R. D’anna, M. Mirisola,
G. De Leo, A. Flugy, A. Salerno, G. Ayala, A. Kouvatsi, R.
Willems, V. Romano, Continental and subcontinental dis-
tributions of mtDNA control region types, Int. J. Legal Med.
116 (2002) 99–108.
[9] K. Imaizumi, T. Parsons, M. Yoshino, M. Holland, A new
database of mitochondrial DNA hypervariable regions I and II
sequences from 162 Japanese individuals, Int. J. Legal Med.
116 (2002) 68–73.
[10] W. Brown, M. George, A. Wilson, Rapid evolution of animal
mitochondrial DNA, Proc. Natl. Acad. Sci. U.S.A. 76 (1979)
1967–1971.
[11] W. Brown, E. Prager, A. Wang, A. Wilson, Mitochondrial
DNA sequences in primates: tempo and mode of evolution, J.
Mol. Evol. 18 (1982) 225–239.
[12] T. Miyata, H. Hayashida, R. Kikuno, M. Hasegawa, M.
Kobayashi, K. Koike, Molecular clock of silent substitution:
at least six-fold preponderance of silent changes in mitochon-
drial genes over those in nuclear genes, J. Mol. Evol. 19
(1982) 28–35.
[13] D. Wallace, J. Ye, S. Neckelman, G. Singh, K. Webster, B.
Greenberg, Sequence analysis of cDNAs for human and
bovine ATP synthase beta subunit: mitochondrial DNA genes
sustain seventeen times more mutations, Curr. Genet. 12
(1987) 81–90.
[14] D. Wallace, M. Lott, A. Torroni, J. Shoffner, Report of the
committee on human mitochondrial DNA, Cytogenet. Cell
Genet. 58 (1991) 1103–1123.
[15] T.J. Parsons, D.S. Muniec, K. Sullivan, N. Woodyatt, R.
Alliston-Greiner, M.R. Wilson, D.L. Berry, K.A. Holland,
V.W. Weedn, P. Gill, M.M. Holland, A high observed
substitution rate in the human mitochondrial DNA control
region, Nat. Genet. 15 (1997) 363–368.
[16] K.M. Sullivan, R. Alliston-Greiner, F.I.A. Archampong, R.
Piercy, G. Tully, P. Gill, C. Lloyd-Davies, A single difference
in mtDNA control region sequence observed between hair
shaft and reference samples from a single donor, in:
Proceedings from the Seventh International Symposium on
Human Identification, Scottsdale, AZ, Promega Corporation,
1996, pp. 126–129.
[17] M.R. Wilson, D. Polanskey, J. Repogle, J.A. DiZinno, B.
Budowle, A family exhibiting heteroplasmy in the human
mitochondrial DNA control region reveals both somatic
mosaicism and pronounced segregation of mitotypes, Hum.
Genet. 100 (1997) 167–171.
[18] C.D. Calloway, R.L. Reynolds, G.L. Herrin, W.W. Anderso,
The frequency of heteroplasmy in the HVII region of mtDNA
differs across tissue types and increases with age, Am. J.
Hum. Genet. 66 (2000) 1384–1397.
[19] A. Salas, M. Lareu, A. Carracedo, Heteroplasmy in mtDNA
and the weight of evidence in forensic mtDNA analysis: a
case report, Int. J. Legal Med. 114 (2001) 186–190.
[20] C. Linch, D. Whiting, M. Holland, Human hair histogenesis
for the mitochondrial DNA forensic scientist, J. Forensic Sci.
46 (2001) 844–853.
[21] M.M. Holland, T.J. Parsons, Mitochondrial DNA sequence
analysis—validation and use for forensic casework, Forensic
Sci. Rev. 11 (1999) 21–50.
[22] W. Bar, B. Brinkmann, B. Budowle, A. Carracedo, P. Gill, M.
Holland, P.J. Lincoln, W. Mayr, N. Morling, B. Olaisen, P.M.
Schneider, G. Tully, M. Wilson, DNA commission of the
international society for forensic genetics: guidelines for
mitochondrial DNA typing, Forensic Sci. Int. 110 (2000)
79–85.
[23] G. Tully, W. Bar, B. Brinkmann, A. Carracedo, P. Gill, N.
Morling, W. Parson, P. Schneider, Considerations by the
European DNA Profiling (EDNAP) Group on the working
practices, Forensic Sci. Int. 124 (2001) 83–91.
[24] R. Sparkes, C. Kimpton, S. Watson, N. Oldroyd, T. Clayton,
L. Barnett, J. Arnold, C. Thompson, R. Hale, J. Chapman, A.
Urquhart, P. Gill, The validation of a 7-locus multiplex STR
test for use in forensic casework. I. Mixtures, ageing,
degradation and species studies, Int. J. Legal Med. 109
(1996) 186–194.
[25] L. Tully, T. Parsons, R. Steighner, M. Holland, M. Marino, V.
Prenger, A sensitive DGGE assay reveals a high frequency of
heteroplasmy in hypervariable region one of the human
mitochondrial DNA control region, Am. J. Hum. Genet. 67
(2000) 432–443.
[26] H. Bandelt, L. Quintana-Murci, A. Salas, V. Macaulay, The
fingerprint of phantom mutations in mitochondrial DNA data,
Am. J. Hum. Genet. 71 (2002) 1150–1160.
[27] B. Malyarchuk, I. Rogozin, V. Berikov, M. Derenko, Analysis
of phylogenetically reconstructed mutational spectra in
human mitochondrial DNA control region, Hum. Genet. 111
(2002) 453.
[28] Y. Yao, Q. Kong, H. Bandelt, T. Kivisild, Y. Zhang,
Phylogenetic differentiation of mitochondrial DNA in Han
Chinese, Am. J. Hum. Genet. 70 (2002) 635–651.
[29] J. Huhne, H. Pfeiffer, K. Waterkamp, B. Brinkmann,
Mitochondrial DNA in human hair shafts—evidence of
intra-individual differences? Int. J. Legal Med. 112 (1999)
172–175.
[30] T. Grzybowski, Extremely high levels of heteroplasmy in
single hair roots, Electrophoresis 21 (2000) 548–553.
[31] G. Tully, M. Lareu, Letters to the editor: controversies over
heteroplasmy, Electrophoresis 22 (2001) 180–182.
[32] B. Budowle, M. Allard, M. Wilson, Critique of interpretation
of high levels of heteroplasmy in the human mitochondrial
DNA hypervariable region I from hair roots, Forensic Sci. Int.
126 (2002) 30–33.
[33] S. Anderson, A.T. Bankier, B.G. Barrell, M.H.L. de Bruijin,
A.R. Coulson, J. Drouin, I.C. Eperon, D.P. Nierlich, B. A
Roe, F. Sanger, P.H. Schreier, A.J.H. Smith, R. Staden, I.G.
Young, Sequence and organization of the human mitochon-
drial genome, Nature 290 (1981) 457–463.
G. Tully et al. / Forensic Science International 140 (2004) 1–11 11