Abstract
Mutations in the MPV17 gene are associated with hepatocerebral form of mitochondrial depletion syndrome. The mechanisms through which MPV17 mutations cause respiratory chain dysfunction and mtDNA depletion is still unclear. The MPV17 gene encodes an inner membrane mitochondrial protein that was recently described to function as a non-selective channel. Although its exact function is unknown, it is thought to be important in the maintenance of mitochondrial membrane potential (ΔΨm). To obtain more information about the role of MPV17 in human disease, we investigated the effect of MPV17 knockdown and of selected known MPV17 mutations associated with MPV17 disease in vitro. We used different approaches in order to evaluate the cellular consequences of MPV17 deficiency. We found that lower levels of MPV17 were associated with impaired mitochondrial respiration and with a quiescent energetic metabolic profile. All the mutations studied destabilized the protein, resulting in reduced protein levels. We also demonstrated that different mutations caused different cellular abnormalities, including increased ROS production, decreased oxygen consumption, loss of ΔΨm, and mislocalization of MPV17 protein. Our study provides novel insight into the molecular effects of MPV17 mutations and opens novel possibilities for testing therapeutic strategies for a devastating group of disorders.
Introduction
Mitochondrial DNA depletion syndromes (MDDS) are autosomal recessive disorders characterized by a severe decrease in mitochondrial DNA (mtDNA) copy number in affected tissues. MDDS are one of the most common forms of childhood respiratory chain deficiencies (Sarzi et al., ; Yamazaki et al., ). The clinical presentation can be very heterogeneous, with overlapping phenotypes, but MDDS are usually classified as purely myopathic, hepatocerebral, encephalomyopathic, and neurogastrointestinal (Suomalainen and Isohanni, ; Spinazzola, ).
Several genes are associated with MDDS, affecting mtDNA maintenance by three possible mechanisms: (i) impaired mitochondrial nucleotide (dNTP) synthesis [mutations in TK2 (Saada et al., ), SUCLA2 (Elpeleg et al., ), SUCLG1 (Ostergaard et al., ), RRMB2 (Bourdon et al., ), DGUOK (Mandel et al., ), TYMP (Nishino et al., ), ABAT (Besse et al., ), SLC25A4 (Echaniz-Laguna et al., ), and AGK (Mayr et al., )]; (ii) deficient mtDNA replication [mutations in POLG (Naviaux et al., ), POLG2 (Young et al., ), TWINKLE (Sarzi et al., ), TFAM (Stiles et al., ), RNASEH1 (Reyes et al., ), MGME1 (Kornblum et al., ), and DNA2 (Ronchi et al., )], which leads to defects in mtDNA synthesis and repair; or (iii) defective mitochondrial dynamics [OPA1 (Amati-Bonneau et al., ), MFN2 (Rouzier et al., ), and FBXL4 (Antoun et al., )]. All these mechanisms cause deficient turnover and/or segregation of mtDNA to daughter cells, leading to a decrease in its copy number. Ultimately, mtDNA depletion results in inadequate synthesis of respiratory chain complexes with consequent deficient energy production and organ dysfunction (Sarzi et al., ; Spinazzola et al., ).
MPV17 is one of the genes associated with the hepatocerebral form of MDDS (MTDPS6; NNH; MIM #256810). To date, ~100 affected individuals have been reported with MPV17-related disease, with a total number of 48 MPV17 pathogenic variants described (Karadimas et al., ; Spinazzola et al., , ; Nogueira et al., ; Wong et al., ; Navarro-Sastre et al., ; Kaji et al., ; Parini et al., ; El-Hattab et al., , ; AlSaman et al., ; Merkle et al., ; Mendelsohn et al., ; Uusimaa et al., ; Bijarnia-Mahay et al., ; Choi et al., ; Dalla Rosa et al., ; Kim et al., ).
Despite some clinical heterogeneity, the most common clinical manifestations in patients harboring mutations in MPV17 are infantile-onset hepatic and neurologic manifestations, with cerebral leukoencephalopathy in brain MRI and metabolic abnormalities, including lactic acidosis. Progressive liver failure is frequently the cause of death in early childhood. More rarely (around 4%), MPV17 mutations have also been associated with late-onset neuromyopathy with multiple mtDNA deletions in skeletal muscle with minimal or no hepatic disease (Blakely et al., ).
The human MPV17 gene encodes a 176 amino acid protein of 18 kDa. MPV17 was previously thought to be a peroxisomal protein (Zwacka et al., ) but was later found to localize in the inner mitochondrial membrane (El-Hattab et al., ). Although the function of MPV17 and the mechanisms leading to mtDNA dysfunction are still unclear, MPV17 was described as a non-selective channel that modulates mitochondrial membrane potential (ΔΨm) (Antonenkov et al., ). Several factors regulate the gating properties of this putative channel, which is maintained closed in normally functioning mitochondria (Antonenkov et al., ). In damaged mitochondria, the MPV17 channel is prone to be fully open. Some point mutations in MPV17 are also predicted to produce a leaky channel, unable to maintain a high (normal) ΔΨm. Furthermore, loss of MPV17 function causes a total loss of channel activity, increasing ΔΨm (Antonenkov et al., ). In both cases there is excessive ROS production (Binder et al., ; Zorov et al., ). Loss of MPV17 function results in failure of oxidative phosphorylation (OXPHOS) and in mtDNA depletion in affected individuals (Karadimas et al., ; Wong et al., ; Uusimaa et al., ). In MPV17−/− mice (Viscomi et al., ) a reduced rate of mtDNA replication has been observed, due to dNTP insufficiency (Dalla Rosa et al., ).
To better understand the role of MPV17 in human disease, we investigated the effect of MPV17 knockdown in human cell lines and characterized the effect of selected MPV17 mutations associated with MDDS on cellular and organelle localization, and on mitochondrial bioenergetics.
Materials and Methods
Plasmids
Human MPV17 cDNA (a kind gift from Prof. M. Zevianni, Milano, Italy) was cloned with an HA C-terminal tag, in pcDNA3.1 (Invitrogen) as previously described (El-Hattab et al., ). The MPV17 mutations were generated by site-directed mutagenesis using the QuickChange® Kit (Stratagene). The sequences of the primers used for mutagenesis are shown in Supplementary Table 1. Mutagenesis was confirmed by Sanger sequencing (Supplementary Material 2). MPV17-specific shRNA-pLKO plasmid (clone ID TRCN0000129921; Santa Cruz Biotechnology) was used to knockdown the expression of MPV17 in cell lines.
Cell Culture and Transfections
Cell lines were cultured in DMEM (Dulbecco's modified Eagles's medium) (Invitrogen) supplemented with 10% (v/v) FBS (fetal bovine serum), 10 U/mL penicillin and 100 μg/mL streptomycin at 37°C, in a 5% CO2 humidified incubator. HEK293T knockdown cell line was supplemented with 2.5 mM pyruvate and 0.2 mM uridine. Cells were transiently transfected with the different constructs using FuGENE®6 transfection reagent (Roche), according to the manufacturer's instructions.
Subcellular Fractionation
Standard methods were used to prepare total cell lysates, mitochondrial, and post-mitochondrial fractions from cultured cells. In brief, cells were harvested, washed, and resuspended in MIB (Mitochondrial Isolation Buffer [225 mM mannitol, 75 mM sucrose, 10 mM MOPS, 1 mM EGTA, 0,5% BSA, pH 7.2]. Cells were homogenized with a glass Pyrex homogenizer (type B pestle, 30 strokes). Unbroken cells, large plasma membrane and nuclei were removed by two centrifugations at 1,000 × g for 10 min. The resultant supernatant was again centrifuged at 8,000 × g for 15 min at 4°C to separate the mitochondrial fraction (pellet) from the cytosolic fraction (supernatant). Mitoplasts were isolated from mitochondria by digitonin treatment (Gallet et al., ).
Western Blot Analysis
Protein extracts were separated in 15% SDS-polyacrylamide gel and transferred to nitrocellulose membranes. Equal amounts of protein were loaded - 80 μg for whole cell homogenate and for endogenous MPV17, 50 μg for mitochondria, 50 μg mitoplasts. Protein concentration was determined using BCA Kit (Pierce). Primary antibodies used were: goat polyclonal anti-MPV17 antibody (1:250, SantaCruz Biotechnology), anti-HA (1:1000, SantaCruz Biotechnology), and anti-Complex IV (1:5000, Invitrogen). The secondary antibodies were developed using ECL kit (Bio-Rad) and band intensities were quantified using ImageJ software.
Viral Production
The shRNA was co-transfected with the packaging plasmid pCMV-dR8.9 (Addgene) and envelope plasmid pCMV-VSV-G (Addgene) using FuGENE®6 transfection reagent (Roche). Non-specific shRNA-pLKO plasmid was used as control. Infection of 293T cells using lentivirus was performed by exposing cells to the virus-containing medium and incubating for 4 h. Cells were maintained under selection conditions with puromycin. Stable knockdown cell lines were periodically checked by standard Western blotting, using MPV17 antibody. Rescue studies in MPV17 knockdown cells were carried out by overexpression of HA-tagged constructs of wild-type MPV17 and mutant MPV17 cDNA obtained by mutagenesis.
Real-Time Quantitative PCR Analysis of Mitochondrial DNA Content
The relative mtDNA copy number was evaluated using quantitative real-time PCR (qRT-PCR). In brief, DNA samples were diluted to 0.4 ng/μL using Tris-EDTA buffer solution. For each sample, two primer pairs were used to amplify the tRNA Leu(UUR) gene - mtDNA and ß-2-microglobulin (ß2M) gene - nuclear DNA, in separate wells. The primers sequences and qRT-PCR amplification conditions were previously described (Venegas et al., ). The amplification for each sample was performed in a final volume of 10 μL, using the 2 × SYBR SuperMix (iTaq SYBR Green Supermix with ROX, BioRad) and the primer concentration used was 5 μM. All samples were run twice, in triplicate for both mitochondrial and nuclear genes using 7500 Fast Real-time PCR system (qRTPCR; PE7500 real-time PCR instrument; Applied Biosystems, Foster City, CA, USA). A negative and a positive control were also included in each run in order to verify possible contaminations and to act as intern calibrator of values, respectively. Standard deviations for the cycle of threshold value were accepted at 0.50. The results were analyzed with the 7500® v2.0.4 software (Applied Biosystems).
Confocal Microscopy
Cells were grown to around 60% confluency on glass coverslips coated with poly-D-lysine (Sigma). Cells were then sequentially incubated with medium containing 200 nM MitoTracker Red CMXRos (Molecular Probes), fixed in 4% paraformaldehyde (PFA), permeabilized with 0.1% Triton-100, blocked in 10% normal goat serum and incubated with primary anti-HA antibody (1:100, SantaCruz Biotechnology). Next, cells were incubated with conjugated secondary antibody anti-rabbit Alexa Fluor 488 (1:1000, Invitrogen) and DAPI. Images were captured in a Zeiss LSM 510 META confocal microscope (Carl Zeiss) with 63x/1.4 oil immersion objective. Acquisitions were made using the MetaMorph software (Universal Imaging). The images were analyzed in ImageJ.
Cellular Metabolic Activity Studies
Cellular metabolic activity was used as a measure of viability and was determined using resazurin (Sigma-Aldrich, St. Louis, Missouri, USA) assay in 96-well plates. Fluorescence intensity was measured (excitation 560 nm, emission 590 nm) in a microplate reader Tecan Infinite® 200 (Tecan, MÄnnedorf, Switzerland). The release of lactate dehydrogenase (LDH) into the culture media was used as a measure of cytotoxicity, and was performed according to the manufacturers' instructions (Sigma-Aldrich).
Mitochondrial Respiration Studies
Wildtype (control) and MPV17-knockdown HEK293T cells were grown in standard conditions. KD cell line was supplemented with 2.5 mM pyruvate and 0.2 mM uridine. The assay was performed in the Seahorse Bioscience XF24. Oxygen Consumption Rate (OCR) and Extra-Cellular Acidification Rate (ECAR) measurements were measured under basal conditions and after the sequential injection of 1 μg/ml oligomycin, 0.5 μM of FCCP and 1 μM rotenone/antimycin A, allowing for the estimation of basal respiration, proton leak, maximal respiration rate, spare respiratory capacity and ATP production (Tan et al., ). ATP production is calculated from the decrease in respiration rate when ATPase activity is blocked, meaning the difference between de OCR before and after oligomycin injection. Proton leak is the respiration that remains after ATPase activity inhibition by oligomycin minus the non-mitochondrial respiration. Maximal respiration rate (MRR) is the difference between maximal mitochondrial respiration measured as OCR in mitochondrial uncoupled by FCCP and the non-mitochondrial OCR measured in mitochondrial exposed to rotenone/antimycin A. Spare respiratory capacity (SRC), an indicator of the bioenergetics capacity, is the difference between OCR-F and OCR-B. The respiratory control ratio (RCR) is an index of mitochondrial coupling and is calculated as the ratio between the OCR-F and OCR-O.
ROS Measurement
The total ROS content of cells was assessed using 2',7'-dichlorodihydrofluorescein diacetate (H2DCF-DA). Fluorescence intensity was measured (excitation 495 nm, emission 527 nm). Mitochondrial superoxide production was monitored using the superoxide-sensitive probe MitoSOX™ Red (Molecular Probes). Fluorescence intensity was measured (excitation 510 nm, emission 580 nm). A microplate reader Tecan Infinite® 200 (Tecan, MÄnnedorf, Switzerland) was used.
Measurement of Mitochondrial Membrane Potential (ΔΨm)
Mitochondrial membrane potential was measured using JC-1 probe (Invitrogen), according to manufacturer's protocol. Fluorescence intensity was measured (excitation 560 nm, emission 595 nm for J-aggregates in healthy cells and excitation 485 nm, emission 535 nm for monomers in apoptotic cells) in the microplate reader Tecan Infinite® 200 (Tecan, MÄnnedorf, Switzerland). The ratio of fluorescent intensity of J-aggregates to fluorescence intensity of monomers was used as an indicator of cell health.
Lactate and Pyruvate Measurement
L-lactate and pyruvate were measured in the culture medium with L-Lactate Assay Kit (Abcam, #ab65330) and Pyruvate Assay Kit (Abcam, ab65342), respectively, in 96-well plates, according to the manufacturer's instructions.
Statistics
Statistical analyses were performed with GraphPad Prism (Version 7.0). One-way ANOVA test with Bonferroni post-hoc test were performed as shown. Values are plotted as the mean ± SD. p-Values are based on the following pattern: ns (not significant; p > 0.05), * (significant; p < 0.05), ** (significant; p < 0.01), and *** (significant; p < 0.001).
Results
Reduction of MPV17 Protein Levels Affect Cell Metabolic Activity
To investigate the molecular effects of MPV17 loss of function, we developed a cellular model to analyse protein content, subcellular localization and mitochondrial bioenergetics. We assessed the endogenous levels of MPV17 in three different cell lines by immunoblotting analyses. The protein was expressed at similar levels in HEK293T, H4, and HepG2 (Figure 1A). To generate an MPV17 knockdown (MPV17KD) cell line, we used a lentiviral vector encoded shRNA. In parallel, we generated a control cell line using a scrambled shRNA sequence. In the KD cells, the levels of MPV17 were efficiently reduced when compared to the control cells (Figure 1B). We restored MPV17 levels in the MPV17KD cell line by transfecting wild-type MPV17 (Figure 1C; Supplementary Figure 3). Next, we assessed the metabolic activity of cells knocked down for MPV17 using the resazurin reduction assay. Under basal conditions, MPV17KD cells showed decreased resazurin reduction by 38.5% when compared to control. As resazurin reduction can be used as a measure of cell viability but also as an indicator of overall mitochondrial function, this indicates MPV17KD reduced cell viability and/or deficient mitochondrial oxidative capacity, as differences in cell proliferation were not detected by cell count (Figure 1D). Mitochondrial DNA (mtDNA) determination, using RT-qPCR, in MPV17KD cells showed no significant difference when compared to cells expressing wild-type levels of MPV17 (Figure 1E).
Figure 1
Reduction of MPV17 Protein Levels Decreases Mitochondrial Bioenergetics
To understand how MPV17 knockdown affects mitochondrial function, we analyzed the cellular respiratory profile using a Seahorse Flux Analyser, through the successive addition of oligomycin, FCCP and rotenone/antimycin A to block mitochondrial electron transport chain. In control conditions, oxygen consumption rate (OCR-B) is inhibited when ATP synthesis is blocked by oligomycin (OCR-O) and is stimulated when oxidation is uncoupled from phosphorylation by FCCP (OCR-F). OCR-B, -O, and –F were significantly lower (p < 0.001) in the MPV17KD cells when compared with the control (Figures 2A,B). Values of OCR-B, -O, and -F were used to calculate other bioenergetics parameters, i.e., proton leak, maximal respiration rate (MRR), spare respiratory capacity (SRC), and estimate ATP production (Figure 2C). MRR, SRC, and ATP production were all significantly lower in the mutant cells when compared with control. Overall, these results indicate decreased mitochondrial bioenergetic capacity in MPV17KD cells. Despite the slightly lower ECAR in basal conditions (Figures 2D,E), in stressful conditions (by the addition of oligomycin and FCCP) MPV17KD cells showed lower ECAR when compared with WT (CTRL) cells. OCR-B/ECAR ratio was decreased in MPV17KD cells (Figure 2F), and the bioenergetic map that plots OCR against ECAR reveal a limited glycolytic capacity of KD cells when compared with control. Overall, these results suggest that loss of MPV17 causes significant impairment in mitochondrial respiratory capacity and limited metabolic capacity to meet high energy demand (Figure 2G).
Figure 2
MPV17 Mutants Maintain Mitochondrial Localization but Are Less Stable Than WT Protein
To investigate the effect of selected disease-associated mutations, we expressed MPV17 mutants in MPV17KD cells. Considering the clinical criteria, and in silico analysis of the predicted effects of the amino acid substitutions, we selected five MDDS-related mutations for in depth molecular studies: arginine to glutamine substitution at position 50 (MPV17R50Q) which is the most frequent mutation; arginine to tryptophan substitution at position 50 (MPV17R50W) which usually associates with Navajo neurohepatopathy with a severe clinical prognosis; deletion of three amino acids encoding for Gly79_Thr81 (MPV1779−81del) which are predicted to delete a putative kinase C phosphorylation site; substitution of glycine for an arginine at position 94 (MPV17G94R) in the first amino acid of exon 5 which has a possible effect on splicing; and a substitution of a serine for phenylalanine at position 170 (MPV17S170F), localized at the C-terminal of the protein and possibly abolishes a phosphorylation site (Supplementary Material 2) Wild-type HA-tagged MPV17 (MPV17WT) was used as control.
To determine whether MPV17 mutations affect the subcellular localization of the protein, HEK293T knockdown cells were transfected with the various constructs and imaged using confocal microscopy. As expected, HA-tagged MPV17WT displayed a punctate distribution in the cytoplasm that colocalized with the mitochondrial marker. Additionally, we observed mitochondrial localization for all mutant proteins, except for the MPV17G94R mutant. In cells expressing the MPV17G94R mutant protein, the mitochondria network was ill-defined or not staining with MitoTracker Red CMXRos, a ΔΨm-sensitive dye (Tan et al., ). These mitochondria devoid of Mitotracker staining coexisted side by side with non-transfected cells with normal mitochondria staining. These results suggested that the MPV17G94R mutation induced the dissipation of ΔΨm in the transfected cells (Figure 3A).
Figure 3
Although the mutations did not affect mitochondrial targeting of the protein, we hypothesized that they might affect the localization of the protein inside the mitochondria. To test this, we performed subcellular localization studies in cells expressing the different MPV17 variants. Both WT and mutant protein were enriched in the mitochondrial fraction (Figure 3B). However, MPV17 levels were markedly decreased for all mutants, especially for MPV1779−81del, MPV17G94R, and MPV17S170F (<10% of levels found in MPV17WT) (Figure 3C). Mitoplasts were obtained by removing the mitochondrial outer membrane from isolated mitochondria. As expected, wild-type protein (MPV17WT) was found in mitoplasts, consistent with its inner mitochondrial membrane localization. In the presence of MPV17G94R mutation, the mutant protein was not found in the mitoplast fraction (Figure 3D). All other mutations had no effect on the localization of the protein. To investigate if the reduced levels of mutant MPV17 protein were a consequence of decreased transcription, mRNA levels were assessed by RT-qPCR. When compared with WT, mutants mRNA levels are lower particularly for G94R but also for R50W and S170F mutants (Figure 3E), suggesting a compromised transcription for these mutants. The G94R mutant is mislocalized as well as transcribed at low levels when compared with the WT protein.
MPV17 Mutations Affect Mitochondrial Function and the Metabolic Capacity of the Cells
To further understand how mutations in MPV17 result in mitochondrial dysfunction, MPV17KD cells were transfected with either wild-type (WT) or mutant forms of MPV17. Several parameters of mitochondrial function were assessed. We first determined whether MPV17 mutations affected cell viability/proliferation. Although we cannot completely exclude that this might be explained by differences in cell proliferation (which we could not detect by comparing cell numbers), cytotoxicity measured by LDH release did not account for the differences observed (Figure 4A). Cells expressing R50W, G79-T81del, and G94R mutations displayed decreased resazurin reduction capacity (Figure 4B). As mentioned before, this suggests either decreased cell viability and/or deficient mitochondrial oxidative capacity. Mitochondrial membrane potential (ΔΨm) is a key parameter of mitochondrial function and a useful indicator of cell viability. As suggested by microscopy experiments, the MPV17G94R mutation appears to cause the loss of mitochondrial membrane potential. In steady-state cells, only the MPV17S170F mutation seemed to cause significant loss of ΔΨm, as seen by decreased J-aggregates/monomers (red/green fluorescence) ratio when compared with cells expressing wild-type protein (Figure 4C). To investigate the effect of the various MPV17 forms on glycolytic activity, we measured the levels of lactate, pyruvate and the lactate/pyruvate ratio. Mutations in MPV17 were associated with increased dependence on glycolytic ATP production (Figure 4D) and showed decreased mitochondrial competency, as evaluated by the resazurin assay. Cells expressing the MPV17R50W mutation showed an 8-fold increase in lactate levels when compared with control cells. These differences in lactate levels were not accompanied by significant differences in pyruvate. LDH activity did not explain the differences in lactate levels between the samples (Figure 4A). Mitochondria are the major source of ROS in the cell, as a byproduct of the oxidative phosphorylation (OXPHOS) process. Thus, we assessed the possible effect of MPV17 mutations on intracellular ROS levels. We used the DCFH-DA dye in order to measure ROS levels in cells, and MitoSOX to evaluate mitochondrial superoxide () levels. Cells expressing most MPV17 mutants showed a tendency toward elevated levels of intracellular ROS when compared with control cells. In the case of mitochondrial superoxide, this difference was detectable both in basal conditions (with the exception of cells expressing MPV17R50W and MPV17S170F) and in H2O2-induced stress conditions (Figure 4E).
Figure 4
Discussion
MPV17 is a non-selective channel of the inner mitochondrial membrane. Its exact function as well as the mechanism leading to mtDNA instability is unknown. To gain insight into the molecular mechanisms associated with disease in patients carrying MPV17 mutations, we evaluated the cellular consequences of MPV17 deficiency.
Lack of functional MPV17 protein in mitochondria may initiate tissue-specific cell-death and degeneration (Meyer zum Gottesberge and Felix, ). Our experiments indicate that MPV17KD cells exhibited reduced viability while proliferation was not affected (Figure 1D) (O'Brien et al., ).
In the absence of MPV17, cells exhibit decreased oxygen consumption, indicating respiratory chain and OXPHOS impairment. Extracellular acidification rate (ECAR) is a surrogate measure of the production of lactic acid and, therefore, an indicator of anaerobic glycolysis. The OCR-B/ECAR ratio indicates the balance between mitochondrial respiration (OXPHOS) and glycolysis.
Frequently, defects in the former are associated with compensatory increase in the latter. As expected, OCR-B/ECAR ratio was decreased in MPV17KD cells, reflecting impaired mitochondrial respiratory capacity and limited glycolytic capacity (Figures 2F,G). However, and contrary to our initial expectations, stress-induced ECAR in MPV17 deficient cells were lower than in control cells, reflecting an energetic quiescent metabolic profile. These results are consistent with those already observed in MPV17−/− murine podocytes (Casalena et al., ) that showed ≈40% reduced ECAR when compared with wild-type. This is compatible with a role of MPV17 in mitochondrial homeostasis and may indicate that MPV17 is important for the activation of compensatory mechanisms in situations of impaired OXPHOS.
To date, almost 50 pathogenic mutations have been associated with MPV17-related MDDS. Here, we describe metabolic parameters in five cell lines expressing MPV17 mutations, which enable the study of protein levels, localization, and effects on overall bioenergetics. Although no clear genotype-phenotype correlation currently exists, some mutations carry a better prognosis, implying that different mutations may cause distinct consequences in MPV17. Our findings demonstrate that all mutations affect the levels of the MPV17 protein (Figure 3B). Using two different techniques, we show that the MPV17 mutations studied do not affect proper targeting of the protein to the mitochondria (Figures 3A,C). Nevertheless, MPV17G94R behaved differently, as this mutant protein is found in mitochondrial enriched extract but it is not detected in the mitoplast (Figure 3D), reflecting that G94R mutation abolishes protein localization in the inner mitochondrial membrane. Although further studies are needed, one possible explanation is that the mutant protein probably is retained in the intermembrane space (IMS) or is loosely attached to the outer mitochondrial membrane (OMM).
In normally functioning mitochondria, the MPV17 channel is closed (Antonenkov et al., ). Nevertheless, the gating properties of the channel are modified by membrane potential, phosphorylation status and redox state. The regulation of the ΔΨm is vital to sustain mitochondrial homeostasis, as many functions like mitochondrial biogenesis, dynamics, ROS production and mitophagy depend on it. The MPV17G94R mutation affects ΔΨm as seen in microscopy studies using the Mitotracker Red dye. The fact that this was not obvious with the use of JC-1, another probe for measuring ΔΨm, is not clear, but might be due to different sensitivity of this probe, especially in the context of a heterogeneous population of transfected cells. The new findings regarding the effect of the G94R mutation suggests that mislocalization of MPV17 may induce MPV17 dysfunction. Moreover, together with loss of membrane potential we reasoned that this affects the protein function and is likely to play a major role in MPV17 disease in patients affected by this particular mutation.
Mutations in phosphorylation sites are responsible for several neurodegenerative disorders, including MDDS cases (Uusimaa et al., ). The amino acid sequence of the MPV17 protein has three predicted phosphorylation sites (Thr53, Thr80, and Ser170), and mutations in all three were found to cause MDDS. Studies employing mutants that mimic phosphorylated and unphosphorylated forms of the protein indicate that phosphorylation of MPV17 may lead to a complete closing of the channel, particularly in reducing conditions (Antonenkov et al., ). We demonstrate that MPV17S170F causes a lower steady-state ΔΨm. However, MPV1779−81del, a mutation that also abolishes a phosphorylation site, does not appear the affect ΔΨm (Figure 4C). Together, this suggests that mutations on phosphorylation sites do not have a uniform pathogenic effect on the protein. Second, cells expressing the MPV17S170F display lower ROS levels when compared with other mutants (Figure 4E). We hypothesize that the mutation creates a “leaky” MPV17 channel, unable to stay completely closed to preserve a normally high ΔΨm. At the same time, this may also reflect some intact compensatory mechanisms to prevent excessive ROS accumulation.
Mitochondrial DNA and proteins are particularly sensitive to ROS as they are located in close proximity to the respiratory chain. Except for MPV17S170F, all mutants have increased levels of mitochondrial ROS (superoxide) production. However, this increase is only significant in MPV1779−81del and MPV17G94R mutants. We interpret this as a common end-stage effect of different mutations, independently of the exact mechanisms by which each mutation leads to protein dysfunction. A cascade of negative effects of ROS on mtDNA stability and OXPHOS function is expected and, therefore, a common clinical phenotype of the MPV17-associated MDDS.
Although the molecular implications of our findings in simple cell systems may not be straightforward, the cellular abnormalities identified inform on possible cellular pathologies taking place. We hypothesize the effects of the mutations might be even larger in such tissues, given their high energy demand, increased number of mitochondria and dependence of OXPHOS for ATP supply.
In conclusion, we demonstrate that MPV17 mutations can cause respiratory chain impairment by different mechanisms, including increased ROS production, decreased oxygen consumption, potential dissipation and sub-mitochondrial mislocalization. For each mutation, isolated or several simultaneous mechanisms may contribute to a cumulative deleterious effect that leads to disease. Future studies, based on our findings, will inform on the best strategy for targeting the effects of MPV17 mutations in MDDS.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
SJ, TC-O, and TO designed the experimental work. SJ performed all the experiments, data analysis except Seahorse mitochondrial respiration studies, mtDNA copy number, and drafted the manuscript. PG performed and analyzed Seahorse mitochondrial respiration studies. MS and MG performed and interpreted mtDNA copy number. PG and RM trained SJ in performing experimental procedures and supervised the experimental work. TO and AR supervised the study and critically reviewed the manuscript. All authors reviewed, commented, and approved the manuscript.
Funding
SJ was funded by Fundação Calouste Gulbenkian. TO was supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany's Excellence Strategy - EXC 2067/1- 390729940). AR and TC-O were funded by Fundação para a Ciência e a Tecnologia and European Regional Development Fund (CENTRO-01-0145-FEDER-000012- HealthyAging2020), (POCI-010145-FEDER-007440), and (POCI-01-145-FEDER-29297). MS and MG were funded by UID/NEU/04539/2019 (Strategic Plan CNC.IBILI).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fncel.2021.641264/full#supplementary-material
References
1
AlSamanA.TomoumH.InvernizziF.ZevianiM. (2012). Hepatocerebral form of mitochondrial DNA depletion syndrome due to mutation in MPV17 gene. Saudi J. Gastroenterol.18, 285–289. 10.4103/1319-3767.98439
2
Amati-BonneauP.ValentinoM. L.ReynierP.GallardoM. E.BornsteinB.BoissièreA.et al. (2008). OPA1 mutations induce mitochondrial DNA instability and optic atrophy “plus” phenotypes. Brain131, 338–351. 10.1093/brain/awm298
3
AntonenkovV. D.IsomursuA.MennerichD.VapolaM. H.WeiherH.KietzmannT.et al. (2015). The human mitochondrial DNA depletion syndrome gene MPV17 encodes a non-selective channel that modulates membrane potential. J. Biol. Chem.290, 13840–13861. 10.1074/jbc.M114.608083
4
AntounG.McBrideS.VanstoneJ. R.LinesM. (2015). Detailed biochemical and bioenergetic characterization of FBXL4-related encephalomyopathic mitochondrial DNA depletion. JIMD Rep.27, 1–9. 10.1007/8904_2015_491
5
BesseA.WuP.BruniF.DontiT.GrahamB. H.CraigenW. J.et al. (2015). The GABA transaminase, ABAT, is essential for mitochondrial nucleoside metabolism. Cell Metab.21, 417–427. 10.1016/j.cmet.2015.02.008
6
Bijarnia-MahayS.MohanN.GoyalD.VermaI. C. (2014). Mitochondrial DNA depletion syndrome causing liver failure. Indian Pediatr.51, 1–2. 10.1007/s13312-014-0475-z
7
BinderC. J.WeiherH.ExnerM.KerjaschkiD. (2011). Glomerular overproduction of oxygen radicals in Mpv17 gene-inactivated mice causes podocyte foot process flattening and proteinuria. Am. J. Pathol.154, 1067–1075. 10.1016/S0002-9440(10)65359-X
8
BlakelyE. L.ButterworthA.HaddenR. D.BodiI.HeL.McFarlandR.et al. (2012). MPV17 mutation causes neuropathy and leukoencephalopathy with multiple mtDNA deletions in muscle. Neuromuscul. Disord.22, 587–591. 10.1016/j.nmd.2012.03.006
9
BourdonA.MinaiL.SerreV.JaisJ. P.SarziE.AubertS.et al. (2007). Mutation of RRM2B, encoding p53-controlled ribonucleotide reductase (p53R2), causes severe mitochondrial DNA depletion. Nat. Genet.39, 776–780. 10.1038/ng2040
10
CasalenaG.KrickS.DaehnI.YuL.JuW.ShiS.et al. (2014). Mpv17 in mitochondria protects podocytes against mitochondrial dysfunction and apoptosis in vivo and in vitro. AJP Ren. Physiol.306, F1372–F1380. 10.1152/ajprenal.00608.2013
11
ChoiY. R.HongY. B.JungS. C.LeeJ. H.KimY. J.ParkH. J.et al. (2015). A novel homozygous MPV17 mutation in two families with axonal sensorimotor polyneuropathy. BMC Neurol.15:179. 10.1186/s12883-015-0430-1
12
Dalla RosaI.CámaraY.DurigonR.MossC. F.VidoniS.AkmanG.et al. (2016). MPV17 loss causes deoxynucleotide insufficiency and slow DNA replication in mitochondria. PLoS Genet.12:e1005779. 10.1371/journal.pgen.1005779
13
Echaniz-LagunaA.ChassagneM.CeresuelaJ.RouvetI.PadetS.AcquavivaC.et al. (2012). Complete loss of expression of the ANT1 gene causing cardiomyopathy and myopathy. J. Med. Genet.49, 146–150. 10.1136/jmedgenet-2011-100504
14
El-HattabA. W.LiF. Y.SchmittE.ZhangS.CraigenW. J.WongL. J. (2010). MPV17-associated hepatocerebral mitochondrial DNA depletion syndrome: new patients and novel mutations. Mol. Genet. Metab.99, 300–308. 10.1016/j.ymgme.2009.10.003
15
El-HattabA. W.WangJ.DaiH.AlmannaiM.StaufnerC.AlfadhelM.et al. (2018). MPV17-related mitochondrial DNA maintenance defect: new cases and review of clinical, biochemical, and molecular aspects. Hum. Mutat.39, 1–29. 10.1002/humu.23387
16
ElpelegO.MillerC.HershkovitzE.Bitner-GlindziczM.Bondi-RubinsteinG.RahmanS.et al. (2005). Deficiency of the ADP-forming succinyl-CoA synthase activity is associated with encephalomyopathy and mitochondrial DNA depletion. Am. J. Hum. Genet.76, 1081–1086. 10.1086/430843
17
GalletP. F.ZachowskiA.JulienR.FellmannP.DevauxP. F.MaftahA. (1999). Transbilayer movement and distribution of spin-labelled phospholipids in the inner mitochondrial membrane. Biochim. Biophys. Acta1418, 61–70. 10.1016/S0005-2736(99)00022-X
18
KajiS.MurayamaK.NagataI.NagasakaH.TakayanagiM.OhtakeA.et al. (2009). Fluctuating liver functions in siblings with MPV17 mutations and possible improvement associated with dietary and pharmaceutical treatments targeting respiratory chain complex II. Mol. Genet. Metab.97, 292–296. 10.1016/j.ymgme.2009.04.014
19
KaradimasC. L.VuT. H.HolveS. A.ChronopoulouP.QuinziiC.JohnsenS. D.et al. (2006). Navajo neurohepatopathy is caused by a mutation in the MPV17 gene. Am. J. Hum. Genet.79, 544–548. 10.1086/506913
20
KimJ.KangE.KimY.KimJ. M.LeeB. H.MurayamaK.et al. (2016). MPV17 mutations in patients with hepatocerebral mitochondrial DNA depletion syndrome. Mol. Genet. Metab. Rep.8, 74–76. 10.1016/j.ymgmr.2016.06.006
21
KornblumC.NichollsT. J.HaackT. B.SchölerS.PeevaV.DanhauserK.et al. (2013). Loss-of-function mutations in MGME1 impair mtDNA replication and cause multisystemic mitochondrial disease. Nat. Genet.45, 214–219. 10.1038/ng.2501
22
MandelH.SzargelR.LabayV.ElpelegO.SaadaA.ShalataA.et al. (2001). The deoxyguanosine kinase gene is mutated in individuals with depleted hepatocerebral mitochondrial DNA. Nat. Genet.29, 337–341. 10.1038/ng746
23
MayrJ. A.HaackT. B.GrafE.ZimmermannF. A.WielandT.HaberbergerB.et al. (2012). Lack of the mitochondrial protein acylglycerol kinase causes sengers syndrome. Am. J. Hum. Genet.90, 314–320. 10.1016/j.ajhg.2011.12.005
24
MendelsohnB. A.MehtaN.HameedB.PekmezciM.PackmanS.RalphJ. (2013). Adult-onset fatal neurohepatopathy in a woman caused by MPV17 mutation. JIMD Rep.21, 37–41. 10.1007/8904_2013_267
25
MerkleA. N.NasceneD. R.McKinneyA. M. (2012). MR imaging findings in the reticular formation in siblings with MPV17-related mitochondrial depletion syndrome. Am. J. Neuroradiol.33, E34–E35. 10.3174/ajnr.A2455
26
Meyer zum GottesbergeA. M.FelixH. (2005). Abnormal basement membrane in the inner ear and the kidney of the Mpv17−/− mouse strain: ultrastructural and immunohistochemical investigations. Histochem. Cell Biol.124, 507–516. 10.1007/s00418-005-0027-7
27
Navarro-SastreA.Martín-HernándezE.CamposY.QuintanaE.MedinaE.de Las HerasR. S.et al. (2008). Lethal hepatopathy and leukodystrophy caused by a novel mutation in MPV17 gene: description of an alternative MPV17 spliced form. Mol. Genet. Metab.94, 234–239. 10.1016/j.ymgme.2008.01.012
28
NaviauxR. K.NyhanW. L.BarshopB. A.PoultonJ.MarkusicD.KarpinskiN. C.et al. (1999). Mitochondrial DNA polymerase gamma deficiency and mtDNA depletion in a child with Alpers' syndrome. Ann. Neurol.45, 54–58. 10.1002/1531-8249(199901)45:1<54::AID-ART10>3.0.CO;2-B
29
NishinoI.SpinazzolaA.HiranoM. (1999). Thymidine phosphorylase gene mutations in MNGIE, a human mitochondrial disorder. Science283, 689–692. 10.1126/science.283.5402.689
30
NogueiraC.NunesJ.EvangelistaT.FattoriF.TessaA.PereiraC.et al. (2007). A new mtDNA-tRNA(Glu) mutation (14728T>C) presenting a late-onset mitochondrial encephalomyopathy. Mitochondrion7, 396–398. 10.1016/j.mito.2007.08.002
31
O'BrienJ.WilsonI.OrtonT.PognanF. (2000). Investigation of the Alamar Blue (resazurin) fluorescent dye for the assessment of mammalian cell cytotoxicity. Eur. J. Biochem.267, 5421–5426. 10.1046/j.1432-1327.2000.01606.x
32
OstergaardE.ChristensenE.KristensenE.MogensenB.DunoM.ShoubridgeE. A.et al. (2007). Deficiency of the α subunit of succinate–coenzyme A ligase causes fatal infantile lactic acidosis with mitochondrial DNA depletion. Am. J. Hum. Genet.81, 383–387. 10.1086/519222
33
PariniR.FurlanF.NotarangeloL.SpinazzolaA.UzielG.StrisciuglioP.et al. (2009). Glucose metabolism and diet-based prevention of liver dysfunction in MPV17 mutant patients. J. Hepatol.50, 215–221. 10.1016/j.jhep.2008.08.019
34
ReyesA.MelchiondaL.NascaA.CarraraF.LamanteaE.ZanoliniA.et al. (2015). RNASEH1 mutations impair mtDNA replication and cause adult-onset mitochondrial encephalomyopathy. Am. J. Hum. Genet.97, 186–193. 10.1016/j.ajhg.2015.05.013
35
RonchiD.Di FonzoA.LinW.BordoniA.LiuC.FassoneE.et al. (2013). Mutations in DNA2 link progressive myopathy to mitochondrial DNA instability. Am. J. Hum. Genet.92, 293–300. 10.1016/j.ajhg.2012.12.014
36
RouzierC.BannwarthS.ChaussenotA.ChevrollierA.VerschuerenA.Bonello-PalotN.et al. (2012). The MFN2 gene is responsible for mitochondrial DNA instability and optic atrophy ‘plus’ phenotype. Brain135, 23–34. 10.1093/brain/awr323
37
SaadaA.ShaagA.MandelH.NevoY.ErikssonS.ElpelegO. (2001). Mutant mitochondrial thymidine kinase in mitochondrial DNA depletion myopathy. Nat. Genet.29, 342–344. 10.1038/ng751
38
SarziE.BourdonA.ChrétienD.ZarhrateM.CorcosJ.SlamaA.et al. (2007a). Mitochondrial DNA depletion is a prevalent cause of multiple respiratory chain deficiency in childhood. J. Pediatr.150, 531–534, 534.e1–534.e6. 10.1016/j.jpeds.2007.01.044
39
SarziE.GoffartS.SerreV.ChrétienD.SlamaA.MunnichA.et al. (2007b). Twinkle helicase (PEO1) gene mutation causes mitochondrial DNA depletion. Ann. Neurol.62, 579–587. 10.1002/ana.21207
40
SpinazzolaA. (2011). Mitochondrial DNA mutations and depletion in pediatric medicine. Semin. Fetal Neonatal Med.16, 190–196. 10.1016/j.siny.2011.04.011
41
SpinazzolaA.InvernizziF.CarraraF.LamanteaE.DonatiA.DiroccoM.et al. (2009). Clinical and molecular features of mitochondrial DNA depletion syndromes. J. Inherit. Metab. Dis.32, 143–158. 10.1007/s10545-008-1038-z
42
SpinazzolaA.SanterR.AkmanO. H.TsiakasK.SchaeferH.DingX.et al. (2008). Hepatocerebral form of mitochondrial DNA depletion syndrome. Arch. Neurol.65, 1108–1113. 10.1001/archneur.65.8.1108
43
SpinazzolaA.ViscomiC.Fernandez-VizarraE.CarraraF.D'AdamoP.CalvoS.et al. (2006). MPV17 encodes an inner mitochondrial membrane protein and is mutated in infantile hepatic mitochondrial DNA depletion. Nat. Genet.38, 570–575. 10.1038/ng1765
44
StilesA. R.SimonM. T.StoverA.EftekharianS.KhanlouN.WangH. L.et al. (2016). Mutations in TFAM, encoding mitochondrial transcription factor A, cause neonatal liver failure associated with mtDNA depletion. Mol. Genet. Metab.119, 91–99. 10.1016/j.ymgme.2016.07.001
45
SuomalainenA.IsohanniP. (2010). Mitochondrial DNA depletion syndromes–many genes, common mechanisms. Neuromuscul. Disord.20, 429–437. 10.1016/j.nmd.2010.03.017
46
TanM. G. K.OoiL. L. P. J.AwS. E.HuiK. M. (2004). Cloning and identification of hepatocellular carcinoma down-regulated mitochondrial carrier protein, a novel liver-specific uncoupling protein. J. Biol. Chem.279, 45235–45244. 10.1074/jbc.M403683200
47
UusimaaJ.EvansJ.SmithC.ButterworthA.CraigK.AshleyN.et al. (2013). Clinical, biochemical, cellular and molecular characterization of mitochondrial DNA depletion syndrome due to novel mutations in the MPV17 gene. Eur. J. Hum. Genet.22, 184–191. 10.1038/ejhg.2013.112
48
VenegasV.WangJ.DimmockD.WongL. J. (2011). Real-time quantitative PCR analysis of mitochondrial DNA content. Curr. Protoc. Hum. Genet. Chapter 19:Unit 19.7. 10.1002/0471142905.hg1907s68
49
ViscomiC.SpinazzolaA.MaggioniM.Fernandez-VizarraE.MassaV.PaganoC.et al. (2009). Early-onset liver mtDNA depletion and late-onset proteinuric nephropathy in Mpv17 knockout mice. Hum. Mol. Genet.18, 12–26. 10.1093/hmg/ddn309
50
WongL. J.Brunetti-PierriN.ZhangQ.YazigiN.BoveK. E.DahmsB. B.et al. (2007). Mutations in the MPV17 gene are responsible for rapidly progressive liver failure in infancy. Hepatology46, 1218–1227. 10.1002/hep.21799
51
YamazakiT.MurayamaK.ComptonA. G.SugianaC.HarashimaH.AmemiyaS.et al. (2014). Molecular diagnosis of mitochondrial respiratory chain disorders in Japan: focusing on mitochondrial DNA depletion syndrome. Pediatr. Int.56, 180–187. 10.1111/ped.12249
52
YoungM. J.HumbleM. M.DeBalsiK. L.SunK. Y.CopelandW. C. (2015). POLG2 disease variants: analyses reveal a dominant negative heterodimer, altered mitochondrial localization and impaired respiratory capacity. Hum. Mol. Genet.24, 5184–5197. 10.1093/hmg/ddv240
53
ZorovD. B.JuhaszovaM.SollottS. J. (2014). Mitochondrial reactive oxygen species (ROS) and ROS-induced ROS release. Physiol. Rev.94, 909–950. 10.1152/physrev.00026.2013
54
ZwackaR. M.ReuterA.PfaffE.MollJ.GorgasK.KarasawaM.et al. (1994). The glomerulosclerosis gene Mpv17 encodes a peroxisomal protein producing reactive oxygen species. EMBO J.13, 5129–5134. 10.1002/j.1460-2075.1994.tb06842.x
Summary
Keywords
Mpv17 mutations, mitochondrial depletion syndrome, mitochondrial dysfunction, protein mislocation, neurode generation
Citation
Jacinto S, Guerreiro P, de Oliveira RM, Cunha-Oliveira T, Santos MJ, Grazina M, Rego AC and Outeiro TF (2021) MPV17 Mutations Are Associated With a Quiescent Energetic Metabolic Profile. Front. Cell. Neurosci. 15:641264. doi: 10.3389/fncel.2021.641264
Received
13 December 2020
Accepted
08 February 2021
Published
17 March 2021
Volume
15 - 2021
Edited by
Riccardo Fesce, Humanitas University, Italy
Reviewed by
Nicholas D. Mazarakis, Imperial College London, United Kingdom; Massimo Zeviani, University of Padua, Italy
Updates
Copyright
© 2021 Jacinto, Guerreiro, de Oliveira, Cunha-Oliveira, Santos, Grazina, Rego and Outeiro.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Tiago F. Outeiro touteir@gwdg.de
This article was submitted to Cellular Neurophysiology, a section of the journal Frontiers in Cellular Neuroscience
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.