Thymidine Kinase 2 Deficiency-Induced mtDNA ...

2 downloads 272 Views 5MB Size Report
Mar 7, 2013 - 1 Division of Clinical Microbiology F-68, Karolinska Institutet, ... Sweden, 5 Division of Clinical Chemistry, Department of Laboratory Medicine, Lund University, Lund, ... Thymidine kinase 2 (TK2) deficiency in humans causes ...
Thymidine Kinase 2 Deficiency-Induced mtDNA Depletion in Mouse Liver Leads to Defect b-Oxidation Xiaoshan Zhou1., Kristina Kannisto2., Sophie Curbo1*., Ulrika von Do¨beln3, Kjell Hultenby4, Sindra Isetun3, Mats Ga˚fvels2,5, Anna Karlsson1 1 Division of Clinical Microbiology F-68, Karolinska Institutet, Department of Laboratory Medicine, Karolinska University Hospital, Huddinge, Sweden, 2 Division of Clinical Chemistry, C1-72, Karolinska Institutet, Department of Laboratory Medicine, Karolinska University Hospital, Huddinge, Sweden, 3 Division of Metabolic Diseases, Karolinska Institutet, Department of Laboratory Medicine, Karolinska University Hospital, Huddinge, Sweden, 4 Division of Clinical Research Center, Karolinska Institutet, Department of Laboratory Medicine, Karolinska University Hospital, Huddinge, Sweden, 5 Division of Clinical Chemistry, Department of Laboratory Medicine, Lund University, Lund, Sweden

Abstract Thymidine kinase 2 (TK2) deficiency in humans causes mitochondrial DNA (mtDNA) depletion syndrome. To study the molecular mechanisms underlying the disease and search for treatment options, we previously generated and described a TK2 deficient mouse strain (TK22/2) that progressively loses its mtDNA. The TK22/2 mouse model displays symptoms similar to humans harboring TK2 deficient infantile fatal encephalomyopathy. Here, we have studied the TK22/2 mouse model to clarify the pathological role of progressive mtDNA depletion in liver for the severe outcome of TK2 deficiency. We observed that a gradual depletion of mtDNA in the liver of the TK22/2 mice was accompanied by increasingly hypertrophic mitochondria and accumulation of fat vesicles in the liver cells. The levels of cholesterol and nonesterified fatty acids were elevated and there was accumulation of long chain acylcarnitines in plasma of the TK22/2 mice. In mice with hepatic mtDNA levels below 20%, the blood sugar and the ketone levels dropped. These mice also exhibited reduced mitochondrial b-oxidation due to decreased transport of long chain acylcarnitines into the mitochondria. The gradual loss of mtDNA in the liver of the TK22/2 mice causes impaired mitochondrial function that leads to defect b-oxidation and, as a result, insufficient production of ketone bodies and glucose. This study provides insight into the mechanism of encephalomyopathy caused by TK2 deficiency-induced mtDNA depletion that may be used to explore novel therapeutic strategies. Citation: Zhou X, Kannisto K, Curbo S, von Do¨beln U, Hultenby K, et al. (2013) Thymidine Kinase 2 Deficiency-Induced mtDNA Depletion in Mouse Liver Leads to Defect b-Oxidation. PLoS ONE 8(3): e58843. doi:10.1371/journal.pone.0058843 Editor: David C. Samuels, Vanderbilt University Medical Center, United States of America Received November 30, 2012; Accepted February 7, 2013; Published March 7, 2013 Copyright: ß 2013 Zhou et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Funding: This work was supported by grants from the Swedish Research Council [A.K. grant number 2010–2828] and the Medical Faculty of Karolinska Institutet [A.K., S.C.]. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing Interests: The authors have declared that no competing interests exist. * E-mail: [email protected] . These authors contributed equally to this work.

associated mtDNA depletion in brain and liver in addition to muscle, confirming that TK2 mutations can cause multi-organ involvement [9]. The different phenotypes of TK2 and DGUOK are probably due to lack of enzymes with overlapping substrate specificities in the affected tissues or poor residual activity in tissues with high demand [14,15,16]. So far, there has been no human case of complete lack of functional TK2 protein, suggesting that residual TK2 activity might be required for human life. The phenotype in patients seems to correspond to the severity of the TK2 deficiency, where partial reductions of TK2 activity (14–45% of normal) lead to myopathy and severe reductions (,10% of normal) cause encephalomyophathy [9,10,17]. Currently there is no cure available for MDS and to increase the understanding of the mechanisms behind the syndrome we previously generated and described a mouse strain deficient in TK2 (TK22/2) that we use as a model for human TK2 deficiency, and human TK2 deficient infantile fatal encephalomyopathy in particular [15]. Since the TK22/2 mice completely lack TK2 activity and display multi-organ involvement, the model is similar to patients with TK2 mutations that involve brain and liver as well as muscle. The TK22/2 mice were born with mtDNA depletion

Introduction Mitochondrial DNA (mtDNA) depletion syndrome (MDS) is a group of rare and severe autosomal recessive disorders that are characterized by quantitative reduction of mtDNA in affected tissues [1]. MDS presents in early childhood or infancy with diverse symptoms such as progressive myopathy, hepatopathy and encephalopathy and is often fatal before the age of 3 years [2,3]. The different clinical symptoms have been associated with different mutations in nuclear genes encoding proteins involved in the mtDNA metabolism [4]. Several mutations in the thymidine kinase 2 (TK2), deoxyguanosine kinase (DGUOK), polymerase gamma (POLG) and MpV17 mitochondrial inner membrane protein (MPV17) genes have been linked to inheritable MDS in humans [5,6,7,8]. Humans carrying TK2 deficiency primarily present with severe myopathy and display neurological phenotypes whereas DGUOK, POLG and MPV17 deficiencies have been associated with liver failure and encephalomyopathy [5,6]. However, there are also reports of patients harbouring TK2 deficiency that have liver and brain involvement [9,10,11,12,13]. Patients with a certain TK2 mutation (R172W/R172W) have been shown to harbor TK2-

PLOS ONE | www.plosone.org

1

March 2013 | Volume 8 | Issue 3 | e58843

TK2 Deficient Mice Have an Altered Lipid Phenotype

Figure 1. Transmission electron microscopy pictures of mitochondria in livers from TK2+/+ and TK22/2 F6 mice. Bars = 2 mm. Sections shown are representative of the samples. doi:10.1371/journal.pone.0058843.g001

defect b-oxidation and, as a result, insufficient production of ketone bodies and glucose.

in adipose tissue, brain and skeletal muscle. Although the mice behaved and grew normally the first 7 days they demonstrated progressive loss of mtDNA content in several organs and at postnatal day 12 the mtDNA content was also decreased in liver. Although the liver structure was similar in the TK2+/+ and TK22/ 2 mice, the latter showed an increased number of intracellular lipid vesicles. The TK22/2 mice were also progressively hypothermic and displayed abnormal brown adipose tissue and an almost complete loss of hypodermal fat at 14 days of age when this fat layer is normally well developed. In addition, the adipocytes in the 14 days old TK22/2 mice showed heterogeneity in size and accumulation of lipid vesicles in contrast to a homogeneous pattern seen in the TK2+/+ mice [15]. The attained data suggested an involvement of the lipid metabolism in the progressive loss of essential body functions in the TK22/2 mice. Here we analysed the TK22/2 mouse model to clarify the importance of mtDNA for liver function by studies of the lipid phenotype and functionality of the hepatic mitochondria as well as alterations in plasma. The obtained data show that the TK2deficiency induced gradual loss of mtDNA in the liver leads to

