Oxidative Medicine and Cellular Longevity

Volume 2012 (2012), Article ID 824275, 13 pages

doi:10.1155/2012/824275

Research Article

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

Tsu-Kung Lin,1,2 Hung-Yu Lin,1,2,3 Shang-Der Chen,1,2 Yao-Chung Chuang,1,2,3 Jiin-Haur Chuang,2,4 Pei-Wen Wang,2,5 Sheng-Teng Huang,2,6 Mao-Meng Tiao,2,7 Jin-Bor Chen,2,8 and Chia-Wei Liou1,2

¹Department of Neurology, Kaohsiung Chang Gung Memorial Hospital and Chang Gung University College of Medicine, Kaohsiung 833, Taiwan

²The Mitochondrial Research Unit, Kaohsiung Chang Gung Memorial Hospital and Chang Gung University College of Medicine, Kaohsiung 833, Taiwan

³Department of Biological Sciences, National Sun Yat-Sen University, Kaohsiung 804, Taiwan

⁴Division of Pediatric Surgery, Kaohsiung Chang Gung Memorial Hospital and Chang Gung University College of Medicine, Kaohsiung 833, Taiwan

⁵Department of Internal Medicine, Kaohsiung Chang Gung Memorial Hospital and Chang Gung University College of Medicine, Kaohsiung 833, Taiwan

⁶Department of Chinese Medicine, Kaohsiung Chang Gung Memorial Hospital and Chang Gung University College of Medicine, Kaohsiung 833, Taiwan

⁷Department of Pediatrics, Kaohsiung Chang Gung Memorial Hospital and Chang Gung University College of Medicine, Kaohsiung 833, Taiwan

⁸Department of Nephrology, Kaohsiung Chang Gung Memorial Hospital and Chang Gung University College of Medicine, Kaohsiung 833, Taiwan

Received 28 June 2012; Revised 15 September 2012; Accepted 17 October 2012

Academic Editor: Jon Jay Ramsey

Copyright © 2012 Tsu-Kung Lin et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Mitochondrial DNA (mtDNA) haplogroups may contribute to the development of aging-related diseases. A reliable in vitro cellular system for investigating the physiologic significance of mtDNA haplogroups is essential. This study aims to construct and characterize a series of cybrid cell lines harboring variant mtDNA haplogroups collected from healthy Taiwanese volunteers. Cybrid cells harboring different mtDNA haplogroups like B4a, B4b, B4c, B4d, B5, R, F1a, F2, D4e, D4a, D5b, D5a, E, M8, C, and N9a were prepared. Luminex 1000 and full-length mtDNA sequencing were used to confirm that mtDNA haplogroups of transmitochondrial cybrids were identical to their original donors. Cybrid B4b had a significantly lower oxygen consumption rate and higher mitochondrial membrane potential compared to F1a, B5, D5a, D4a, and N9a but had more susceptibility to H₂O₂-induced oxidative stress than cybrid F1a, D4a, and N9a. Cybrid N9a had better oxygen consumption and H₂O₂-challenged viability compared to B4b, F1a, B5, D5a, and D4a. A series of cybrid cells harboring the main haplogroups of the Taiwanese population with ethnic Chinese background has been developed in vitro. With this mtDNA haplogroup population, the underlying mechanisms of aging-related diseases may be better understood, and therapeutic interventions can be accelerated.

1. Introduction

A double-membrane structure and a circular mitochondrial genome with its own transcription, translation, and protein assembly systems support the endosymbiotic origin of the mitochondria [1]. While most mitochondrial proteins are encoded by nuclear DNA and imported into the organelle, 13 are encoded by mitochondrial DNA. Human mitochondrial DNA (mtDNA) is a 16,568 bp circular double-stranded DNA molecule coding for 13 polypeptide components of the mitochondrial oxidative phosphorylation (OXPHOS) machinery, 2 ribosomal RNA molecules, and a set of 22 transfer RNA molecules [1–4]. More than 99% of mitochondrial proteins are encoded by the nuclear genome, translated on cytoplasmic ribosomes, and selectively imported into the appropriate mitochondrial compartment [5].

Mitochondrial DNA is inherited maternally, and each mitochondrion contains 2–10 mtDNA molecules [6]. Mitochondrial DNA has a very high mutation rate and when a mutation occurs, normal and mutant mitochondrial DNA can coexist within the same cell, known as heteroplasmy [5, 7]. Mitochondrial DNA defects are linked with a wide range of diseases known as mitochondrial encephalomyopathies. Over the last two decades, more than 100 pathologic defects of mtDNA have been characterized, many due to gross structural rearrangements (i.e., single or multiple deletions or duplications) or from point mutations in the mtDNA [7–9]. Large-scale deletions of mtDNA in cancer cells imply that deficits in the mitochondrial genome can change cellular behavior similar to cross-talking with nuclear genome [10].

Phylogeographic studies on human mtDNA have revealed a remarkable correlation between mtDNA lineages and the geographic origins of indigenous populations. Regional mtDNA lineages are groups of related individual mtDNA sequences (haplotypes) known as haplogroups [11]. Unique sets of common SNPs define groups of mtDNA (mtDNA haplogroups) that have evolved from the same ancestors [12]. Cann et al. found the mtDNA evolution originated from a common African ancestor, named “mitochondrial Eve” [11]. Haplogroups L1, L2, and L3 are specific for the African population, whereas all non-African mtDNA haplogroups are descended from L3 and belong to either the M or the N superclades. European populations belong to nine haplogroups (i.e., H, I, J, K, T, U, V, W, and X) from the N superclades [13], while the haplogroups A, B, C, D, and E are specific to Asian populations.

Recently, mtDNA epidemiologic studies have shown that mtDNA haplogroups are associated with various phenotypes. There are significant epidemiologic associations between mtDNA haplogroup and longevity among the Japanese and Chinese [14–16], neurodegenerative [17–19] psychiatric [20, 21], cardiovascular [22, 23], and metabolic diseases [24, 25]. An association between human mtDNA variants and diseases has also been reported [26–28].

To further confirm the molecular pathogenesis and screen effective treatment for mtDNA-associated diseases, optimal cell models for reproducing the progress of diseases are valuable. A cybrid model that carries the same nuclear background and harbors variant mtDNA can exclude variations from the nuclear genome. Within analogous nuclear background, the physiologic function of mitochondrial genome can be reproduced in vitro, enabling the validation of epidemiologic observations. Thus, a complete set of cybrid cells carrying different mtDNA haplogroups are warranted to exclusively investigate the impact of mtDNA common variations on mitochondrial function. The present study successfully developed the most commonly occurring cybrids harboring mtDNA haplogroups from the Taiwanese population with ethnic Chinese background.

