Next Article in Journal
Genetic Diversity Analysis of Banana Cultivars (Musa sp.) in Saudi Arabia Based on AFLP Marker
Next Article in Special Issue
Revealing the History and Mystery of RNA-Seq
Previous Article in Journal
Peroxiredoxins and Hypoxia-Inducible Factor-1α in Duodenal Tissue: Emerging Factors in the Pathophysiology of Pediatric Celiac Disease Patients
Previous Article in Special Issue
Cooperative Binding of SRSF3 to Structured 3’ss-α Exon RNA during α Exon Inclusion in the ZO-1 mRNA
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
CIMB
Volume 45
Issue 3
10.3390/cimb45030115
cimb-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of the m.13513G>A Variant on the Functions of the OXPHOS System and Cell Retrograde Signaling

by
Dita Kidere
1,
Pawel Zayakin
1,
Diana Livcane
1,
Marina Makrecka-Kuka
2,
Janis Stavusis
1,
Baiba Lace
1,3,
Tsu-Kung Lin
4,5,
Chia-Wei Liou
4,5 and
Inna Inashkina
1,*
1
Latvian Biomedical Research and Study Centre, LV-1067 Riga, Latvia
2
Latvian Institute of Organic Synthesis, LV-1006 Riga, Latvia
3
Children’s Clinical University Hospital, LV-1004 Riga, Latvia
4
Department of Neurology, Kaohsiung Chang Gung Memorial Hospital and Chang Gung University College of Medicine, Kaohsiung 83305, Taiwan
5
Center for Mitochondrial Research and Medicine, Kaohsiung Chang Gung Memorial Hospital and Chang Gung University College of Medicine, Kaohsiung 83301, Taiwan
*
Author to whom correspondence should be addressed.
Curr. Issues Mol. Biol. 2023, 45(3), 1794-1809; https://doi.org/10.3390/cimb45030115
Submission received: 23 December 2022 / Revised: 8 February 2023 / Accepted: 16 February 2023 / Published: 22 February 2023
(This article belongs to the Special Issue Studying the Function of RNAs Using Omics Approaches)
Browse Figures
Versions Notes

Abstract

:
Mitochondria are involved in many vital functions in living cells, including the synthesis of ATP by oxidative phosphorylation (OXPHOS) and regulation of nuclear gene expression through retrograde signaling. Leigh syndrome is a heterogeneous neurological disorder resulting from an isolated complex I deficiency that causes damage to mitochondrial energy production. The pathogenic mitochondrial DNA (mtDNA) variant m.13513G>A has been associated with Leigh syndrome. The present study investigated the effects of this mtDNA variant on the OXPHOS system and cell retrograde signaling. Transmitochondrial cytoplasmic hybrid (cybrid) cell lines harboring 50% and 70% of the m.13513G>A variant were generated and tested along with wild-type (WT) cells. The functionality of the OXPHOS system was evaluated by spectrophotometric assessment of enzyme activity and high-resolution respirometry. Nuclear gene expression was investigated by RNA sequencing and droplet digital PCR. Increasing levels of heteroplasmy were associated with reduced OXPHOS system complex I, IV, and I + III activities, and high-resolution respirometry also showed a complex I defect. Profound changes in transcription levels of nuclear genes were observed in the cell lines harboring the pathogenic mtDNA variant, indicating the physiological processes associated with defective mitochondria.

1. Introduction

Mitochondria are dynamic cell structures that possess their own intrinsic complement of circular DNA (mtDNA), ribosomes, and proteins [1]. Mitochondria play essential roles in eukaryotic aerobic metabolism, including in ATP synthesis and fatty acid oxidation, as well as being involved in the biosynthesis of intermediary metabolites, the regulation of calcium homeostasis, citric acid cycle metabolism, heat generation, and apoptosis [2,3,4,5,6,7,8,9]. Mitochondria are also involved in retrograde signaling, a communication system in which mitochondria generate signals that regulate nuclear gene expression [10]. The retrograde response may be stimulated by the ATP/ADT ratio, disruption of the mitochondrial membrane potential, reactive oxygen species (ROS) production, or general cellular stress [11,12]. Mitochondrial dysfunction in mammalian cells results in reductions in mitochondrial membrane potential, incrementally altering intracellular calcium levels [10]. These alterations result in the activation of calcium-dependent kinases and phosphatases (e.g., calcineurin, PKC, CamKIV, JNK, and MAPK), culminating in the activation of various transcription factors, such as ATF2, CEBP/δ, NFAT, CREB, Egr-1, CHOP, and NFκB [10,13,14,15,16].
Mutations in human mtDNA have been associated with various multisystem disorders known as mitochondrial diseases, with the severity of these disorders depending on the mutant mtDNA load in cells [17]. Furthermore, pathogenic variants in mtDNA, reduced mtDNA copy number, or disruption of mitochondrial respiratory chain (RC) complexes activate mitochondria-to-nucleus retrograde signaling, which induces a global change in nuclear gene expression, ultimately contributing to various pathologic conditions. These mitochondrial changes cause transcriptional reprogramming of nuclear genes at the epigenetic [18], transcriptional [19,20,21], and post-transcriptional [22,23] levels, although the exact mechanism and components of this cross-talk remain unclear.
Transmitochondrial cytoplasmic hybrid systems (cybrids) are created by fusing mitochondria-enriched platelets with mitochondria-depleted cells [24,25]. These cybrids allow the expression of selected mtDNA sequences against a fixed nuclear DNA (nDNA) background, providing a method for studying the effects of known mtDNA variants on mitochondrial biochemistry and mtDNA-nDNA crosstalk. RNA seq has been used to evaluate the transcriptional profiles of cybrids with different levels of heteroplasmy of the MT-TL1 m.3243A>G variant, showing that quantitative changes in mtDNA mutant levels can qualitatively regulate nuclear gene expression [26].
Isolated complex I deficiency represents the most common mitochondrial respiratory chain defect, involved in mitochondrial disorders, and it has been associated with heterogeneous manifestations, including phenotypic overlaps of mitochondrial encephalomyopathy with lactic acidosis and stroke-like episodes (MELAS), Leigh syndrome, and Leber’s hereditary optic neuropathy (LHON) [27]. A previous report from our laboratory described a patient with Leigh syndrome harboring the mtDNA variant m.13513G>A in the MT-ND5 gene at 60% heteroplasmy level [28]. The present study was designed to investigate the effect of this m.13513G>A variant on the oxidative phosphorylation (OXPHOS) system and cell retrograde signaling. Cybrid cell models differing in heteroplasmy levels were generated, and mitochondrial OXPHOS enzyme activities, oxygen consumption, and transcriptional profiles were analyzed in these cybrid cell lines differing in mutant mtDNA load in comparison to cells harboring the wild-type (WT) variant.

2. Materials and Methods

2.1. Recruitment of Participants and Ethics Statement

A patient with mitochondrial disease and healthy volunteers were recruited for this study. The study protocols were approved by the Central Committee of Medical Ethics of Latvia (protocol No.2016-5, chapter 2, approved 24 November 2016, and protocol No.2020-6 II, chapter-3, approved 22 May 2020), which covers all consent- and data-handling-related issues for genetic research on human subjects. All participants and/or their legal guardians provided written informed consent prior to enrollment. The study conformed to the guidelines of the Declaration of Helsinki.

2.2. Cell Culture, Generation of ρ0 Cells, and Construction of Cybrids

Human 143B osteosarcoma cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM, high glucose, Gibco, Carlsbad, CA, USA) supplemented with 10% heat-inactivated fetal bovine serum (FBS; Gibco) and antibiotic–antimycotic solution (Gibco) at 37 °C in an atmosphere containing 5% CO2. Human 143B ρ0 cells were generated by long-term culture of 143B cells with 50 ng/mL ethidium bromide (Merck KGaA, Darmstadt, Germany), 1 mM pyruvate (Gibco), and 50 μg/mL uridine (Sigma, St. Louis, MO, USA) [29]. The cells were subcultured every 2–3 days and kept at 50–70% confluence. The absence of mtDNA was determined by real-time PCR [30].
Transmitochondrial cybrids were generated by fusing 143B ρ0 cells with human platelets in the absence of pyruvate and uridine, as previously reported [24] and modified [25]. Platelets containing the heteroplasmic m.13513G>A variant were isolated and fused with ρ0 cells in the presence of polyethylene glycol 1500 (50% w/v; Sigma). The mixture of fused cells was allowed to recover for 1 week in ρ0 medium and cultured in pyruvate/uridine-free DMEM, a medium, in which even cybrids with very low Cox-II expression had been shown to grow. On days 14–30 after fusion, successful cybrid clones growing in this medium were isolated. The presence of exogenous mtDNA or specific mtDNA haplogroups in each cybrid was confirmed by Sanger sequencing of isolated DNA [28].

2.3. Mitochondria Isolation

