Next Article in Journal
Influence of Cardiovascular Risk Factors, Comorbidities, Medication Use and Procedural Variables on Remote Ischemic Conditioning Efficacy in Patients with ST-Segment Elevation Myocardial Infarction
Next Article in Special Issue
Deletion of OGG1 Results in a Differential Signature of Oxidized Purine Base Damage in mtDNA Regions
Previous Article in Journal
LMM24 Encodes Receptor-Like Cytoplasmic Kinase 109, Which Regulates Cell Death and Defense Responses in Rice
Previous Article in Special Issue
Genetic Defects in Mitochondrial Dynamics in Caenorhabditis elegans Impact Ultraviolet C Radiation- and 6-hydroxydopamine-Induced Neurodegeneration
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Analysis of Mitochondrial DNA Polymorphisms in the Human Cell Lines HepaRG and SJCRH30

by
Matthew J. Young
1,*,
Anitha D. Jayaprakash
2 and
Carolyn K. J. Young
1
1
Department of Biochemistry and Molecular Biology, Southern Illinois University School of Medicine, Carbondale, Illinois 62901, USA
2
Girihlet Inc., Oakland, CA 94609, USA
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(13), 3245; https://doi.org/10.3390/ijms20133245
Submission received: 31 May 2019 / Revised: 26 June 2019 / Accepted: 1 July 2019 / Published: 2 July 2019
(This article belongs to the Special Issue mtDNA and Mitochondrial Stress Signaling in Human Diseases)

Abstract

:
The mitochondrial DNA (mtDNA) sequences of two commonly used human cell lines, HepaRG and SJCRH30, were determined. HepaRG originates from a liver tumor obtained from a patient with hepatocarcinoma and hepatitis C while SJCRH30 originates from a rhabdomyosarcoma patient tumor. In comparison to the revised Cambridge Reference Sequence, HepaRG and SJCRH30 mtDNA each contain 14 nucleotide variations. In addition to an insertion of a cytosine at position 315 (315insC), the mtDNA sequences from both cell types share six common polymorphisms. Heteroplasmic variants were identified in both cell types and included the identification of the 315insC mtDNA variant at 42 and 75% heteroplasmy in HepaRG and SJCRH30, respectively. Additionally, a novel heteroplasmic G13633A substitution in the HepaRG ND5 gene was detected at 33%. Previously reported cancer-associated mtDNA variants T195C and T16519C were identified in SJCRH30, both at homoplasmy (100%), while HepaRG mtDNA harbors a known prostate cancer-associated T6253C substitution at near homoplasmy, 95%. Based on our sequencing analysis, HepaRG mtDNA is predicted to lie within haplogroup branch H15a1 while SJCRH30 mtDNA is predicted to localize to H27c. The catalog of polymorphisms and heteroplasmy reported here should prove useful for future investigations of mtDNA maintenance in HepaRG and SJCRH30 cell lines.

1. Introduction

More than 1000 human mitochondrial proteins are encoded by the nuclear genome and must be imported into mitochondria following translation on cytoplasmic ribosomes [1]. The mitochondrial DNA genome (mtDNA) is an ~16.6 kilobase pair (kbp) covalently closed circular molecule that contains 13 genes for polypeptides, 2 genes for rRNAs, and 22 genes for tRNAs [2,3,4]. Our maternally inherited mtDNA is critical to cellular viability as exemplified by the numerous disease mutations associated with it and by observations that knocking out mtDNA maintenance genes results in embryonic lethality in various mouse models [5,6,7]. Maintenance of the mitochondrial genome is also required to avoid apoptosis induced by mtDNA damage [8,9]. A single human cell can contain several thousand copies of mtDNA that are distributed within hundreds of individual mitochondria or throughout an elaborate mitochondrial reticular network [10,11,12]. Homoplasmy is the existence of clonal copies of mtDNA within a cell; however, mtDNA can also exist in heteroplasmy or a mixed population of different mtDNA genotypes within a cell [13,14]. Mitochondrial fusion and fission allow for complementation of heteroplasmic mtDNA genes in trans (e.g., a wild-type gene product complementing a mutant gene product); however, as the proportion of wild-type to mutant mtDNA decreases, the ability of a cell to produce energy via oxidative phosphorylation (OXPHOS) can decline [15]. Interestingly, enrichment of mtDNA heteroplasmic variants have been reported in certain cancers [16,17], and heteroplasmic variants have been identified in human peripheral blood mononuclear cells, 501T fibroblast cell lines, and cancer-derived cell lines [14,18].
Human cancer-derived cell lines are useful preclinical tools for studies of mitochondrial metabolism and cancer cell biology, as well as high-throughput drug toxicity and drug discovery experiments. SJCRH30 is a cell line that was derived from the tumor of a 17-year-old white male with rhabdomyosarcoma. SJCRH30 cells harbor attenuated sarcomere structures resembling those found in primitive rhabdomyoblasts [19]. SJCRH30 has been used to evaluate the cytotoxicity of chemotherapeutic drugs, such as cisplatin, doxorubicin, topotecan, and others [20,21]. Additionally, SJCRH30 has been used as a model of human myoblasts to study the regulation of mitochondrial biogenesis and cellular oxygen consumption rates [22,23].
HepaRG was originally derived from a liver tumor obtained from a female patient suffering from hepatitis C infection and hepatocarcinoma [24]. The HepaRG cell line expresses various mature hepatocyte-like functions, including cytochrome P450s associated with xenobiotic metabolism [25]. Importantly, HepaRG displays sensitivity to hepatotoxic compounds such as acetaminophen and aflatoxin B1 [25,26]. Recent studies support that HepaRG is a suitable model to test drug-induced mitochondrial toxicity and mtDNA homeostasis [3,26,27,28]. Here we report the mtDNA genome polymorphisms and heteroplasmy of SJCRH30 and HepaRG.