2. Materials and Methods

2.1. Cell Culture and Generation of Cells

Human 143B osteosarcoma (BCRC 60439) (Bioresource Collection and Research Center, Taiwan, Republic of China) were grown in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10% heat-inactivated fetal bovine serum (growth medium) at 37°C in 5% CO₂. The FBS was thawed at 37°C and then heated to 56°C for 30 min for heat inactivation. Production and culturing of ρ ⁰ cells required growth medium containing 1 mM pyruvate and 50 μg/mL uridine (ρ ⁰ medium) to support growth [29]. The ρ ⁰ cells were produced by culturing 143B osteosarcoma cells in the presence of ethidium bromide (EtBr, 50 ng/mL) for 8 weeks. The cells were subcultured every 2-3 days and kept in 50–70% of confluence.

2.2. Isolation and Transfusion of Platelets to ρ ⁰ Cells

Healthy Taiwanese volunteers donated their peripheral blood. All of them provided written informed consent in accordance with the protocols approved by the institutional review board of the contributing center. The studies were conducted according to the guidelines of the Declaration of Helsinki and adopted the methods previously reported by Chomyn [30]. Briefly, platelets were isolated and fused with ρ ⁰ cells in the presence of polyethylene glycol 1500 (50% w/v) (Roche, Nutley, NJ, USA). After 24 h of incubation, fresh ρ ⁰ medium was introduced. The cells were allowed to recover for 1 week in ρ ⁰ medium, with medium changed every 2 days. After 1 week of recovery in ρ ⁰ medium, the cells were placed on a 96-well plate in a 1 : 2 limited dilution in DMEM-based selection medium composed of 10% dialyzed FBS (Gemini), which lacks uridine and pyruvate. After 10–14 days, the cybrid populations forming single colonies were screened and transferred to a 35 mm dish. With continuing growth, the cybrids were transferred to 100 mm dish for further 20 generations of culture in order to ensure that nuclear-encoded mitochondrial proteins were incorporated into the newly introduced mitochondria.

2.3. MtDNA Copy Number Analysis for ρ ⁰ and Cybrid Cells

Total DNA was isolated from peripheral blood leucocytes using the PUREGENE DNA Purification kit (Qiagen, Chatsworth, CA, USA). The extracted DNA samples were frozen at −20°C until the assay without repeated freeze-thawing cycles. The relative mtDNA copy numbers were measured using real-time PCR and corrected by simultaneous measurement of nuclear DNA. The forward and reverse primers for nuclear gene that were complementary to chromosome 1 genomic contig (accession no. NT_004487.19) were 5′-GGCTCTGTGAGGGATATAAAGACA-3′ and 5′-CAAACCACCCGAGCAACTAATCT-3′. The forward and reverse primers for mtDNA that were complementary to the sequence of the ND2 gene (accession no. NC_012920.1) were 5′-CACAGAAGCTGCCATCAAGTA-3′ and 5′-CCGGAGAGTATATTGTTGAAGAG-3′. The PCR was performed by a LightCycler 480 Real-Time PCR System (Roche Diagnostics; F. Hoffmann-La Roche Ltd.) using the SYBR GREEN I PCR MASTER MIX (Roche Diagnostics; F. Hoffmann-La Roche Ltd.). The DNA (10 ng) was mixed with 10 μL SYBR GREEN I PCR MASTER MIX containing 400 nmol of forward and reverse primers in a final volume of 20 μL. The PCR conditions were 10 min at 95°C, followed by 40 cycles of denaturation at 95°C for 15 s, annealing at 60°C for 20 s and primer extension at 72°C for 15 s. The amplified products were denatured and reannealed at different temperatures to detect their specific melting points. Samples showing primer dimmers or unspecific fragments were excluded.

The threshold cycle number (Ct) values of the chromosome 1 sequence and the mitochondrial ND2 gene were determined for each individual in the same quantitative PCR run. Each measurement was carried out at least three times and normalized in each experiment against a serial dilution of control DNA sample. There was good reproducibility within and between runs. The intraassay coefficients of variation of Ct values were 2.2% and 3.8% for ND2 and chromosome 1, respectively. The Ct values were used as a measure of the input copy number, and Ct value differences were used to quantify the mtDNA copy number relative to the chromosome 1 using the following equation:

where results were expressed as the copy number of mtDNA/cell.

2.4. Aerobic Metabolism

Cultured cells used glucose as substrates for glycolysis to generate ATP. By replacing galactose with glucose in the culture medium, energy production shunted glycolysis to oxidation phosphorylation depending on the mitochondria [31, 32]. To examine the aerobic capacity of mitochondria, cells were grown in DMEM containing 10 mM galactose and 10% FBS in the absence of glucose.

2.5. Oxygen Consumption Analysis for and Cybrid Cells

Detection of O₂ consumption in permeabilized cells was as described previously with some modifications [33–35]. O₂ consumption was monitored by a Clark electrode (Mitocell S200 micro respirometry system; Strathkelvin Instruments, Motherwell, UK). Cells (100 μL at 5 × 10⁶ cells/mL) in KCl medium (100 mM KCl, 3 mM MgCl₂, 20 mM HEPES, 1 mM EDTA, and 5 mM KH₂PO₄; pH 7.4) were permeabilized by digitonin (optimal concentration 32.5 μg/mL as determined by trypan blue staining) and loaded into a 200 μL MT200 Respirometer Chamber, suspended by a fixed-speed solid-state magnetic stirrer inside the chamber and maintained at 37°C by a circulating water bath. Glutamate/malate (10 mM) and succinate (10 mM), combined with ADP (0.1 mM), were used as substrates to test the function of complex I- and complex II-linked respiration. Changes in O₂ levels were calculated by respirometry software. For control recordings, 100 μL of buffer was injected into the chamber.

2.6. Measurement of Mitochondrial Membrane Potential

Cells (5 × 10⁵) were placed in six-well plates and allowed to attach for 16–18 h and then incubated in 100 nM TMRE (Sigma) for 30 min at 37°C. They were then washed twice with PBS, harvested by trypsinization, and resuspended in 500 μL PBS containing 2% FBS. Cytofluorimetric analysis was conducted using a fluorescence-activated cell scanner machine (BD Biosciences FACScan system).

2.7. Cellular Viability Assay