Mitochondria were isolated from cytoplasmic hybrid cells as previously described [31], with modifications. Briefly, cells were suspended in 10 mM Tris (pH 7.4) and homogenized with a glass–glass Dounce homogenizer on ice. Sucrose was added immediately to a concentration of 0.25 M, and the suspension was centrifuged at 600× g for 10 min at 4 °C. The postnuclear supernatant was decanted and again centrifuged at 10,000× g for 10 min at 4 °C. Mitochondria-enriched pellets were washed with a buffer containing 250 mM sucrose, 2 mM EGTA, 40 mM KCl, and 20 mM Tris (pH 7.4); suspended in the same buffer; and stored at −70 °C.

2.4. Enzyme Measurements

Enzymatic activities of citrate synthase (CS), OXPHOS complexes I–IV, and complexes I + III and II + III in isolated mitochondria were measured spectrophotometrically (Uvicon 922, Kontron, Italy) at 37 °C, as previously described [32,33,34], with modifications. Briefly, CS activity was measured by adding 10 μL of isolated mitochondria suspension to a reaction mixture containing 50 mM Tris/HCl (pH 8.1), 0.1% bovine serum albumin (BSA), 0.1% TritonX-100, 0.2 mM 5,5′-dithio-bis(2-nitrobenzoic acid) (DTNB), and 0.15 mM acetyl-CoA. Thiolase activity was recorded for 2–3 min, 0.5 mM oxaloacetate was added, and CS activity was recorded for 8 min at 412 nm.
To measure complex I activity, 40 µL of isolated mitochondria were first disrupted by resuspension in water and incubated for 3 min at 37 °C. Then, a reaction mixture containing 50 mM Tris (pH 8.1), 0.1% BSA, 0.3 mM KCN, and 50 µM decylubiquinone was added, and the reaction was initiated by adding 0.2 mM NADH. The NADH oxidation reaction was allowed to proceed for 10 min, and complex I activity was measured at 340 nm, after which 3 µM rotenone, a complex I inhibitor, was added, and nonspecific NADH oxidation was recorded.
To measure complex II activity, 10 µL of isolated mitochondria was added to a reaction mixture containing 50 mM potassium phosphate buffer (pH 7.8), 2 mM ethylenediaminetetraacetic acid (EDTA), 0.1% BSA, 3 µM rotenone, 80 µM 2,6-dichlorophenolindophenol, 50 µM decylubiquinone, 1 µM antimycin A, 0.2 mM ATP, and 0.3 mM KCN. After preincubation for 10 min, the reaction was started by adding 10 mM succinate and allowed to proceed for 5 min, during which complex II activity was measured at 600 nm. Subsequently, 10 mM malonate was added to stop the reaction and determine nonspecific activity.
Complex III activity was measured by adding 5 µL of isolated mitochondria to a reaction mixture containing 50 mM potassium phosphate buffer (pH 7.8), 0.08% Tween-20, 2 mM EDTA, 0.3 mM KCN, 100 µM cytochrome c, and 200 µM reduced decylubiquinol. Cytochrome c reduction was measured at 550 nm for 4 min, after which 1 µM antimycin A was added to inhibit the reaction and determine the nonspecific activity.
The activity of complexes I + III was measured by incubating 10 µL of isolated mitochondria in water for 3 min at 37 °C, followed by addition of a reaction mixture containing 50 mM Tris (pH 8.1), 0.1% BSA, 0.3 mM KCN, and 50 µM cytochrome c. The reaction was started by adding 0.2 mM NADH and allowed to proceed for 3–4 min, with activity measured at 550 nm. The reaction was stopped by adding 3 µM rotenone, and nonspecific activity was determined.
The activity of complexes II + III was determined by incubating 10 µL of isolated mitochondria with a reaction mixture containing 50 mM potassium phosphate buffer (pH 7.8), 2 mM EDTA, 0.1% BSA, 3 µM rotenone, 0.2 mM ATP, and 0.3 mM KCN at 37 °C. The reactions were started by adding 50 µM cytochrome c and 10 mM succinate. Activity was measured for 3 min, and 10 mM malonate was added to stop the reaction and determine nonspecific activity.
Complex IV activity was measured as described [35], with modifications. The reaction was started by adding 60 µM reduced cytochrome c to a reaction mixture containing 50 mM potassium phosphate buffer (pH 7.0) and 0.1% BSA, and 10 µL of isolated mitochondria. Activity was measured at 550 nm for 10 min. Subsequently, 1 µL of bovine heart mitochondria was added to reach the endpoint of the reaction, and the final absorbance at 550 nm was recorded after 30 min.

2.5. High-Resolution Respirometry

Mitochondrial functionality was determined by high-resolution respirometry using an OxygraphO2k (O2k; Oroboros Instruments, Innsbruck, Austria), according to the manufacturer’s instructions. All measurements were performed in Mir05 media (110 mM sucrose, 60 mM K-lactobionate, 0.5 mM EGTA, 3 mM MgCl2, 20 mM taurine, 10 mM KH2PO4, 20 mM HEPES, and 0.1% BSA essentially free of fatty acids) at 37 °C. Mitochondrial respiration was initially measured in intact cells (ROUTINE state), followed by the addition of 5 µg/mL digitonin to permeabilize the cells. Complex I-linked oxygen consumption was measured using 5 mM pyruvate plus 2 mM malate (LEAK state), followed by subsequent additions of 5 mM ADP (OXPHOS state) and 10 mM glutamate. Complex II-linked mitochondrial respiration was evaluated by adding 10 mM succinate and 0.5 µM rotenone. Residual oxygen consumption (ROX) was determined by adding 2.5 µM antimycin A, and complex IV-linked oxygen consumption was measured by adding 0.5 mM N,N,N′,N′-tetramethyl-p-phenylenediamine dihydrochloride (TMPD) and 2 mM ascorbate. Complex IV activity was inhibited, and chemical background determined by the addition of 100 mM sodium azide. The contribution of each substrate to the respiration rate was determined by calculating the flux control factor using the formula: Flux control factor = 1 − (respiratory rate before the addition of substrate)/(respiratory rate after the addition of substrate).

2.6. RNA Sample Preparation and Next Generation Sequencing

Cultivated cybrid cells containing 50% and 70% heteroplasmy of the m.13513G>A variant and cells containing WT m.13513G were collected and washed with PBS. Total RNA was isolated using RNeasy Mini Kits (QIAGEN Science, Germantown, MD, USA), according to the manufacturer’s instructions. The integrity of the extracted RNA was evaluated by determining the RNA integrity number (RIN) using the Agilent 2100 Bioanalyzer system (Agilent, Santa Clara, CA, USA). Ribosomal RNA was removed from 500 ng aliquots of total RNA from each sample using the Low Input RiboMinus Eukaryote System v2 (Thermo Fisher Scientific, Waltham, MA, USA). cDNA libraries were prepared using the Ion Total RNA-Seq Kits v2 (Thermo Fisher Scientific, USA) and sequenced using the Ion Proton System and Ion PI Chip (Thermo Fisher Scientific, USA), according to the manufacturer’s instructions. All experiments were performed on two biological replicates.

2.7. Bioinformatic Analysis

The adapters were trimmed from sequencing reads with Cutadapt [36]. For inclusion in subsequent analyses, sequencing reads after trimming were required to be of a minimum length of 42 bp, with reads corresponding to rRNA removed with Sortmerna [37]. Sequencing reads were mapped against the human reference genome GRCh38.p12 with STAR 2.7.1a [38], and read counts were calculated with HTSEQ-count [39]. The obtained read counts were normalized, principal component analysis was performed by package PCAtools [40], and differentially expressed genes (DEGs) were estimated using the Bioconductor package DESeq2 [41] and EnhancedVolcano [42] in R. Multiple testing correction was implemented using the Benjamini–Hochberg procedure, with significant DEGs defined as those having an absolute log fold change (FC) > 2 and a false discovery rate (FDR) < 0.05. GO and enrichment analysis were performed with R packages GOstat [43], rentrez [44], GO.db [45], and ShinyGO [46].

2.8. Selfie-Digital RT-PCR

Cell culture media were aspirated from the culture flasks, and each cell layer was lysed in 100ST DNA/RNA/Protein Solubilization Reagent (DireCtQuant, Lleida, Spain) at a concentration of 250,000 cells/mL, to preserve the nucleic acids in the samples. The lysates were incubated at 90 °C for 3 min with shaking and centrifuged at 5000× g for 1 min, and the supernatants were used directly. Strand-specific absolute levels of gene transcription were analyzed using Selfie-digital PCR, as described previously [47]. The reverse transcription step was performed using the reverse primer listed in Table S1 and Maxima H Minus Reverse Transcriptase (Thermo Fisher Scientific, USA). Digital PCR was performed using EvaGreen SuperMix (Bio-Rad, Hercules, CA, USA) and the QX200 ddPCR platform. The number of RNA transcripts was calculated by subtracting the number of amplicons measured in the reaction without reverse transcriptase (RT−) from the number measured in the reaction with reverse transcriptase (RT+) and dividing by (RT-). The results are expressed as the number of transcripts per gene. These experiments were repeated twice, with each containing two technical replicates.