2. Results

2.1. HepaRG mtDNA Nucleotide Variants and Heteroplasmy

The HepaRG mtDNA genome contains a single C at positions 3106–3107, which is in agreement with the revised Cambridge Reference Sequence, rCRS [29]. In contrast to the rCRS, HepaRG mtDNA contains 14 nucleotide variations, Table 1. With exception to the A16T transversion substitution and the insertion of C at position 315 (315insC), the remaining variants are transitions with A to G, G to A, T to C, or C to T changes. The MITOMASTER human mtDNA sequence analysis tool predicts HepaRG mtDNA localizes within haplogroup branch H15a1 [30]. H15a1 is of European origin [31].
Of the 14 variants, six have been previously reported as polymorphisms, A263G (non-coding), A750G (RNR1, 12 S ribosomal RNA), A1438G (RNR1), A4769G (synonymous substitution in the NADH dehydrogenase subunit 2 gene, ND2), A8860G (Thr112Ala missense in ATP synthase F0 subunit 6 gene, ATP6), and A15326G/Thr194Ala missense in the cytochrome b gene, CYB [18]. Additionally, the previously reported prostate cancer-associated T6253C missense substitution was detected in the HepaRG mtDNA genome [32]. The T6253C/Met117Thr substitution alters an evolutionarily conserved amino acid residue in cytochrome c oxidase subunit 1 (COX1) that could affect mitochondrial function as predicted by MutationAssessor [33]. MutationAssessor predicts the potential deleterious impact of a DNA mutation that changes a protein’s amino acid residue. To make this prediction, multiple sequence alignments of homologous proteins are grouped into families and subfamilies. Using the conservation pattern information generated from the alignments, MutationAssessor generates a functional impact (FI) score to rate a mutation’s impact as either high, medium, low, or neutral. The COX1 Met117Thr variant has an FI score of 2.62 and is predicted to have a medium impact on protein function.
Three of the identified variants are not frequently associated with branch H15a1, A16T, 315insC, and G13633A. The non-coding A16T substitution is listed in the Single Nucleotide Polymorphism Database (dbSNP) as a single nucleotide variation (refSNP id: rs1556422363) and occurs at near homoplasmic levels, 94%. The 315insC insertion is located in the mtDNA control region and was detected at 42% heteroplasmy. To the best of our knowledge, the 33% heteroplasmic non-synonymous G13633A/Gly433Ser substitution localizing to the NADH dehydrogenase subunit 5 (ND5) gene has not been identified to date. The ND5 G13633A/Gly433Ser mutation has an FI score of 2.55 and is predicted to have a medium impact on protein function. Therefore, the Gly433Ser variant could negatively affect mitochondrial function. Maintenance of HepaRG mtDNA heteroplasmy during growth in tissue culture was confirmed by sequencing at passages 11 (41% 315insC and 37% G13633A) and 16 (42% 315insC and 33% G13633A, Table 1). With the exception of the 315insC and the G13633A heteroplasmic variants, the remaining substitutions exist at ≥88% within the population of mtDNA molecules and are therefore near homoplasmic. The remaining mtDNA variants include two silent substitutions, C14953T and T11410C, as well as two non-coding substitutions, T55C and T57C.

2.2. SJCRH30 mtDNA Nucleotide Variants and Heteroplasmy

Similar to the rCRS and to the HepaRG mtDNA sequences, the SJCRH30 mtDNA genome contains a single C at positions 3106–3107. SJCRH30 mtDNA contains 14 nucleotide variations relative to the rCRS, and MITOMASTER predicts localization within haplogroup branch H27c, Table 2.
The 315insC variant found in HepaRG mtDNA is also present in the SJCRH30 mitochondrial genome at 75% heteroplasmy. With the exception of 315insC, the remaining variants are transitions with A to G, G to A, or T to C changes. Six of the previously mentioned polymorphisms identified in HepaRG were detected in SJCRH30 mtDNA (A263G, A750G, A1438G, A4769G, A8860G, and A15326G). The common G11719A variant was also identified in SJCRH30 mtDNA [18]. Additionally, three of the nucleotide changes detected do not conventionally occur in H27c, namely T195C (non-coding), 315insC (non-coding), and the missense T14634C/Met14Val NADH dehydrogenase subunit 6 (ND6) variant. The ND6 Met14Val substitution has an FI score of −0.075 and is predicted to be neutral via MutationAssessor analysis. The mtDNA T195C is a previously reported European Caucasian melanoma-associated substitution while T16519C has been reported to increase a woman’s risk of developing breast cancer or to be in linkage disequilibrium with a functional SNP that increases a woman’s risk [36,38]. With the exception of the heteroplasmic 315insC variant, the remaining mtDNA substitutions were homoplasmic or near homoplasmic. Three of the near homoplasmic mtDNA variants include the synonymous T4838C (91%) and the non-coding A16316G (90%) and G16129A (95%) variants.