Cells (8 × 10³) were seeded into 96-well plates 24 h prior to treatment. Cell survival was determined using cell counting kit 8 (Wako) with WST-8 [2-(2-methoxy-4-nitrophenyl)-3-(4-nitrophenyl)-5-(2,4-disulfophenyl)-2-H-tetrazolium, monosodium salt] as substrate. The WST-8 was reduced by dehydrogenase to form formazan that was proportional to living cells. The amount of formazan was assayed by spectrophotometry at OD450.

2.8. Western Blot

Cells were lysed in RIPA lysis buffer (50 mM Tris, pH 7.4; 150 mM NaCl; 1 mM PMSF (phenylmethanesulfonyl fluoride); 1 mM EDTA (ethylenediaminetetraacetic acid); 1% Triton X-100; 1% sodium deoxycholate; 0.1% SDS) with the addition of Protease Inhibitor Cocktail (Roche) and Phosphatase Inhibitor Cocktail I (Sigma). The antibodies used included anti-COX II antibody (Abcam), peroxidase-labeled anti-rabbit IgG (H + L) secondary antibody (Abcam), and anti-β-actin (Cell Signaling, Danvers, MA). The signals were developed by ECL plus (GE Healthcare Bio-Sciences AB, Uppsala, Sweden) using X-ray films.

2.9. Genotyping of Mitochondrial Haplogroups

The selected polymorphic sites were based on Nishigaki’s and Kong’s work [36] (Table 1). Venous blood (7 mL) was collected from each subject into tubes containing 50 mmol/L EDTA (disodium salt). Genomic DNA was isolated using a commercial kit (Qiagen, Chatsworth, CA, USA). For amplifying mtDNA fragments, 24 pair primers for the mitochondrial DNA were used for gene amplification by PCR using primer sequences described previously [28]. Amplicon size ranged from 190 to 300 bp. The PCR amplification was conducted in a final reaction volume of 20 μL, containing 10 ng human genomic DNA, 10 mM Tris-HCl (pH 8.3), 50 mM KCl, 1.5 mM MgCl₂, 0.2 mM dNTP (i.e., dATP, dTTP, dCTP, and dGTP), 0.5 μM of each primer, and 1 unit of ExTaq DNA polymerase (Takara Bio Inc.). The PCR conditions were as follows: an initial denaturation step at 94°C for 10 min, followed by 40 cycles of denaturation at 94°C for 35 s, annealing at 60°C for 45 s, and extension at 72°C for 45 s, with a final extension of 10 min at 72°C and kept at 25°C using an Eppendorf Mastercycler gradient.

Table 1: Polymorphism sites characteristic [B4a, B4b, B4c, B5, R9, F1a, F2, D4a, D4e, D5a, D5b, E, M8, C, and N9a].

Haplogroup
Polymorphism

B4a
8272 (9-bp deletion), G9123A

B4b
8272 (9-bp deletion), C15535T, G4820A

B4c
8272 (9-bp deletion), G15346A

B4d
8272 (9-bp deletion), C15535T, A13942G

B5
8272 (9-bp deletion), T9950C, A10398G

R
G13928C, C3970T

F1a
C3970T, G10310A, G13928C, G12406A, C4086T

F2
G13928C, C3970T, G10310A, T12338C

D4a
G3010A, C5178A, C4883T, T14979C

D4e
G3010A, C5178A, C4883T, A15924G

D5a
C5178A, C4883T, A10397G

D5b
T681C, C5178A, C4883T, A10397G

E
G4491A, T3027C

M8
A4715G, C8684T, G6179A, A15487T

C
A4715G, A15487T, A13263G

N9a
G5231A, G12372A, A12358G

After PCR amplification, 60 probes were used for mitochondrial haplogroup definition. Probes for haplotyping were as described previously [28]. Briefly, synthesized oligonucleotide probes were modified at the 5′-end with a terminal anion group and covalently bound to carboxylated fluorescent microbeads using ethylene dichloride (EDC), following the manufacturer’s instructions. Oligonucleotide-labeled microbeads (oligobeads) were mixed together to make a 5000/5 μL oligobead mixture for hybridization. The 5′-biotin-labeled PCR amplicons were hybridized to the oligobeads in a total volume of 50 μL per well in a 96-well plate by adding 5 μL of the PCR amplicon to 40 μL of hybridization buffer (3.75 M TMAC, 62.5 mM TB [pH 8.0], 0.5 mM EDTA, and 0.125% N-lauroylsarcosine). This reaction mixture was denatured at 95°C for 2 min and then hybridized at 55°C for 30 min using the Eppendorf Mastercycler gradient.

After hybridization, the oligobeads were washed with 1× PBS (pH 7.5, containing 0.01% Tween 20) and centrifuged at 2200 ×g for 2 min. The supernatant was discarded (washed twice), and the pelleted oligobeads were reacted with 70 μL aliquot of a 1000× diluted solution of SA-PE (Prozyme). The hybridized amplicons were labeled with SA-PE at 52°C for 15 min using the Eppendorf Mastercycler gradient. Then the reactions were measured by the Luminex¹⁰⁰ flow cytometer.

2.10. Determination of mtDNA Full-Length Sequence

2.10.1. First PCR

The entire mitochondrial genome was amplified as six fragments, each approximately 3.0 kb in length, by a symmetric PCR method with primer pairs as stated previously [28]. The PCR amplification was conducted in a final reaction volume of 20 μL, containing 50 ng human genomic DNA, 10 mM Tris-HCl (pH 8.3), 50 mM KCl, 1.5 mM MgCl₂, 0.2 mM concentration of each dNTP, 1.0 μM concentration of each primer, and 1 unit of pfu DNA polymerase (Takara Bio Inc.). The PCR conditions used were as follows: an initial denaturation step at 94°C for 5 min, followed by 40 cycles of denaturation at 94°C for 15 s, annealing at 58°C for 35 s, and extension at 72°C for 3 min, with a final extension of 10 min at 72°C. The amplified fragments were analyzed by electrophoresis on a 1% agarose gel and visualized by ethidium bromide staining.

2.10.2. Second PCR

The first PCR DNA templates for sequence analysis of the entire mitochondrial genome were amplified as 32 overlapping segments, each approximately 500–1100 bp as shown previously [28]. The PCR amplifications were conducted in a final reaction volume of 20 μL, containing 200 ng of the first PCR product, 10 mM Tris-HCl (pH 8.3), 50 mM KCl, 1.5 mM MgCl₂, 0.2 mM concentration of each dNTP, 1.0 μM concentration of each primer, and 0.5 units of pfu DNA polymerase (Takara Bio Inc.). The PCR conditions used were the same as those for the first PCR except that the annealing temperature used was 60°C. The second PCR fragments were analyzed by electrophoresis on a 1.5% agarose gel and visualized by ethidium bromide staining.