3. Results

The effects of the genetic variant m.13513G>A in the MT-ND5 gene on the OXPHOS system and nuclear gene expression in cells were investigated by constructing cybrid cell lines. Platelets from a patient harboring 60% of mutant mtDNA in blood were fused with mtDNA-free ρ0 cells. Heteroplasmic cybrid cell lines containing 50% and 70% of mutant mtDNA were generated, along with a cybrid cell line harboring only wild-type (WT) mtDNA.

3.1. Functionality of the OXPHOS System

To determine whether different heteroplasmy levels of the m.13513G>A variant impaired the performance of the OXPHOS system and mitochondrial respiration, enzyme activities were measured, and high-resolution respirometry was performed. All mitochondrial OXPHOS complex activities were expressed relative to CS activity, as CS is a marker of mitochondrial content [48], as a ratio called “complex activity” throughout the text. Measurements of the activity of enzymes in the OXPHOS system revealed that increasing levels of heteroplasmy caused dose-dependent reductions in complex I activity, with cell lines containing 50% and 70% mutated m.13513G>A showing 38% and 77% reductions, respectively, in complex I activity compared with WT cells (Figure 1A). Similarly, these mutant cell lines showed 54% and 47% reductions, respectively, in complex III activity (Figure 1B). Furthermore, the cell line harboring 70% heteroplasmy showed 67% and 20% reductions in complex IV and complexes I + III activities relative to WT (Figure 1A). All these reductions in OXPHOS complex activities were statistically significant (p < 0.05).
Evaluation of mitochondrial respiration in intact cells showed that ROUTINE respiration was 72% lower in cells containing 70% of the m.13513G>A variant than in cells containing WT (p < 0.05) (Figure 2). Assessments of mitochondrial functionality in permeabilized cells showed that complex I-dependent oxygen consumption was lower in cell lines with 50% and 70% heteroplasmy. Respiration rates in the OXPHOS state using pyruvate and malate as substrates were 40% and 86% lower in cell lines containing 50% and 70% mutant mtDNA, respectively, than in WT cells (Figure 2). Moreover, respiration rates of both heteroplasmic cell lines did not increase after addition of glutamate (Figure 2 and Figure 3), indicating that the reductions in respiration were directly related to complex I function and not substrate specific metabolism. Further evaluation of the cell line harboring 70% of the m.13513G>A variant showed that rotenone had a significantly reduced effect, whereas succinate had a significantly increased effect on respiration rate (Figure 3). In contrast, the cell line harboring 50% of the m.13513G>A variant showed rotenone and succinate effects similar to the WT cells (Figure 3).

3.2. Effect of the Pathogenic mtDNA Variant on Differential Gene Expression

The effects of the pathogenic mtDNA variants on differential gene expression were assessed by performing transcriptome analyses of cybrid cells containing 50% and 70% of the m.13513G>A variant and cells containing WT mtDNA. Principal component analysis revealed that most variability is explained by the different mutational loads (50% versus 70%) within the heteroplasmic cells (Figure 4A). Furthermore, a second source of variability arises from the difference between WT cells and both cell lines harboring the m.13513G>A variant (Figure 4B).
RNA sequencing analysis of both heteroplasmic cell lines versus WT cells showed that the levels of transcription of 84 nuclear genes were significantly affected (adj. p < 0.05) by the m.13513G>A variant (ArrayExpress accession No. E-MTAB-11035). The corresponding volcano plot is shown in Figure 4C. Evaluation of DEGs, defined as those with absolute logFC >2 and FDR < 0.05 cut-offs, showed that the levels of expression of 31 genes were altered more than fourfold, with 25 of those genes being downregulated (log2 FC > 2) and six being upregulated (log2 FC < –2) (Table S2). All of these genes were nuclear.
To determine the affected molecular pathways and biological functions, the DEGs were subjected to Gene Ontology (GO) enrichment analysis (Figure 5 and Figure 6, Tables S3 and S4). Cells harboring the m.13513G>A variant (both 50% and 70% heteroplasmy) showed enrichment of genes associated with both active and passive transmembrane transporter activity (GRID1, NIPAL4, SLC24A4, FABP3, RYR3, ABCG1), particularly calcium release activity (RYR3, SLC24A4)(adj. p = 0.093). The HDKC1 gene was also highlighted in the GO enrichment analysis (adj. p = 0.093); it encodes kinase, which is involved in carbohydrate (sucrose, fructose, and mannose) metabolism, and is downregulated according to RNAseq results. Cluster analysis of the molecular functions of DEGs, based on GO terms, showed a cluster of genes related to calcium transport/release (Table S3). These findings are in agreement with the GO terms of cellular components, which showed enrichment of integral components of the membrane (Table S3), although it did not reach statistical significance (adj. p = 0.430).
The RNAseq analysis of cells harboring 50% versus 70% of the m.13513G>A variant identified 1247 nuclear genes whose transcription levels were significantly affected (adj. p < 0.05; Figure 4D). 437 of these genes had levels of expression that were altered more than four-fold (Table S2).
Ten (DYDC2, TCEA3, GRID1, HKDC1, SEL1L3, FABP3, RYR3, GRM2, DYPSL4, and ABCG1) of the thirty-one most pronounced DEGs were selected for further evaluation. These genes were chosen because of their putative capacity to be involved in cellular physiological and pathological processes. Genes involved in cancer physiology, however, were excluded, as changes in their expression might be associated with the cancerous nature of ρ0 cells used for the production of cybrid cells. The results of RNA expression analysis of the selected genes were validated using Selfie-digital PCR in freshly propagated cybrid cells (Figure S1). Only four genes yielded results resembling the general trend demonstrated by the RNA seq analysis. Compared with WT cells, the TCEA3, SEL1L3, and ABCG1 genes were upregulated, whereas the FABP3 gene was downregulated, in cells bearing the m.13513G>A variant. The GRM2 gene was upregulated in the cells harboring 50% of the variant, but downregulated in the cells harboring 70% heteroplasmy of the variant, as well as being downregulated in the RNA seq analysis.

4. Discussion