3. Discussion

Reduced mtDNA copy number and the presence of mtDNA mutations that alter OXPHOS have been reported to be common in cancer; however, functional mitochondria and mtDNA are necessary for cancer cell growth and tumorigenesis [15,17,39,40]. Human cancer cells harbor both homoplasmic and heteroplasmic mtDNA mutations, and cell culture experiments utilizing cancer cell lines have demonstrated that heteroplasmy can be stably maintained after many passages [14,16]. For these reasons, mtDNA polymorphisms and heteroplasmy were evaluated in two commonly used cancer-derived cell lines, HepaRG and SJCRH30. Previously reported cancer-associated mtDNA variants were found in both cell lines. HepaRG harbors the prostate cancer-associated T6253C substitution at near homoplasmy (95%) while SJCRH30 contains both the melanoma-associated T195C and the T16519C breast cancer-associated variants at homoplasmy [32,36,38]. According to the MITOMASTER database, both T16519C (SJCRH30, H27c) and T6253C (HepaRG, H15a1) occur at 100% in their respective haplogroup branches. To date, the T195C variant is not associated with the two representative H27c sequences found in the MITOMASTER database, Table 2.
In HepaRG and SJCRH30 mtDNA, six previously reported polymorphisms were found to be homoplasmic (fixed) or nearly homoplasmic, occurring at ≥88% within the population of mtDNA molecules (Table 1 and Table 2). These polymorphisms appear at frequencies of ≥96.8% and of 100% in the H15a1 and H27c branches, respectively. The HepaRG mtDNA variants are maintained at near homoplasmy and occur at ≥88%, with the exception of 315insC (42% heteroplasmy) and G13633A (33% heteroplasmy), Table 1. Similarly, SJCRH30 mtDNA variants occur at ≥90% with the exception of the 315insC variant (75% heteroplasmy), Table 2. The 315insC insertion is found in the mtDNA control region while the G13633A substitution alters an evolutionarily conserved codon that may affect mitochondrial function. We hypothesize that the 315insC and the G13633A variants are under negative selection and that wild-type mtDNAs are being sustained to maintain functional mitochondria. In mitochondrial genomes from both cell types, many of the substitutions were near homoplasmy, ranging from 88% to 100% occurrence within the population of mtDNAs. We predict that the presence of near homoplasmy variants, e.g., synonymous substitutions, may expand to 100% homoplasmy during further passaging in cell culture. Moreover, a non-coding substitution or missense variant that does not negatively impact OXPHOS could expand to 100% homoplasmy. As the majority of the mtDNA nucleotide variants identified in both HepaRG and SJCRH30 are transitions, these substitutions likely arose from erroneous incorporation by the replicative mtDNA polymerase gamma or from deamination events as previously proposed [41,42,43]. The SJCRH30 T195C transition substitution is atypical for haplogroup branch H27c and localizes to the mtDNA control region heavy-strand origin of replication, which is important for initiation of mtDNA replication. T195C may have been subject to selection in melanoma and rhabdomyosarcoma to modify mtDNA replication and, by extension, mitochondrial metabolism. Numerous modifications in mitochondrial function have been implicated in cancer biology including shifting energy production, disrupting apoptosis signaling, increasing mutation of mtDNA, and altering antioxidant activity and reactive oxygen species production [39]. Interestingly, T195C was found to be significantly higher in melanoma patients in comparison to unrelated control individuals [36]. The HepaRG A16T transversion substitution lies within the non-coding 7S DNA control region and is not typical of the H15a1 haplogroup branch. One hypothesis is that 7S DNA can be utilized as a primer to initiate heavy-strand replication [44]. Perhaps the A16T 7S DNA modification was subject to selection in hepatocarcinoma to modify mtDNA replication initiation in cancer cells. Future next-generation sequencing studies of mtDNA from tumors and unaffected control tissues could determine whether T195C and A16T are rhabdomyosarcoma- and hepatocarcinoma-specific changes, respectively.

4. Materials and Methods

4.1. Cell Culture

