Research article | Open | Open Peer Review | Published:
Mitochondrial DNA variant m.15218A > G in Finnish epilepsy patients who have maternal relatives with epilepsy, sensorineural hearing impairment or diabetes mellitus
BMC Medical Geneticsvolume 14, Article number: 73 (2013)
Mitochondrial diseases caused by mutations in mitochondrial DNA (mtDNA) affect tissues with high energy demand. Epilepsy is one of the manifestations of mitochondrial dysfunction when the brain is affected. We have studied here 79 Finnish patients with epilepsy and who have maternal first- or second-degree relatives with epilepsy, sensorineural hearing impairment or diabetes mellitus.
The entire mtDNA was studied by using conformation sensitive gel electrophoresis and PCR fragments that differed in mobility were directly sequenced.
We found a common nonsynonymous variant m.15218A > G (p.T158A, MTCYB) that occurs in haplogroup U5a1 to be more frequent in patients with epilepsy. The m.15218A > G variant was present in five patients with epilepsy and in four out of 403 population controls (p = 0.0077). This variant was present in two branches in the phylogenetic network constructed on the basis of mtDNA variation among the patients. Three algorithms predicted that m.15218A > G is damaging in effect.
We suggest that the m.15218A > G variant is mildly deleterious and that mtDNA involvement should be considered in patients with epilepsy and who have a maternal history of epilepsy, sensorineural hearing impairment or diabetes mellitus.
Epilepsies are a heterogeneous group of complex seizure disorders, which can be classified according to underlying pathology or clinical manifestation. Main classes of epilepsy in adults are localization related (focal, partial), generalized, undetermined and epilepsies associated with special syndromes. Causes for epilepsy include structural damage of brain tissue such as trauma, stroke and tumours, central nervous system infections and genetic causes [1, 2]. Some epilepsy cases remain etiologically undetermined.
Epilepsy is a frequent manifestation of mitochondrial disorders, as the brain is often severely affected in these diseases . Seizures are commonly seen in Alpers-Huttenlocher syndrome, Leigh syndrome, MELAS syndrome (mitochondrial encephalomyopathy, lactic acidosis and stroke-like episodes), MEMSA syndrome (myoclonic epilepsy, myopathy, sensory ataxia) and in ataxia neuropathy spectrum [4, 5]. The most common mitochondrial DNA (mtDNA) mutation causing epilepsy is m.8344A > G in the MTTK gene. It causes MERRF syndrome (myoclonus epilepsy associated with ragged red fibers) . Mutations in the MTTF gene have also been reported to cause severe epilepsy . Pathogenic mtDNA mutations are often heteroplasmic so that both the mutant variant and the wild type variant are present in the mitochondria. The proportion of the mutant variant must exceed a certain threshold before symptoms of mitochondrial disease become manifest . Other mitochondrial diseases which can present with seizure symptoms are neuropathy, ataxia, retinitis pigmentosa syndrome and Kearns-Sayre syndrome . Furthermore, mitochondrial oxidative stress may play a role in epileptogenesis by virtue of affecting neuronal excitability .
It has been suggested that mildly deleterious mtDNA variants are involved in the pathogenesis of many degenerative diseases [10–12]. Such variants are present in the population without being subjected to selection , but are maintained more by genetic drift . The increased risk for diseases with mitochondrial contribution in pathogenesis is most likely determined by a balance between beneficial and damaging mtDNA mutations within haplogroups. In the case of some haplogroups, the balance leans more towards the deleterious side. Also, deleterious mtDNA mutations have been suggested to be population specific so that a given variant may be associated with a disease in one population, but not in the other .
Deleterious mtDNA mutations are the cause of the MELAS syndrome  and the MERRF syndrome . In addition to the full-blown encephalomyopathic MELAS syndrome, patients with the m.3243A > G mutation often present with diabetes mellitus or sensorineural hearing impairment (SNHI) . These two phenotypes as well as epilepsy are common in the population and, therefore, patients with these diseases are good subjects in the search of mildly deleterious mtDNA variants. In order to assess the role of such variants and mtDNA haplogroups we determined the entire mtDNA sequence of 79 Finnish patients with epilepsy and with a maternal history of epilepsy, diabetes mellitus or SNHI. The pathogenic potential of all nonsynonymous mutations was also analysed.
Subjects and samples
Most adult patients with epilepsy in the province of Northern Ostrobothnia make regular visits to the outpatient Neurology Department of the Oulu University Hospital. During a one year period, a physician involved in the study reviewed the charts of the visiting patients and once the diagnosis of epilepsy was confirmed, the patient was requested to fill out a family history questionnaire. We did not make distinctions between the etiologies or types of epilepsy. We identified 223 patients with epilepsy that had any combination of epilepsy, sensorineural hearing impairment or diabetes mellitus in their first- or second-degree maternal relatives. Blood samples were obtained from 165 patients and then 79 patients were selected if the number of maternal relatives with epilepsy, diabetes or sensorineural hearing impairment was ≥ 2 (38 patients) or if a family history score was ≥ 0.1 in families with one maternal relative with epilepsy or sensorineural hearing impairment (41 patients). The family history score was calculated using the formula Naffected/Ntotal. Patients were tested not to harbour m.3243A > G or m.8344A > G mutation using restriction fragment analysis prior to selecting samples for this study.
Samples from healthy blood donors (N = 403) obtained at local Finnish Red Cross offices served as population controls . It was required that the donors did not report neurological ailments, diabetes and sensorineural hearing impairment of their own or of their mothers. It was also required that the donor and the mother had been born in the same province. All the patients signed an informed consent to participate in the study. The Ethics Committee of the Medical Faculty, University of Oulu, and that of Finnish Red Cross have approved the study protocol.
Conformation Sensitive Gel Electrophoresis (CSGE)
Total DNA was extracted from blood using the QIAmp Blood Kit (Qiagen, Valencia, CA, U.S.A.). Mitochondrial DNA haplogroups were determined by restriction fragment analysis . CSGE was performed as described earlier [21, 22]. MtDNA coding region (nucleotides m.577-16090) was amplified in 63 partially overlapping fragments covering the entire mtDNA coding region. The primers were ~22 bp in length and they were designed as previously reported . The mean PCR fragment size was 350 bp and each amplified fragment overlapped by about 80 nucleotides with its neighbouring fragments [see Additional file 1]. PCR fragments were amplified in a total volume of 30 μl in 30 cycles: denaturation at 94˚C for 1 min, annealing at a primer-specific temperature and extension at 72˚C for 1 min and a final extension for 10 min. Touchdown-PCR protocol was used in parallel and yielded similar results. 3-10 μl of the PCR product was used for heteroduplex formation. Each amplified fragment was mixed with a reference sample and denatured at 95˚C for 5 min and annealed at 68˚C for 30 min for heteroduplex formation. Samples were electrophoresed on 15% polyacrylamide gel overnight at a constant voltage of 400 V at room temperature. The gel was stained in 150 μg/l of ethidium bromide for 5 min and destained in water after which it was transferred to UV transluminator and photographed (Grab-IT Annotating Grabber 2.04.7; UVP Inc. Upland, CA, U.S.A.). Previous data have suggested that the limit of detection of heteroplasmy in various heteroduplex screening assays is < 10% .
PCR fragments that differed in mobility in CSGE were purified with exonuclease I and shrimp alkaline phosphatase  and sequenced (ABI PRISM™ 377 and ABI PRISM™ 3100 Sequencers using DYEnamic ET Terminator Cycle Sequencing Kit; Amersham Pharmacia Biotech Inc. Piscataway, NJ, U.S.A.). Same primers were used for sequencing and amplification of the 63 PCR fragments. CSGE has previously been shown to be a sensitive (98.8%) and specific (100%) method for detecting mtDNA variants . The D-loop (nucleotides m.15975 – 725) was amplified in a unique PCR fragment and directly sequenced.
Analysis of substitutions
Sequences were compared to the revised Cambridge reference sequence NC_012920  and to mtDNA sequences available in the GiiB-JST mtSNP database [http://mtsnp.tmig.or.jp/mtsnp/index_e.shtml] , mtDB Human Mitochondrial Genome database [http://www.mtdb.igp.uu.se] , HmtDB database [http://www.hmtdb.uniba.it]  and Mitomap [http://www.mitomap.org]  accessed in June 2012. Nonsynonymous polymorphisms were considered rare if a total of three or less sequences had been reported in the above mentioned databases. Sequences were also compared to complete mtDNA sequences of 192 population controls obtained from the same area as the patients [19, 31].
Novel substitutions were confirmed by restriction fragment analysis and/or sequencing in both directions at least twice from separate PCR products. Variant was regarded as a previously reported pathogenic mutation if it was listed as such in the Mitomap. Base conservation was determined using the GiiB-JST mtSAP evaluation http://mtsnp.tmig.or.jp/cgi-bin/mtsnp/specAlign/ctrlSpecAlignE.cgi for nonsynonymous variants  and Mamit-tRNA for tRNA variants http://mamit-trna.u-strasbg.fr/. Phylogenetic networks based on mtDNA sequence variation in the coding region and the D-loop were constructed using the median algorithm .
Analysis of pathogenic potential of nonsynonymous mtDNA variants
The pathogenic potential of nonsynonymous mtDNA variants detected in epilepsy patients was analyzed with PolyPhen-2, version 2.2.2 [http://genetics.bwh.harvard.edu/pph2/] , SIFT BLink [http://sift.jcvi.org/www/SIFT_BLink_submit.html]  and PMut [http://mmb.pcb.ub.es] . PolyPhen-2 uses a Bayesian algorithm to calculate the pathogenic potential of a nonsynonymous mutation using both evolutionary conservation data and the structural similarities/differences of amino acids. The mutations are classified as benign (< 50% chance of pathogenic effect), possibly deleterious (> 50% chance of pathogenic effect) and probably deleterious (> 90% chance of pathogenic effect). SIFT BLink uses evolutionary conservation data to assess whether a nonsynonymous variant is pathogenic. It predicts a variant to be either functionally important (0 being the most pathogenic) or unimportant (1 being the most benign). PMut is based on evolutionary information derived from multiple sequence alignments. Variants are classified as pathogenic or benign and a reliability score is given (0, lowest reliable; 9, very reliable).
Nonsynonymous variants predicted to be damaging
CSGE and sequencing were used to determine complete mtDNA sequences of 79 patients with epilepsy. MtDNA sequences were then used to construct a phylogenetic network separately for the coding region (Figure 1) and the D-loop (Figure 2). Exact test of population differentiation  revealed no differences in the frequencies of mtDNA haplogroups between patients and controls (Table 1).
We found 52 different nonsynonymous mtDNA variants in the 79 patients. Evaluation of the pathogenic potential of all the nonsynonymous variants was then carried out by using PolyPhen-2, SIFT BLink and PMut. PolyPhen-2 predicted the lowest and SIFT BLink the highest number of pathogenic mutations (Figure 3). All three algorithms concordantly predicted five nonsynonymous variants to be damaging (Table 2). One of these variants was m.15218A > G (p.T158A, MTCYB), which was present in two mtDNA haplotypes in five patients (Figure 4). Four out of 403 population controls belonged to the major haplotype while none belonged to the minor haplotype (p = 0.0077 for difference between cases and controls, Fisher’s exact two-tailed test). The presence of m.15218A > G in the five patients and four controls and the absence of the variant in samples 40 and 48 was confirmed by sequencing the fragment amplified in two separate reactions in both L- and H-directions. Patients 37 and 75 with m.15218A > G had identical mtDNA sequences, but family information did not suggest immediate relatedness. The clinical features of the five patients were variable (Table 3). Another variant predicted to be damaging was m.9903T > C (p.F233L, MTCO3) that occurred in two patients but not in the controls (p = 0.08, Fisher’s exact two-tailed test). The remaining three variants did not differ in frequency between patients and controls.
Nonsynonymous variants with previously reported disease associations
We observed five nonsynonymous mtDNA variants with previously reported disease associations. Three of them were situated in the MTCO1 gene. The m.6489C > A transversion (p.L196I, MTCO1) has been reported in a patient with therapy-resistant epilepsy  and in patients with matrilineal diabetes . The second variant m.6253T > C (p.M117T, MTCO1) has been reported with an association to prostate cancer  and the third variant m.7444G > A (p.*514K, MTCO1) in association with Leber’s hereditary optic neuropathy or deafness [44–46]. The remaining two variants were located in MTND5 including m.12811T > C (p.Y159H) and m.13637A > G (p.Q434R). These variants have been postulated to affect the phenotype of Leber’s hereditary optic neuropathy [47–49]. Prediction algorithms did not support a pathogenic role for these variants.
Rare mtDNA variants
Four nonsynonymous polymorphisms found in the patients were considered rare (Table 4). Three of them have been previously reported in Finnish individuals. Two of the rare nonsynonymous mutations, m.8923A > G and m.15725C > T, were predicted to be damaging with SIFT BLink and PMut, but not with PolyPhen-2. In addition, a novel m.11392A > G synonymous change in MTND4 was also discovered in a sample belonging to haplogroup K.
The m.15218A > G (p.T158A, MTCYB) variant was found in five patients with epilepsy and with a maternal history of epilepsy, SNHI or diabetes. The patients belonged to haplogroup U5a. Haplogroup U5 has previously been suggested to have possible disease associations, as it has been found to be more frequent in patients with occipital stroke and migrainous stroke . Stroke-like episodes in mitochondrial diseases can resemble occipital epileptic seizures  and mitochondrial diseases often manifest themselves only in epileptic seizures . Haplogroup U5a has been associated with a rapid progression to AIDS and death among HIV-1 infected patients, especially those harbouring the m.15218A > G variant .
Haplogroup U5 is more common in Finland than the rest of Europe. Its frequency is ~ 30% in the Finnish population compared to 18 - 22% in Western and Eastern Europe . Haplogroup U5 is especially common in Northern Finland and its variation is low there. The reason for this is suspected to be a founder effect of a relatively small settler group that colonized northern Finland after the 16th century . Furthermore, haplogroup U constitutes 32-52% of mtDNA in the indigenous Saami people of Northern Scandinavia and Northwestern Russia. Admixture of Finns and the Saami has been reported in the North of Finland [19, 53]. In order to avoid bias caused by geographical differences in haplogroup frequencies we collected our patient and control samples from the same area.
We used three different methods for predicting the pathogenic potential of nonsynonymous mutations. PolyPhen-2 predicted the highest number of benign polymorphisms, whilst SIFT BLink predicted the highest number of pathogenic mutations. PolyPhen-2 has been reported to be sensitive in predicting neutral polymorphisms  with a 5% false negative rate for benign variants. On the other hand, SIFT BLink has been reported to be more reliable in detecting truly pathogenic variants. Both PMut and SIFT BLink make use of evolutionary conservation data and yield comparable results, while PolyPhen-2 makes use of physical and chemical properties of amino acids in addition to evolutionary conservation data. We considered the variants that all three algorithms predicted pathogenic to be the best candidates for pathogenic mutations.
The m.15218A > G mutation in MTCYB was predicted to be pathogenic by the three algorithms. This mutation leads to substitution of a conserved hydrophilic threonine by a hydrophobic alanine at position 159 of the cytochrome b subunit of complex III. Out of the 60 species compared, 53 species harbour threonine in this position  and only Barbary macaque (Macaca sylvanus) harbours alanine in this position of the subunit. The cytochrome b subunit of complex III forms the catalytic core of cytochrome bc1 complex along with cytochrome c1 and Rieske iron-sulphur protein. Changes in this subunit have been shown to cause changes to the catalytic function of complex III, thus leading to complex III defects .
The m.15218A > G variant is present in haplogroups M7, M10, HV, H13 and U5a1  indicating that m.15218A > G is a homoplasic variant. We found this variant in two haplotypes within haplogroup U5a1, the major haplotype being identical to that detected previously in Finnish population controls  and the minor haplotype being identical to that in a Finnish patient with diabetes mellitus . Evolutionary conservation of the amino acid position, results from the prediction algorithms and the homoplasic nature of the variant suggest that m.15218A > G is mildly deleterious rather than a neutral variant.
The m.14198G > A variant (p.T159M, MTND6) found in one patient was also predicted to be damaging. This variant has been found to exist in haplogroups G and H [28, 29], although in our patient data it was present in haplogroup U5a1. PMut predicted it to be pathogenic but the reliability score was low and, therefore, we cannot make conclusions of its true pathogenic potential. It is noteworthy though that it was also found in haplogroup U5a1.
Other mutations predicted to be pathogenic by the three algorithms were m.9903T > C (p.F233L, MTCO3) and m.9055G > A (p.A177T, MTATP6). The m.9033T > C variant occurred in two patients belonging to haplogroup K but was not present in the controls. It is present in one haplogroup L2a sequence [GenBank: EF657335], in two haplogroup K1c1 sequences [GenBank: EU262720, GenBank: EU753433] and in a haplogroup X2 sequence  in the HmtDB and mtDB databases [28, 29]. We suggest that m.9033T > C is a rare homoplasic variant occurring in haplogroups K, X and L and that its pathogenic nature is uncertain. The five nonsynonymous mutations with previous disease association were not predicted to be deleterious but, interestingly, m.6489C > A (p.L196I, MTCO1) has been reported to be associated with therapy-resistant epilepsy .
The m.9055G > A (p.A177T, MTATP6) mutation is the defining variant for haplogroup K. It was predicted to be damaging by the three algorithms, but was present both in the patients and in the controls. It has been suspected to be deleterious, as an excess of haplogroup K patients has been found to convey a greater risk to lipodystrophy among HIV-patients on antiretroviral therapy . It is interesting that both m.9055G > A and m.9903T > C are associated with haplogroup K. This haplogroup might harbour several mtDNA variants that could impair mitochondrial function.
In this study, we found that the nonsynonymous mtDNA variant m.15218A > G (p.T158A, MTCYB) is more common in epilepsy patients with maternal family members with epilepsy, SNHI or diabetes mellitus. The m.15218A > G variant has previously been suggested to be mildly deleterious, which was supported by our results. We suggest that mildly deleterious mtDNA mutations, such as m.15218A > G, may be considered a risk factor for epilepsy in patients with maternal relatives with SNHI, epilepsy or diabetes mellitus.
Sensorineural hearing impairment
Noncoding displacement loop in the mitochondrial DNA
Conformation sensitive gel electrophoresis.
Patel M: Mitochondrial dysfunction and oxidative stress: cause and consequence of epileptic seizures. Free Radic Biol Med. 2004, 37 (12): 1951-1962. 10.1016/j.freeradbiomed.2004.08.021.
Turnbull J, Lohi H, Kearney JA, Rouleau GA, Delgado-Escueta AV, Meisler MH, Cossette P, Minassian BA: Sacred disease secrets revealed: the genetics of human epilepsy. Hum Mol Genet. 2005, 14 Spec No. 2: 2491-2500.
Federico A, Cardaioli E, Da Pozzo P, Formichi P, Gallus GN, Radi E: Mitochondria, oxidative stress and neurodegeneration. J Neurol Sci. 2012, 322 (1–2): 254-262.
Finsterer J, Zarrouk Mahjoub S: Epilepsy in mitochondrial disorders. Seizure. 2012, 21 (5): 316-321. 10.1016/j.seizure.2012.03.003.
Rahman S: Mitochondrial disease and epilepsy. Dev Med Child Neurol. 2012, 54 (5): 397-406. 10.1111/j.1469-8749.2011.04214.x.
Shoffner JM, Lott MT, Lezza AM, Seibel P, Ballinger SW, Wallace DC: Myoclonic epilepsy and ragged-red fiber disease (MERRF) is associated with a mitochondrial DNA tRNA(Lys) mutation. Cell. 1990, 61 (6): 931-937. 10.1016/0092-8674(90)90059-N.
Zsurka G, Kunz WS: Mitochondrial dysfunction in neurological disorders with epileptic phenotypes. J Bioenerg Biomembr. 2010, 42 (6): 443-448. 10.1007/s10863-010-9314-7.
Chinnery PF, Elliott HR, Hudson G, Samuels DC, Relton CL: Epigenetics, epidemiology and mitochondrial DNA diseases. Int J Epidemiol. 2012, 41 (1): 177-187. 10.1093/ije/dyr232.
Waldbaum S, Patel M: Mitochondria, oxidative stress, and temporal lobe epilepsy. Epilepsy Res. 2010, 88 (1): 23-45. 10.1016/j.eplepsyres.2009.09.020.
Lehtonen MS, Moilanen JS, Majamaa K: Increased variation in mtDNA in patients with familial sensorineural hearing impairment. Hum Genet. 2003, 113 (3): 220-227. 10.1007/s00439-003-0966-9.
Coskun P, Wyrembak J, Schriner SE, Chen HW, Marciniack C, Laferla F, Wallace DC: A mitochondrial etiology of Alzheimer and Parkinson disease. Biochim Biophys Acta. 2012, 1820 (5): 553-564. 10.1016/j.bbagen.2011.08.008.
Poulton J, Luan J, Macaulay V, Hennings S, Mitchell J, Wareham NJ: Type 2 diabetes is associated with a common mitochondrial variant: evidence from a population-based case-control study. Hum Mol Genet. 2002, 11 (13): 1581-1583. 10.1093/hmg/11.13.1581.
Moilanen JS, Majamaa K: Phylogenetic network and physicochemical properties of nonsynonymous mutations in the protein-coding genes of human mitochondrial DNA. Mol Biol Evol. 2003, 20 (8): 1195-1210. 10.1093/molbev/msg121.
Breen MS, Kondrashov FA: Mitochondrial pathogenic mutations are population-specific. Biol Direct. 2010, 5: 68-10.1186/1745-6150-5-68.
Majamaa K, Moilanen JS, Uimonen S, Remes AM, Salmela PI, Kärppä 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.
Berkovic SF, Shoubridge EA, Andermann F, Andermann E, Carpenter S, Karpati G: Clinical spectrum of mitochondrial DNA mutation at base pair 8344. Lancet. 1991, 338 (8764): 457-
Sproule DM, Kaufmann P: Mitochondrial encephalopathy, lactic acidosis, and strokelike episodes: basic concepts, clinical phenotype, and therapeutic management of MELAS syndrome. Ann N Y Acad Sci. 2008, 1142: 133-158. 10.1196/annals.1444.011.
van Esch A, Steyerberg EW, Berger MY, Offringa M, Derksen-Lubsen G, Habbema JD: Family history and recurrence of febrile seizures. Arch Dis Child. 1994, 70 (5): 395-399. 10.1136/adc.70.5.395.
Meinilä M, Finnilä S, Majamaa K: Evidence for mtDNA admixture between the Finns and the Saami. Hum Hered. 2001, 52 (3): 160-170. 10.1159/000053372.
Torroni A, Huoponen K, Francalacci P, Petrozzi M, Morelli L, Scozzari R, Obinu D, Savontaus ML, Wallace DC: Classification of European mtDNAs from an analysis of three European populations. Genetics. 1996, 144 (4): 1835-1850.
Finnilä S, Hassinen IE, Ala-Kokko L, Majamaa K: Phylogenetic network of the mtDNA haplogroup U in Northern Finland based on sequence analysis of the complete coding region by conformation-sensitive gel electrophoresis. Am J Hum Genet. 2000, 66 (3): 1017-1026. 10.1086/302802.
Körkkö J, Annunen S, Pihlajamaa T, Prockop DJ, Ala-Kokko L: Conformation sensitive gel electrophoresis for simple and accurate detection of mutations: comparison with denaturing gradient gel electrophoresis and nucleotide sequencing. Proc Natl Acad Sci USA. 1998, 95 (4): 1681-1685. 10.1073/pnas.95.4.1681.
Wong LJ, Boles RG: Mitochondrial DNA analysis in clinical laboratory diagnostics. Clin Chim Acta. 2005, 354 (1–2): 1-20.
Werle E, Schneider C, Renner M, Volker M, Fiehn W: Convenient single-step, one tube purification of PCR products for direct sequencing. Nucleic Acids Res. 1994, 22 (20): 4354-4355. 10.1093/nar/22.20.4354.
Hinttala R, Smeets R, Moilanen JS, Ugalde C, Uusimaa J, Smeitink JA, Majamaa K: Analysis of mitochondrial DNA sequences in patients with isolated or combined oxidative phosphorylation system deficiency. J Med Genet. 2006, 43 (11): 881-886. 10.1136/jmg.2006.042168.
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 (2): 147-10.1038/13779.
Tanaka M, Cabrera VM, Gonzalez AM, Larruga JM, Takeyasu T, Fuku N, Guo LJ, Hirose R, Fujita Y, Kurata M, Shinoda K, Umetsu K, Yamada Y, Oshida Y, Sato Y, Hattori N, Mizuno Y, Arai Y, Hirose N, Ohta S, Ogawa O, Tanaka Y, Kawamori R, Shamoto-Nagai M, Maruyama W, Shimokata H, Suzuki R, Shimodaira H: Mitochondrial genome variation in eastern Asia and the peopling of Japan. Genome Res. 2004, 14 (10A): 1832-1850. 10.1101/gr.2286304.
Ingman M, Gyllensten U: mtDB: Human Mitochondrial Genome Database, a resource for population genetics and medical sciences. Nucleic Acids Res. 2006, 34 (Database issue): D749-D751.
Rubino F, Piredda R, Calabrese FM, Simone D, Lang M, Calabrese C, Petruzzella V, Tommaseo-Ponzetta M, Gasparre G, Attimonelli M: HmtDB, a genomic resource for mitochondrion-based human variability studies. Nucleic Acids Res. 2012, 40 (Database issue): D1150-D1159.
Ruiz-Pesini E, Lott MT, Procaccio V, Poole JC, Brandon MC, Mishmar D, Yi C, Kreuziger J, Baldi P, Wallace DC: An enhanced MITOMAP with a global mtDNA mutational phylogeny. Nucleic Acids Res. 2007, 35 (Database issue): D823-D828.
Finnilä S, Lehtonen MS, Majamaa K: Phylogenetic network for European mtDNA. Am J Hum Genet. 2001, 68 (6): 1475-1484. 10.1086/320591.
Putz J, Dupuis B, Sissler M, Florentz C: Mamit-tRNA, a database of mammalian mitochondrial tRNA primary and secondary structures. RNA. 2007, 13 (8): 1184-1190. 10.1261/rna.588407.
Bandelt HJ, Forster P, Sykes BC, Richards MB: Mitochondrial portraits of human populations using median networks. Genetics. 1995, 141 (2): 743-753.
Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, Kondrashov AS, Sunyaev SR: A method and server for predicting damaging missense mutations. Nat Methods. 2010, 7 (4): 248-249. 10.1038/nmeth0410-248.
Ng PC, Henikoff S: Predicting deleterious amino acid substitutions. Genome Res. 2001, 11 (5): 863-874. 10.1101/gr.176601.
Ferrer-Costa C, Orozco M, de la Cruz X: Sequence-based prediction of pathological mutations. Proteins. 2004, 57 (4): 811-819. 10.1002/prot.20252.
Excoffier L, Laval G, Schneider S: Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evol Bioinform Online. 2007, 1: 47-50.
Gasparre G, Porcelli AM, Bonora E, Pennisi LF, Toller M, Iommarini L, Ghelli A, Moretti M, Betts CM, Martinelli GN, Ceroni AR, Curcio F, Carelli V, Rugolo M, Tallini G, Romeo G: Disruptive mitochondrial DNA mutations in complex I subunits are markers of oncocytic phenotype in thyroid tumors. Proc Natl Acad Sci USA. 2007, 104 (21): 9001-9006. 10.1073/pnas.0703056104.
Seoane M, Mosquera-Miguel A, Gonzalez T, Fraga M, Salas A, Costoya JA: The mitochondrial genome is a "genetic sanctuary" during the oncogenic process. PLoS One. 2011, 6 (8): e23327-10.1371/journal.pone.0023327.
Berg AT, Berkovic SF, Brodie MJ, Buchhalter J, Cross JH, van Emde BW, Engel J, French J, Glauser TA, Mathern GW, Moshe SL, Nordli D, Plouin P, Scheffer IE: Revised terminology and concepts for organization of seizures and epilepsies: report of the ILAE Commission on Classification and Terminology, 2005-2009. Epilepsia. 2010, 51 (4): 676-685. 10.1111/j.1528-1167.2010.02522.x.
Varlamov DA, Kudin AP, Vielhaber S, Schroder R, Sassen R, Becker A, Kunz D, Haug K, Rebstock J, Heils A, Elger CE, Kunz WS: Metabolic consequences of a novel missense mutation of the mtDNA CO I gene. Hum Mol Genet. 2002, 11 (16): 1797-1805. 10.1093/hmg/11.16.1797.
Soini HK, Moilanen JS, Finnilä S, Majamaa K: Mitochondrial DNA sequence variation in Finnish patients with matrilineal diabetes mellitus. BMC Res Notes. 2012, 5 (1): 350-10.1186/1756-0500-5-350.
Petros JA, Baumann AK, Ruiz-Pesini E, Amin MB, Sun CQ, Hall J, Lim S, Issa MM, Flanders WD, Hosseini SH, Marshall FF, Wallace DC: mtDNA mutations increase tumorigenicity in prostate cancer. Proc Natl Acad Sci USA. 2005, 102 (3): 719-724. 10.1073/pnas.0408894102.
Yuan H, Qian Y, Xu Y, Cao J, Bai L, Shen W, Ji F, Zhang X, Kang D, Mo JQ, Greinwald JH, Han D, Zhai S, Young WY, Guan MX: Cosegregation of the G7444A mutation in the mitochondrial COI/tRNA(Ser(UCN)) genes with the 12S rRNA A1555G mutation in a Chinese family with aminoglycoside-induced and nonsyndromic hearing loss. Am J Med Genet A. 2005, 138A (2): 133-140. 10.1002/ajmg.a.30952.
Yuan H, Chen J, Liu X, Cheng J, Wang X, Yang L, Yang S, Cao J, Kang D, Dai P, Zhai S, Han D, Young WY, Guan MX: Coexistence of mitochondrial 12S rRNA C1494T and CO1/tRNA(Ser(UCN)) G7444A mutations in two Han Chinese pedigrees with aminoglycoside-induced and non-syndromic hearing loss. Biochem Biophys Res Commun. 2007, 362 (1): 94-100. 10.1016/j.bbrc.2007.07.161.
Kokotas H, Grigoriadou M, Yang L, Lodahl M, Rendtorff ND, Gyftodimou Y, Korres GS, Ferekidou E, Kandiloros D, Korres S, Tranebjaerg L, Guan MX, Petersen MB: Homoplasmy of the G7444A mtDNA and heterozygosity of the GJB2 c.35delG mutations in a family with hearing loss. Int J Pediatr Otorhinolaryngol. 2011, 75 (1): 89-94. 10.1016/j.ijporl.2010.10.016.
Cai W, Fu Q, Zhou X, Qu J, Tong Y, Guan MX: Mitochondrial variants may influence the phenotypic manifestation of Leber's hereditary optic neuropathy-associated ND4 G11778A mutation. J Genet Genomics. 2008, 35 (11): 649-655. 10.1016/S1673-8527(08)60086-7.
Huoponen K, Lamminen T, Juvonen V, Aula P, Nikoskelainen E, Savontaus ML: The spectrum of mitochondrial DNA mutations in families with Leber hereditary optic neuroretinopathy. Hum Genet. 1993, 92 (4): 379-384. 10.1007/BF01247339.
Yang J, Zhu Y, Chen L, Zhang H, Tong Y, Huang D, Zhang Z, Chen S, Han X, Ma X: Novel A14841G mutation is associated with high penetrance of LHON/C4171A family. Biochem Biophys Res Commun. 2009, 386 (4): 693-696. 10.1016/j.bbrc.2009.06.102.
Finnilä S, Hassinen IE, Majamaa K: Phylogenetic analysis of mitochondrial DNA in patients with an occipital stroke. Evaluation of mutations by using sequence data on the entire coding region. Mutat Res. 2001, 458 (1–2): 31-39.
Cesaroni E, Scarpelli M, Zamponi N, Polonara G, Zeviani M: Mitochondrial encephalomyopathy lactic acidosis and strokelike episodes mimicking occipital idiopathic epilepsy. Pediatr Neurol. 2009, 41 (2): 131-134. 10.1016/j.pediatrneurol.2009.02.018.
Hendrickson SL, Hutcheson HB, Ruiz-Pesini E, Poole JC, Lautenberger J, Sezgin E, Kingsley L, Goedert JJ, Vlahov D, Donfield S, Wallace DC, O'Brien SJ: Mitochondrial DNA haplogroups influence AIDS progression. AIDS. 2008, 22 (18): 2429-2439. 10.1097/QAD.0b013e32831940bb.
Sajantila A, Lahermo P, Anttinen T, Lukka M, Sistonen P, Savontaus ML, Aula P, Beckman L, Tranebjaerg L, Gedde-Dahl T, Issel-Tarver L, DiRienzo A, Pääbo S: Genes and languages in Europe: an analysis of mitochondrial lineages. Genome Res. 1995, 5 (1): 42-52. 10.1101/gr.5.1.42.
Wang J, Schmitt ES, Landsverk ML, Zhang VW, Li FY, Graham BH, Craigen WJ, Wong LJ: An integrated approach for classifying mitochondrial DNA variants: one clinical diagnostic laboratory's experience. Genet Med. 2012, 14 (6): 620-626. 10.1038/gim.2012.4.
Legros F, Chatzoglou E, Frachon P, Ogier De Baulny H, Laforet P, Jardel C, Godinot C, Lombes A: Functional characterization of novel mutations in the human cytochrome b gene. Eur J Hum Genet. 2001, 9 (7): 510-518. 10.1038/sj.ejhg.5200678.
van Oven M, Kayser M: Updated comprehensive phylogenetic tree of global human mitochondrial DNA variation. Hum Mutat. 2009, 30 (2): E386-E394. 10.1002/humu.20921.
De Luca A, Nasi M, Di Giambenedetto S, Cozzi-Lepri A, Pinti M, Marzocchetti A, Mussini C, Fabbiani M, Bracciale L, Cauda R, Cossarizza A: Mitochondrial DNA haplogroups and incidence of lipodystrophy in HIV-infected patients on long-term antiretroviral therapy. J Acquir Immune Defic Syndr. 2012, 59 (2): 113-120. 10.1097/QAI.0b013e31823daff3..
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/14/73/prepub
We wish to thank Ms. Pirjo Keränen, Ms. Anja Heikkinen and Ms. Irma Vuoti for their expert technical assistance. This work was supported by grants from the Academy of Finland (project numbers 107174 and 127764), the Sigrid Juselius Foundation and the Finnish Graduate School for Population Genetics.
The authors declare that they have no competing interests.
HKS carried out the laboratory experiments, result analysis and was responsible for writing the manuscript. JSM collected and edited the mtDNA sequence data into sequence files, submitted mtDNA sequences, contributed to the study plan and data analysis. TVK carried out part of the CSGE analysis. SF contributed to the study plan and data analysis. KM planned the study protocol, selected and collected the epilepsy patients and reviewed the manuscript critically. All authors read and approved the final manuscript.
Electronic supplementary material
About this article
- Mitochondrial DNA
- Mitochondrial haplogroups
- Sequence variation
- Nonsynonymous mutation
- Maternal inheritance
- Haplogroup U5a1