PLOS ONE | www.plosone.org

Table 1. Relative volume of mitochondria in TK2+/+ and TK22/2 mice.

Age (days)

Relative Volume (%) TK2+/+

TK22/2

7

20.563.2

24.061.5

12

27.860.2

45.564.5*

A stereological method was used to determine the volume density of the mitochondria in the cytoplasm from digital electron images taken by a Morada camera. Three different animals were analysed for each genotype and age. For each animal 9 randomly taken images from two different sections (18 images/ animal) were counted making a total of 54 images/group. The P-value for statistical comparison (two-tailed unpaired Student’s t-test) between TK2+/+ and TK22/2 of the same age is shown, *p#0.05. doi:10.1371/journal.pone.0058843.t001

2

March 2013 | Volume 8 | Issue 3 | e58843

TK2 Deficient Mice Have an Altered Lipid Phenotype

High Capacity Reverse Transcriptase Kit (Applied Biosystems). The cDNA was diluted 1 to 10 before real-time PCR was performed on the Applied Biosystems PRISM 7000 using the SYBRgreen reagent (Applied Biosystems). Reactions were performed in triplicates and hypoxanthine-guanine phosphoribosyltransferase (Hprt) was selected as endogenous control. Data was collected on individual cDNA samples. To identify additional genes involved in fatty acid oxidation that might be altered in TK22/2 mice, a PCR screening array for 84 genes in the mouse fatty acid oxidation pathway (catalog no. PAMM-007, Qiagen) was conducted with liver tissue from 14 days old mice (three of each genotype). Total RNA was extracted using the RNA easy mini kit (Qiagen). The RNA integrity was assessed spectrophotometrically, and all samples had 260/280 ratios between 1.8–2.0 and 260/230 ratios above 1.7. Reverse transcription was done on 0.5 mg total RNA per sample using the RT2 First strand kit (Qiagen) according to the manufacturer’s protocol. SYBR green RT-PCR analysis was performed on an ABI 7500 fast cycler. The array evaluated the expression of 84 genes involved in fatty acid oxidation. The acquired data were analyzed for statistical differences using Qiagen’s software RT2 Profiler PCR Array Data Analysis version 3.5. Briefly the PCR array data were calculated by the comparative cycle threshold method, normalized against housekeeping genes, and expressed as mean fold change in TK22/2 samples relative to TK2+/+ control samples. Genes of TK22/2 mice that differed by .2-fold compared with TK2+/+ controls are included and highlighted in bold when p#0.05.

Materials and Methods Animals Official Swedish regulations for the use and care of laboratory animals were followed throughout and the experimental protocol was approved by the local ethical committee (Permit numbers: S104-09 and S135-11). TK2+/2 and TK22/2 mice were developed and maintained as previously described [15]. F8 TK22/2 mice were compared to TK2+/+ litter mates of the same generation, unless otherwise specified. Livers were immediately frozen in liquid nitrogen. Whole blood was collected through heart puncture and centrifuged at 800 g for 10 min. Serum and livers were stored at 280uC until use.

Transmission Electron Microscopy Electron images of mitochondria were obtained as previously described [18]. 40–50 nm ultra-thin sections from longitude parts were cut and examined in a Tecnai 10 transmission electron microscope (Fei, The Netherlands) at 80 kV. Digital images were randomly taken by using a Morada camera (Olympus Soft Imaging Solutions, Germany). A stereological method was used to determine the volume density of the mitochondria in the cytoplasm of the hepatocytes. The areas of mitochondria profiles on printed images and the total area of the cytoplasm was measured by point counting using 2 cm2 lattices according to Weibel [19]. The ratio between the mitochondria areas and cytoplasm area was calculated. A pilot study was performed to determine the number of blocks and images needed for an appropriate sample using cumulative mean plot for evaluation [19]. Thus, 9 randomly taken images from two different sections were counted per animal (18 images/animal) and three different animals were used per group (54 images/group). A total of four groups were studied (7 and 12 days old TK2+/+ and TK22/2). To determine the numbers of mitochondria per cell, estimations were done based on the two-dimensional structures in the images and the assumption that the size of mitochondria in the tissue correlates with the size of the area in the image.

Protein Expression Liver tissues from 7 and 12 days old mice (three of each genotype and age) were ground in liquid nitrogen, dissolved in 400 ml RIPA buffer, incubated on ice for two hours, centrifuged at 16000 g for 20 min at 4uC and the supernatants retrieved. Western blot analyses to determine adipophilin, cytochrome c oxidase subunit II (COX II) and voltage-dependent anion channel (VDAC) protein levels in liver protein extracts were performed as previously described [20]. Bio-RAD Quantity one software was used to determine the intensity of the bands. The expression of adipophilin and COXII was related to the expression of VDAC for each individual sample and the mean of the expression ratios were calculated. The adipophilin antibody was purchased from Progen biotechnik and the COX II antibody was a generous gift of professor Nils-Go¨ran Larsson [20].

mRNA Expression Total RNA from 7 and 12 days old mice (four of each genotype and age) was extracted by TRIzolH (Life Technologies) and quantified with the NanoDropTM 2000 (Thermo Scientific). RNA quality assessment was performed using the ExperionTM RNA StdSens analysis kit (Bio-Rad Laboratories). Oligo-dT-primed cDNA synthesis was performed with 1 mg total RNA using the

Figure 2. Hepatic expression of adipophilin, voltage dependent anion channel (VDAC) and cytochrome oxidase subunit II (COXII). Protein expression in livers of A) 7 days old and B) 12 days old TK2+/+ and TK22/2 mice detected with Western blot. Protein samples are from three individuals of each genotype and age. Bio-RAD Quantity one software was used to determine the intensity of the bands and statistical comparisons (two-tailed unpaired Student’s t-test) between TK2+/+ and TK22/2 were done. VDAC was used to normalize the data. Only 12 days old TK22/2 mice exhibit 150% higher expression of adipophilin compared to TK2+/+ mice (p#0.05), no other difference was detected. doi:10.1371/journal.pone.0058843.g002

PLOS ONE | www.plosone.org

3

March 2013 | Volume 8 | Issue 3 | e58843

TK2 Deficient Mice Have an Altered Lipid Phenotype

Figure 3. Hepatic gene expression and lipid content. A) Hepatic cholesteryl esters in 7 and 12 days old TK2+/+ and TK22/2 mice determined with GC-MS. B) Relative hepatic mRNA expression of 3-hydroxy-3-methyl-glutaryl-Coenzyme A reductase (Hmgcr) and 3-hydroxy-3-methyl-glutarylCoenzyme A synthase (Hmgcs) in 7 and 12 days old TK2+/+ and TK22/2 mice determined with real-time PCR. C) Hepatic triglycerides in 7, 10 and 12 days old TK2+/+ and TK22/2 mice determined with an enzymatic assay. D) Relative hepatic mRNA expression of fatty acid synthase (Fas) and sterol regulatory element binding protein 1c (Srebp1c) in 7 and 12 days old TK2+/+ and TK22/2 mice determined with real-time PCR. Data presented as mean 6 SEM. Statistically significant difference (two-tailed unpaired Student’s t-test) compared to TK2+/+ of the same age, *p#0.05. doi:10.1371/journal.pone.0058843.g003