Proliferating HepaRG (Biopredic International Saint-Grégoire, France) were grown until passage 16, as previously described [24,26]. SJCRH30 cells (RC13, RMS 13, SJRH30, ATCC® CRL­2061TM) were grown until passage 5 according to ATCC recommendations. Human cell lines were (1) cultured for the recommended number of passages to avoid marked phenotypic and morphological changes and not more than 20 passages if recommendations were not available, (2) screened for correct morphology and bacterial contamination via microscopy on days of feeding and passaging to ensure correct growth rates and appearances of cells, and (3) routinely screened for potential mycoplasma infection using in-house comparative PCR analysis with primers and controls kindly provided by Dr. Uphoff [45]. Dulbecco’s phosphate-buffered saline (DPBS)-washed cell pellets were frozen at −80 °C and then processed as described below.

4.2. mtDNA Next-Generation Sequencing and Data Analysis

The Mseek method of sample processing and deep sequencing of mtDNA, as well as the procedure for data analysis, was conducted as previously reported [14]. Briefly, Mseek consists of (i) isolating total DNA from a thawed DPBS-washed cell pellet, (ii) digesting linear nuclear DNA (nDNA) with Exonuclease V, (iii) purifying the products using Ampure beads to remove short fragments, (iv) testing the results of the digestion with PCR primers specific for mtDNA and nDNA using 1 μl of the digested sample, (v) fragmenting the remaining sample using Covaris and end-repair, (vi) ligating barcoded adapters compatible with the sequencing platform to the fragments, (vii) amplifying the library utilizing universal adapters, and (viii) loading samples onto the Illumina NextSeq 500 platform.
The sequencing data were generated as fastq files, as previously described [14]. Briefly, the sequences were filtered for quality (sequences with >10 consecutive nucleotides with Q < 20 were eliminated) and mapped to the revised Cambridge Reference Sequence (rCRS), accession NC_012920. Identical reads were identified as being clonal and were considered only once, irrespective of the number of copies, toward variant calling. A variant call was made only if there were at least three non-clonal reads carrying the variant, and a minimum coverage of 10 was required at the variant. Variants occurring on reads on one strand (with a skew greater than 0.1 or 10%) of the mtDNA were excluded to further reduce errors. The error rate in NextSeq reads are usually <1 in 1000 (phred score Q > 30) and requiring at least three non-clonal reads reduces the error rate to well under one in a million. Nuclear contamination was estimated using sequences that map to repeat elements such as long interspersed nuclear elements (LINEs) and short interspersed nuclear elements (SINEs), which only occur in nDNA. This enables a reliable estimation of the level of nDNA contamination. The fastq files have been submitted to the NCBI Sequence Read Archive (SRA) BioProject ID: PRJNA545541.

5. Conclusions

We determined the mtDNA sequence polymorphisms and heteroplasmy of two commonly used cell lines derived from patients with hepatocarcinoma (HepaRG) and rhabdomyosarcoma (SJCRH30). We expect that the catalog of polymorphisms and heteroplasmy reported here will prove useful for future investigations of mtDNA maintenance utilizing these cell lines.

Author Contributions

Conceptualization, M.J.Y.; Formal analysis, M.J.Y. and A.D.J.; Funding acquisition, M.J.Y; Methodology, A.D.J., M.J.Y., and C.K.J.Y.; Project administration, M.J.Y. and C.K.J.Y.; Resources, M.J.Y.; Supervision, M.J.Y.; Writing—original draft, M.J.Y.; Writing—review and editing, M.J.Y., A.D.J., and C.K.J.Y.

Funding

This research was supported by a National Institute of Environmental Health Sciences Pathway to Independence Award to M.J.Y. (5R00ES022638-04) and by start-up funds from SIU School of Medicine to M.J.Y.

Acknowledgments

The authors acknowledge the work of Christiane Guguen-Guillouzo, Philippe Gripon, and Christian Trepo for the isolation and the characterization of HepaRG. We also thank Joel Wheeler and Israa Hassan for technical assistance.

Conflicts of Interest

A.D.J. is a co-founder of Girihlet Inc. The other authors have no financial conflicts of interest.

Abbreviations

mtDNAmitochondrial DNA
rCRSrevised Cambridge Reference Sequence
OXPHOSoxidative phosphorylation
kbpkilobase pair
dbSNPThe Single Nucleotide Polymorphism Database
ND5NADH dehydrogenase subunit 5 gene