2.11. Sequence Analysis and mtSNP Identification

The second PCR template and primers were sent to the Genome Research Center, Sequencing Center of National Yang-Ming University, which performed the sequence using Applied Biosystems. The primers used for sequencing were as described previously [28].

3. Results

3.1. Production of ρ ⁰ Cells from 143B Osteosarcoma Cells

EtBr was bound to the double-helix DNA and inhibited replication. MtDNA was more susceptible to EtBr than nuclear DNA in an indicated concentration of EtBr. Low dosages of EtBr were used for long-term exposure to fully deplete mtDNA while maintaining cellular growth. The reversion of mtDNA in ρ ⁰ cell in the absence of EtBr was reported to be minimized by prolonged culturing [37]. Thus, instead of shorter durations of EtBr exposure [29, 38], 143B cells were grown in EtBr for 8 weeks to reduce the possibility of mtDNA reversion. In order to obtain stable clones, individual clones were further isolated. The reversion of mtDNA was examined in the absence of EtBr.

The mtDNA of 143B cells decreased to 0.93 copy/cell after 8 weeks (Table 2). However, mtDNA reversed to the level of the parental 143B cells after EtBr removal for 14 days. The possibility of mtDNA reversion was eliminated by isolating single clones, and six of 16 isolated clones maintained mtDNA; a level of less than one copy per cell after EtBr withdrawal for 50 days. Clone no. 38 was selected for another 60 days of culture without EtBr to confirm the irreversibility of 143B-derived ρ ⁰ cells (Table 2). The ρ ⁰ cells died within 4 weeks without pyruvate and uridine (data not shown).

Table 2: mtDNA copy number of 143B osteosarcoma following EtBr exposure.

Cell type
Days without EtBr
mtDNA (copy/cell)

143B
—
1060.11

08w
0
0.93

8w revertant
14
1009.9

Clone #23
50
0.98

Clone #24
50
0.86

Clone #31
50
0.98

Clone #35
50
0.85

Clone #36
50
0.97

Clone #38
50
0.97

Clone #38
110
0.69

143B: parental osteosarcoma; 8w: 143B following 8 weeks of EtBr exposure; 8w revertant: recovered mtDNA level in the absence of EtBr for 14 days.

Clones 23, 24, 31, 35, 36, and 38 were selected clones through limiting dilution.

3.2. Loss of Mitochondrial Oxygen Consumption and Defects in Respiratory Enzyme Complexes

Since subunit 2 of cytochrome c oxidase (COX2) was mtDNA encoded, COX2 was examined for ρ ⁰ cells. Immunoblotting demonstrated undetectable COX2 of ρ ⁰ cells (Figure 1). In terms of oxygen consumption ability of the ρ ⁰ cells, mitochondrial complex I- and II-linked respiration monitored in the presence of glutamate/malate and succinate, respectively, were deficient in ρ ⁰ cells (Figure 3).

(a)

(b)

Figure 1: (a) The absence of mtDNA-encoded protein in ρ ⁰ cells. The expression of mtDNA-encoded COX2 and nuclear-encoded β-actin in the parental 143B osteosarcoma and ρ ⁰ cell was analyzed by Western blotting. (b) The ρ ⁰ cells failed to survive in the galactose medium over 48 h, in contrast to the normal growth of the parental 143B cells.

3.3. Production of Cybrids Harboring Different Haplogroups from Human Platelet Mitochondria

Sixteen different haplogroups were selected from healthy volunteers for mitochondrial genotyping. Platelets from the 16 known mtDNA haplogroups were fused with ρ ⁰ cells, followed by selection. All of the transmitochondrial experiments contained negative control groups in which ρ ⁰ cells underwent the same fusion process in the absence of human platelets. All of the negative control cells died after selection, while cells fused with platelets survived (data not shown), suggesting that fusion was successful, and the mtDNA of ρ ⁰ cells were unable to revert.

3.4. Mitochondrial Respiratory Function and Mitochondrial DNA Content Reversed after Transmitochondrial Fusion

The mtDNA copy number of parental 143B, ρ ⁰ cells and cybrids harboring human mtDNA suggested that the depleted mtDNA content in ρ ⁰ cells was reversed following transmitochondrial fusion (Figure 2). Because most cellular metabolisms relied on normal mitochondria, the respiratory function of the cybrids was examined. Parental 143B cells, ρ ⁰ cells, and cybrids harboring the B5 haplogroup assayed for oxygen consumption rate indicated that ρ ⁰ cells with depleted mtDNA demonstrated loss of respiratory function, whereas transmitochondrial fusion successfully restored the lost respiration (Figure 3).

Figure 2: The mtDNA content of 143B osteosarcoma, ρ ⁰ cell, and the cybrid harboring the B5 haplogroup was analyzed by real-time PCR and shown as mtDNA copy number. Data was expressed as mean ± SEM for three independent experiments.

(a)

(b)

(c)

Figure 3: Oxygen consumption of 143B osteosarcoma, ρ ⁰ cell, and the cybrid harboring the mtDNA haplogroup of B5. (a) Complex I-linked respiration was performed by adding 10 mM glutamate and 10 mM malate in the presence of 0.1 mM ADP. The complex I inhibitor, rotenone, terminated the respiration. (b) Complex II-linked respiration preincubated with rotenone was performed by adding succinate in the presence of 0.1 mM ADP. The complex III inhibitor, antimycin A, terminated the respiration. (c) Quantitative data of oxygen consumption expressed in slope were shown as mean ± SEM for three independent experiments.

3.5. Mitochondrial Polymorphisms Generated Cybrids Harboring Identical Mitochondrial Genomarkers as the Original Platelet Donors

To verify the mtDNA haplogroup of the 16 cybrids, polymorphism sites were reexamined by Luminex 1000. All 16 cybrids had similar mtDNA haplogroups as their original donors (data not shown). To further verify the mtDNA haplogroup of these constructed cybrids, full-length sequencing for mtDNA was performed. Nonsynonymous sites in mtDNA genome caused amino acid substitution in the 16 cybrids (Tables 3 and 4), showing characterizations of the genetic aspect and corresponding features of the amino acid.

Table 3: Nonsynonymous substitution of 16 cybrid cell lines.

Gene
Substitution
AA change
B4a
B4b
B4c
B4d
B5
R
F1a
F2
D4a
D4e
D5a
D5b
E
M8
C4
N9a
143B

G4491A
VI
—
—
—
—
—
—
—
—
—
—
—
—
A
—
—
—
—