The genetic variant m.13513G>A, first described by Santorelli et al. in a patient with MELAS [49], has been reported to be the most common pathogenic variant in the MT-ND5 gene [50]. This variant has later been identified in many individuals with heterogeneous clinical manifestations, including Leigh syndrome [51], Wolff Parkinson-White (WPW) syndrome [52], and LHON, as well as overlapping phenotypes [49]. Later, this spectrum was expanded by the addition of nephropathy [53,54,55], sensorineural deafness, and subcortical and cerebellar atrophy, features commonly seen in MELAS [56]. Usually, the heteroplasmic load of the mutant variant shows a good correlation with disease severity and the spectrum of clinical phenotypes, ranging from fatal neonatal forms to asymptomatic cases [57]. Although most pathogenic mtDNA variants must be present at high heteroplasmy levels (>90%) to cause disease, mutated m.13513G>A loads as low as 30–50% can cause complex I deficiency [58]. Moreover, even lower levels of heteroplasmy (~20–30%) have been reported in patients with isolated LHON [59,60,61].
The present study analyzed the effects of the m.13513G>A variant on the OXPHOS system and alterations in cell transcriptional profiles by comparing transmitochondrial cybrid cell lines harboring 50% and 70% of mutated DNA with WT cells. Increased levels of heteroplasmy were found to result in reduced complex I functionality. As expected, complex I activity was reduced in both mutant cell lines, as the MT-ND5 gene encodes a complex I subunit essential for complex I activity [62]. Furthermore, both mutant cell lines showed reduced complex III activity, and cells with 70% heteroplasmy showed decreases in complex IV and complex I + III activities. Complexes I, III, and IV form respiratory supercomplexes [63]. Although complex I dysfunction does not usually lead to disturbed complex III and complex IV functionality, mutations in nuclear gene NDUFS4, which encodes a structural subunit of complex I, were found to cause combined complex I and III deficiencies [64,65,66]. We suggest that the m.13513 G>A mutation in the MT-ND5 gene may also cause secondary complex III and complex IV deficiencies due to impaired supercomplex formation.
High-resolution respirometry confirmed the complex I defect in both mutant cell lines, with both showing reductions in complex I-dependent oxygen consumption, and cells with 70% heteroplasmy showing little response to rotenone addition. Our obtained data were consistent with previously published results of high-resolution respirometry measurements in fibroblasts harboring the MT-ND5 m.13513G>A variant, where complex I contribution to respiration was reduced, and the total mitochondrial mass was diminished, leading to a defect in complexes I, III + V, and IV [67]. Furthermore, induced pluripotent stem cells (iPSC) generated from these fibroblasts had reduced basal respiration and a prominent defect in the activity of complexes I and III. Further confirming this, the induced neuron cells generated from these iPSCs also showed complex I deficiency [67].
Moreover, cells with 70% heteroplasmy demonstrated a significantly elevated input of succinate, suggesting a compensatory mechanism of the mitochondria for the production of ATP, as both complex I and complex II are entry points for electrons to the electron transport chain. Complex I deficiencies can be bypassed by supplying the cells with substrates for complex II, such as cell-permeable succinate prodrugs [68]. Interestingly, in contrast to cells with 70% heteroplasmy, cells with 50% heteroplasmy showed sensitivity to rotenone in a similar manner to WT cells, and did not have increased succinate input. This would indicate that, despite diminished complex I-dependent respiration, cells harboring 50% of the m.13513G>A variant are still able to produce enough energy through complex I-dependent oxygen consumption. The considerable difference in effect after succinate input seen in both heteroplasmic cell lines may be explained by the well-known fact that symptoms of mitochondrial diseases, and therefore also the biochemical manifestations, become more severe as the heteroplasmy level grows [17,26,69]. We hypothesize there might be a threshold between 50% and 70% heteroplasmy for the m.13513G>A variant that affects the ability to produce energy through complex I-dependent oxygen consumption and the need for compensation through complex II-dependent respiration. This would be supported by the fact that the affected patient, described by us in detail previously, carries 60% heteroplasmy for the m.13513G>A variant [28].
Little is known about changes in gene expression levels, associated with Leigh syndrome, and in particular with the m.13513G>A variant in the MT-ND5 gene. Wahedi et al. performed transcriptome analysis of 113 genes, currently associated with Leigh syndrome, and showed neuronal specificity of target genes, with significant enrichment in hippocampal and somatosensory pyramidal neurons as well as interneurons of the brain [70]. A previously described, transcriptomic profile analysis of fibroblasts harboring pathogenic variants in the LRPPRC gene, associated with Leigh syndrome, French-Canadian type (MIM #220111), obtained a set of DEGs and identified NDUFA4L2 gene as a target for further functional studies [71].
The RNA sequencing analysis of cells harboring 50% versus 70% of the m.13513G>A variant identified 1247 nuclear genes whose transcription levels were significantly affected, with 437 of them having expression levels altered more than four-fold.
The patient, who participated in our research as a donor of altered mtDNA, which was used for creation of cybrid cell lines, had a heteroplasmy level of 60% in the blood and exhibited features of Leigh syndrome. Confirming this, our studies on OXPHOS system functionality revealed that cells with both 50% and 70% heteroplasmy levels showed reduced complex I and III activities as well as complex I related defect in oxygen consumption. Additionally, it was previously described that, cells harboring 20–30% of the m.3243A>G variant were found to have mutually similar transcriptional profiles [26]. Furthermore, cells harboring 50–90% of the same variant also had mutually similar transcriptional profiles, although these were different from those in cells harboring 20–30% of this variant. This suggests that gradual increases in mtDNA heteroplasmy levels are associated with phased alterations in nDNA and mtDNA gene expression profiles [26]. Taking into consideration all the above-mentioned arguments, we suggest that more weight is to be given to the results of combined RNA sequencing analysis from cells with 50% and 70% heteroplasmy, compared with WT cells. Transcriptome analysis revealed that, compared with WT cells, the levels of expression of 84 nuclear genes were significantly altered in cells bearing the m.13513G>A variant. GO enrichment analysis and clustering of selected DEGs showed the involvement of genes related to calcium transport and release activity. This was not unexpected, as mitochondrial dysfunction was shown to alter levels of intracellular calcium [72]. Cybrids harboring 50%–90% of the 3243A>G variant showed increased expression of glycolytic genes, indicating bioenergetic adaptation by these cells [26]. In contrast, the expression levels of genes involved in glycolysis were not significantly altered in cells bearing the m.13513G>A variant (data not shown).
Selfie-digital PCR analysis of five genes (TCEA3, SEL1L3, ABCG1, FABP3, and GRM2) showed a trend in expression consistent with the results of transcriptome analysis. The expression of the TCEA3 gene, which encodes transcription elongation factor, A3 which releases RNA polymerase II from transcriptional arrest [73], is tissue-specific and highly expressed in skeletal muscle, where it promotes the transcription elongation of muscle specific genes by binding to myogenic regulatory transcription factors (MRFs) MyoD and myogenin. It is a required co-factor for MRF-driven gene expression during myogenesis [74]. In proliferating cells, TCEA3 is expressed at low levels, however, upon differentiation, TCEA3 is upregulated [74]. Because exercise intolerance and myopathy are often associated with the m.13513G>A variant, we can speculate that upregulation of TCEA3 expression in cells harboring the variant may serve as a compensatory mechanism for defects in muscle functions. The SEL1L3 gene, which encodes protein sel-1 homolog 3, suppressor of lin-12-like protein 3, is induced in response to mitochondrial damage [75]. SEL1L3 was identified as one of ten common mitochondrial defect (CMD) genes, and was shown to contribute to CMD-mediated regulatory mechanism, resulting in mitochondrial defect and subsequent retrograde signals for transcriptional reprogramming during hepatocellular carcinoma progression [75]. The ABCG1 gene, which encodes ATP binding cassette subfamily G member 1, plays a prominent role in apolipoprotein mediated cholesterol transport and suppresses β-amyloid (Aβ) accumulation. Dysregulated Abcg1 expression is associated with intracellular cholesterol accumulation in the brains of the APP/PS1 mouse model [76]. The upregulation of the ABCG1 gene in cells harboring the m.13513G>A variant could be considered a compensatory mechanism, possibly unspecific, for Aβ induced early mitochondria-associated ER membrane dysregulation. The FABP3 gene, which encodes fatty acid binding protein 3, is critical for the uptake of alpha-synuclein and its aggregation into intracellular filamentous-shaped inclusion bodies in a mouse cell model, as well as being involved in MPP+-induced neuronal retraction, mitochondrial activity, and reactive oxygen species (ROS) formation [77]. Fabp3 knockout mice showed decreased alpha-synuclein oligomerization and neuronal degeneration of tyrosine hydroxylase positive neurons in vivo. Published data indicate that FABP3 is critical for preventing synucleinopathies; therefore, its downregulation in cells harboring the m.13513G>A variant could lead to impaired alpha-synuclein uptake and formation of aggregates [77]. This in turn could disturb the normal functioning of neurons and could be associated with neuronal dysfunction observed in Leigh syndrome. The GRM2 gene, which encodes glutamate metabotropic receptor 2, plays an important role in normal brain function, and is a potential target for the development of drugs for neurological and psychiatric diseases [78]. In a recent study, it was shown that GRM2 is negatively regulated by miR-342-p, one of miRNA expressed in the epileptic hippocampus. Upon pharmacological induction of tonic–clonic seizures in rats, the downregulation of miR-342-p and upregulation of its target gene GRM2 were observed [79]. These results suggested a possible involvement of GMR2 in epileptogenesis. Seizures are one of the clinical manifestations associated with the m.13513G>A variant. Therefore, we propose that upregulation of GRM2 may play a role in their development.
In summary, mutant cell lines harboring the m.13513G>A variant of the MT-ND5 gene showed defects in the OXPHOS system, which had a more pronounced effect in cells bearing 70% than 50% heteroplasmy. These findings are in agreement with results showing that higher levels of heteroplasmy cause more severe symptoms. The pathogenic mtDNA variant had pronounced effects on nuclear gene expression in cybrid cell lines, with DEGs in cells harboring the m.13513G>A mtDNA variant generally reflecting the physiological processes associated with damaged mitochondria.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cimb45030115/s1. Figure S1: The dynamics in transcription level of selected genes evaluated by Selfie-digital PCR in cybrid cells bearing m.13513G>A mtDNA variant. Graph shows representative experiment. Table S1: The primer used for Selfie-digital PCR. Table S2: List of differentially expressed genes. Selected genes that correspond to chosen criteria (padj < 0.05 and absolute logFold change > 2) are shown in red. Table S3: Ontology terms determined for the selected 31 DEGs differing significantly in cybrid cell lines harboring 50% and 70% of the m.13513G>A variant in the MT-ND5 gene compared with WT cells, ranked by statistical significance (p < 0.05 and absolute logFold change > 2). Table S4: Gene Ontology terms determined for the selected 437 DEGs that differed significantly between cybrid cell lines harboring 50% and 70% of the m.13513G>A variant in MT-ND5 gene, ranked by statistical significance (p < 0.05 and absolute logFold change > 2).

Author Contributions

Conceptualization, D.K. and I.I.; methodology, D.K., M.M.-K. and I.I.; software, P.Z.; validation, D.K. and I.I.; formal analysis, P.Z.; investigation, D.K., D.L., M.M.-K., I.I. and J.S.; resources, M.M.-K., B.L. and I.I.; data curation, P.Z.; writing—original draft preparation, D.K., J.S. and I.I.; writing—review & editing, D.K., M.M.-K., I.I. and J.S.; visualization, D.K., P.Z., M.M.-K. and I.I.; supervision, T.-K.L., C.-W.L. and I.I.; project administration, B.L., T.-K.L., C.-W.L. and I.I.; funding acquisition, B.L., C.-W.L. and I.I. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Latvian Council of Science, Riga, Latvia Taiwan–Latvian–Lithuanian Collaboration Project “Functional model for the mitochondrial disease evaluation and biomarker development” and the Latvian Council of Science, Riga, Latvia [grant No: lzp-2018/1-0180].

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved the Central Committee of Medical Ethics of Latvia (protocol No.2016-5, chapter 2, approved 24 November 2016, and protocol No.2020-6 II, chapter-3, approved 22 May 2020).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study or their legal guardians.

Data Availability Statement

Publicly available datasets were analyzed in this study. These data can be found at [ArrayExpress/ E-MTAB-11035].

Acknowledgments

We express our deepest gratitude to all the patients, their families, and healthy volunteers. We acknowledge Petar Podlesny (BioRad Laboratories, Inc.) for help with primer design for Selfie-digital PCR analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. John, P.; Whatley, F.R. Paracoccus Denitrificans and the Evolutionary Origin of the Mitochondrion. Nature 1975, 254, 495–498. [Google Scholar] [CrossRef]
  2. Deluca, H.F.; Engstrom, G.W. Calcium Uptake by Rat Kidney Mitochondria. Proc. Natl. Acad. Sci. USA 1961, 47, 1744–1750. [Google Scholar] [CrossRef] [Green Version]
  3. Sancak, Y.; Markhard, A.L.; Kitami, T.; Kovács-Bogdán, E.; Kamer, K.J.; Udeshi, N.D.; Carr, S.A.; Chaudhuri, D.; Clapham, D.E.; Li, A.A.; et al. EMRE is an Essential Component of the Mitochondrial Calcium Uniporter Complex. Science 2013, 342, 1379–1382. [Google Scholar] [CrossRef] [Green Version]
  4. McBride, H.M.; Neuspiel, M.; Wasiak, S. Mitochondria: More than Just a Powerhouse. Curr. Biol. CB 2006, 16, R551–R560. [Google Scholar] [CrossRef] [Green Version]
  5. Lax, N.Z.; Turnbull, D.M.; Reeve, A.K. Mitochondrial Mutations: Newly Discovered Players in Neuronal Degeneration. Neuroscientist 2011, 17, 645–658. [Google Scholar] [CrossRef]
  6. Koopman, W.J.H.; Willems, P.H.G.M.; Smeitink, J.A.M. Monogenic Mitochondrial Disorders. N. Engl. J. Med. 2012, 366, 1132–1141. [Google Scholar] [CrossRef] [Green Version]
  7. Kroemer, G.; Galluzzi, L.; Brenner, C. Mitochondrial Membrane Permeabilization in Cell Death. Physiol. Rev. 2007, 87, 99–163. [Google Scholar] [CrossRef]
  8. Hockenbery, D.; Nuñez, G.; Milliman, C.; Schreiber, R.D.; Korsmeyer, S.J. Bcl-2 is an Inner Mitochondrial Membrane Protein That Blocks Programmed Cell Death. Nature 1990, 348, 334–336. [Google Scholar] [CrossRef]
  9. Edinger, A.L.; Thompson, C.B. Death by Design: Apoptosis, Necrosis and Autophagy. Curr. Opin. Cell Biol. 2004, 16, 663–669. [Google Scholar] [CrossRef]
  10. Butow, R.A.; Avadhani, N.G. Mitochondrial Signaling: The Retrograde Response. Mol. Cell 2004, 14, 1–15. [Google Scholar] [CrossRef]
  11. da Cunha, F.M.; Torelli, N.Q.; Kowaltowski, A.J. Mitochondrial Retrograde Signaling: Triggers, Pathways, and Outcomes. Oxid. Med. Cell. Longev. 2015, 2015, 482582. [Google Scholar] [CrossRef] [Green Version]
  12. Walker, B.R.; Moraes, C.T. Nuclear-Mitochondrial Interactions. Biomolecules 2022, 12, 427. [Google Scholar] [CrossRef]
  13. Srinivasan, V.; Kriete, A.; Sacan, A.; Jazwinski, S.M. Comparing the Yeast Retrograde Response and NF-ΚB Stress Responses: Implications for Aging. Aging Cell 2010, 9, 933–941. [Google Scholar] [CrossRef] [Green Version]
  14. Jazwinski, S.M.; Kriete, A. The Yeast Retrograde Response as a Model of Intracellular Signaling of Mitochondrial Dysfunction. Front. Physiol. 2012, 3, 139. [Google Scholar] [CrossRef] [Green Version]
  15. Cagin, U.; Enriquez, J.A. The Complex Crosstalk between Mitochondria and the Nucleus: What Goes in Between? Int. J. Biochem. Cell Biol. 2015, 63, 10–15. [Google Scholar] [CrossRef]
  16. Biswas, G.; Adebanjo, O.A.; Freedman, B.D.; Anandatheerthavarada, H.K.; Vijayasarathy, C.; Zaidi, M.; Kotlikoff, M.; Avadhani, N.G. Retrograde Ca2+ Signaling in C2C12 Skeletal Myocytes in Response to Mitochondrial Genetic and Metabolic Stress: A Novel Mode of Inter-Organelle Crosstalk. EMBO J. 1999, 18, 522–533. [Google Scholar] [CrossRef] [Green Version]
  17. Tuppen, H.A.L.; Blakely, E.L.; Turnbull, D.M.; Taylor, R.W. Mitochondrial DNA Mutations and Human Disease. Biochim. Biophys. Acta BBA Bioenerg. 2010, 1797, 113–128. [Google Scholar] [CrossRef] [Green Version]
  18. Guha, M.; Srinivasan, S.; Guja, K.; Mejia, E.; Garcia-Diaz, M.; Johnson, F.B.; Ruthel, G.; Kaufman, B.A.; Rappaport, E.F.; Glineburg, M.R.; et al. HnRNPA2 is a Novel Histone Acetyltransferase That Mediates Mitochondrial Stress-Induced Nuclear Gene Expression. Cell Discov. 2016, 2, 16045. [Google Scholar] [CrossRef] [Green Version]
  19. Rooijers, K.; Loayza-Puch, F.; Nijtmans, L.G.; Agami, R. Ribosome Profiling Reveals Features of Normal and Disease-Associated Mitochondrial Translation. Nat. Commun. 2013, 4, 2886. [Google Scholar] [CrossRef] [Green Version]
  20. Chae, S.; Ahn, B.Y.; Byun, K.; Cho, Y.M.; Yu, M.-H.; Lee, B.; Hwang, D.; Park, K.S. A Systems Approach for Decoding Mitochondrial Retrograde Signaling Pathways. Sci. Signal. 2013, 6, rs4. [Google Scholar] [CrossRef]
  21. Kandel, J.; Picard, M.; Wallace, D.C.; Eckmann, D.M. Mitochondrial DNA 3243A>G Heteroplasmy is Associated with Changes in Cytoskeletal Protein Expression and Cell Mechanics. J. R. Soc. Interface 2017, 14, 20170071. [Google Scholar] [CrossRef] [Green Version]
  22. Pignataro, D.; Francia, S.; Zanetta, F.; Brenna, G.; Brandini, S.; Olivieri, A.; Torroni, A.; Biamonti, G.; Montecucco, A. A Missense MT-ND5 Mutation in Differentiated Parkinson Disease Cytoplasmic Hybrid Induces ROS-Dependent DNA Damage Response Amplified by DROSHA. Sci. Rep. 2017, 7, 9528. [Google Scholar] [CrossRef] [Green Version]
  23. Meseguer, S.; Martínez-Zamora, A.; García-Arumí, E.; Andreu, A.L.; Armengod, M.-E. The ROS-Sensitive MicroRNA-9/9* Controls the Expression of Mitochondrial TRNA-Modifying Enzymes and is Involved in the Molecular Mechanism of MELAS Syndrome. Hum. Mol. Genet. 2015, 24, 167–184. [Google Scholar] [CrossRef] [Green Version]
  24. Chomyn, A. Platelet-Mediated Transformation of Human Mitochondrial DNA-Less Cells. Methods Enzymol. 1996, 264, 334–339. [Google Scholar] [CrossRef]
  25. Lin, T.-K.; Lin, H.-Y.; Chen, S.-D.; Chuang, Y.-C.; Chuang, J.-H.; Wang, P.-W.; Huang, S.-T.; Tiao, M.-M.; Chen, J.-B.; Liou, C.-W. The Creation of Cybrids Harboring Mitochondrial Haplogroups in the Taiwanese Population of Ethnic Chinese Background: An Extensive in Vitro Tool for the Study of Mitochondrial Genomic Variations. Oxid. Med. Cell. Longev. 2012, 2012, 824275. [Google Scholar] [CrossRef] [Green Version]
  26. Picard, M.; Zhang, J.; Hancock, S.; Derbeneva, O.; Golhar, R.; Golik, P.; O’Hearn, S.; Levy, S.; Potluri, P.; Lvova, M.; et al. Progressive Increase in MtDNA 3243A>G Heteroplasmy Causes Abrupt Transcriptional Reprogramming. Proc. Natl. Acad. Sci. USA 2014, 111, E4033–E4042. [Google Scholar] [CrossRef] [Green Version]
  27. Rodenburg, R.J. Mitochondrial Complex I-Linked Disease. Biochim. Biophys. Acta 2016, 1857, 938–945. [Google Scholar] [CrossRef]
  28. Pelnena, D.; Burnyte, B.; Jankevics, E.; Lace, B.; Dagyte, E.; Grigalioniene, K.; Utkus, A.; Krumina, Z.; Rozentale, J.; Adomaitiene, I.; et al. Complete MtDNA Sequencing Reveals Mutations m.9185T>C and m.13513G>A in Three Patients with Leigh Syndrome. Mitochondrial DNA Part A DNA Mapp. Seq. Anal. 2018, 29, 1115–1120. [Google Scholar] [CrossRef]
  29. King, M.P.; Attardi, G. Human Cells Lacking MtDNA: Repopulation with Exogenous Mitochondria by Complementation. Science 1989, 246, 500–503. [Google Scholar] [CrossRef]
  30. Chang, Y.-C.; Liao, C.-B.; Hsieh, P.-Y.C.; Liou, M.-L.; Liu, Y.-C. Expression of Tumor Suppressor P53 Facilitates DNA Repair but Not UV-Induced G2/M Arrest or Apoptosis in Chinese Hamster Ovary CHO-K1 Cells. J. Cell. Biochem. 2008, 103, 528–537. [Google Scholar] [CrossRef]
  31. Mayr, J.A.; Havlícková, V.; Zimmermann, F.; Magler, I.; Kaplanová, V.; Jesina, P.; Pecinová, A.; Nusková, H.; Koch, J.; Sperl, W.; et al. Mitochondrial ATP Synthase Deficiency Due to a Mutation in the ATP5E Gene for the F1 Epsilon Subunit. Hum. Mol. Genet. 2010, 19, 3430–3439. [Google Scholar] [CrossRef] [Green Version]
  32. Feichtinger, R.G.; Zimmermann, F.; Mayr, J.A.; Neureiter, D.; Hauser-Kronberger, C.; Schilling, F.H.; Jones, N.; Sperl, W.; Kofler, B. Low Aerobic Mitochondrial Energy Metabolism in Poorly- or Undifferentiated Neuroblastoma. BMC Cancer 2010, 10, 149. [Google Scholar] [CrossRef] [Green Version]
  33. Aminzadeh-Gohari, S.; Weber, D.D.; Catalano, L.; Feichtinger, R.G.; Kofler, B.; Lang, R. Targeting Mitochondria in Melanoma. Biomolecules 2020, 10, 1395. [Google Scholar] [CrossRef]
  34. Rustin, P.; Chretien, D.; Bourgeron, T.; Gérard, B.; Rötig, A.; Saudubray, J.M.; Munnich, A. Biochemical and Molecular Investigations in Respiratory Chain Deficiencies. Clin. Chim. Acta Int. J. Clin. Chem. 1994, 228, 35–51. [Google Scholar] [CrossRef]
  35. Berger, A.; Mayr, J.A.; Meierhofer, D.; Fötschl, U.; Bittner, R.; Budka, H.; Grethen, C.; Huemer, M.; Kofler, B.; Sperl, W. Severe Depletion of Mitochondrial DNA in Spinal Muscular Atrophy. Acta Neuropathol. 2003, 105, 245–251. [Google Scholar] [CrossRef] [Green Version]
  36. Martin, M. Cutadapt Removes Adapter Sequences from High-Throughput Sequencing Reads. EMBnet. J. 2011, 17, 10–12. [Google Scholar] [CrossRef]
  37. Kopylova, E.; Noé, L.; Touzet, H. SortMeRNA: Fast and Accurate Filtering of Ribosomal RNAs in Metatranscriptomic Data. Bioinformatics 2012, 28, 3211–3217. [Google Scholar] [CrossRef] [Green Version]
  38. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast Universal RNA-Seq Aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef]
  39. Anders, S.; Pyl, P.T.; Huber, W. HTSeq—A Python Framework to Work with High-Throughput Sequencing Data. Bioinformatics 2015, 31, 166–169. [Google Scholar] [CrossRef] [Green Version]
  40. Blighe, K. Haplotype Classification Using Copy Number Variation and Principal Components Analysis. Open Bioinforma. J. 2013, 7, 19–24. [Google Scholar] [CrossRef] [Green Version]
  41. Love, M.I.; Huber, W.; Anders, S. Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [Green Version]
  42. Blighe, K. EnhancedVolcano: Publication-Ready Volcano Plots with Enhanced Colouring and Labeling 2023. Available online: https://github.com/kevinblighe/EnhancedVolcano (accessed on 15 January 2023).
  43. Falcon, S.; Gentleman, R. Using GOstats to Test Gene Lists for GO Term Association. Bioinformatics 2007, 23, 257–258. [Google Scholar] [CrossRef] [Green Version]
  44. Winter, D.J. Rentrez: An R Package for the NCBI EUtils API. R J. 2017, 9, 520–526. [Google Scholar] [CrossRef]
  45. GO.Db. Available online: http://bioconductor.org/packages/GO.db/ (accessed on 20 October 2022).
  46. Ge, S.X.; Jung, D.; Yao, R. ShinyGO: A Graphical Gene-Set Enrichment Tool for Animals and Plants. Bioinformatics 2020, 36, 2628–2629. [Google Scholar] [CrossRef]
  47. Podlesniy, P.; Trullas, R. Absolute Measurement of Gene Transcripts with Selfie-Digital PCR. Sci. Rep. 2017, 7, 8328. [Google Scholar] [CrossRef] [Green Version]
  48. Larsen, S.; Nielsen, J.; Hansen, C.N.; Nielsen, L.B.; Wibrand, F.; Stride, N.; Schroder, H.D.; Boushel, R.; Helge, J.W.; Dela, F.; et al. Biomarkers of Mitochondrial Content in Skeletal Muscle of Healthy Young Human Subjects. J. Physiol. 2012, 590, 3349–3360. [Google Scholar] [CrossRef]
  49. Santorelli, F.M.; Tanji, K.; Kulikova, R.; Shanske, S.; Vilarinho, L.; Hays, A.P.; DiMauro, S. Identification of a Novel Mutation in the MtDNA ND5 Gene Associated with MELAS. Biochem. Biophys. Res. Commun. 1997, 238, 326–328. [Google Scholar] [CrossRef]
  50. Danhelovska, T.; Kolarova, H.; Zeman, J.; Hansikova, H.; Vaneckova, M.; Lambert, L.; Kucerova-Vidrova, V.; Berankova, K.; Honzik, T.; Tesarova, M. Multisystem Mitochondrial Diseases Due to Mutations in MtDNA-Encoded Subunits of Complex I. BMC Pediatr. 2020, 20, 41. [Google Scholar] [CrossRef]
  51. Shanske, S.; Coku, J.; Lu, J.; Ganesh, J.; Krishna, S.; Tanji, K.; Bonilla, E.; Naini, A.B.; Hirano, M.; DiMauro, S. The G13513A Mutation in the ND5 Gene of Mitochondrial DNA as a Common Cause of MELAS or Leigh Syndrome: Evidence from 12 Cases. Arch. Neurol. 2008, 65, 368–372. [Google Scholar] [CrossRef] [Green Version]
  52. Wang, S.-B.; Weng, W.-C.; Lee, N.-C.; Hwu, W.-L.; Fan, P.-C.; Lee, W.-T. Mutation of Mitochondrial DNA G13513A Presenting with Leigh Syndrome, Wolff-Parkinson-White Syndrome and Cardiomyopathy. Pediatr. Neonatol. 2008, 49, 145–149. [Google Scholar] [CrossRef] [Green Version]
  53. Ng, P.S.; Pinto, M.V.; Neff, J.L.; Hasadsri, L.; Highsmith, E.W.; Fidler, M.E.; Gavrilova, R.H.; Klein, C.J. Mitochondrial Cerebellar Ataxia, Renal Failure, Neuropathy, and Encephalopathy (MCARNE). Neurol. Genet. 2019, 5, e314. [Google Scholar] [CrossRef] [Green Version]
  54. Bakis, H.; Trimouille, A.; Vermorel, A.; Redonnet, I.; Goizet, C.; Boulestreau, R.; Lacombe, D.; Combe, C.; Martin-Négrier, M.-L.; Rigothier, C. Adult Onset Tubulo-Interstitial Nephropathy in MT-ND5-Related Phenotypes. Clin. Genet. 2020, 97, 628–633. [Google Scholar] [CrossRef]
  55. Barone, V.; La Morgia, C.; Caporali, L.; Fiorini, C.; Carbonelli, M.; Gramegna, L.L.; Bartiromo, F.; Tonon, C.; Morandi, L.; Liguori, R.; et al. Case Report: Optic Atrophy and Nephropathy With m.13513G>A/MT-ND5 MtDNA Pathogenic Variant. Front. Genet. 2022, 13, 887696. [Google Scholar] [CrossRef]
  56. Gramegna, L.L.; Evangelisti, S.; Di Vito, L.; La Morgia, C.; Maresca, A.; Caporali, L.; Amore, G.; Talozzi, L.; Bianchini, C.; Testa, C.; et al. Brain MRS Correlates with Mitochondrial Dysfunction Biomarkers in MELAS-Associated MtDNA Mutations. Ann. Clin. Transl. Neurol. 2021, 8, 1200–1211. [Google Scholar] [CrossRef]
  57. Ng, Y.S.; Lax, N.Z.; Maddison, P.; Alston, C.L.; Blakely, E.L.; Hepplewhite, P.D.; Riordan, G.; Meldau, S.; Chinnery, P.F.; Pierre, G.; et al. MT-ND5 Mutation Exhibits Highly Variable Neurological Manifestations at Low Mutant Load. EBioMedicine 2018, 30, 86–93. [Google Scholar] [CrossRef] [Green Version]
  58. Kirby, D.M.; Boneh, A.; Chow, C.W.; Ohtake, A.; Ryan, M.T.; Thyagarajan, D.; Thorburn, D.R. Low Mutant Load of Mitochondrial DNA G13513A Mutation Can Cause Leigh’s Disease. Ann. Neurol. 2003, 54, 473–478. [Google Scholar] [CrossRef]
  59. Krylova, T.D.; Sheremet, N.L.; Tabakov, V.Y.; Lyamzaev, K.G.; Itkis, Y.S.; Tsygankova, P.G.; Andreeva, N.A.; Shmelkova, M.S.; Nevinitsyna, T.A.; Kadyshev, V.V.; et al. Three Rare Pathogenic MtDNA Substitutions in LHON Patients with Low Heteroplasmy. Mitochondrion 2020, 50, 139–144. [Google Scholar] [CrossRef]
  60. Sun, C.-B.; Bai, H.-X.; Xu, D.-N.; Xiao, Q.; Liu, Z. Mitochondrial 13513G>A Mutation with Low Mutant Load Presenting as Isolated Leber’s Hereditary Optic Neuropathy Assessed by Next Generation Sequencing. Front. Neurol. 2021, 12, 601307. [Google Scholar] [CrossRef]
  61. Vázquez-Justes, D.; Carreño-Gago, L.; García-Arumi, E.; Traveset, A.; Montoya, J.; Ruiz-Pesini, E.; López, R.; Brieva, L. Mitochondrial m.13513G>A Point Mutation in ND5 in a 16-Year-Old Man with Leber Hereditary Optic Neuropathy Detected by Next-Generation Sequencing. J. Pediatr. Genet. 2019, 8, 231–234. [Google Scholar] [CrossRef]
  62. Chomyn, A. Mitochondrial Genetic Control of Assembly and Function of Complex I in Mammalian Cells. J. Bioenerg. Biomembr. 2001, 33, 251–257. [Google Scholar] [CrossRef]
  63. Zhu, J.; Vinothkumar, K.R.; Hirst, J. Structure of Mammalian Respiratory Complex I. Nature 2016, 536, 354–358. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Moreno-Lastres, D.; Fontanesi, F.; García-Consuegra, I.; Martín, M.A.; Arenas, J.; Barrientos, A.; Ugalde, C. Mitochondrial Complex I Plays an Essential Role in Human Respirasome Assembly. Cell Metab. 2012, 15, 324–335. [Google Scholar] [CrossRef] [Green Version]
  65. Ugalde, C.; Janssen, R.J.R.J.; van den Heuvel, L.P.; Smeitink, J.A.M.; Nijtmans, L.G.J. Differences in Assembly or Stability of Complex I and Other Mitochondrial OXPHOS Complexes in Inherited Complex I Deficiency. Hum. Mol. Genet. 2004, 13, 659–667. [Google Scholar] [CrossRef]
  66. Budde, S.M.; van den Heuvel, L.P.; Janssen, A.J.; Smeets, R.J.; Buskens, C.A.; DeMeirleir, L.; Van Coster, R.; Baethmann, M.; Voit, T.; Trijbels, J.M.; et al. Combined Enzymatic Complex I and III Deficiency Associated with Mutations in the Nuclear Encoded NDUFS4 Gene. Biochem. Biophys. Res. Commun. 2000, 275, 63–68. [Google Scholar] [CrossRef]
  67. Galera-Monge, T.; Zurita-Díaz, F.; Canals, I.; Hansen, M.G.; Rufián-Vázquez, L.; Ehinger, J.K.; Elmér, E.; Martin, M.A.; Garesse, R.; Ahlenius, H.; et al. Mitochondrial Dysfunction and Calcium Dysregulation in Leigh Syndrome Induced Pluripotent Stem Cell Derived Neurons. Int. J. Mol. Sci. 2020, 21, 3191. [Google Scholar] [CrossRef]
  68. Ehinger, J.K.; Piel, S.; Ford, R.; Karlsson, M.; Sjövall, F.; Frostner, E.Å.; Morota, S.; Taylor, R.W.; Turnbull, D.M.; Cornell, C.; et al. Cell-Permeable Succinate Prodrugs Bypass Mitochondrial Complex I Deficiency. Nat. Commun. 2016, 7, 12317. [Google Scholar] [CrossRef] [Green Version]
  69. Stewart, J.B.; Chinnery, P.F. The Dynamics of Mitochondrial DNA Heteroplasmy: Implications for Human Health and Disease. Nat. Rev. Genet. 2015, 16, 530–542. [Google Scholar] [CrossRef]
  70. Wahedi, A.; Soondram, C.; Murphy, A.E.; Skene, N.; Rahman, S. Transcriptomic Analyses Reveal Neuronal Specificity of Leigh Syndrome Associated Genes. J. Inherit. Metab. Dis. 2023, 1–18. [Google Scholar] [CrossRef]
  71. Bchetnia, M.; Tardif, J.; Morin, C.; Laprise, C. Expression Signature of the Leigh Syndrome French-Canadian Type. Mol. Genet. Metab. Rep. 2022, 30, 100847. [Google Scholar] [CrossRef]
  72. Karakaidos, P.; Rampias, T. Mitonuclear Interactions in the Maintenance of Mitochondrial Integrity. Life 2020, 10, 173. [Google Scholar] [CrossRef]
  73. Labhart, P.; Morgan, G.T. Identification of Novel Genes Encoding Transcription Elongation Factor TFIIS (TCEA) in Vertebrates: Conservation of Three Distinct TFIIS Isoforms in Frog, Mouse, and Human. Genomics 1998, 52, 278–288. [Google Scholar] [CrossRef]
  74. Kazim, N.; Adhikari, A.; Davie, J. The Transcription Elongation Factor TCEA3 Promotes the Activity of the Myogenic Regulatory Factors. PLoS ONE 2019, 14, e0217680. [Google Scholar] [CrossRef] [PubMed]
  75. Lee, Y.-K.; Jee, B.A.; Kwon, S.M.; Yoon, Y.-S.; Xu, W.G.; Wang, H.-J.; Wang, X.W.; Thorgeirsson, S.S.; Lee, J.-S.; Woo, H.G.; et al. Identification of a Mitochondrial Defect Gene Signature Reveals NUPR1 as a Key Regulator of Liver Cancer Progression. Hepatology 2015, 62, 1174–1189. [Google Scholar] [CrossRef]
  76. Völgyi, K.; Badics, K.; Sialana, F.J.; Gulyássy, P.; Udvari, E.B.; Kis, V.; Drahos, L.; Lubec, G.; Kékesi, K.A.; Juhász, G. Early Presymptomatic Changes in the Proteome of Mitochondria-Associated Membrane in the APP/PS1 Mouse Model of Alzheimer’s Disease. Mol. Neurobiol. 2018, 55, 7839–7857. [Google Scholar] [CrossRef]
  77. Kawahata, I.; Bousset, L.; Melki, R.; Fukunaga, K. Fatty Acid-Binding Protein 3 Is Critical for α-Synuclein Uptake and MPP+-Induced Mitochondrial Dysfunction in Cultured Dopaminergic Neurons. Int. J. Mol. Sci. 2019, 20, 5358. [Google Scholar] [CrossRef] [Green Version]
  78. Wei, L.; He, F.; Zhang, W.; Chen, W.; Yu, B. Bioinformatics Analysis of Microarray Data to Reveal the Pathogenesis of Diffuse Intrinsic Pontine Glioma. Biol. Res. 2018, 51, 26. [Google Scholar] [CrossRef] [Green Version]
  79. Pala, M.; Meral, I.; Pala Acikgoz, N.; Gorucu Yilmaz, S.; Okur, S.K.; Acar, S.; Polat, Y.; Akbas, F. Downregulatory Effect of MiR-342-3p on Epileptogenesis in the PTZ-Kindling Model. Mol. Biol. Rep. 2022, 49, 11997–12006. [Google Scholar] [CrossRef]