References

  1. Calvo, S.E.; Clauser, K.R.; Mootha, V.K. MitoCarta2.0: An updated inventory of mammalian mitochondrial proteins. Nucleic Acids Res. 2016, 44, D1251–D1257. [Google Scholar] [CrossRef] [PubMed]
  2. Young, M.J.; Copeland, W.C. Human mitochondrial DNA replication machinery and disease. Curr. Opin. Genet. Dev. 2016, 38, 52–62. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Wheeler, J.H.; Young, C.K.J.; Young, M.J. Analysis of Human Mitochondrial DNA Content by Southern Blotting and Nonradioactive Probe Hybridization. Curr. Protoc Toxicol. 2019, 80, e75. [Google Scholar] [CrossRef] [PubMed]
  4. Wallace, K.B. Drug-induced mitochondrial neuropathy in children: A conceptual framework for critical windows of development. J. Child Neurol. 2014, 29, 1241–1248. [Google Scholar] [CrossRef] [PubMed]
  5. Park, C.B.; Larsson, N.G. Mitochondrial DNA mutations in disease and aging. J. Cell Biol. 2011, 193, 809–818. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Humble, M.M.; Young, M.J.; Foley, J.F.; Pandiri, A.R.; Travlos, G.S.; Copeland, W.C. Polg2 is essential for mammalian embryogenesis and is required for mtDNA maintenance. Hum. Mol. Genet. 2013, 22, 1017–1025. [Google Scholar] [CrossRef] [PubMed]
  7. Wallace, D.C.; Chalkia, D. Mitochondrial DNA genetics and the heteroplasmy conundrum in evolution and disease. Cold Spring Harb. Perspect Biol. 2013, 5, a021220. [Google Scholar] [CrossRef] [PubMed]
  8. Tann, A.W.; Boldogh, I.; Meiss, G.; Qian, W.; Van Houten, B.; Mitra, S.; Szczesny, B. Apoptosis induced by persistent single-strand breaks in mitochondrial genome: critical role of EXOG (5’-EXO/endonuclease) in their repair. J. Biol. Chem. 2011, 286, 31975–31983. [Google Scholar] [CrossRef] [PubMed]
  9. Santos, J.H.; Hunakova, L.; Chen, Y.; Bortner, C.; Van Houten, B. Cell sorting experiments link persistent mitochondrial DNA damage with loss of mitochondrial membrane potential and apoptotic cell death. J. Biol. Chem. 2003, 278, 1728–1734. [Google Scholar] [CrossRef]
  10. Young, M.J.; Humble, M.M.; DeBalsi, K.L.; Sun, K.Y.; Copeland, W.C. POLG2 disease variants: Analyses reveal a dominant negative heterodimer, altered mitochondrial localization and impaired respiratory capacity. Hum. Mol. Genet. 2015, 24, 5184–5197. [Google Scholar] [CrossRef] [PubMed]
  11. Young, M.J. Off-Target Effects of Drugs that Disrupt Human Mitochondrial DNA Maintenance. Front Mol. Biosci. 2017, 4, 74. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Spelbrink, J.N. Functional organization of mammalian mitochondrial DNA in nucleoids: history, recent developments, and future challenges. Iubmb Life 2010, 62, 19–32. [Google Scholar] [CrossRef] [PubMed]
  13. Freyer, C.; Cree, L.M.; Mourier, A.; Stewart, J.B.; Koolmeister, C.; Milenkovic, D.; Wai, T.; Floros, V.I.; Hagstrom, E.; Chatzidaki, E.E.; et al. Variation in germline mtDNA heteroplasmy is determined prenatally but modified during subsequent transmission. Nat. Genet. 2012, 44, 1282–1285. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Jayaprakash, A.D.; Benson, E.K.; Gone, S.; Liang, R.; Shim, J.; Lambertini, L.; Toloue, M.M.; Wigler, M.; Aaronson, S.A.; Sachidanandam, R. Stable heteroplasmy at the single-cell level is facilitated by intercellular exchange of mtDNA. Nucleic Acids Res. 2015, 43, 2177–2187. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Wallace, D.C. Mitochondria and cancer. Nat. Rev. Cancer 2012, 12, 685–698. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. He, Y.; Wu, J.; Dressman, D.C.; Iacobuzio-Donahue, C.; Markowitz, S.D.; Velculescu, V.E.; Diaz, L.A., Jr.; Kinzler, K.W.; Vogelstein, B.; Papadopoulos, N. Heteroplasmic mitochondrial DNA mutations in normal and tumour cells. Nature 2010, 464, 610–614. [Google Scholar] [CrossRef] [Green Version]
  17. Larman, T.C.; DePalma, S.R.; Hadjipanayis, A.G.; The Cancer Genome Atlas Research Network; Protopopov, A.; Zhang, J.; Gabriel, S.B.; Chin, L.; Seidman, C.E.; Kucherlapati, R.; et al. Spectrum of somatic mitochondrial mutations in five cancers. Proc. Natl. Acad. Sci. USA 2012, 109, 14087–14091. [Google Scholar] [CrossRef] [Green Version]
  18. Herrnstadt, C.; Preston, G.; Andrews, R.; Chinnery, P.; Lightowlers, R.N.; Turnbull, D.M.; Kubacka, I.; Howell, N. A high frequency of mtDNA polymorphisms in HeLa cell sublines. Mutat. Res. 2002, 501, 19–28. [Google Scholar] [CrossRef]
  19. Roberts, W.M.; Douglass, E.C.; Peiper, S.C.; Houghton, P.J.; Look, A.T. Amplification of the gli gene in childhood sarcomas. Cancer Res. 1989, 49, 5407–5413. [Google Scholar]
  20. Ohnstad, H.O.; Paulsen, E.B.; Noordhuis, P.; Berg, M.; Lothe, R.A.; Vassilev, L.T.; Myklebost, O. MDM2 antagonist Nutlin-3a potentiates antitumour activity of cytotoxic drugs in sarcoma cell lines. BMC Cancer 2011, 11, 211. [Google Scholar] [CrossRef]
  21. Barretina, J.; Caponigro, G.; Stransky, N.; Venkatesan, K.; Margolin, A.A.; Kim, S.; Wilson, C.J.; Lehar, J.; Kryukov, G.V.; Sonkin, D.; et al. The Cancer Cell Line Encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature 2012, 483, 603–607. [Google Scholar] [CrossRef] [PubMed]
  22. Morino, K.; Petersen, K.F.; Sono, S.; Choi, C.S.; Samuel, V.T.; Lin, A.; Gallo, A.; Zhao, H.; Kashiwagi, A.; Goldberg, I.J.; et al. Regulation of mitochondrial biogenesis by lipoprotein lipase in muscle of insulin-resistant offspring of parents with type 2 diabetes. Diabetes 2012, 61, 877–887. [Google Scholar] [CrossRef] [PubMed]
  23. Civitarese, A.E.; MacLean, P.S.; Carling, S.; Kerr-Bayles, L.; McMillan, R.P.; Pierce, A.; Becker, T.C.; Moro, C.; Finlayson, J.; Lefort, N.; et al. Regulation of skeletal muscle oxidative capacity and insulin signaling by the mitochondrial rhomboid protease PARL. Cell Metab. 2010, 11, 412–426. [Google Scholar] [CrossRef] [PubMed]
  24. Gripon, P.; Rumin, S.; Urban, S.; Le Seyec, J.; Glaise, D.; Cannie, I.; Guyomard, C.; Lucas, J.; Trepo, C.; Guguen-Guillouzo, C. Infection of a human hepatoma cell line by hepatitis B virus. Proc. Natl. Acad. Sci. USA 2002, 99, 15655–15660. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Aninat, C.; Piton, A.; Glaise, D.; Le Charpentier, T.; Langouët, S.; Morel, F.; Guguen-Guillouzo, C.; Guillouzo, A. Expression of cytochromes P450, conjugating enzymes and nuclear receptors in human hepatoma HepaRG cells. Drug Metab. Dispos. 2006, 34, 75–83. [Google Scholar] [CrossRef] [PubMed]
  26. Young, C.K.J.; Young, M.J. Comparison of HepaRG cells following growth in proliferative and differentiated culture conditions reveals distinct bioenergetic profiles. Cell Cycle (Georget. Tex) 2019, 18, 476–499. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Kamalian, L.; Douglas, O.; Jolly, C.E.; Snoeys, J.; Simic, D.; Monshouwer, M.; Williams, D.P.; Kevin Park, B.; Chadwick, A.E. The utility of HepaRG cells for bioenergetic investigation and detection of drug-induced mitochondrial toxicity. Toxicol Vitr. 2018, 53, 136–147. [Google Scholar] [CrossRef] [PubMed]
  28. Le Guillou, D.; Bucher, S.; Begriche, K.; Hoet, D.; Lombes, A.; Labbe, G.; Fromenty, B. Drug-Induced Alterations of Mitochondrial DNA Homeostasis in Steatotic and Nonsteatotic HepaRG Cells. J. Pharm. Exp. 2018, 365, 711–726. [Google Scholar] [CrossRef]
  29. Andrews, R.M.; Kubacka, I.; Chinnery, P.F.; Lightowlers, R.N.; Turnbull, D.M.; Howell, N. Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA. Nat. Genet. 1999, 23, 147. [Google Scholar] [CrossRef]
  30. Lott, M.T.; Leipzig, J.N.; Derbeneva, O.; Xie, H.M.; Chalkia, D.; Sarmady, M.; Procaccio, V.; Wallace, D.C. mtDNA Variation and Analysis Using Mitomap and Mitomaster. Curr. Protoc. Bioinform. 2013, 44, 1.23.1–1.23.26. [Google Scholar] [CrossRef] [Green Version]
  31. Heupink, T.H.; Subramanian, S.; Wright, J.L.; Endicott, P.; Westaway, M.C.; Huynen, L.; Parson, W.; Millar, C.D.; Willerslev, E.; Lambert, D.M. Ancient mtDNA sequences from the First Australians revisited. Proc. Natl. Acad. Sci. USA 2016, 113, 6892–6897. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Petros, J.A.; Baumann, A.K.; Ruiz-Pesini, E.; Amin, M.B.; Sun, C.Q.; Hall, J.; Lim, S.; Issa, M.M.; Flanders, W.D.; Hosseini, S.H.; et al. mtDNA mutations increase tumorigenicity in prostate cancer. Proc. Natl. Acad. Sci. USA 2005, 102, 719–724. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Reva, B.; Antipin, Y.; Sander, C. Predicting the functional impact of protein mutations: Application to cancer genomics. Nucleic Acids Res. 2011, 39, e118. [Google Scholar] [CrossRef] [PubMed]
  34. Fish, J.; Raule, N.; Attardi, G. Discovery of a major D-loop replication origin reveals two modes of human mtDNA synthesis. Science 2004, 306, 2098–2101. [Google Scholar] [CrossRef] [PubMed]
  35. Collins, D.W.; Gudiseva, H.V.; Trachtman, B.; Bowman, A.S.; Sagaser, A.; Sankar, P.; Miller-Ellis, E.; Lehman, A.; Addis, V.; O’Brien, J.M. Association of primary open-angle glaucoma with mitochondrial variants and haplogroups common in African Americans. Mol. Vis. 2016, 22, 454–471. [Google Scholar] [PubMed]
  36. Ebner, S.; Lang, R.; Mueller, E.E.; Eder, W.; Oeller, M.; Moser, A.; Koller, J.; Paulweber, B.; Mayr, J.A.; Sperl, W.; et al. Mitochondrial haplogroups, control region polymorphisms and malignant melanoma: A study in middle European Caucasians. PLoS ONE 2011, 6, e27192. [Google Scholar] [CrossRef]
  37. Rollins, B.; Martin, M.V.; Sequeira, P.A.; Moon, E.A.; Morgan, L.Z.; Watson, S.J.; Schatzberg, A.; Akil, H.; Myers, R.M.; Jones, E.G.; et al. Mitochondrial variants in schizophrenia, bipolar disorder, and major depressive disorder. PLoS ONE 2009, 4, e4913. [Google Scholar] [CrossRef]
  38. Bai, R.K.; Leal, S.M.; Covarrubias, D.; Liu, A.; Wong, L.J. Mitochondrial genetic background modifies breast cancer risk. Cancer Res. 2007, 67, 4687–4694. [Google Scholar] [CrossRef]
  39. Singh, K.K.; Modica-Napolitano, J.S. Special Issue: Mitochondria in Cancer. Semin Cancer Biol. 2017, 47, iv–vi. [Google Scholar] [CrossRef]
  40. Singh, B.; Modica-Napolitano, J.S.; Singh, K.K. Defining the momiome: Promiscuous information transfer by mobile mitochondria and the mitochondrial genome. Semin Cancer Biol. 2017, 47, 1–17. [Google Scholar] [CrossRef]
  41. Schmitt, M.W.; Kennedy, S.R.; Salk, J.J.; Fox, E.J.; Hiatt, J.B.; Loeb, L.A. Detection of ultra-rare mutations by next-generation sequencing. Proc. Natl. Acad. Sci. USA 2012, 109, 14508–14513. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Song, S.; Pursell, Z.F.; Copeland, W.C.; Longley, M.J.; Kunkel, T.A.; Mathews, C.K. DNA percursor asymmetries in mammalian tissue mitochondrial and possible contribution to mitochondrial mutagenesis through reduced replication fidleity. Proc. Natl. Acad. Sci. USA 2005, 102, 4990–4995. [Google Scholar] [CrossRef] [PubMed]
  43. Modrich, P. DNA mismatch correction. Ann. Rev. Biochem. 1987, 56, 435–466. [Google Scholar] [CrossRef] [PubMed]
  44. Nicholls, T.J.; Minczuk, M. In D-loop: 40 years of mitochondrial 7S DNA. Exp. Gerontol. 2014, 56, 175–181. [Google Scholar] [CrossRef] [PubMed]
  45. Uphoff, C.C.; Drexler, H.G. Comparative PCR analysis for detection of mycoplasma infections in continuous cell lines. Vitr. Cell Dev. Biol. Anim 2002, 38, 79–85. [Google Scholar] [CrossRef]