Tecan) and analysed using the MagellanTM Data Analysis Software (Tecan).

Lipid Parameters in Liver Liver tissues were extracted for lipids as described by Folch [21]. The delipidated liver tissues were dried and dissolved in 1 M NaOH for protein determination by Bio-Rad DC Protein Assay (Bio-Rad Laboratories). All hepatic lipids were related to the amount of protein in the liver tissues. For quantification of hepatic free and esterified cholesterol, d6-cholesterol was added as an internal standard to parts of the lipid extracts from 7 and 12 days old mice (four of each genotype and age). Cholesterol samples were hydrolyzed (for total cholesterol) and purified on ISOLUTEH MFC18 SPE columns (Biotage). The purified samples were trimethylsialylated and dissolved in hexane before analysis on GC-MS, and the amounts of cholesterol were quantified by isotope dilution mass spectrometry. Parts of extracts from 7 days old (eight TK2+/+ and nine TK22/2), 10 days old (two of each genotype), and 12 days old (five of each genotype) mice were also used for quantification of hepatic triglyceride content using the TG GPO-PAP reagent (Roche). Samples were run in duplicates and the measurement was repeated once with new extracts. Absorption was measured at 492 nm (InfiniteH F500 plate reader,

PLOS ONE | www.plosone.org

Serum Parameters Size exclusion chromatography of plasma lipoproteins in samples from 7 and 12 days old mice (four of each genotype and age) was performed as described earlier [22] and their lipid content was measured online with cholesterol CHOD-PAP and TG GPO-PAP reagents (Roche). For quantitation of non-esterified fatty acids (NEFA), plasma was collected and immediately assayed according to the ACS-ACOD method by using the NEFA kit (WAKO, USA). Samples were run in duplicates or triplicates and the values compared to a standard curve performed on the same plate as the samples. Absorption was measured at 550 nm and the background at 750 nm (InfiniteH F500 plate reader, Tecan).The assay was performed three times with samples from 12 days old mice (eight of each genotype), and twice for 7 days old mice (seven of each genotype). Blood glucose levels were determined from whole blood by using the Accu-Chek (Roche). Blood ketone levels were determined in whole blood by using the Precision Exceed (Abbott). Six 4

March 2013 | Volume 8 | Issue 3 | e58843

TK2 Deficient Mice Have an Altered Lipid Phenotype

Table 2. Genes differentially expressed in liver tissue of 14 days old TK2+/+ and TK22/2 mice identified by a PCR-focused array. Fold change (TK22/2/TK2+/+)

Gene name Acetyl-Coenzyme A acyltransferase 1A

24.3

Acetyl-Coenzyme A acyltransferase 2 (mitochondrial 3-oxoacyl-Coenzyme A thiolase)

24.0

Acyl-Coenzyme A dehydrogenase family, member 10

22.6

Acyl-Coenzyme A dehydrogenase family, member 9

24.1

Acyl-Coenzyme A dehydrogenase, long-chain

22.2

Acetyl-Coenzyme A acetyltransferase 1

24.8

Acetyl-Coenzyme A acetyltransferase 2

24.4

Acyl-CoA thioesterase 3

27.5

Acyl-CoA synthetase long-chain family member 1

24.0

Acyl-CoA synthetase long-chain family member 3

24.5

Acyl-CoA synthetase long-chain family member 4

2.2

Acyl-CoA synthetase medium-chain family member 5

26.8

3-hydroxybutyrate dehydrogenase, type 1

234.8

3-hydroxybutyrate dehydrogenase, type 2

27.4

Carnitine palmitoyltransferase 1a, liver

22.9

Carnitine palmitoyltransferase 1b, muscle

2.2

Carnitine palmitoyltransferase 2

23.1

Carnitine O-octanoyltransferase

213.4

2.4-dienoyl CoA reductase 1, mitochondrial

25.1

Fatty acid binding protein 1, liver

22.7

Fatty acid binding protein 2, intestinal

244.3

Fatty acid binding protein 3, muscle and heart

23.5

Fatty acid binding protein 5, epidermal

22.5

Fatty acid binding protein 6, ileal (gastrotropin)

23.9

Glycerol kinase 2

22.7

Glycerol phosphate dehydrogenase 2, mitochondrial

24.4

Glycerol kinase

23.3

Hydroxyacyl-Coenzyme A dehydrogenase/3-ketoacyl-Coenzyme A thiolase/enoyl-Coenzyme A hydratase (trifunctional protein), alpha subunit

24.2

3-hydroxy-3-methylglutaryl-Coenzyme A lyase

22.8

3-hydroxy-3-methylglutaryl-Coenzyme A synthase 1

23.2

3-hydroxy-3-methylglutaryl-Coenzyme A synthase 2

23.7

Lipase, hormone sensitive

22.4

3-oxoacid CoA transferase 2A

23.4

Peroxisomal trans-2-enoyl-CoA reductase

25.4

Pyrophosphatase (inorganic) 1

22.4

Protein kinase, AMP-activated, alpha 1 catalytic subunit

22.7

Protein kinase, AMP-activated, beta 1 non-catalytic subunit

22.2

Protein kinase, AMP-activated, beta 2 non-catalytic subunit

22.7

Protein kinase, cAMP dependent, catalytic, alpha

22.4

Protein kinase, cAMP dependent, catalytic, beta

22.4

Protein kinase, AMP-activated, gamma 2 non-catalytic subunit

22.2

Protein kinase, AMP-activated, gamma 3 non-catatlytic subunit

2.3

Solute carrier family 27 (fatty acid transporter), member 3 2/2

24.6 +/+

samples relative to TK2 control samples. The data analysis, including statistical comparisons, were done using Data are expressed as mean fold change in TK2 Qiagen’s software RT2 Profiler PCR Array Data Analysis version 3.5. Genes of TK22/2 mice (n = 3) that differed by .2-fold compared with TK2+/+ controls (n = 3) are included and highlighted in bold when p#0.05. doi:10.1371/journal.pone.0058843.t002

PLOS ONE | www.plosone.org

5

March 2013 | Volume 8 | Issue 3 | e58843

TK2 Deficient Mice Have an Altered Lipid Phenotype

were analyzed in separate experiments. Comparisons of mitochondrial b-oxidation were made within each run.

Table 3. Serum levels of cholesterol and triglycerides in ApoB-containing, ApoA1-containing and total lipoprotein particles, as well as total serum nonesterified fatty acids (NEFA).

Serum lipids

Genotype

TK2

+/+

Days of age 7

7

12

1.0560.10

1.1260.05

1.3160.17

1.3260.05*

ApoA1-cont 1.9260.04

2.3960.02

2.0360.09

Total Triglycerides ApoB-cont

2.9760.09

3.5060.06

Mitochondria were isolated from liver tissue samples of 14 days old mice (three of each genotype). Mitochondrial ATP production rates (MAPR) and citrate synthase (CS) activity were determined as previously described [26]. Briefly, CS activity was measured at 35uC using a Spectronic GENESYS 10 UV-Vis spectrophotometer (Thermo Fisher Scientific). ATP production rates were determined at 25uC using a firefly luciferase and a VarioskanH Flash Spectral Scanning Multimode reader (Thermo electron corporation) in the presence of the following different substrate combinations: glutamate+malate (G+M), glutamate+succinate (G+S), pyruvate+malate (P+M), palmitoyl-l-carnitine+malate (PC+M), succinate (S) and succinate+rotenone (S+R). Activities were expressed as units per unit of CS activity in the mitochondrial suspension.

2/2

12

Cholesterol ApoB-cont (mmol/l)

TK2

ATP Production Rate and Citrate Synthase Activity in Liver

3.3460.12

2.5860.04** *

3.9060.06**

0.8460.02

0.6960.09

0.8960.04

0.8360.11

(mmol/l)

ApoA1-cont 0.0460.01

0.0560.01

0.0460.00

0.0460.01

Total

0.8860.02

0.7460.08

0.9360.04

0.8760.11

NEFA (mmol/l)

Total

0.7960.10

0.6060.02

0.9160.08

0.7660.02**

For each time point data from at least four different mice have been used (n$4) and are presented as mean 6 SEM. The P-value for statistical comparison (twotailed unpaired Student’s t-test) between TK2+/+ and TK22/2 of the same age group is shown, *p#0.05, **p#0.01. doi:10.1371/journal.pone.0058843.t003

Carnitine Palmitoyl Transferase (CPT) Activity in Liver Total CPT activity was measured in isolated liver mitochondria from 14 days old mice (four of each genotype) according to a previously described method [27]. Briefly, each reaction contained 20 mg mitochondrial protein mixed in a 200 ml reaction buffer containing 20 mM HEPES, 1 mM EGTA, 220 mM sucrose, 40 mM KCl, 100 mM 5,5-dithio-bis(2-nitrobenzoic acid)(DTNB), 1.3 mg/ml BSA, and 40 mM palmitoyl-CoA, pH 7.4. The reaction was started by adding 1 mM carnitine (omitted in baseline) and monitored at 412 nm for 4 min at 25uC using a thermospecific spectrophotometer (Genesys 10 UV, VWR).

mice of each genotype were analyzed to determine the average ketone and glucose levels on postnatal day 14.

Acylcarnitine Metabolites in Plasma The determination of acylcarnitines was performed in 12 days old mice (three of each genotype) by a modification of previously described methods [23,24] as follows: Sample preparation: 50 ml plasma was diluted with a corresponding amount of isotope labelled internal standard mix (NSK-B, Cambridge Isotope Laboratories Inc.) and deproteinised with 500 ml acetonitril containing 0.3% v/v formic acid. The sample was vortexed vigorously and centrifuged. The supernatant was transferred to a clean vial and evaporated, and derivatised with 100 ml 3 M HCl in n-butanol (60uC, 15 min.). The sample was evaporated again and reconstituted with 100 ml methanol and 50 ml 10 mM formic acid. LC-MS/MS procedure: 5 ml of the prepared sample was injected onto an Acquity UPLCTM BEH C18 2.1650 mm, 1.7 mm chromatographic column (Waters). A gradient elution with water (mobile phase A) and methanol (mobile phase B), both containing 0.1% v/v formic acid, was set up to separate the analytes. The acylcarnitines in the resulting eluate were ionized in positive ion mode electrospray and determined with a tandem mass spectrometer (XevoTM TQ MS, Waters) in MRM mode. Quantification of each compound was performed by internal standard methodology using the TargetLynx software (Waters). Plasma from 12 days old mice (three of each genotype) were analyzed in separate experiments.

Statistics Two tailed, unpaired Student’s t-test was used to test for statistical significance of difference between TK2+/+ and TK22/2 mice of the same age. Significance was set at p#0.05.

Results Altered Hepatic Mitochondria Transmission electron microscopy was performed on liver samples from the TK2+/+ and the TK22/2 mice. Sections shown are representative of the samples (Fig.1). The volume of the mitochondria in TK22/2 mice was larger compared to TK2+/+ mice of the same age. In 7 days old TK22/2 mice, mitochondria constituted 24.0% of the total cell volume compared to 20.5% in TK2+/+ mice. An even larger difference was observed in 12 days old mice, where the mitochondria constituted 45.5% of the cell volume in TK22/2 mice, compared to 27.8% in TK2+/+ mice (Table 1). To determine if the increased volume was due to hypertrophy or increased number of mitochondria the number of mitochondria profiles were estimated (2189 for TK2+/+ and 2128 for TK22/2) and correlated to the total areas of the cytoplasm giving a mean number of mitochondria/mm2 cytoplasm of 0.5260.02 for TK2+/+ and 0.5460.08 for TK22/2. However, the mean area/mitochondria was 0.54060.010 mm2 for TK2+/+ and 0.86660.046 mm2 for TK22/2, all together showing that the mitochondria in the 12 days old TK22/2 mice were hypertrophic. It was previously shown that the mtDNA content in livers of 12 days old TK22/2 mice constitutes approximately 70% compared to the TK2+/+ [15,28]. To investigate whether the increased mitochondrial volume and mtDNA depletion in the TK22/2 mice affected the protein expression in the hepatic mitochondria the expression levels of the mitochondrial protein cytochrome c

Betaoxidation in Liver Mitochondrial b-oxidation was determined in liver homogenates according to a previously described method [25]. Briefly, 100–200 mg mitochondrial protein was used in each reaction. The conversion of [1-14C] palmitic acid (Sigma) into acid-soluble metabolites was studied in triplicate samples in the presence or absence of 50 mM antimycin A. The peroxisomal b-oxidation (antimycin A insensitive) was subtracted from the total boxidation. Livers of 14 days old mice (three of each genotype)

PLOS ONE | www.plosone.org

6

March 2013 | Volume 8 | Issue 3 | e58843

TK2 Deficient Mice Have an Altered Lipid Phenotype

Table 4. Plasma acylcarnitine metabolites in 12 days old TK2+/+ and TK22/2 mice.

Acylcarnitines

Carnitine

TK2+/+

TK22/2

t-test

[mM]

[mM]

p-value

15.362.8

23.567.1

0.340

Short-Chain Acetyl (C2) Propionyl (C3)

0.131±0.013

0.057±0.011

0.012

Butyryl (C4)

0.53960.114

0.24160.033

0.066

Methylmalonoyl (C4-DC)

0.025760.0079

0.013360.0017

0.200

Tiglyl (C5:1)

0.000760.0007

0.003760.0023

0.284

Isovaleryl (C5)

0.041060.0410

0.031760.0317

0.866

Medium-chain Hexanoyl (C6)

0.085360.0354

0.045760.0122

0.349

Octanoyl (C8)

0.035060.0055

0.020760.0003

0.060

Decenoyl (C10:1)

0.001360.0013

0.001760.0017

0.883

Decanoyl (C10)

0.047060.0145

0.024360.0038

0.206

Dodecenoyl (C12:1)

0.012760.0009

0.061360.0473

0.362

Lauryl (C12)

0.084060.0405

0.068360.0259

0.761

0.331

Long-chain Myristoleyl (C14:1)

0.037760.0182

0.075060.0285

Myristoyl (C14)

0.234±0.009

0.418±0.023

0.002

3-OH myristoyl (C14-OH)

0.0090±0.0015

0.0163±0.0007

0.012

Palmitoleyl (C16:1)

0.034±0.010

0.115±0.008

0.004

3-OH palmitoleyl (C16:1-OH)

0.001060.0010

0.009360.0047

0.156

Palmitoyl (C16)

0.524±0.016

1.433±0.106

0.001

3-OH palmitoyl (C16-OH)

0.0030±0.0020

0.0193±0.0033

0.013

Octadecdienoyl (C18:2)

0.056360.0274

0.109760.0452

0.370

3-OH linoleyl (C18:2-OH)

0

0.00160.001

0.374

Octadecenoyl (C18:1)

0.149±0.033

0.401±0.032

0.005

3-OH oleyl (C18:1-OH)

0.0017±0.0017

0.0123±0.0030

0.035

3-OH stearoyl (C18-OH)

0.004760.0047

0.012760.0127

0.585

Steaoryl (C18)

0.0567±0.0119

0.1267±0.0058

0.006

For each time point data from three different mice have been used (n = 3) and are presented as mean 6 SEM. Chains highlighted in bold are significantly different (comparisons were made with two-tailed unpaired Student’s t-test) between TK2+/+ and TK22/2 mice and exhibit at least p#0.05. doi:10.1371/journal.pone.0058843.t004

same age, but no significant differences were observed between TK22/2 and TK2+/+ mice of older age (Fig. 3C). The hepatic mRNA expression of fatty acid synthase (Fas), a key enzyme in fatty acid synthesis, was decreased in the 12 days old TK22/2 mice compared to TK2+/+ mice of the same age and the hepatic mRNA expression of the sterol regulatory element binding protein 1c (Srebp1c), which is responsible for the regulation of genes required for de novo lipogenesis, such as Fas, also displayed a similar expression pattern (Fig. 3D). To further investigate the expression of genes involved in the fatty acid oxidation pathway in livers severely depleted of mtDNA, liver samples of 14 days old TK22/2 mice containing ,20% mtDNA, were analysed with PCR screening arrays and compared to TK2+/+. The expression levels of several genes were altered in the TK22/2 mice (Table 2). However, only Hmgcs-2 encoding the mitochondrial protein that catalyzes the first reaction of ketogenesis, was statistically significantly down regulated when analyzed using Qiagen’s software RT2 Profiler PCR Array Data Analysis.

oxidase subunit 2 (COXII) were measured in 7 and 12 days old mice. For both time points the COXII expression in the TK22/2 mice was similar to the expression in the TK2+/+ mice (Fig. 2).

Altered Hepatic Gene Expression and Lipid Content The expression levels of adipophilin, a cellular indicator of the amount of stored lipids, were significantly increased in livers of 12 days old TK22/2 mice (150%) compared to TK2+/+ mice (Fig. 2). Increased amounts of cholesteryl esters were detected in livers of 7 days old TK22/2 mice compared to TK2+/+ mice of the same age. In livers of 12 days old mice no statistically significant difference was detected, although the mean value of the TK22/2 mice was higher (p = 0.052) (Fig. 3A). The hepatic mRNA expression of 3-hydroxy-3-methyl-glutaryl-coenzyme A reductase (Hmgcr), the rate limiting enzyme in cholesterol synthesis, was similar in both 7 and 12 days old TK22/2 and TK2+/+ mice, but 3-hydroxy3-methyl-glutaryl-coenzyme A synthase (Hmgcs) -1 displayed a trend towards decreasing expression levels with increasing age in the TK22/2 mice (Fig. 3B). Hepatic triglyceride levels were elevated in 7 days old TK22/2 mice compared to TK2+/+ mice of the

PLOS ONE | www.plosone.org

7

March 2013 | Volume 8 | Issue 3 | e58843

TK2 Deficient Mice Have an Altered Lipid Phenotype

Table 5. Mitochondrial ATP production rates standardized with citrate synthase activity.

genotype Substrate combination (units/unit CS) G+M TK2

+/+

TK22/2

G+S

PC+M

P+M

S

S+R

1.9460.83 3.8162.28 1.3860.54 3.3861.79 3.6362.48 1.5360.2 1.6961.03 2.8661.50 1.1360.51 2.2761.31 1.8761.54 1.3861.25

CS: citrate synthase; G: glutamate; M: malate; S: succinate; P: pyruvate; PC: palmitoyl-L-carnitine; R: rotenone. Mitochondrial CS activity was similar in the genotypes. For each time point data from 14 days old mice, three from each genotype, have been used (n = 3) and are presented as mean 6 SEM. Statistical comparisons (two-tailed unpaired Student’s t-test) between TK2+/+ and TK22/2 were made but no statistical differences were detected. doi:10.1371/journal.pone.0058843.t005

Figure 4. Mitochondrial palmitate oxidation rate in liver homogenates of 14 days old TK2+/+ and TK22/2 mice. Three independent measurements were performed (TK2+/+ n = 3, TK22/2 n = 3) and comparisons between the genotypes were made within each measurement. Data presented as per cent activity (mean 6 SEM) compared to TK2+/+. The P-value for statistical comparison (two-tailed unpaired Student’s t-test) between TK2+/+ and TK22/2 is shown, *p#0.05. doi:10.1371/journal.pone.0058843.g004

Altered Hepatic Mitochondrial b-oxidation The accumulated data pointed to impaired b-oxidation and for that reason the fatty acid oxidation was determined in liver homogenates from 14 days old mice (Fig. 4). The TK22/2 mice demonstrated considerably reduced mitochondrial palmitate oxidation rate within each measurement, exhibiting only 9–40% of the oxidation rate observed in the TK2+/+ mice.

Elevated Levels of Cholesterol and NEFAs in Plasma

Normal MAPR

Plasma total cholesterol levels were elevated in both 7 and 12 days old TK22/2 mice compared to TK2+/+ (Table 3). Separation of the lipoprotein fractions revealed that the cholesterol levels in both ApoA1-containing (HDL) and ApoB-containing (Chylomicron/VLDL/LDL) lipoprotein particles tended to be higher in the TK22/2 mice of both ages, but the elevation was statistically significant only in 12 days old TK22/2 mice. No changes in plasma triglyceride levels were detected between the genotypes (Table 3). The highly hydrophobic NEFAs are generally toxic to cells and therefore kept at low concentrations in plasma. The levels of NEFAs were similar in 7 days old TK22/2 and TK2+/+ mice, but significantly increased in serum of 12 days old TK22/2 mice (Table 3). The blood glucose levels were determined in several mice. However, the blood sugar is increased as a response to food intake and all mice in this study were suckling pups with constant access to food, making it impossible to determine a fasting blood sugar value. Nevertheless, in some of the sacrificed 14 days old TK22/2 mice blood sugar values as low as 2 mmol/l were observed. Such low values were not detected in the TK2+/+ mice. The ketone levels were also significantly reduced in 14 days old TK22/2 mice, on average 1.6 mmol/l in TK22/2 compared to 2.4 mmol/l in TK2+/+ mice.

It is known that mtDNA depletion may interfere with respiratory chain ATP production [27]. In 14 days old TK22/ 2 mice the mtDNA levels had drastically dropped below 20% compared to TK2+/+, however, the mtDNA depletion did not affect the MAPR at this point (Table 5).

Altered CPT Activity The observed MAPR in 14 days old mice was essentially normal also when the substrates entered through fatty acid oxidation (palmitoyl-L-carnitine plus malate). A key regulatory site controlling the flux of long-chain fatty acids through b-oxidation is the mitochondrial carnitine system that catalyzes the transport of longchain fatty acids into the mitochondrial matrix. Since the gene expression level of carnitine palmitoyltransferase 1A (CPT1A), which is responsible for the formation of long-chain acylcarnitines from activated fatty acids and free carnitine, was decreased (Table 2),

Accumulation of Acylcarnitines in Plasma In mitochondrial fatty acid oxidation defects, metabolic blocks exist in the catabolic pathways of mitochondrial coenzyme A esters. Screening for disorders of fatty acid b-oxidation is commonly performed by determination of plasma acylcarnitine levels [29]. To determine if the plasma acylcarnitine metabolites differed between the mice the acylcarnitine profiles in plasma were determined by tandem mass spectrometry (Table 4). The TK22/2 mice exhibited elevated levels of acyl- and 3-OH acylcarnitines with long chain lengths compared to the levels in TK2+/+ mice. In addition, the level of propionylcarnitine (C3) was significantly reduced and the levels of methylmalonylcarnitine (C4-DC) and buturylcarnitine (C4) were also decreased, although not significantly, in the TK22/2 mice. PLOS ONE | www.plosone.org

Figure 5. Carnitine palmitoyltransferase activity in hepatic mitochondria from TK2+/+ and TK22/2 mice 14 days of age. The results represent data from four individuals of each genotype. Data presented as per cent activity (mean 6 SEM) compared to TK2+/+. The Pvalue for statistical comparison (two-tailed unpaired Student’s t-test) between TK2+/+ and TK22/2 is shown, **p#0.01. doi:10.1371/journal.pone.0058843.g005

8

March 2013 | Volume 8 | Issue 3 | e58843

TK2 Deficient Mice Have an Altered Lipid Phenotype

levels in plasma have not yet indicated any difference between the TK22/2 and TK2+/+ mice, although it should be mentioned that these measurements were performed on non-fasted mice (data not shown). Since the accumulated data suggested a defect in the fatty acid oxidation in the TK22/2 mice plasma of 12 days old mice was screened for fatty acid oxidation disorders. Long-chain acylcarnitines and 3-OH acylcarnitines were found to accumulate in TK22/2 mice suggesting an impaired b-oxidation. Therefore the palmitate oxidation rate was analyzed in liver mitochondria and found to be decreased by 60–90% in the TK22/2 mice. However, due to limited sample size and inability to fast the mice the absolute palmitate oxidation rates could not be determined. The severe drop in mtDNA levels in the TK22/2 mice occurred rapidly between postnatal day 12 and 14, and it is likely that residual activities together with activities of newly synthesized mtDNA encoded proteins, were enough to maintain the respiratory chain activity at this point. The increased mitochondrial volume in the TK22/2 liver cells may also contribute to further activate the respiratory chain. Other compensatory mechanisms such as increased mitochondrial transcript or protein half-lives, which have been suggested to occur in muscle tissue harboring severe mtDNA depletion due to TK2 deficiency [31], could also be involved. The observed MAPR was normal also for substrates entering through the fatty acid oxidation as palmitoyl-L-carnitine plus malate, suggesting that the alteration in the b-oxidation would be prior to the formation of acylcarnitines from fatty acids and carnitine. In the analysis of genes involved in the fatty acid oxidation pathway the expression of carnitine palmitoyltransferase CPT1A, which is responsible for conjugation of carnitine to long chain acyl-CoAs, was reduced. To test if the CPT1A activity was altered in 14 days old TK22/2 mice a method determining total CPT activity was used and the activity was found to be reduced to 65% compared to TK2+/+. Since CPT1A is the overall ratelimiting step in the b-oxidation of long chain fatty acids it is likely that the reduced CPT1A activity is participating in the altered boxidation in the TK22/2 mice. CPT1A is regulated both by malony-CoA, which inhibits its activity, and by the curvature and lipid characteristics of the outer mitochondrial membrane [35]. The methylmalonyl-CoA levels were unaffected in the TK22/2 mice but the lipid profile of the TK22/2 mice was altered and it is possible that it contributes to a decrease in CPT activity. In 12 days old TK22/2 mice the fatty acid oxidation was impaired but the observed ketone and glucose levels were normal. However, in 14 days old TK22/2 mice the b-oxidation was reduced in the liver and the observed ketone and glucose levels were dramatically decreased. Since the brain cannot utilize fatty acids directly it is dependent on the liver to produce ketone bodies through fatty acid b-oxidation. These mechanisms are of vital importance for normal brain function when the organism is on a carbohydrate restricted diet. The brain has high ATP requirement and low regenerative capacity and is highly sensitive to mtDNA depletion. We have previously shown that lack of TK2 activity causes a decline in mtDNA levels in neurons in vivo [36]. The depletion of mtDNA caused impaired mitochondrial bioenergetic function and degeneration of selected neuronal types. The symptoms in patients with TK2 mutations that cause mtDNA depletion in brain and liver are similar to what is observed in the TK22/2 mouse model. Besides the major reduction of mtDNA in muscle, liver and brain they also show an early disease onset, accumulation of fat in the liver in form of macrovesicular steatosis, and seizures as a late manifestation similar to the mice [9]. However, the TK22/2 mouse completely lacks TK2 activity and displays more severe neurological features compared to patients

the CPT activity in the livers of 14 days old mice was determined (Fig. 5). The method used determined the activity of both carnitine palmitoyltransferases, CPT1A and CPT2. The livers of TK22/2 mice displayed approximately 65% CPT activity compared to TK2+/+ mice.

Discussion In this study liver and plasma from both symptom-free (7 days old) and symptomatic (12 and 14 days old) TK22/2 mice were used to clarify how hepatic mitochondria are involved in the development of the encephalopathic phenotype of the TK22/2 mice. The hepatic cells of TK22/2 mice progressively lose the mtDNA and contain approximately 90% at postnatal day 7, 70% at postnatal day 12 and less than 20% at postnatal day 14, compared to the TK2+/+ mice. The TK22/2 mice retained normal expression of the mtDNA encoded protein COXII in the hepatic cells both at day 7 and 12, showing that they had enough mtDNA to maintain normal protein levels. It has been shown before that partially mtDNA depleted cells can retain COXII protein levels [30,31]. In a study of brain and heart of 13 days old TK2 knock-in mice several different mitochondrial initiation transcription factors were unaltered although the mtDNA content was reduced [28]. The loss of mtDNA content in hepatocytes was accompanied by an increase in mitochondrial volume, which was almost doubled in the TK22/2 mice at day 12. Mitochondrial volume homeostasis is a housekeeping function that is essential to maintain the structural integrity of the mitochondrion. Increased matrix volume has also been suggested to activate the respiratory chain and thereby increase ATP production [32,33]. The identified increase in matrix volume could, thus, be interpreted as a response to energy depletion. Furthermore, microscopy data showed altered mitochondrial structure and accumulation of lipid vesicles in the liver cells of both 7 and 12 days old TK22/2 mice. In addition, hepatic adipophilin protein expression, a cellular indicator of stored lipids, was significantly increased in 12 days old TK22/2 mice. Macrovesicular steatosis has also been seen periportally in liver of human patients harboring the multi-organ involvement form of TK2 mutations in which mtDNA depletion has been detected in brain and liver as well as muscle [9]. However, measurements of hepatic cholesteryl esters and triglycerides in mice older than 10 days did not show any significant alterations between TK22/2 and TK2+/+, suggesting an alteration in the components of lipid vesicles in these mice which needs further investigation. Lipids, such as cholesterol and NEFAs were elevated in 12 days old TK22/2 mice. The increased blood lipids may stem from a number of factors, such as intake from food and hepatic and peripheral clearance from and output to the bloodstream. However, high levels of circulating lipids may be an indicator of fat being released from the adipose tissue, contributing to the accumulation of lipid droplets in the liver cells of the TK22/2 mice. Consistent with an altered lipid metabolism in the TK22/2 mice abnormal adipose tissue development together with a moderate impairment of mitochondrial respiratory function was also found in TK2 knock-in mice [34]. Increased NEFAs accompanied by an exponential increase in total ketone bodies in plasma are also seen in patients suffering from starvation. However, the gastrointestinal tracts of the sacrificed TK22/2 mice were filled with milk and the ketone-bodies were decreased in the TK22/2 mice, not reflecting a typical profile for starvation. It is possible that hormonal imbalances e.g. caused by loss of pancreatic functions are also contributing to the deteriorating condition. However, our own preliminary investigations of insulin PLOS ONE | www.plosone.org

9

March 2013 | Volume 8 | Issue 3 | e58843

TK2 Deficient Mice Have an Altered Lipid Phenotype

with TK2 mutations [36]. These differences may reflect different tissue-specific dNTP pool regulatory pathways between the species. It has also been suggested that the mouse brain is more dependent on the mitochondrial nucleotide salvage pathway than other tissues. The obtained data of the TK22/2 mice resembles symptoms seen in patients with progressive liver failure due to DGUOK deficiency that commonly present with hypoketotic hypoglycemia in combination with lactic acidosis and mild neurological involvement [37]. These patients have profound mtDNA depletions in liver which has been observed already at the prenatal stage in one patient [37,38]. The cause of hypoglycemia in these patients has been suggested to be due to an islet cell hyperplasia causing hyperinsulinism [37]. A patient with a MPV17 mutation experienced improved liver function upon continuous glucose infusions and regular feeding [39]. Even though reduced b-oxidation or altered ketogenesis have not been described in patients harbouring TK2 deficiency with liver and brain involvement it might not have been studied and warrants further attention in this patient group. Altered ketogenesis has previously been described in patients with other hepatocerebral mtDNA depletion syndromes and it would be relevant to investigate if these patients suffer from reduced b-oxidation or CPT-activity as well. The mtDNA content in the brain of the TK22/2 mice was lower already at birth, and was somewhat further decreased at 14 days of age. However, the TK22/2 mice were born with normal levels of mtDNA in the liver and developed severe mtDNA

depletion by postnatal day 14 (,20% of normal mice). The severe mtDNA depletion coincided with severely reduced mitochondrial b-oxidation in the liver most likely due to loss of CPT1A activity. We propose that the gradually increasing failure of the liver to produce ketone bodies and glucose for the brain diminishes the remaining brain function and leads to loss of consciousness and ultimately death of the TK22/2 mice. A better understanding of the primary and secondary factors that influence TK2 progression, and liver damage due to mtDNA depletion, may contribute to better prognosis and adequate treatments. In the absence of a cure for MDS the achieved results provide additional insight in the pathogenesis that can be used to explore novel therapeutic strategies.

Acknowledgments The authors wish to thank M.Sc. Marjan Amiri for excellent technical assistance and Professor Nils-Go¨ran Larsson for the kind gift of the COX II antibody.

Author Contributions Conceived and designed the experiments: XZ KK SC MG AK. Performed the experiments: XZ KK SC KH SI. Analyzed the data: XZ KK SC UD KH. Contributed reagents/materials/analysis tools: UD KH MG AK. Wrote the paper: XZ KK SC UD KH SI MG AK.

References 17. Wang L, Limongelli A, Vila MR, Carrara F, Zeviani M, et al. (2005) Molecular insight into mitochondrial DNA depletion syndrome in two patients with novel mutations in the deoxyguanosine kinase and thymidine kinase 2 genes. Mol Genet Metab 84: 75–82. 18. Park CB, Asin-Cayuela J, Ca´mara Y, Shi Y, Pellegrini M, et al. (2007) MTERF3 is a negative regulator of mammalian mtDNA transcription. Cell. 130: 273–285. 19. Weibel E (1979) Stereological methods, Practical methods for biological morphometry. Volume 1, Academic Press, London, 91–100. 20. Wredenberg A, Wibom R, Wilhelmsson H, Graff C, Wiener HH, et al. (2002) Increased mitochondrial mass in mitochondrial myopathy mice. Proc Natl Acad Sci USA 99: 5066–15071. 21. Folch J, Lees M, Sloane Stanley GH (1957) A simple method for the isolation and purification of total lipides from animal tissues. J Biol Chem 226: 497–509. 22. Parini P, Johansson L, Broijersen A, Angelin B, Ruding M (2006) Lipoprotein profiles in plasma and interstitial fluid analyzed with an automated gel-filtration system. Eur J Clin Invest 36: 98–104. 23. Chase DH, DiPerna JC, Mitchell BL, Sgroi B, Hofman LF, et al. (2001) Electrospray tandem mass spectrometry for analysis of acylcarnitines in dried postmortem blood specimens collected at autopsy from infants with unexplained cause of death. Clin Chem 47: 1166–1182. 24. Ghoshal AK, Tiedong G, Soukhova N, Soldin SJ (2005) Rapid measurement of plasma acylcarnitines by liquid chromatography-tandem mass spectrometry without derivatization. Clin Chim Acta 358: 104–112. 25. Hirschey MD, Shimazu T, Goetzman E, Jing E, Schwer B, et al. (2010) SIRT3 regulates fatty acid oxidation via reversible enzyme deacetylation. Nature 464: 121–125. 26. Wibom R, Hagenfeldt L, von Do¨beln U (2002) Measurement of ATP production and respiratory chain enzyme activities in mitochondria isolated from small muscle biopsy samples. Anal Biochem 311: 139–151. 27. Boudina S, Sena S, O’Neill B T, Tathireddy P, Young M E, et al. (2005) Reduced mitochondrial oxidative capacity and increased mitochondrial uncoupling impair myocardial energetics in obesity. Circ 112: 2686–2695. 28. Dorado B, Area E, Akman H O, Hirano M (2011) Onset and organ specificity of Tk2 deficiency depends on Tk1 down-regulation and transcriptional compensation. Hum Mol Genet 20: 155–164. 29. Chace DH, Kalas TA, Naylor EW (2003) Use of tandem mass spectrometry for multianalyte screening of dried blood specimens from newborns. Clin Chem 49: 1797–1817. 30. Durham SE, Bonilla E, Samuels DC, DiMauro S, Chinnery PF (2005) Mitochondrial DNA copy number threshold in mtDNA depletion myopathy. Neurology. 65: 453–455. 31. Vila` MR, Villarroya J, Garcı´a-Arumı´ E, Castellote A, Meseguer A, et al. (2008) Selective muscle fiber loss and molecular compensation in mitochondrial myopathy due to TK2 deficiency. J Neurol Sci 267: 137–141.

1. Moraes C T, Shanske S, Tritschler H-J, Aprille JR, Andreetta F, et al. (1991) mtDNA depletion with variable tissue expression: a novel genetic abnormality in mitochondrial diseases. Am J Hum Genet 48: 492–501. 2. DiMauro S, Moraes CT (1993) Mitochondrial Encephalomyopathies. Arch Neurol 50: 1197–1208. 3. McFarland R, Taylor RW, Turnbull DM (2002) The neurology of mitochondrial DNA disease. Lancet Neurol 1: 343–351. 4. Taylor RW, Turnbull DM (2005) Mitochondrial DNA mutations in human disease. Nat Rev Genet 6: 389–402. 5. Saada A, Shaag A, Mandel H, Nevo Y, Eriksson etal. (2001) Mutant mitochondrial thymidine kinase in mitochondrial DNA depletion myopathy. Nat Genet 29: 342–344. 6. Mandel H, Szargel R, Labay V, Elpeleg O, Saada A, et al. (2001)The deoxyguanosine kinase gene is mutated in individuals with depleted hepatocerebral mitochondrial DNA. Nat Genet 29: 337–341. 7. Van Goethem G, Dermaut B, Lo¨fgren A, Martin J-J, Van Broeckhoven C (2001) Mutation of POLG is associated with progressive external ophthalmoplegia characterized by mtDNA deletions. Nat Genet 28: 211–212. 8. Spinazzola A, Viscomi C, Fernandez-Vizarra E, Carrara F, D’Adamo P, et al. (2006) MPV17 encodes an inner mitochondrial membrane protein and is mutated in infantile hepatic mitochondrial DNA depletion. Nat Genet 38: 570– 575. 9. Go¨tz A, Isohanni P, Pihko H, Paetau A, Herva R, et al. (2008) Thymidine kinase 2 defects can cause multi-tissue mtDNA depletion syndrome. Brain 131: 2841– 2850. 10. Lesko N, Naess K, Wibom R, Solaroli N, Nennesmo I, et al. (2010) Two novel mutations in thymidine kinase-2 cause early onset fatal encephalomyopathy and severe mtDNA depletion. Neuromuscul Disord 20: 198–203. 11. Mancuso M, Salviati L, Sacconi S, Otaegui D, Camano P, et al. (2002) Mitochondrial DNA depletion: mutations in thymidine kinase gene with myopathy and SMA. Neurol 59: 1197–1202. 12. Oskoui M, Davidzon G, Pascual J, Erazo R, Gurgel-Giannetti J, et al. (2006) Clinical spectrum of mitochondrial DNA depletion due to mutations in the thymidine kinase 2 gene. Arch Neurol 63: 1122–1126. 13. Carrozzo R, Bornstein B, Lucioli S, Campos Y, de la Pena P, et al. (Mar 19,2003) Mutation analysis in 16 patients with mtDNA depletion. Hum Mutat 10.1002/humu.9135. 14. Akman HO, Dorado B, Lopez LC, Garcı´a-Cazorla A, Vila MR, et al. (2008) Thymidine kinase 2 (H126N) knock-in mice show the essential role of balanced deoxynucleotide pools for mitochondrial DNA maintenance. Hum Mol Genet 17: 2433–2440. 15. Zhou X, Solaroli N, Bjerke M, Stewart JB, Rozell B, et al. (2008) Progressive loss of mitochondrial DNA in thymidine kinase 2-deficient mice. Hum Mol Genet 17: 2329–2335. 16. Saada A (2009) Fishing in the (deoxyribonucleotide) pool. Biochem J 422: e3–6.

PLOS ONE | www.plosone.org

10

March 2013 | Volume 8 | Issue 3 | e58843

TK2 Deficient Mice Have an Altered Lipid Phenotype

36. Bartesaghi S, Betts-Henderson J, Cain K, Dinsdale D, Zhou X, et al. (2010) Loss of thymidine kinase 2 alters neuronal bioenergetics and leads to neurodegeneration. Hum Mol Genet 19: 1669–1677. 37. Pronicka E, We˛glewska-Jurkiewicz A, Taybert J, Pronicki M, Szyman´skaDe˛bin´ska T, et al. (Nov 16, 2010) Post mortem identification of deoxyguanosine kinase (DGUOK) gene mutations combined with impaired glucose homeostasis and iron overload features in four infants with severe progressive liver failure. J Appl Genet 10.1007/s13353–010–0008-y. 38. Mu¨ller-Ho¨cker J, Muntau A, Scha¨fer S, Jaksch M, Staudt F, et al. (2002) Depletion of mitochondrial DNA in the liver of an infant with neonatal giant cell hepatitis. Hum Pathol 33: 247–253. 39. Parini R, Furlan F, Notarangelo L, Spinazzola A, Uziel G, et al. (2009) Glucose metabolism and diet-based prevention of liver dysfunction in MPV17 mutant patients. J Hepatol 50: 215–221.

32. Halestrap AP (1989) The regulation of the matrix volume of mammalian mitochondria in vivo and in vitro and its role in the control of mitochondrial metabolism. Biochim Biophys Acta 973: 355–382. 33. Lim KH, Javadov SA, Das M, Clarke SJ, Suleiman MS, et al. (2002) The effects of ischaemic preconditioning, diazoxide and 5-hydroxydecanoate on rat heart mitochondrial volume and respiration. J Physiol 545: 961–974. 34. Villarroya J, Dorado B, Vila` MR, Garcia-Arumi E, Domingo P, et al. (Dec 27, 2011) Thymidine kinase 2 deficiency-induced mitochondrial DNA depletion causes abnormal development of adipose tissues and adipokine levels in mice. PLoS ONE 10.1371/journal.pone.0029691. 35. Rao JN, Warren GZL, Estolt-Povedano S, Zammit VA, Ulmer TS (2011) An environment-dependent structural switch underlies the regulation of carnitine palmitoyltransferase 1A. J Biol Chem 286: 42545–42554.

PLOS ONE | www.plosone.org

11

March 2013 | Volume 8 | Issue 3 | e58843