Cimb 45 00115 g001 550
Figure 1. Ratios of OXPHOS system complex activities to citrate synthase (CS) activity in cell lines harboring 50% and 70% of the m.13513G>A variant compared with wild-type (WT) cells. (A) Relative activities of the OXPHOS system complexes I, II, IV, I + III, and II + III. (B) Relative activity of the OXPHOS system complex III. Data are shown as the mean ± standard deviation. * p < 0.05 compared with WT by two-sample t-tests.
Figure 1. Ratios of OXPHOS system complex activities to citrate synthase (CS) activity in cell lines harboring 50% and 70% of the m.13513G>A variant compared with wild-type (WT) cells. (A) Relative activities of the OXPHOS system complexes I, II, IV, I + III, and II + III. (B) Relative activity of the OXPHOS system complex III. Data are shown as the mean ± standard deviation. * p < 0.05 compared with WT by two-sample t-tests.
Cimb 45 00115 g001
Cimb 45 00115 g002 550
Figure 2. Comparative respiration rates of cybrid cells harboring 50% and 70% of the m.13513G>A variant and of wild-type (WT) cells. Reduced complex I-dependent oxygen consumption is seen in cells harboring the m.13513G>A variant, with more pronounced effect in cells with 70% heteroplasmy. Data are shown as the mean ± standard deviation. * p < 0.05 compared with WT cells; # p < 0.05 compared with cells having 70% heteroplasmy by two-sample t-tests. Abbreviations: PM, pyruvate + malate; ADP, adenosine diphosphate; G, glutamate; S, succinate; Rot, rotenone; AsTm, ascorbate + TMPD; CI, complex I; CII, complex II; CIV, complex IV; N, NADH; ROUTINE, physiological coupling state; LEAK, non-phosphorylating substrate-dependent state; OXPHOS, oxidative phosphorylation-dependent state.
Figure 2. Comparative respiration rates of cybrid cells harboring 50% and 70% of the m.13513G>A variant and of wild-type (WT) cells. Reduced complex I-dependent oxygen consumption is seen in cells harboring the m.13513G>A variant, with more pronounced effect in cells with 70% heteroplasmy. Data are shown as the mean ± standard deviation. * p < 0.05 compared with WT cells; # p < 0.05 compared with cells having 70% heteroplasmy by two-sample t-tests. Abbreviations: PM, pyruvate + malate; ADP, adenosine diphosphate; G, glutamate; S, succinate; Rot, rotenone; AsTm, ascorbate + TMPD; CI, complex I; CII, complex II; CIV, complex IV; N, NADH; ROUTINE, physiological coupling state; LEAK, non-phosphorylating substrate-dependent state; OXPHOS, oxidative phosphorylation-dependent state.
Cimb 45 00115 g002
Cimb 45 00115 g003 550
Figure 3. Effect of glutamate (A), succinate (B), and rotenone (C) on respiration rates of cell lines harboring 50% and 70% of the m.13513G>A variant and wild-type (WT) cells. Data are shown as the mean ± standard deviation. * p < 0.05 compared with WT cells; # p < 0.05 compared with cells having 70% heteroplasmy by two-sample t-tests.
Figure 3. Effect of glutamate (A), succinate (B), and rotenone (C) on respiration rates of cell lines harboring 50% and 70% of the m.13513G>A variant and wild-type (WT) cells. Data are shown as the mean ± standard deviation. * p < 0.05 compared with WT cells; # p < 0.05 compared with cells having 70% heteroplasmy by two-sample t-tests.
Cimb 45 00115 g003
Cimb 45 00115 g004 550
Figure 4. Differential gene expression in cybrid cells harboring 50% and 70% of the m.13513G>A variant and cells containing WT mtDNA. (A,B) Principal component analysis showing variability between (A) heteroplasmic cells containing m.13513G>A variant (50% and 70%) and WT mtDNA and (B) cells with 50% heteroplasmy vs. cells with 70% heteroplasmy. (C,D) Volcano plots for expressed genes (C) in the cells with 50% and 70% heteroplasmy vs. cells containing WT mtDNA and (D) in the cells with 50% heteroplasmy vs. cells with 70% heteroplasmy.
Figure 4. Differential gene expression in cybrid cells harboring 50% and 70% of the m.13513G>A variant and cells containing WT mtDNA. (A,B) Principal component analysis showing variability between (A) heteroplasmic cells containing m.13513G>A variant (50% and 70%) and WT mtDNA and (B) cells with 50% heteroplasmy vs. cells with 70% heteroplasmy. (C,D) Volcano plots for expressed genes (C) in the cells with 50% and 70% heteroplasmy vs. cells containing WT mtDNA and (D) in the cells with 50% heteroplasmy vs. cells with 70% heteroplasmy.
Cimb 45 00115 g004
Cimb 45 00115 g005 550
Figure 5. Clustered molecular function Gene Ontology terms of the differentially expressed genes (DEGs) differing significantly in cybrid cell lines harboring 50% and 70% of the m.13513G>A variant in the MT-ND5 gene compared with WT cells. DEGs were defined as those with absolute logFC > 2 and false discovery rate (FDR) < 0.05; the data were plotted using R package ShinyGO. Top 30 of Molecular Function GO terms are shown (FDR cutoff 0.1). Hierarchical clustering tree is built on similarity of gene sets listed in GO terms.
Figure 5. Clustered molecular function Gene Ontology terms of the differentially expressed genes (DEGs) differing significantly in cybrid cell lines harboring 50% and 70% of the m.13513G>A variant in the MT-ND5 gene compared with WT cells. DEGs were defined as those with absolute logFC > 2 and false discovery rate (FDR) < 0.05; the data were plotted using R package ShinyGO. Top 30 of Molecular Function GO terms are shown (FDR cutoff 0.1). Hierarchical clustering tree is built on similarity of gene sets listed in GO terms.
Cimb 45 00115 g005
Cimb 45 00115 g006 550
Figure 6. Clustered molecular functions, based on Gene Ontology terms of the differentially expressed genes (DEGs), which differed significantly between cybrid cell lines harboring 50% and 70% of the m.13513G>A variant in MT-ND5 gene. DEGs were defined as genes with absolute logFC > 2 and false discovery rate (FDR) < 0.05; the data were plotted using R package ShinyGO. Top 30 of Molecular Function GO terms are shown (FDR cutoff 0.1). Hierarchical clustering tree is built on similarity of gene sets listed in GO terms.
Figure 6. Clustered molecular functions, based on Gene Ontology terms of the differentially expressed genes (DEGs), which differed significantly between cybrid cell lines harboring 50% and 70% of the m.13513G>A variant in MT-ND5 gene. DEGs were defined as genes with absolute logFC > 2 and false discovery rate (FDR) < 0.05; the data were plotted using R package ShinyGO. Top 30 of Molecular Function GO terms are shown (FDR cutoff 0.1). Hierarchical clustering tree is built on similarity of gene sets listed in GO terms.
Cimb 45 00115 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kidere, D.; Zayakin, P.; Livcane, D.; Makrecka-Kuka, M.; Stavusis, J.; Lace, B.; Lin, T.-K.; Liou, C.-W.; Inashkina, I. Impact of the m.13513G>A Variant on the Functions of the OXPHOS System and Cell Retrograde Signaling. Curr. Issues Mol. Biol. 2023, 45, 1794-1809. https://doi.org/10.3390/cimb45030115