G4596A
VI
—
—
—
—
—
—
—
—
—
—
—
—
—
—
A
—
—

ND2
C5178A
LM
—
—
—
—
—
—
—
—
A
A
A
A
—
—
—
—
—

A5301G
IV
—
—
—
—
—
—
—
—
—
—
G
G
—
—
—
—
—

A5466G
TA
—
—
—
—
—
—
—
—
G
—
—
—
—
—
—
—
—

T6253C
MT
—
—
—
—
—
—
—
—
—
—
—
C
—
—
—
—
—

COI
G6267A
AT
—
—
—
—
—
—
—
—
—
—
—
—
—
—
—
—
A

T7266G
SA
G
G
G
—
—
—
—
—
—
—
—
—
—
—
—
—
—

COII
G7598A
AT
—
—
—
—
—
—
—
—
—
—
—
—
A
—
—
—
—

A7934G
IV
—
—
—
—
—
—
—
—
—
—
—
—
G
—
—
—
—

ATPase8
C8414T
LI
—
—
—
—
—
—
—
—
T
T
—
—
—
—
—
—
—

G8584A
AT
—
—
—
—
A
—
—
—
—
—
—
—
—
A
A
—
—

C8684T
TI
—
—
—
—
—
—
—
—
—
—
—
—
—
T
—
—
—

A8701G
TA
—
—
—
—
—
—
—
—
G
G
G
G
G
G
G
—
—

ATPase6
A8860G
TA
G
G
G
G
G
G
G
G
G
G
G
G
G
G
G
G
G

G9053A
SN
—
—
A
—
—
—
A
—
—
—
—
—
—
—
—
—
—

C9099T
IM
—
—
—
—
—
—
—
T
—
—
—
—
—
—
—
—
—

T9128C
IT
—
—
—
—
—
—
—
—
—
—
—
C
—
—
—
—
—

T10084C
IT
—
C
—
—
—
—
—
—
—
—
—
—
—
—
—
—
—

ND3
A10398G
TA
—
—
—
—
G
—
—
—
G
G
G
G
G
G
G
—
—

C10400T
TA
—
—
—
—
—
—
—
—
T
T
T
T
T
T
T
—
—

ND4L
T10609C
MT
—
—
—
—
—
—
C
—
—
—
—
—
—
—
—
—
—

11138T
(Insertion)

T
T

ND4
G11969A
AT
—
—
—
—
—
—
—
—
—
—
—
—
—
—
A
—
—

A12026G
IV
—
—
—
—
—
—
—
—
—
—
G
—
—
—
—
—
—

T12338C
MT
—
—
—
—
—
—
—
C
—
—
—
—
—
—
—
—
—

A12358G
TA
—
—
—
—
—
—
—
—
—
—
—
—
—
—
—
G
—

G12406A
VI
—
—
—
—
—
—
A
—
—
—
—
—
—
—
—
—
—

G13135A
AT
—
—
—
—
—
—
—
—
—
A
—
—
—
—
—
—
A

G13708A
AT
—
—
—
A
—
—
—
A
—
—
—
—
—
—
—
—
—

ND5
G13759A
AT
—
—
—
—
—
—
A
—
—
—
—
—
—
—
—
—
—

A13834G
TA
—
—
—
—
—
—
—
—
G
—
—
—
—
—
—
—
—

G13928C
ST
—
—
—
—
—
C
C
C
—
—
—
—
—
—
—
—
—

A13942G
TA
—
—
—
G
—
—
—
—
—
—
—
—
—
—
—
—
—

13944C
(Insertion)

C

A13966G
TA
—
—
—
—
—
—
—
—
—
—
—
—
—
—
—
—
G

ND6
T14318C
NS
—
—
—
—
—
—
—
—
—
—
—
—
—
—
C
—
—

T14577C
IV
—
—
—
—
—
—
—
—
—
—
—
—
C
—
—
—
—

C14751T
TI
T
—
—
—
—
—
—
—
—
—
—
—
—
—
—
—
—

C14766T
IT
—
—
—
—
—
—
—
—
—
—
—
T
—
—
—
—
—

Cytb
T14979C
IT
—
—
—
—
—
—
—
—
C
—
—
—
—
—
—
—
—

T15204C
IT
—
—
—
—
—
—
—
—
—
—
—
—
—
—
C
—
—

A15326G
TA
G
G
G
G
G
G
G
G
G
G
G
G
G
G
G
G
G

A15758G
IV
—
G
—
—
—
—
—
—
—
—
—
—
—
—
—
—
—

