Open Access
Open Peer Review

This article has Open Peer Review reports available.

How does Open Peer Review work?

Minisequencing mitochondrial DNA pathogenic mutations

  • Vanesa Álvarez-Iglesias1, 2,
  • Francisco Barros2,
  • Ángel Carracedo1, 2 and
  • Antonio Salas1Email author
BMC Medical Genetics20089:26

https://doi.org/10.1186/1471-2350-9-26

Received: 26 September 2007

Accepted: 10 April 2008

Published: 10 April 2008

Abstract

Background

There are a number of well-known mutations responsible of common mitochondrial DNA (mtDNA) diseases. In order to overcome technical problems related to the analysis of complete mtDNA genomes, a variety of different techniques have been proposed that allow the screening of coding region pathogenic mutations.

Methods

We here propose a minisequencing assay for the analysis of mtDNA mutations. In a single reaction, we interrogate a total of 25 pathogenic mutations distributed all around the whole mtDNA genome in a sample of patients suspected for mtDNA disease.

Results

We have detected 11 causal homoplasmic mutations in patients suspected for Leber disease, which were further confirmed by standard automatic sequencing. Mutations m.11778G>A and m.14484T>C occur at higher frequency than expected by change in the Galician (northwest Spain) patients carrying haplogroup J lineages (Fisher's Exact test, P-value < 0.01). The assay performs well in mixture experiments of wild:mutant DNAs that emulate heteroplasmic conditions in mtDNA diseases.

Conclusion

We here developed a minisequencing genotyping method for the screening of the most common pathogenic mtDNA mutations which is simple, fast, and low-cost. The technique is robust and reproducible and can easily be implemented in standard clinical laboratories.

Background

There are over 100 point mutations putatively associated with human mtDNA diseases [1]; however, only a small percentage of these mutations were properly confirmed http://www.mitomap.org/rimtab2.html. In fact, the pathogenicity of some mutations is still under question due to the fact that the conclusions claimed in a substantial number of studies do not rest on solid grounds or are partially or completely flawed [25].

During this last two decades, many different techniques were designed for the screening of mtDNA pathogenic mutations, including restriction fragment length polymorphism (RFLP) analysis, heteroduplex analysis (HDA), single strand conformation polymorphisms (SSCP), etc [611]. The main disadvantage of these classical screening methods is that, due to the intrinsic nature of the techniques, many mutations usually pass unnoticed. Novel and more efficient strategies have recently been proposed that allow the screening of entire mtDNA genomes; for instance, the mismatch-specific DNA endonuclease "Surveyor™ Nuclease" seems to be particularly useful when the most common pathogenic mutations have been ruled out previously (e.g. [12]). Ideally, the gold standard of sequencing the whole genome would solve the problem; nevertheless, this strategy is labor-intensive and therefore unfeasible for routine molecular diagnosis (taking into account the capacities of standard laboratories). Screening methods are also difficult to standardize and some of them are not suitable for replication assays due to their dependence on for instance electrophoretic conditions.

Minisequencing techniques have demonstrated to be efficient in other fields of research, such as forensic and population genetics [1318], but also in clinical research [19]. This approach, also known as single nucleotide primer extension, allows the determination of the polymorphic position by DNA polymerase addition of the ddNTP complementary to the base interrogated. The primer anneals therefore to its target DNA immediately adjacent to the SNP under scrutiny. Here we present a minisequencing multiplex design which is a rapid, sensitive and low cost assay for screening mtDNA pathogenic mutations in patients suffering mtDNA disease.

Methods

Subject and SNP selection

Blood sample of 15 patients was sent to the Fundación Pública Galega de Medicina Xenómica located in the Hospital Clínico of Santiago de Compostela (Galicia, Spain). Six of these patients are from Galicia (northwest Spain; [20]), and three other patients are from southeast Spain (in the Mediterranean coast); the regional Spanish provenance of the other two patients is unknown. Four of these patients presented clinical features of an undetermined neuropathy, while the rest of the patients were analyzed under suspicion of Leber disease. Written informed consent was obtained from all patients and the study was approved by the institutional review board of the University of Santiago de Compostela (Spain). All the samples were genotyped for a set of 25 variants, 21 are confirmed pathogenic mutations related to different mtDNA diseases such as MELAS, LHON, Leigh Disease, and NARP. These mutations are distributed all around the mtDNA molecule and 'hit' different mtDNA genes (Table 1). All the selected mutations, with the only exception of m.3243A>G, consist of non-synonymous substitutions (both transitions and transversion, but note that the minisequencing technique could be also suitable for interrogating indels [17]).
Table 1

Characteristics of the mutations incorporated in our multiplex design

POSITIONS

Base change

Amino Acid change

Locus

Disease1

3243

A>G

tRNA Leu (UUR)

MT-TL1

MELAS, DM/DMDT, CPEO, MM

3460

G>A

Ala>Thr

MT-ND1

LHON

3697

G>A

Gly>Arg

MT-ND1

MELAS

3946

G>A

Glu>Lys

MT-ND1

MELAS

3949

T>C

Tyr>Lys

MT-ND1

MELAS

7445

A>G

Ter>Ter

MT-CO1

SNHL

7445

A>C

Ter>Ser

MT-CO1

DEAF

8993

T>G

Leu>Arg

MT-ATP6

NARP

8993

T>C

Leu>Pro

MT-ATP6

NARP, Leigh Disease

9176

T>C

Leu>Pro

MT-ATP6

FBSN, Leigh Disease

9176

T>G

Leu>Arg

MT-ATP6

Leigh Disease

10158

T>C

Arg>Pro

MT-ND3

Leigh Disease

10191

T>C

Arg>Pro

MT-ND3

ESOC, Leigh-like Disease

10663

T>C

Val>Ala

MT-ND4

LHON

11777

C>A

Arg>Ser

MT-ND4

Leigh Disease

11778

G>A

Arg>His

MT-ND4

LHON

11832

G>A

Trp>Ter

MT-ND4

Exercise Intolerance

12706

T>C

Phe>Leu

MT-ND5

Leigh Disease

13513

G>A

Asp>Asn

MT-ND5

MELAS, Leigh Disease

13514

A>G

Asp>Gly

MT-ND5

MELAS

14459

G>A

Ala>Val

MT-ND6

LYDT, Leigh Disease

14482

C>A

Met>Ile

MT-ND6

LHON

14482

C>G

Met>Ile

MT-ND6

LHON

14484

T>C

Met>Val

MT-ND6

LHON

14487

T>C

Met>Val

MT-ND6

Dystonia, Leigh Disease

1 List of some diseases where these mutations are frequently observed: CPEO: Chronic Progressive External Ophthalmoplegia; DEAF Maternally inherited DEAFness or aminoglycoside-induced DEAFness DM: Diabetes Mellitus; LDYT: Leber's hereditary optic neuropathy and DysTonia; LHON: Leber Hereditary Optic Neuropathy; MELAS: Mitochondrial Encephalomyopathy, Lactic Acidosis, and Stroke-like episodes; MM: Mitochondrial Myopathy; NARP: Neurogenic muscle weakness, Ataxia, and Retinitis Pigmentosa; alternate phenotype at this locus is reported as Leigh Disease; SNHL: SensoriNeural Hearing Loss

Primer design

The primers both for PCR amplification (Table 2) and minisequencing reaction (Table 3) were designed to have an annealing temperature around 60°C using Primer3 software http://frodo.wi.mit.edu/cgi-bin/primer3. The sequence databases at the National Centre for Biotechnology Information (NCBI; http://www.ncbi.nlm.nih.gov) were interrogated using the online BLAST tool to test the primers against possible repetitive sequences and sequence homologies in the autosomal genome. Each primer pair for PCR amplification and each single base extension primer were selected independently and AutoDimer http://www.cstl.nist.gov/biotech/strbase/AutoDimerHomepage/AutoDimerProgramHomepage.htm was used to test for potential hairpin structures and primer-dimer problems.
Table 2

Amplification primers

POSITIONS

FORWARD

REVERSE

SIZE (bp)

Final Concentration (μM)

3243

tatacccacacccacccaag

ggccatgggtatgttgttaag

118

0.2

3460

ccgaacgaaaaattctaggc

gcggtgatgtagagggtgat

153

0.2

3697

gcctagccgtttactcaatcc

tgagattgtttggggctactgc

94

0.15

3946/39491

tagcagagaccaaccgaacc

gaagattgtagtggtgagggtgt

157

0.2

7445

ccctaccacacattcgaagaa

tggcttgaaaccagctttg

89

0.2

8993

aatgccctagcccacttctt

aggtggcctgcagtaatgtt

140

0.15

9176

aaatcgctgtcgccttaatc

tcattaggagggctgagagg

154

0.3

10158/10191

tcaacaccctcctagcctta

gggtaaaaggagggcaattt

196

0.3

10663

acacccactccctcttagcc

ggccatatgtgttggagattg

110

0.3

11777/11778/11832

cacgggcttacatcctcatt

gggggtaaggcgaggttag

157

0.2

127061

tgtagcattgttcgttacatgg

agttggaataggttgttagcgg

146

0.2

13513/13514

attggcagcctagcattagc

cagggaggtagcgatgagag

131

0.2

14459/14482/14484/14487

ctccatcgctaaccccacta

ttctgaattttgggggaggt

170

0.4

1These primers were reported in [16]

Table 3

Minisequencing primers

POSITIONS

Extension Primer1

Length

Base change

Strand

Final Concentration (μM)3

10158

acaactcaacggctacatagaaaaa

25

T>C

L

0.2

3946

Cgaactagtctcaggcttcaacatc

25

G>A

L

0.2

10663

GACTGcaatattgtgcctattgccatactag

31

T>C

L

0.2

3460

(GACT)2Ggctactacaacccttcgctgac

31

G>A

L

0.2

10191

(GACT)4agtgcggcttcgaccctata

36

T>C

L

0.3

14459

(GACT)2GACctcaggatactcctcaatagccatc

36

G>A

L

0.3

144842

(GACT)3atcgctgtagtatatccaaagacaacYa

40

T>C

L

0.5

3243

(GACT)4GAacagggtttgttaagatggcag

40

A>G

L

0.2

11777

(GACT)5Gcaaactacgaacgcactcacagt

44

C>A

L

0.2

7445

(GACT)5tcgaagaacccgtatacataaaatctag

48

A>G/C

L

0.2

12706

(GACT)9gcggtaactaagattagtatggtaattagga

52

T>C

H

0.2

39492

(GACT)11GAagtctcaggcttcaacatcRaa

52

T>C

L

0.2

144872

(GACT)7GctgtagtatatccaaagacaaccaSca

56

T>C

L

0.6

13513

(GACT)8ttcctcacaggtttctactccaaa

56

G>A

L

0.1

13514

(GACT)10Ttgcggtttcgatgatgtgg

60

A>G

H

0.2

3697

(GACT)9CTaaactcaaactacgccctgatc

60

G>A

L

0.2

11778

(GACT)9GAgaagtccttgagagaggattatgatg

64

G>A

H

0.3

11832

(GACT)8 GACtcaaactctactcccactaatagcttttt

64

G>A

L

0.1

9176

(GACT)11Gatccaagcctacgttttcacacttc

70

T>C/G

L

0.2

14482

(GACT)11gccatcgtcgctgtagtatatccaaagacaac

73

C>A/G

L

0.6

8993

(GACT)13Gcctactcattcaaccaatagccc

76

T>G/C

L

0.2

1 Capital letters indicate the segment of the primer belonging to the tail

2 These primers are 'degenerated' at the indicated position in the second column (see Material and Methods) using the IUB code

Although the simultaneous occurrence of more than one pathogenic mutation in a single individual is infrequent, it could be possible that the presence of m.14484T>C could interfere with the genotyping of m.14487T>G (due to their physical proximity). Therefore, in order to avoid potential artifacts due to deficient annealing, we have designed a degenerate extension primer containing a mixture of nucleotide C and T at site 14484. The same rationale applies to variation at other two positions. According to our design, m.14482C>G mutation could interfere with the genotyping at position 14484; and the presence of m.3946G>A could alter the genotyping of m.3949T>C; we therefore designed another extension degenerate primer containing C and G at site 14482 and another additional extension degenerated primer carrying G and A at site 3946 (Table 3).

PCR multiplex amplification

The SNPs were PCR amplified in 13 amplicons with sizes ranging from 89 and 196 bp. The amplicons are deliberately designed to be small in order to facilitate the analysis of samples that are highly degraded or with low quantity of DNA. Some amplicons encompass several SNPs (Table 1). We performed multiplex reaction using 5 ng of DNA template and PCR master mix of QIAGEN Kit Multiplex PCR (Qiagen, Hilden, Germany), amplification primers and their final concentrations are in Table 2. Amplification was carried out in a 9700 Thermocycler (Applied Biosystems, Foster City, CA, USA). After a 95°C pre-incubation step for 15 minutes, PCR was performed in a total of 30 cycles using the following conditions: 94°C denaturation for 30 seconds, annealing at 60°C for 90 seconds and extension at 72°C during 90 seconds, followed by a 15 minutes of final extension at 72°C and 4°C until removed from thermocycler. PCR products were checked by polyacrylamide gel electrophoresis (T9, C5) visualized by silver staining.

Minisequencing reaction

Previous to minisequening reaction, PCR products are treated with ExoSAP-IT (Amershan Bioscences, Uppsala, Sweden) to remove excess primers and un-incorporated dNTPs: 3 μl of PCR product was incubated with 1.5 μl for 15 minutes at 37°C followed by 15 minutes at 80°C for enzyme inactivation. Minisequencing reaction is performed using SNaPshot™ Kit (AB). We modified the length of the primers (between 25 and 76 bp) by the addition of non-homologous tails, poly(dGACT) added at the 5'- end (Table 3). The minisequencing reaction was carried out in a total volume of 10 μl comprising 3 μl of the SNaPshot™ Kit (AB), 1.5 μl PCR product, 1.5 μl of extension primers mix (final concentrations are between 0.1 and 0.6 μM) (Table 3), and water up to 10 μl. The reaction was performed in a 9700 Thermocycler (AB) following the recommendations of the manufacturer: 25 cycles of denaturation at 96°C for 10 seconds, annealing at 50°C for 5 seconds and extension at 60°C during 30 seconds. Un-incorporated ddNTPs are eliminated with a treatment with SAP (Amershan Bioscences). The final volume (10 μl) was treated with 1 μl of SAP for 60 minutes at 37°C followed by 15 minutes at 80°C for enzyme inactivation.

The minisequencing products (1.5 μl) were mixed with 10 μl of HiDi™ formamide (AB) and 0.25 μl of GeneScan -120 LIZ (AB) and capillary electrophoresis was undertaken on an ABI PRISM 3130 ×l Genetic Analyzer (AB). The data was analyzed using GeneMapper™ 3.7 Software (AB).

Sequencing reaction

All the 13 amplicons were sequenced in an ABI PRISM 3130 ×l Genetic Analyzer (AB) for all the patients with the aim of corroborating the mutations observed with the minisequencing assay. Technical details regarding sequencing reaction are given in [17].

Results

Only 13 amplicons are needed to interrogate the selected 25 variants. This is because (i) 'neighbouring' variants are genotyped using different minisequencing probes annealing at the same amplicon, and (ii) different minisequencing probes were designed to genotype different pathogenic substitutions occurring at the same mtDNA site. Therefore, for each of the sites 7445, 9176, and 14482, there are two common mtDNA mutations described in the literature as pathogenic, the most common one for each of these three sites was already confirmed (m.7445A>G, m.9176T>C, and m.14482C>G), while the alternative one (m.7445A>C, m.9176T>G, and m.14482C>A) still bears the status of 'provisional' in MITOMAP. In addition, position 8993 has two already confirmed variants, namely, T to C and T to G substitutions.

We did not observe mutations in four patients presenting neurological symptoms. On the other hand, a pathogenic mutation was found for each of the other 11 patients with clinical suspicion of Leber disease (Table 4). All the mutations were confirmed by sequencing and were observed as homoplasmic using both minisequencing and automatic sequencing (Figure 1). Nine patients carried mutation m.11778G>A whereas two patients carried m.14484T>C.
Figure 1

Electropherograms showing (a) two different SNaPshot profiles carrying mutations m.11778G>A and m.14484T>C; rCRS [50] electropherogram is also shown indicating the whole set of mutations tested with the SNaPshot reaction; and (b) forward and reverse sequence electropherograms for these same mutations.

Table 4

Patient samples analyzed in the present study

  

3

3

3

3

3

7

8

9

1

1

1

1

1

1

1

1

1

1

1

1

1

 
  

2

4

6

9

9

4

9

1

0

0

0

1

1

1

2

3

3

4

4

4

4

 
  

4

6

9

4

4

4

9

7

1

1

6

7

7

8

7

5

5

4

4

4

4

 
  

3

0

7

6

9

5

3

6

5

9

6

7

7

3

0

1

1

5

8

8

8

 
          

8

1

3

7

8

2

6

3

4

9

2

4

7

Variants1

 

rCRS

A

G

G

G

T

A

T

T

T

T

T

C

G

G

T

G

A

G

T

T

T

 

#1

I-111

-

-

-

-

-

-

-

-

-

-

-

-

A

-

-

-

-

-

-

-

-

 

#2

D-112

-

-

-

-

-

-

-

-

-

-

-

-

A

-

-

-

-

-

-

-

-

11812

#3

I-138

-

-

-

-

-

-

-

-

-

-

-

-

A

-

-

-

-

-

-

-

-

10172

#4

F-393

-

-

-

-

-

-

-

-

-

-

-

-

A

-

-

-

-

-

-

-

-

 

#5

C-436

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

 

#6

B-703

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

 

#7

B-900

-

-

-

-

-

-

-

-

-

-

-

-

A

-

-

-

-

-

-

-

-

 

#8

B-901

-

-

-

-

-

-

-

-

-

-

-

-

A

-

-

-

-

-

-

-

-

 

#9

D-992

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

 

#10

1586

-

-

-

-

-

-

-

-

-

-

-

-

A

-

-

-

-

-

-

-

-

 

#11

5699

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

C

-

12696

#12

5869

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

C

-

12696

#13

G-495

-

-

-

-

-

-

-

-

-

-

 

-

A

-

-

-

-

-

-

-

-

 

#14

7218

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

-

 

#15

9615

-

-

-

-

-

-

-

-

-

-

-

-

A

-

-

-

-

-

-

-

-

12705C/T

1 Variants detected by automatic sequencing (putatively unrelated to the disease outcome). Position 12705 is heteroplasmic C/T.

The two carriers of m.14484T>C lived in the same geographic (relatively isolated) village (namely, Burela) in the northern cost of Galicia [20, 21]. These samples were submitted at different times to the laboratory and according to the carriers, their families were unrelated. Both patients share also the transition m.12696T>C (observed by sequencing), indicating that these individuals could actually share some recent common ancestor. The m.12696T>C transition is not a frequent variant in Europe and likely constitutes a diagnostic variant of a minor HV1 sub-lineage (although this variant is also recurrent in at least other eight non-European haplogroups, including haplogroup J, see below): it appears in at least six haplogroup HV1 complete genomes (see for instance [22, 23]). Therefore, the m.14484T>C transition could represent a founder pathogenic mutation responsible of Leber disease in Galicia (as it is the case of a wide spectrum of different non-mtDNA diseases [24, 25].

In order to further investigate this hypothesis, we additionally genotyped all the Galician samples for a set of haplogroup diagnostic mtDNA SNPs following [16] (Table 5). The two carriers of m.14484T>C but also two other Galician samples carrying m.11778G>A could be allocated to haplogroup J (and not to HV1 as it could be inferred by the presence of m.12696T>C). Therefore, the observation of the m.14484T>C mutation in two unrelated patients from the above mentioned small Galician village could just reflect the already reported high incidence of Leber mutations within haplogroup J [26, 27]. In fact, we also observe that the frequency of haplogroup J in healthy unrelated Galician individuals is ~14.7% [20, 21, 28] but it occurs in four out of our six Galician patients (~66.7%); the difference in prevalence is statistically significant (Fisher's Exact test, P-value < 0.006)
Table 5

Coding SNP and haplogroup status of the Galician patients

           

1

1

1

1

1

1

1

1

 
  

3

3

4

4

4

4

4

6

7

0

0

0

0

2

2

3

4

 
  

9

9

2

5

5

7

7

7

0

3

4

4

8

3

7

9

7

 
  

1

9

1

2

8

6

9

7

2

9

0

6

7

0

0

6

6

 
  

5

2

6

9

0

9

3

6

8

8

0

3

3

8

5

6

6

HG

 

rCRS

G

C

T

A

G

A

A

T

C

A

C

T

T

A

C

A

C

 

#3

I-138

-

-

C

-

-

-

-

-

T

G

-

-

-

-

-

-

T

J

#4

F-393

-

-

C

-

-

-

-

-

T

G

-

-

-

-

-

-

T

J

#10

1586

-

-

-

-

-

-

-

G

-

-

-

-

-

-

-

-

-

H3

#11

5699

-

-

C

-

-

-

-

-

T

G

-

-

-

-

-

-

T

J

#12

5869

-

-

C

-

-

-

-

-

T

G

-

-

-

-

-

-

T

J

#13

G-495

-

-

-

-

-

-

-

-

T

-

-

-

-

G

-

-

T

U

HG = haplogroup

Other variants were also observed by sequencing analysis in m.11778G>A carriers. For instance, the transition A11812G is diagnostic of haplogroup T2 and G10172A identifies haplogroup J2b; these variants also appear sporadically in other haplogroup backgrounds. In addition, we observe a heteroplasmy at position 12705; transition C12705T is a well-known diagnostic site that (together with T16223C) leads from macro-haplogroup N to R.

Mixtures of two different DNA bearing different SNP profiles were carried out in order to simulate mtDNA heteroplasmy, a common state in mtDNA disease patients. As shown in Figure 2, the minisequencing assay is perfectly able to detect the mixtures at different proportions. According to [29] "the single most important disadvantage of the SNaPshot is that it seems less accurate than conventional PCR-RFLP analysis for [detecting] high levels of mutant mtDNA". Our results however indicate that the ability of SNaPshot for detecting heteroplasmy is actually equivalent to standard sequencing. Some fluorocromes 'project' more intensity in the electropherogram than other (as it also occurs with some sequencing chemistries), and this should be taking into account when quantifying the real proportions of wild:mutant DNA (heteroplasmy).
Figure 2

Partial electropherograms showing the performance of the SNaPshot assay with mixtures of wild:mutant DNAs. Numbers in the top indicate the percentages of wild:mutant mixture.

Although our primer design take into account the proximity of two pathologic mutations, it is not possible to consider all the potential polymorphisms that could hit within the annealing segments, and therefore alter the PCR efficacy (and in the worst of the cases, produce an artefact related to e.g. the spurious amplification of a NUMT; see below). This is however an universal problem affecting most of the genotyping techniques (including RFLP and sequencing) and therefore, the expert should be always alert under the possibility of false negative results. The experience of previous SNaPshot designs [16, 18, 30] indicates however that the method is robust under the presence of polymorphisms.

Finally, the multiplex reaction was also tested in 102 control individuals. There were not false positives in this control sample and the SNP-minisequencing profiles were phylogenetically concordant [31, 32].

Discussion

MtDNA mutations are responsible for various clinical features, which often make diagnosis a considerable challenge. Since the mitochondria constitute the energetic 'factory' of all nucleated cells, mtDNA diseases affect many tissues with variable clinical outcomes. When searching the Mitomap database for pathogenic mtDNA mutations, many of them harbour the status of 'provisional', while some other stated as 'confirmed'; among the later, some were recently questioned [2, 33, 34]. We here selected a total of 25 mutations, 21 of them were confirmed as pathogenic in common mtDNA diseases in multiple independent studies and considering different criteria, including the score system proposed by [33] (see also [35]) and also the complementary phylogenetic approach [2]. For instance, we included all the MTND gene mutations evaluated with the complex I pathogenicity scoring system of [33] (see their Table 2) and reaching the 16 top values (≥ 30). Apart from the 21 confirmed mutations included in the design, the remaining four mutations were additionally included because they occur at lower frequencies at four of the initially 21 selected positions.

Among the selected mutations, we have the transition m.3243A>G which is a typical cause of MELAS [36] and MIDD [37]; transition m.8993T>G, responsible of a number of patients suffering NARP [38] and MILS [39]; transitions m.3460G>A, m.11778G>A, m.14484T>C, which are common causal mutations in LHON patients [40, 41], etc. (see [42] for a review). It can be tentatively said that the selected SNPs actually cover an ample spectra of the fully confirmed causal mutations in the most common mtDNA diseases. For instance, according to [43], the prevalence of m.3243A>G is about 1 in 6,135 in the general (European) population; see also [44]. It is also worth to mention that it is straightforward to add mutations to our initial multiplex design because SNaPshot minisequencing is very flexible in this regard [17].

Conclusion

The minisequencing reaction presented here detected 11 mutations in an easy and straightforward manner in 11 different patients with clinical suspicion of Leber disease. All the mutations and their homoplasmic status were finally confirmed by automatic sequencing. We did not observed false positives in a substantial number of controls. We also simulated heteroplasmic states by mixing wild:mutant DNAs; the minisequencing technique demonstrates a good performance in detecting at least mixtures up to 1:10.

The multiplex reaction shows various advantages in the clinical field with respect to other techniques (see [17] for more details). For instance, the primers were designed in order to obtain amplicons of size ranging 89–196, a feature of special interest in a clinical context where often the laboratories have to deal with suboptimum samples containing low amounts of DNA or highly degraded DNA (paraphine imbibed samples, biopsies, etc.) as it is also the case with forensic samples [13]. Finally, we estimate that the cost of genotyping 25 SNPs is at least 20 times lower than the cost of RFLP genotyping (considering genotyping of one mutation at the time), and it can be done in less than three hours (including post-PCR purification, minisequencing reaction, purification of minisequencing products, electrophoresis, and documentation) for as many samples as the number of capillaries of the automatic sequencer.

Also important is the fact that minisequencing is robust respect potential artefacts related to false positives at NUMTs (author's unpublished data). Multiplex genotyping also prevents from artificial recombination; the later is a kind of common artefact that more easily arise as more independent fragments are genotyped for the same sample [31, 32, 4549].

Finally, the screening method presented here should not be considered as a substitute of complete genome sequencing. The latter should be performed when the minisequencing screening (or the use of other alternative strategies) fails to identify the causal mutation.

Availability and requirements

There are over 100 point mutations putatively associated with human mtDNA diseases [1]; however, only a small percentage of these mutations were properly confirmed http://www.mitomap.org/rimtab2.html. The primers both for PCR amplification (Table 2) and minisequencing reaction (Table 3) were designed to have an annealing temperature around 60°C using Primer3 software http://frodo.wi.mit.edu/cgi-bin/primer3. The sequence databases at the National Centre for Biotechnology Information (NCBI; http://www.ncbi.nlm.nih.gov) were interrogated using the online BLAST tool to test the primers against possible repetitive sequences and sequence homologies in the autosomal genome. Each primer pair for PCR amplification and each single base extension primer were selected independently and AutoDimer http://www.cstl.nist.gov/biotech/strbase/AutoDimerHomepage/AutoDimerProgramHomepage.htm was used to test for potential hairpin structures and primer-dimer problems.

Declarations

Acknowledgements

The 'Ramón y Cajal' Spanish programme from the Ministerio de Educación y Ciencia (RYC2005-3), the grant from the Xunta de Galicia (PGIDIT06PXIB208079PR), and a grant from the Fundación de Investigación Médica Mutua Madrileña awarded to AS supported this project.

Authors’ Affiliations

(1)
Unidade de Xenética, Instituto de Medicina Legal, Facultad de Medicina, Universidad de Santiago de Compostela
(2)
Fundación Pública Galega de Medicina Xenómica (FPGMX), Hospital Clínico Universitario, Universidad de Santiago de Compostela

References

  1. Schapira AH: Mitochondrial disease. Lancet. 2006, 368 (9529): 70-82. 10.1016/S0140-6736(06)68970-8.View ArticlePubMedGoogle Scholar
  2. Bandelt HJ, Yao YG, Salas A, Kivisild T, Bravi CM: High penetrance of sequencing errors and interpretative shortcomings in mtDNA sequence analysis of LHON patients. Biochem Biophys Res Commun. 2007, 352: 283-291. 10.1016/j.bbrc.2006.10.131.View ArticlePubMedGoogle Scholar
  3. Yao YG, Salas A, Bravi CM, Bandelt HJ: A reappraisal of complete mtDNA variation in East Asian families with hearing impairment. Hum Genet. 2006, 119 (5): 505-515. 10.1007/s00439-006-0154-9.View ArticlePubMedGoogle Scholar
  4. Bandelt HJ, Yao YG, Salas A: The search of ‘novel’ mtDNA mutations in hypertrophic cardiomyopathy: MITOMAPping as a risk factor. Int J Cardiol. 2007.Google Scholar
  5. Kong QP, Bandelt HJ, Sun C, Yao YG, Salas A, Achilli A, Wang CY, Zhong L, Zhu CL, Wu SF, Torroni A, Zhang YP: Updating the East Asian mtDNA phylogeny: a prerequisite for the identification of pathogenic mutations. Hum Mol Genet. 2006, 15 (13): 2076-2086. 10.1093/hmg/ddl130.View ArticlePubMedGoogle Scholar
  6. Wallace DC, Brown MD, Lott MT: Mitochondrial DNA variation in human evolution and disease. Gene. 1999, 238 (1): 211-230. 10.1016/S0378-1119(99)00295-4.View ArticlePubMedGoogle Scholar
  7. Barros F, Lareu MV, Salas A, Carracedo A: Rapid and enhanced detection of mitochondrial DNA variation using single-strand conformation analysis of superposed restriction enzyme fragments from polymerase chain reaction-amplified products. Electrophoresis. 1997, 18 (1): 52-54. 10.1002/elps.1150180110.View ArticlePubMedGoogle Scholar
  8. Salas A, Rasmussen EM, Lareu MV, Morling N, Carracedo A: Fluorescent SSCP of overlapping fragments (FSSCP-OF): a highly sensitive method for the screening of mitochondrial DNA variation. Forensic Sci Int. 2001, 124 (2-3): 97-103.View ArticlePubMedGoogle Scholar
  9. Zeviani M, Gellera C, Antozzi C, Rimoldi M, Morandi L, Villani F, Tiranti V, DiDonato S: Maternally inherited myopathy and cardiomyopathy: association with mutation in mitochondrial DNA tRNA(Leu)(UUR). Lancet. 1991, 338 (8760): 143-147. 10.1016/0140-6736(91)90136-D.View ArticlePubMedGoogle Scholar
  10. DiMauro S: Mitochondrial DNA medicine. Biosci Rep. 2007, 27 (1-3): 5-9. 10.1007/s10540-007-9032-5.View ArticlePubMedGoogle Scholar
  11. Hudson G, Carelli V, Spruijt L, Gerards M, Mowbray C, Achilli A, Pyle A, Elson J, Howell N, La Morgia C, Valentino ML, Huoponen K, Savontaus ML, Nikoskelainen E, Sadun AA, Salomao SR, Belfort R, Griffiths P, Man PY, de Coo RF, Horvath R, Zeviani M, Smeets HJ, Torroni A, Chinnery PF: Clinical Expression of Leber Hereditary Optic Neuropathy Is Affected by the Mitochondrial DNA-Haplogroup Background. Am J Hum Genet. 2007, 81 (2): 228-233. 10.1086/519394.View ArticlePubMedPubMed CentralGoogle Scholar
  12. Bannwarth S, Procaccio V, Paquis-Flucklinger V: Surveyor Nuclease: a new strategy for a rapid identification of heteroplasmic mitochondrial DNA mutations in patients with respiratory chain defects. Hum Mutat. 2005, 25 (6): 575-582. 10.1002/humu.20177.View ArticlePubMedGoogle Scholar
  13. Crespillo M, Paredes MR, Prieto L, Montesino M, Salas A, Albarrán C, Álvarez-Iglesias V, Amorin A, Berniell-Lee G, Brehm A, Carril JC, Corach D, Cuevas N, Di Lonardo AM, Doutremepuich C, Espinheira RM, Espinoza M, Gómez F, González A, Hernández A, Hidalgo M, Jimenez M, Leite FP, López AM, López-Soto M, Lorente JA, Pagano S, Palacio AM, Pestano JJ, Pinheiro MF, Raimondi E, Ramon MM, Tovar F, Vidal-Rioja L, Vide MC, Whittle MR, Yunis JJ, Garcia-Hirschfel J: Results of the 2003-2004 GEP-ISFG collaborative study on mitochondrial DNA: focus on the mtDNA profile of a mixed semen-saliva stain. Forensic Sci Int. 2006, 160 (2-3): 157-167. 10.1016/j.forsciint.2005.09.005.View ArticlePubMedGoogle Scholar
  14. Brandstätter A, Parsons TJ, Parson W: Rapid screening of mtDNA coding region SNPs for the identification of west European Caucasian haplogroups. Int J Legal Med. 2003, 117 (5): 291-298. 10.1007/s00414-003-0395-2.View ArticlePubMedGoogle Scholar
  15. Bandelt HJ, Kivisild T, Parik J, Villems R, Bravi CM, Yao YG, Brandstätter A, Parson W: Lab-specific mutation processes. Human mitochondrial DNA and the evolution of Homo sapiens Series: Nucleic Acids and Molecular Biology. Edited by: H.-J. Bandelt MRVM. 2006, Berlin-Heidelberg , Springer-Verlag, 18: 119–150-Google Scholar
  16. Quintáns B, Álvarez-Iglesias V, Salas A, Phillips C, Lareu MV, Carracedo : Typing of mitochondrial DNA coding region SNPs of forensic and anthropological interest using SNaPshot minisequencing. Forensic Sci Int. 2004, 140 (2-3): 251-257. 10.1016/j.forsciint.2003.12.005.View ArticlePubMedGoogle Scholar
  17. Álvarez-Iglesias V, Jaime JC, Carracedo , Salas A: Coding region mitochondrial DNA SNPs: targeting East Asian and Native American haplogroups. Forensic Sci Int: Genet. 2007, 1: 44-55. 10.1016/j.fsigen.2006.09.001.View ArticleGoogle Scholar
  18. Brandstätter A, Salas A, Niederstätter H, Gassner C, Carracedo , Parson W: Dissection of mitochondrial superhaplogroup H using coding region SNPs. Electrophoresis. 2007, 2541-2550.Google Scholar
  19. Filippini S, Blanco A, Fernandez-Marmiesse A, Alvarez-Iglesias V, Ruiz-Ponte C, Carracedo A, Vega A: Multiplex SNaPshot for detection of BRCA1/2 common mutations in Spanish and Spanish related breast/ovarian cancer families. BMC Med Genet. 2007, 8: 40-10.1186/1471-2350-8-40.View ArticlePubMedPubMed CentralGoogle Scholar
  20. Salas A, Comas D, Lareu MV, Bertranpetit J, Carracedo : mtDNA analysis of the Galician population: a genetic edge of European variation. Eur J Hum Genet. 1998, 6 (4): 365-375. 10.1038/sj.ejhg.5200202.View ArticlePubMedGoogle Scholar
  21. Salas A, Lareu MV, Sánchez-Diz P, Calafell F, Carracedo A: mtDNA hypervariable region II (HVII) sequences in human evolution studies: impact of mutation rate heterogeneity. Progress in Forensic Genetics. 2000, 8: 329-331.Google Scholar
  22. Achilli A, Rengo C, Magri C, Battaglia V, Olivieri A, Scozzari R, Cruciani F, Zeviani M, Briem E, Carelli V, Moral P, Dugoujon JM, Roostalu U, Loogvali EL, Kivisild T, Bandelt HJ, Richards M, Villems R, Santachiara-Benerecetti AS, Semino O, Torroni A: The molecular dissection of mtDNA haplogroup H confirms that the Franco-Cantabrian glacial refuge was a major source for the European gene pool. Am J Hum Genet. 2004, 75 (5): 910-918. 10.1086/425590.View ArticlePubMedPubMed CentralGoogle Scholar
  23. Herrnstadt C, Elson JL, Fahy E, Preston G, Turnbull DM, Anderson C, Ghosh SS, Olefsky JM, Beal MF, Davis RE, Howell N: Reduced-median-network analysis of complete mitochondrial DNA coding-region sequences from the major African, Asian, and European haplogroups. Am J Hum Genet. 2002, 70: 1152-1171. 10.1086/339933.View ArticlePubMedPubMed CentralGoogle Scholar
  24. Fernández-Marmiesse A, Salas A, Vega A, Fernández-Lorenzo JR, Barreiro J, Carracedo : Mutation spectra of ABCC8 gene in Spanish patients with Hyperinsulinism of Infancy (HI). Hum Mutat. 2006, 27 (2): 214-10.1002/humu.9401.View ArticlePubMedGoogle Scholar
  25. Brage A, Tomé S, Garcia A, Carracedo , Salas A: Clinical and molecular characterization of Wilson disease in Spanish patients. Hepatol Res. 2007, 37 (1): 18-26.View ArticlePubMedGoogle Scholar
  26. Brown MD, Sun F, Wallace DC: Clustering of Caucasian Leber hereditary optic neuropathy patients containing the 11778 or 14484 mutations on an mtDNA lineage. Am J Hum Genet. 1997, 60 (2): 381-387.PubMedPubMed CentralGoogle Scholar
  27. Howell N, Oostra RJ, Bolhuis PA, Spruijt L, Clarke LA, Mackey DA, Preston G, Herrnstadt C: Sequence analysis of the mitochondrial genomes from Dutch pedigrees with Leber hereditary optic neuropathy. Am J Hum Genet. 2003, 72 (6): 1460-1469. 10.1086/375537.View ArticlePubMedPubMed CentralGoogle Scholar
  28. González AM, Brehm A, Pérez JA, Maca-Meyer N, Flores C, Cabrera VM: Mitochondrial DNA affinities at the Atlantic fringe of Europe. Am J Phys Anthropol. 2003, 120 (4): 391-404. 10.1002/ajpa.10168.View ArticlePubMedGoogle Scholar
  29. Cassandrini D, Calevo MG, Tessa A, Manfredi G, Fattori F, Meschini MC, Carrozzo R, Tonoli E, Pedemonte M, Minetti C, Zara F, Santorelli FM, Bruno C: A new method for analysis of mitochondrial DNA point mutations and assess levels of heteroplasmy. Biochem Biophys Res Commun. 2006, 342 (2): 387-393. 10.1016/j.bbrc.2006.01.152.View ArticlePubMedGoogle Scholar
  30. Álvarez-Iglesias V, Salas A, Cerezo M, Ramos-Luis E, Jaime JC, Lareu MV, Carracedo : Genotyping coding region mtDNA SNPs for Asian and Native American haplogroup assignation. Int Congress Series. 2006, 11 (1288): 4-6. 10.1016/j.ics.2005.11.036.View ArticleGoogle Scholar
  31. Bandelt HJ, Kong QP, Parson W, Salas A: More evidence for non-maternal inheritance of mitochondrial DNA?. J Med Genet. 2005, 42: 957-960. 10.1136/jmg.2005.033589.View ArticlePubMedPubMed CentralGoogle Scholar
  32. Bandelt HJ, Salas A, Bravi CM: Problems in FBI mtDNA database. Science. 2004, 305 (5689): 1402-1404. 10.1126/science.305.5689.1402b.View ArticlePubMedGoogle Scholar
  33. Mitchell AL, Elson JL, Howell N, Taylor RW, Turnbull DM: Sequence variation in mitochondrial complex I genes: mutation or polymorphism?. J Med Genet. 2006, 43 (2): 175-179. 10.1136/jmg.2005.032474.View ArticlePubMedPubMed CentralGoogle Scholar
  34. McFarland R, Taylor RW, Elson JL, Lightowlers RN, Turnbull DM, Howell N: Proving pathogenicity: when evolution is not enough. Am J Med Genet A. 2004, 131 (1): 107-8; author reply 109-10. 10.1002/ajmg.a.30318.View ArticlePubMedGoogle Scholar
  35. DiMauro S, Schon EA: Mitochondrial DNA mutations in human disease. Am J Med Genet. 2001, 106 (1): 18-26. 10.1002/ajmg.1392.View ArticlePubMedGoogle Scholar
  36. Goto Y, Nonaka I, Horai S: A mutation in the tRNA(Leu)(UUR) gene associated with the MELAS subgroup of mitochondrial encephalomyopathies. Nature. 1990, 348 (6302): 651-653. 10.1038/348651a0.View ArticlePubMedGoogle Scholar
  37. van den Ouweland JM, Lemkes HH, Ruitenbeek W, Sandkuijl LA, de Vijlder MF, Struyvenberg PA, van de Kamp JJ, Maassen JA: Mutation in mitochondrial tRNA(Leu)(UUR) gene in a large pedigree with maternally transmitted type II diabetes mellitus and deafness. Nat Genet. 1992, 1 (5): 368-371. 10.1038/ng0892-368.View ArticlePubMedGoogle Scholar
  38. Holt IJ, Harding AE, Petty RK, Morgan-Hughes JA: A new mitochondrial disease associated with mitochondrial DNA heteroplasmy. Am J Hum Genet. 1990, 46 (3): 428-433.PubMedPubMed CentralGoogle Scholar
  39. de Vries DD, van Engelen BG, Gabreels FJ, Ruitenbeek W, van Oost BA: A second missense mutation in the mitochondrial ATPase 6 gene in Leigh's syndrome. Ann Neurol. 1993, 34 (3): 410-412. 10.1002/ana.410340319.View ArticlePubMedGoogle Scholar
  40. Howell N, Bindoff LA, McCullough DA, Kubacka I, Poulton J, Mackey D, Taylor L, Turnbull DM: Leber hereditary optic neuropathy: identification of the same mitochondrial ND1 mutation in six pedigrees. Am J Hum Genet. 1991, 49 (5): 939-950.PubMedPubMed CentralGoogle Scholar
  41. Johns DR, Neufeld MJ, Park RD: An ND-6 mitochondrial DNA mutation associated with Leber hereditary optic neuropathy. Biochem Biophys Res Commun. 1992, 187 (3): 1551-1557. 10.1016/0006-291X(92)90479-5.View ArticlePubMedGoogle Scholar
  42. Taylor RW, Turnbull DM: Mitochondrial DNA mutations in human disease. Nat Rev Genet. 2005, 6 (5): 389-402. 10.1038/nrg1606.View ArticlePubMedPubMed CentralGoogle Scholar
  43. Majamaa K, Moilanen JS, Uimonen S, Remes AM, Salmela PI, Karppa M, Majamaa-Voltti KA, Rusanen H, Sorri M, Peuhkurinen KJ, Hassinen IE: Epidemiology of A3243G, the mutation for mitochondrial encephalomyopathy, lactic acidosis, and strokelike episodes: prevalence of the mutation in an adult population. Am J Hum Genet. 1998, 63 (2): 447-454. 10.1086/301959.View ArticlePubMedPubMed CentralGoogle Scholar
  44. Schaefer AM, Taylor RW, Turnbull DM, Chinnery PF: The epidemiology of mitochondrial disorders--past, present and future. Biochim Biophys Acta. 2004, 1659 (2-3): 115-120. 10.1016/j.bbabio.2004.09.005.View ArticlePubMedGoogle Scholar
  45. Bandelt HJ, Quintana-Murci L, Salas A, Macaulay V: The fingerprint of phantom mutations in mitochondrial DNA data. Am J Hum Genet. 2002, 71 (5): 1150-1160. 10.1086/344397.View ArticlePubMedPubMed CentralGoogle Scholar
  46. Bandelt HJ, Salas A, Lutz-Bonengel S: Artificial recombination in forensic mtDNA population databases. Int J Legal Med. 2004, 118 (5): 267-273. 10.1007/s00414-004-0455-2.View ArticlePubMedGoogle Scholar
  47. Salas A, Bandelt HJ, Macaulay V, Richards MB: Phylogeographic investigations: The role of trees in forensic genetics. Forensic Sci Int. 2007, 168: 1-13. 10.1016/j.forsciint.2006.05.037.View ArticlePubMedGoogle Scholar
  48. Salas A, Carracedo , Macaulay V, Richards M, Bandelt HJ: A practical guide to mitochondrial DNA error prevention in clinical, forensic, and population genetics. Biochem Biophys Res Commun. 2005, 335 (3): 891-899. 10.1016/j.bbrc.2005.07.161.View ArticlePubMedGoogle Scholar
  49. Salas A, Yao YG, Macaulay V, Vega A, Carracedo , Bandelt HJ: A critical reassessment of the role of mitochondria in tumorigenesis. PLoS Med. 2005, 2 (11): e296-10.1371/journal.pmed.0020296.View ArticlePubMedPubMed CentralGoogle Scholar
  50. Andrews RM, Kubacka I, Chinnery PF, Lightowlers RN, Turnbull DM, Howell N: Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA. Nat Genet. 1999, 23: 147-10.1038/13779.View ArticlePubMedGoogle Scholar
  51. Pre-publication history

    1. The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/9/26/prepub

Copyright

© Álvarez-Iglesias et al; licensee BioMed Central Ltd. 2008

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.