AMA Style

Kidere D, Zayakin P, Livcane D, Makrecka-Kuka M, Stavusis J, Lace B, Lin T-K, Liou C-W, Inashkina I. Impact of the m.13513G>A Variant on the Functions of the OXPHOS System and Cell Retrograde Signaling. Current Issues in Molecular Biology. 2023; 45(3):1794-1809. https://doi.org/10.3390/cimb45030115

Chicago/Turabian Style

Kidere, Dita, Pawel Zayakin, Diana Livcane, Marina Makrecka-Kuka, Janis Stavusis, Baiba Lace, Tsu-Kung Lin, Chia-Wei Liou, and Inna Inashkina. 2023. "Impact of the m.13513G>A Variant on the Functions of the OXPHOS System and Cell Retrograde Signaling" Current Issues in Molecular Biology 45, no. 3: 1794-1809. https://doi.org/10.3390/cimb45030115

APA Style

Kidere, D., Zayakin, P., Livcane, D., Makrecka-Kuka, M., Stavusis, J., Lace, B., Lin, T. -K., Liou, C. -W., & Inashkina, I. (2023). Impact of the m.13513G>A Variant on the Functions of the OXPHOS System and Cell Retrograde Signaling. Current Issues in Molecular Biology, 45(3), 1794-1809. https://doi.org/10.3390/cimb45030115

Article Metrics

Back to TopTop