Table 1. Sequence changes and heteroplasmy identified in HepaRG mitochondrial DNA (mtDNA).
Table 1. Sequence changes and heteroplasmy identified in HepaRG mitochondrial DNA (mtDNA).
MtDNA nt. Change aLocation bCoverage c%Variant d%Freq in H15a1 eRemarks f
A16TCR: 7S37,783940SNV
T55CCR: 7S57,8459080.65SNV
T57CCR: 7S, HVS2, OH5757,8479061.29SNV
A263GCR: HVS2, OH4,8799896.77SNV g
315insCCR:HVS2, OH5,392424.84Insertion of C
A750GRNR148,18088100SNV g
A1438GRNR170,03395100Benign g
A4769GND2/M10014,2239096.77Synonymous variant, ATA > ATG g
T6253CCOX1/M117T58,70495100Missense variant, ATA > ACA; prostate cancer associated h
A8860GATP6/T112A39,6149596.77Missense variant, ACA > GCA g
T11410CND4/P21728,04892100Synonymous variant, CCT > CCC
G13633AND5/G433S43,013330Missense variant, GGT > AGT
C14953TCYB/I6961,4118998.39Synonymous variant, ATC > ATT
A15326GCYB/T194A62,03893100Missense variant, ACA > GCA g
a Nucleotide (nt.) positions are numbered according to the revised Cambridge Reference Sequence (rCRS) light-strand, NC_012920.1. The rCRS was derived from a single person from haplogroup H2a2, see www.mitomap.org for details. The rCRS nucleotides are listed on the left-hand side of the nucleotide positions while detected changes (variants) are listed on the right-hand side. b The non-coding control region (CR) that contains the 7S DNA (7S), the hypervariable segment 2 (HVS2), the heavy-strand origins of replication (OH), and in HeLa, A549, 143B, and TK cancer-derived cell lines and immortalized lymphocytes the major replication origin at position 57, OH57 [34], https://www.mitomap.org/foswiki/bin/view/MITOMAP/GenomeLoci. c Coverage at the indicated position. d The percentage of the variant in the plus strand represents heteroplasmy; 100% represents homoplasmy. e MITOMASTER predicts that HepaRG mtDNA lies within the haplogroup branch H15a1 of European origin [30,31]; the MITOMASTER-predicted frequency of each variant within H15a is indicated and is based on a total of 62 H15a sequences. f SNV, single nucleotide variation as listed in NCBI dbSNP short genetic variations, https://www.ncbi.nlm.nih.gov/snp/; benign as listed in ClinVar, https://www.ncbi.nlm.nih.gov/clinvar/; synonymous and missense variant codon changes are shown. g These polymorphisms occur in most mtDNA genomes except for a small subcluster of haplogroup H that includes the rCRS [18]. h T6253C has previously been found to be associated with prostate cancer [32] and with primary open-angle glaucoma [35].
Table 2. Sequence changes and heteroplasmy identified in SJCRH30 mtDNA.
Table 2. Sequence changes and heteroplasmy identified in SJCRH30 mtDNA.
MtDNA nt. Change aLocation bCoverage c%Variant d%Freq in H27c eRemarks f
T195CCR: HVS2, OH1271000Melanoma-associated g
A263GCR: HVS2, OH60100100SNV h
315insCCR: HVS2, OH24750Insertion of C
A750GRNR1117100100SNV h
A1438GRNR1193100100Benign h
A4769GND2/M10079100100Synonymous variant, ATA > ATG h
T4838CND2/P12314191100Synonymous variant, CCT > CCC
A8860GATP6/T112A149100100Missense variant, ACA > GCA h
G11719AND4/G32012197100SNV, synonymous variant, GGG > GGA
T14634CND6/M14V136900Missense variant, ATG > GTG i
A15326GCYB/T194A12995100Missense variant, ACA > GCA h
G16129ACR: HVS1, TAS2, 7S15495100SNV
A16316GCR: HVS1, 7S12090100SNV
T16519CCR: 7S162100100Associated with breast cancer risk j
a Nucleotide (nt.) changes are labeled as described in the footnote for Table 1; b The non-coding control region (CR) that contains the 7S DNA (7S), the hypervariable segments 1 and 2 (HVS1 and HVS2), the heavy-strand origins of replication (OH), and the extended termination-associated sequence (TAS2), https://www.mitomap.org/foswiki/bin/view/MITOMAP/GenomeLoci. c Coverage at the indicated position. d The percentage of the variant in the plus strand represents heteroplasmy; 100% represents homoplasmy. e MITOMASTER predicts that SJCRH30 mtDNA lies within the haplogroup branch H27c; the MITOMASTER-predicted frequency of each variant within H27c is indicated and is based on a total of two H27c sequences. f SNV, benign, and codon changes are as described in the footnote for Table 1; g T195C has been found to be associated with melanoma in European Caucasians [36] and bipolar disorder [37]. h These polymorphisms occur in most mtDNA genomes except for a small subcluster of haplogroup H that includes the rCRS [18]. i The heavy-strand is the ND6 coding strand (the RNA is transcribed from, and could hybridize to, the light-strand). j T16519C was reported to increase a woman’s risk of developing breast cancer or is in linkage disequilibrium with a functional SNP that increases a woman’s risk [38].