Nucleotide change, corresponding amino acid change, and the reference mitochondrial sequence were based on MITOMAP (http://www.mitomap.org/MITOMAP).

AA: amino acid; ND2: NADH dehydrogenase subunit 2; COI: cytochrome c oxidase subunit I; COII: cytochrome c oxidase subunit II; ATPase8: ATP synthase F0 subunit 8; ATPase6: ATP synthase F0 subunit 6; ND3: NADH dehydrogenase subunit 3; ND4L: NADH dehydrogenase subunit 4L; ND4: NADH dehydrogenase subunit 4; ND5: NADH dehydrogenase subunit 5; ND6: NADH dehydrogenase subunit 6; Cytb: cytochrome b.

Table 4: Summary of changes in each gene in each haplogroup.

B4a
B4b
B4c
B4d
B5
R
F1a
F2
D4a
D4e
D5a
D5b
E
M8
C4
N9a
143B

ND2

C5178A,
C5178A
C5178A,
C5178A,
G4491A

G4596A

A5466G

A5301G
A5301G

ND3

T10084C

A10398G

A10398G,
A10398G,
A10398G,
A10398G,
A10398G,
A10398G,
A10398G,

C10400T
C10400T
C10400T
C10400T
C10400T
C10400T
C10400T

ND4

A12026G

ins11138T
ins11138T,

G11969A

ND4L

T10609C

G13708A,

13928,
G12406A,
T12338C,
A13834G
G13135A

A12358G
G13135A,

ND5

A13942G

insl3944C
G13759A,
G13708A,

A13966G

G13928C
G13928C,

ND6

T14577C

T14318C

Cytb
T14751C,
A15326G,
A15326G
A15326G
A15326G
A15326G
A15326G
A15326G
T14979C,
A15326G
A15326G
C14766T,
A15326G
A15326G
T15204C,
A15326G
A15326G

A15326G
A15758G

A15326G

A15326G

A15326G

COI
T7266G
T7266G
T7266G

T6253C

G6267A

COII

G7598A,

A7934G

ATPase6
A8860G
A8860G
A8860G,
A8860G
G8584A,
A8860G
A8860G,
A8860G,
A8701G,
A8701G,
A8701G,
A8701G,
A8701G,
G8584A,
G8584A,
A8860G
A8860G

G9053A

A8860G

G9053A
C9099T
A8860G
A8860G
A8860G
A8860G,
A8860G
C8684T,
A8701G,

T9128C

A8701G,
A8860G

A8860G

ATPase8

C8414T
C8414T

ND2: NADH dehydrogenase subunit 2; COI: cytochrome c oxidase subunit I; COII: cytochrome c oxidase subunit II; ATPase8: ATP synthase F0 subunit 8; ATPase6: ATP synthase F0 subunit 6; ND3: NADH dehydrogenase subunit 3; ND4L: NADH dehydrogenase subunit 4L; ND4: NADH dehydrogenase subunit 4; ND5: NADH dehydrogenase subunit 5; ND6: NADH dehydrogenase subunit 6; cytb: cytochrome b; ins: insertion.

3.6. Cybrids with Various mtDNA Haplogroups Exhibited Differences in Oxygen Consumption Capacity and Mitochondrial Membrane Potential

Haplogroups N9a and D4a were reported to be protective against diabetes [24, 25] and promoted longevity [14], respectively. Haplogroup B4 was diabetes susceptible, while haplogroup D4 was diabetes-resistant [28]. To clarify if mtDNA genetic variations among these haplogroups influenced mitochondrial respiratory capacity, cellular oxygen consumption rates were examined in cell lines created from these haplogroups. Cybrid-harboring mtDNA haplogroups N9a and D4a had significantly higher oxygen consumption rates compared to cybrid-harboring mtDNA haplogroups B4b, F1a, B5, and D5a in complex I- and II-linked respiration (Figure 4(a)). Moreover, compared to cybrid-harboring mtDNA haplogroup B4b, haplogroups F1a, B5, and D5a also exhibited higher O₂ consumption rates in complex I-linked respiration.

(a)

(b)

Figure 4: Oxygen consumption rate and mitochondrial membrane potential of cybrids harboring B4b, F1a, B5, D5a, D4a and N9a. (a) Complex I-linked respiration was performed by adding 10 mM glutamate and 10 mM malate in the presence of 0.1 mM ADP. Complex II-linked respiration pre-incubated with rotenone was performed by adding succinate in the presence of 0.1 mM ADP. (b) The mean TMRE fluorescence value of B4b was set to 100%. Data were expressed as mean ± SEM for four independent experiments. The statistical significance of differences between groups was determined by one-way analysis of variance followed by Bonferroni post hoc test. compared to (1) B4b, (2) F1a, (3) B5, (4) D5a, and (5) D4a.

Measuring mitochondrial membrane potential (ΔΨm) among these cybrid lines, the ΔΨm of cybrid-harboring mtDNA haplogroups B4b and F1a was significantly higher compared to those of B5, D5a, D4a, and N9a (Figure 4(b)). This was consistent with a previous report that a cybrid with higher oxygen consumption rate had lower mitochondrial membrane potential [39].

3.7. Cybrids Harboring Different mtDNA Haplogroups Exhibited Discriminating Tolerance against Oxidative Stress

To assess whether differences in OXPHOS capacity affected cellular tolerance to oxidative stress, the viability of cybrid cells exposed to various doses of hydrogen peroxide (H₂O₂) was examined. Under higher H₂O₂ concentration, cybrids-harboring mtDNA haplogroups F1a, D4a, and N9a had significant resistance and higher cellular viability compared to B4b (Figure 5), suggesting that cybrid B4b was much more sensitive to oxidative stress induced by H₂O₂. Notably, cybrid-harboring mtDNA haplogroup N9a demonstrated superior viability under oxidative stress compared to those with B4b, F1a, B5, D5a, and D4a.

Figure 5: Cellular viability of cybrids under various concentration of H₂O₂ for 4 h. WST-8 assay for viability was determined by spectrophotometry at 450 nm. Data (mean ± SEM, n = 8) expressed as percentage relative to B4b in the untreated control was determined by one-way analysis of variance followed by Bonferroni post hoc test. compared to B4b and to (1) F1a, (2) B5, (3) D5a, and (4) D4a.

4. Discussion

The present study illustrates the process behind the construction of cybrid lines using samples from a representative sample of the Taiwanese population harboring different mtDNA haplogroups. The study further demonstrates that the mitochondria of recomposed cells are functionally intact. The ρ ⁰ cells defective in mtDNA and lacking in O₂ consumption ability have been successfully constructed from 143B osteosarcoma cells. Their defective mitochondrial function is restored after transmitochondrial fusion with mitochondria donated from volunteer platelet donors of different mtDNA haplogroups. The recomposed cybrid lines generated and tested systemically harbor the sixteen most commonly found mtDNA haplogroups identified in the Taiwanese population.

Associations between degenerative diseases and mtDNA haplogroups have previously been documented. Epidemiologic studies reveal that certain mtDNA haplogroups are prone to disease generation, whereas other haplogroups are resistant. While the relationship between mtDNA variants and disease generation is generally accepted, some studies on the association of clinical diseases and haplogroup variations remain controversial. In terms of the 16189 T to C transition variant of mtDNA and its association with type II diabetes, both positive and negative reports have been published [27]. Thus, reliable and reproducible material is essential for studies concerning the impact of mtDNA haplogroups on the development of these diseases.

Approximately two decades ago, Attardi’s group demonstrated the construction of transmitochondrial cybrids and provided a platform for studying mitochondrial function in cells carrying mtDNA variations in the same nuclear background [29]. Since then, studies using cybrids have shown the importance of mitochondrial functions and genetics for normal physiology. An unpublished finding has previously suggested that Taiwanese harboring the haplogroup B4 have higher prevalence of diabetes mellitus, whereas those harboring haplogroup B5 are resistant to Parkinson’s disease. Other studies involving Eastern Asian subjects also reveal that certain haplogroups exhibit dominant physiologic performance, including resistance to Parkinson’s disease generation [18]. However, to date, there is still no complete set of cybrid cells carrying different mtDNA haplogroups from the general Taiwanese population. This study, therefore, has created cybrid lines carrying the most commonly occurring haplotypes found among ethnic Chinese population for future functional and mechanical studies.

One of the problems with constructing cybrid cell lines is the reversion of the ρ ⁰ cells’ mtDNA, causing contamination of the parental mtDNA. EtBr chelation may distort the double-strand DNA structure and is widely used in producing mtDNA-less or mtDNA-depleted cells. Appropriate EtBr dosage chelating mtDNA with minimal damage to the nucleus is considered ideal for interfering with the replication process of mtDNA without affecting cell survival. Very long-term exposure to EtBr can generate stable ρ ⁰ cells with permanent loss of mtDNA. As such, to harvest stable ρ ⁰ cells, a modest period of exposure followed by isolating single-cell clones has been used in the present study. The data suggest that these ρ ⁰ cells must be stable to prevent mtDNA reversion and hence ensure the mtDNA consistency of subsequent cybrids. Characteristics were examined by genetic testing, translational protein assessment, and functional respiration testing to verify the validity of transmitochondrial fusion. In addition, the identity of mtDNA haplogroups of original donors and cybrids has been reconfirmed by multiple probes, while full-length sequencing has further demonstrated the solidity of cybrids.

Physiologic differences in common European mtDNA haplogroups using the cybrid model have been reported in the works of Pello et al. [40] and of Gómez-Durán et al. [39]. Thus, it is feasible to test the mitochondrial respiratory function of the cybrids generated in this study. A previous work has identified the diabetes-susceptible group (haplogroup B4) to consist of 10398A and diabetes-resistant group (haplogroup D4) to consist of C5178A, A10398G, and T152C [28]. As was expected, the cybrid harboring the mtDNA haplogroup B4b that carries the 10398A variant is more susceptible to developing diabetes and has a significantly lower oxygen consumption capacity compared to other commonly seen haplogroups in Asian populations of ethnic Chinese background (Figure 4(a)). Consistent with a previous clinical genetic study, cybrids harboring the specific mtDNA haplogroup carrying A10398G (B5) or in combination with C5178A and A10398G (D4a and D5a) display superior oxygen consumption capacity compared to cybrids with B4b (Figure 4(a)). Furthermore, the superior oxygen consumption capacity of cybrids with N9a and D4a is also compatible with findings of Tanaka and colleagues, who reported that D4a is related to a longer lifespan, and haplogroups N9a and D5 are resistant against metabolic syndromes [14, 25].

Closely related to oxygen consumption and cellular energy need, the value of ΔΨm is a key indicator of cellular metabolic activity. Compatible with a previous report [39], ΔΨm is inversely related to the decrease in oxygen consumption rate in the cybrid related to disease conditions (B4b) (Figure 4(b)). This suggests that efficient energy production through the OXPHOS machinery is associated with less ΔΨm and related to cellular resistance to stressful environment. Testing cellular response of the cybrid cell lines to H₂O₂-induced oxidative stress, cybrid cells carrying various mitochondrial genomic backgrounds have significant differences in survival rates in response to ROS exposure (Figure 5). Consistent with oxygen consumption data, the cybrid harboring the haplogroup N9a is the most resistant to H₂O₂-induced damage.

In conclusion, the present study has generated cybrid lines harboring the most commonly occurring haplogroups in the Taiwanese population of ethnic Chinese background. Mitochondrial respiratory function, membrane potential, and cellular survival in response to oxidative stress are also different among these common cybrid lines. This set of cybrid cell lines may be useful for future studies regarding the impact of common human mtDNA variances on disease generation and progression.

Authors’ Contribution

Drs. T.-K. Lin and H.-Y. Lin contributed equally to this work.

Acknowledgments

The work was supported by Grants from Chang Gung Memorial Hospital (CMRPG870881-3; CMRPG891101-2; CMRPG891081-2) and the National Science Council, Taiwan, (NMRPG896131-2). The authors thank the Chang Gung Medical Foundation Kaohsiung Chang Gung Memorial Hospital Tissue Bank (CLRPG8B0031) for the technical support, and James Waddell of the Han Mei Institute and ATS Medical Editing and Review Solutions for the English editing and proofreading services.

References

M. W. Gray, G. Burger, and B. F. Lang, “Mitochondrial evolution,” Science, vol. 283, no. 5407, pp. 1476–1481, 1999.

S. Anderson, A. T. Bankier, and B. G. Barrell, “Sequence and organization of the human mitochondrial genome,” Nature, vol. 290, no. 5806, pp. 457–465, 1981.

R. O. Poyton and J. E. McEwen, “Crosstalk between nuclear and mitochondrial genomes,” Annual Review of Biochemistry, vol. 65, pp. 563–607, 1996.

Y. H. Wei, “Oxidative stress and mitochondrial DNA mutations in human aging,” Proceedings of the Society for Experimental Biology and Medicine, vol. 217, no. 1, pp. 53–63, 1998.

D. R. Johns, “Seminars in medicine of the Beth Israel Hospital, Boston. Mitochondrial DNA and disease,” New England Journal of Medicine, vol. 333, no. 10, pp. 638–644, 1995.

I. E. Scheffler, “Mitochondria make a comeback,” Advanced Drug Delivery Reviews, vol. 49, no. 1-2, pp. 3–26, 2001.

D. C. Wallace, “Mitochondrial diseases in man and mouse,” Science, vol. 283, no. 5407, pp. 1482–1488, 1999.

P. F. Chinnery and D. M. Turnbull, “Clinical features, investigation, and management of patients with defects of mitochondrial DNA,” Journal of Neurology Neurosurgery and Psychiatry, vol. 63, no. 5, pp. 559–563, 1997.

E. Martínez-Fernández, A. Gil-Peralta, R. García-Lozano et al., “Mitochondrial disease and stroke,” Stroke, vol. 32, no. 11, pp. 2507–2510, 2001.

T. I. Rogounovitch, V. A. Saenko, Y. Shimizu-Yoshida et al., “Large deletions in mitochondrial DNA in radiation-associated human thyroid tumors,” Cancer Research, vol. 62, no. 23, pp. 7031–7041, 2002.

R. L. Cann, M. Stoneking, and A. C. Wilson, “Mitochondrial DNA and human evolution,” Nature, vol. 325, no. 6099, pp. 31–36, 1987.

D. A. Merriwether, A. G. Clark, S. W. Ballinger, et al., “The structure of human mitochondrial DNA variation,” Journal of Molecular Evolution, vol. 33, no. 6, pp. 543–555, 1991.

V. Macaulay, M. Richards, E. Hickey et al., “The emerging tree of west Eurasian mtDNAs: a synthesis of control-region sequences and RFLPs,” American Journal of Human Genetics, vol. 64, no. 1, pp. 232–249, 1999.

E. Bilal, R. Rabadan, G. Alexe et al., “Mitochondrial DNA haplogroup D4a is a marker for extreme longevity in Japan,” PLoS ONE, vol. 3, no. 6, Article ID e2421, 2008.

W. H. Ren, X. H. Li, and H. G. Zhang, “Mitochondrial DNA haplogroups in a Chinese Uygur population and their potential association with longevity,” Clinical and Experimental Pharmacology and Physiology, vol. 35, no. 12, pp. 1477–1481, 2008.

S. Takasaki, “Mitochondrial SNPs associated with Japanese centenarians, Alzheimer's patients, and Parkinson's patients,” Computational Biology and Chemistry, vol. 32, no. 5, pp. 332–337, 2008.

P. F. Chinnery, G. A. Taylor, N. Howell et al., “Mitochondrial DNA haplogroups and susceptibility to AD and dementia with Lewy bodies,” Neurology, vol. 55, no. 2, pp. 302–304, 2000.

J. M. van der Walt, K. K. Nicodemus, E. R. Martin, et al., “Mitochondrial polymorphisms significantly reduce the risk of Parkinson disease,” American Journal of Human Genetics, vol. 72, no. 4, pp. 804–811.

M. Ban, J. Elson, A. Walton et al., “Investigation of the role of mitochondrial DNA in multiple sclerosis susceptibility,” PLoS ONE, vol. 3, no. 8, Article ID e2891, 2008.

C. Magri, R. Gardella, S. D. Barlati, P. Valsecchi, E. Sacchetti, and S. Barlati, “Mitochondrial DNA haplogroups and age at onset of schizophrenia,” American Journal of Medical Genetics, Part B, vol. 144, no. 4, pp. 496–501, 2007.

B. Rollins, M. V. Martin, P. A. Sequeira et al., “Mitochondrial variants in schizophrenia, bipolar disorder, and major depressive disorder,” PLoS ONE, vol. 4, no. 3, Article ID e4913, 2009.

Y. Nishigaki, Y. Yamada, N. Fuku et al., “Mitochondrial haplogroup A is a genetic risk factor for atherothrombotic cerebral infarction in Japanese females,” Mitochondrion, vol. 7, no. 1-2, pp. 72–79, 2007.

Y. Nishigaki, Y. Yamada, N. Fuku et al., “Mitochondrial haplogroup N9b is protective against myocardial infarction in Japanese males,” Human Genetics, vol. 120, no. 6, pp. 827–836, 2007.

N. Fuku, K. S. Park, Y. Yamada et al., “Mitochondrial haplogroup N9a confers resistance against type 2 diabetes in Asians,” American Journal of Human Genetics, vol. 80, no. 3, pp. 407–415, 2007.

M. Tanaka, N. Fuku, Y. Nishigaki et al., “Women with mitochondrial haplogroup N9a are protected against metabolic syndrome,” Diabetes, vol. 56, no. 2, pp. 518–521, 2007.

P. W. Wang, T. K. Lin, S. W. Weng, and C. W. Liou, “Mitochondrial DNA variants in the pathogenesis of type 2 diabetes—relevance of asian population studies,” Review of Diabetic Studies, vol. 6, no. 4, pp. 237–246, 2009.

C. W. Liou, T. K. Lin, J. B. Chen et al., “Association between a common mitochondrial DNA D-loop polycytosine variant and alteration of mitochondrial copy number in human peripheral blood cells,” Journal of Medical Genetics, vol. 47, no. 11, pp. 723–728, 2010.

C. W. Liou, J. B. Chen, and M. M. Tiao, “Mitochondrial DNA coding and control region variants as genetic risk factors for type 2 diabetes mellitus,” Diabetes, vol. 61, no. 10, pp. 2642–2651, 2012.

M. P. King and G. Attardi, “Human cells lacking mtDNA: repopulation with exogenous mitochondria by complementation,” Science, vol. 246, no. 4929, pp. 500–503, 1989.

A. Chomyn, “Platelet-mediated transformation of human mitochondrial DNA-less cells,” Methods in Enzymology, vol. 264, pp. 334–338, 1996.

L. J. Reitzer, B. M. Wice, and D. Kennell, “Evidence that glutamine, not sugar, is the major energy source for cultured HeLa cells,” Journal of Biological Chemistry, vol. 254, no. 8, pp. 2669–2676, 1979.

V. M. Gohil, S. A. Sheth, R. Nilsson et al., “Nutrient-sensitized screening for drugs that shift energy metabolism from mitochondrial respiration to glycolysis,” Nature Biotechnology, vol. 28, no. 3, pp. 249–255, 2010.

S. Djafarzadeh, M. Vuda, J. Takala, M. Ochs, and S. M. Jakob, “Toll-like receptor-3-induced mitochondrial dysfunction in cultured human hepatocytes,” Mitochondrion, vol. 11, no. 1, pp. 83–88, 2011.

D. L. Granger and A. L. Lehninger, “Sites of inhibition of mitochondrial electron transport in macrophage-injured neoplastic cells,” Journal of Cell Biology, vol. 95, no. 2, pp. 527–535, 1982.

M. Vrbacký, J. Krijt, Z. Drahota, and Z. Mělková, “Inhibitory effects of Bcl-2 on mitochondrial respiration,” Physiological Research, vol. 52, no. 5, pp. 545–554, 2003.

Q. P. Kong, H. J. Bandelt, C. Sun et al., “Updating the East Asian mtDNA phylogeny: a prerequisite for the identification of pathogenic mutations,” Human Molecular Genetics, vol. 15, no. 13, pp. 2076–2086, 2006.

S. W. Miller, P. A. Trimmer, W. D. Parker Jr., and R. E. Davis, “Creation and characterization of mitochondrial DNA-depleted cell lines with “neuronal-like” properties,” Journal of Neurochemistry, vol. 67, no. 5, pp. 1897–1907, 1996.

I. Trounce, S. Neill, and D. C. Wallace, “Cytoplasmic transfer of the mtDNA nt 8993 T → G (ATP6) point mutation associated with Leigh syndrome into mtDNA-less cells demonstrates cosegregation with a decrease in state III respiration and ADP/O ratio,” Proceedings of the National Academy of Sciences of the United States of America, vol. 91, no. 18, pp. 8334–8338, 1994.

A. Gómez-Durán, D. Pacheu-Grau, E. López-Gallardo et al., “Unmasking the causes of multifactorial disorders: OXPHOS differences between mitochondrial haplogroups,” Human Molecular Genetics, vol. 19, no. 17, pp. 3343–3353, 2010.

R. Pello, M. A. Martín, V. Carelli et al., “Mitochondrial DNA background modulates the assembly kinetics of OXPHOS complexes in a cellular model of mitochondrial disease,” Human Molecular Genetics, vol. 17, no. 24, pp. 4001–4011, 2008.