Share and Cite

MDPI and ACS Style

Young, M.J.; Jayaprakash, A.D.; Young, C.K.J. Analysis of Mitochondrial DNA Polymorphisms in the Human Cell Lines HepaRG and SJCRH30. Int. J. Mol. Sci. 2019, 20, 3245. https://doi.org/10.3390/ijms20133245

AMA Style

Young MJ, Jayaprakash AD, Young CKJ. Analysis of Mitochondrial DNA Polymorphisms in the Human Cell Lines HepaRG and SJCRH30. International Journal of Molecular Sciences. 2019; 20(13):3245. https://doi.org/10.3390/ijms20133245

Chicago/Turabian Style

Young, Matthew J., Anitha D. Jayaprakash, and Carolyn K. J. Young. 2019. "Analysis of Mitochondrial DNA Polymorphisms in the Human Cell Lines HepaRG and SJCRH30" International Journal of Molecular Sciences 20, no. 13: 3245. https://doi.org/10.3390/ijms20133245

APA Style

Young, M. J., Jayaprakash, A. D., & Young, C. K. J. (2019). Analysis of Mitochondrial DNA Polymorphisms in the Human Cell Lines HepaRG and SJCRH30. International Journal of Molecular Sciences, 20(13), 3245. https://doi.org/10.3390/ijms20133245

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop