- Research article
- Open Access
- Open Peer Review
Mitochondrial DNA mutations in oxyphilic and chief cell parathyroid adenomas
BMC Endocrine Disordersvolume 7, Article number: 8 (2007)
The potential pathogenetic significance of mitochondrial DNA (mtDNA) mutations in tumorigenesis is controversial. We hypothesized that benign tumorigenesis of a slowly replicating tissue like the human parathyroid might constitute an especially fertile ground on which a selective advantage conferred by mtDNA mutation could be manifested and might contribute to the oxyphilic phenotype observed in a subset of parathyroid tumors.
We sought acquired mitochondrial DNA mutations by sequencing the entire 16.6 kb mitochondrial genome of each of thirty sporadic parathyroid adenomas (18 chief cell and 12 oxyphil cell), eight independent, polyclonal, parathyroid primary chief cell hyperplasias plus corresponding normal control samples, five normal parathyroid glands, and one normal thyroid gland.
Twenty-seven somatic mutations were identified in 15 of 30 (9 of 12 oxyphil adenomas, 6 of 18 chief cell) parathyroid adenomas studied. No somatic mutations were observed in the hyperplastic parathyroid glands.
Features of the somatic mutations suggest that they may confer a selective advantage and contribute to the molecular pathogenesis of parathyroid adenomas. Importantly, the statistically significant differences in mutation prevalence in oxyphil vs. chief cell adenomas also suggest that mtDNA mutations may contribute to the oxyphil phenotype.
Single parathyroid adenomas, which cause approximately 85% of the cases of primary hyperparathyroidism  are well-differentiated, benign, clonal tumors, which produce hypercalcemia through excessive secretion of parathyroid hormone (PTH). Activation of the cyclin D1 oncogene and inactivation of the MEN1 tumor suppressor gene are established pathogenetic contributors, but understanding of molecular pathogenesis in this disease remains incomplete. The mitochondrial genome has been identified as a possible target for somatic mutations that may promote tumorigenesis [2–10]. Mitochondrial abnormalities, including changes in structure, number, and respiratory enzyme components and transport systems, have been observed in many cancers. Homoplasmic somatic mutations of the mitochondrial genome have also been identified in many tumor types [2–10]. However, the pathogenetic significance of the mutations detected in those tissues is unclear.
The human parathyroid is an intrinsically low-turnover tissue , as such, benign tumors arising from this tissue have relatively low proliferative rates. The mitochondrial genome is an especially attractive potential target for mutations that could drive tumorigenesis in such a tissue. Mitochondrial DNA replicates frequently and independently of the nuclear genome , and over many years the potential for enhanced DNA damage from locally generated reactive oxygen species (ROS) may further enhance the accumulation of somatic mutations in the mitochondrial genome.
In many endocrine organs, especially the parathyroid glands, with age, increasing numbers of oxyphil cells appear [13, 14]. The oxyphil cells of normal parathyroid glands have a characteristic eosinophilic granular cytoplasm (Figure 1) which on electron microscopic examination is filled with abnormally shaped mitochondria [15, 16]. In contrast to the chief cells, the oxyphil cells in normal glands lack the organelles associated with protein synthesis and secretion, and lack immunohistochemically identifiable chromogranin and PTH . The majority of parathyroid adenomas are composed predominately of chief cells (Figure 1), however there is a subset of hyperfunctioning parathyroid adenomas that contain greater than 90% oxyphil cells and are classified as oxyphil adenomas These tumors, in addition to having a cytoplasm filled with abnormal mitochondria, contain the organelles of protein synthesis and secretion and contain histochemically identifiable PTH. While the molecular basis for this phenotypic change is unknown, it has been hypothesized that mitochondria may accumulate in oxyphil cells due to abnormalities of the mitochondrial genome, supported by the finding of increased numbers of oxyphil cells in normal glands with increasing age, analogous to increased mtDNA damage that occurs with aging .
We therefore sought acquired somatic mutations of the mitochondrial genome in both typical chief cell parathyroid adenomas and oxyphil cell parathyroid adenomas using an especially rigorous approach, and compared these findings with the mitochondrial genomes of normal and hyperplastic parathyroid glands.
Patient and tumor samples
Tissue samples were obtained, with Human Studies Institutional Review Board approval, from thirty-four patients who had undergone parathyroidectomy for management of primary hyperparathyroidism. Thirty samples were surgically and pathologically proven to be single parathyroid adenomas with no malignant features; of these, eighteen were determined to be typical chief cell adenomas and twelve were determined to be oxyphil adenomas, defined as containing at least 90% oxyphil cells. Eight samples from four patients were consistent with primary chief cell hyperplasia, involving multiple glands. None of the patients from whom tissue was obtained had a known history of head or neck irradiation, or a history of familial hyperparathyroidism or multiple endocrine neoplasia. One patient with primary chief cell hyperplasia had a history of lithium use. After surgical removal, tissues were dissected and quickly frozen in liquid nitrogen before being stored at -80°C. Peripheral blood leukocyte samples were obtained from the same patients, to serve as non-parathyroid germline controls. Also included in this study were five normal parathyroid glands and one normal thyroid gland from five patients with primary hyperparathyroidism due to a single adenoma. One normal parathyroid gland and a normal thyroid gland were obtained on autopsy. The remaining normal parathyroid glands were obtained from patients who had the glands removed incidental to an operation for thyroid disease.
Mitochondrial Genome Sequencing
Genomic DNA was extracted from each sample using either proteinase K digestion for surgical samples or sucrose gradient centrifugation for blood samples, followed by phenol-chloroform extraction and ethanol precipitation.
The entire 16.6 kb mitochondrial genome from each parathyroid sample was amplified and sequenced. To avoid amplification of nuclear mitochondrial pseudogenes, the mitochondrial genome was amplified as per Polyak et al , as 8 overlapping 1–3 kb PCR products using a step down PCR protocol or as 11 overlapping PCR products described as follows: PCR reactions were performed in 20 μl reaction volume including 25 ng DNA, 20 pmol of each primer, 200 μM each dNTP, 1.25 U of Amplitaq Gold DNA Polymerase (Applied Biosystems, Foster City, CA) and 1.5 mM magnesium chloride; thermocycling consisted of an initial denaturation step of 95°C for 10 minutes, 35 cycles of 95°C for 30 seconds, 55°C for 30 seconds, and 72°C for 3 minutes, and a final extension step at 72°C for 10 minutes. For difficult to amplify templates, additional PCR primers were designed and are available on request. PCR products were either gel purified using the Qiagen gel extraction kit (Qiagen, Valencia, CA), or 10 U Exonuclease I and 1 U Shrimp Alkaline Phosphatase (Amersham Pharmacia Biotech, Piscataway, NJ). Purified PCR products were subjected to automated sequencing using the Big Dye Terminator Cycle Sequencing Ready Reaction Kit (Applied Biosystems). We also made and utilized additional sequencing primers, sequences available on request. Sequencing products were purified through Sephadex G-50 columns (Sigma-Aldrich, St. Louis, MO), and electrophoresed on 4.75% Long Ranger gels (BioWhittaker Molecular Applications, Rockland, ME) on an ABI Prism 377 Sequencer (Applied Biosystems).
To determine if sequence changes identified in patients with hyperparathyroidism were somatic and/or clonal, the corresponding mitochondrial genomic regions from the patient's normal control peripheral blood leukocytes were also sequenced as above. Resulting sequence data files were analyzed and compared with the standard Anderson mtDNA sequence using the sequence analysis programs Sequencing Analysis and Auto Assembler (Applied Biosystems). Deviations from the standard sequence were compared with the on-line Mitomap database of previously reported mtDNA mutations and polymorphisms . All apparent somatic clonal mutations were confirmed by sequencing products of independent PCR reactions, and by sequencing the opposite DNA strand.
A total of 393 homoplasmic sequence variants were identified (8–45 variants per sample) as compared to the standard Anderson mtDNA sequence . Of the 393 variants, 375 were single base pair substitutions, 13 were single base pair insertions or deletions (4 deletions and 9 insertions), 2 were 2 base pair deletions, 2 were 2 base pair insertions and 1 was an 8 base pair deletion. Evaluation utilizing the mitochondrial database, Mitomap , indicated that of the 393 variants, 117 were in non-coding regions, 244 were in protein coding genes and 32 were predicted to affect rRNA or tRNA. Novel sequence variants are shown in Table 1.
Eleven somatic mutations were identified in 6 of 18 chief cell adenomas (Table 2 and Figure 2). Two tumors contained mutations that were predicted to result in early stop codons in either the ND4 or ND6 genes, both subunits of mitochondrial Complex I. Three tumors contained mutations that were predicted to result in non-synonymous changes in subunits of either Complex I or Complex III. Two tumors contained mutations only in the highly polymorphic, but non-coding D-loop region.
Sixteen somatic mutations were identified in 9 of 12 oxyphil adenomas. Four tumors contained mutations that were predicted to result in early stop codons in either the ND1, ND4 or ND5 genes, all subunits of mitochondrial Complex I. Two tumors contained non-synonymous changes in Complex I subunits, one tumor contained a mutation in CO1 not predicted to change the amino acid sequence and one tumor contained a mutation affecting a tRNA gene. One tumor contained a D-loop mutation only. All mutations appeared to be homoplasmic, with the exception of one, which appeared to be heteroplasmic, a G to A transition at position 559, within the D-loop region. Homoplasmy refers to the presence of a particular sequence variant or mutation in virtually every copy of the mitochondrial genome within a single tumor cell sample (and within the limits of sensitivity of the sequencing technology), while heteroplasmy refers to the presence of more than one such detectable allele . Heteroplasmy within a tumor would be indistinguishable from homoplasmy in a sub-population of the tumor's cells, since each tumor's cells are being analyzed collectively rather than on an individual cell-by-cell basis.
Mutations affecting Complex I genes, which represent 38.2% of the mitochondrial genome, accounted for 55.6% of all mutations and 22.2% of mutations occurred in the highly variable D-loop region.
To determine if the observed trend towards increased numbers of likely or potentially functionally significant mtDNA mutations in oxyphil adenomas as compared to chief cell adenomas was statistically significant, Fisher's Exact test was performed yielding a significant p-value of 0.02.
No somatic mutations occurred in any of the hyperplastic parathyroid glands. The one normal parathyroid for which paired thyroid tissue was available also contained no somatic mutations. No clustering of germline sequence variants occurred in any of the samples studied.
Mitochondrial abnormalities have been observed in many cancers and previous studies have demonstrated altered energy metabolism in cancer cells [2–10]. Homoplasmic somatic mutations of the mitochondrial genome have been identified in various tumor types. The mitochondrial genome is an especially attractive potential target for mutations that could drive tumorigenesis in a low turnover tissue such as the human parathyroid, given the frequent replication of mitochondrial DNA, independent of the nuclear genome and the enhanced potential for mitochondrial DNA damage by ROS. One might expect an increased potential for accumulation of somatic mutations in the mitochondrial genome in such low replicative tissues as the parathyroid over many years. Our data provide initial support for these hypotheses, and suggest that other slow-growing tumors arising from low-replication tissues should be similarly analyzed.
Mitochondrial DNA mutations in parathyroid adenomas tend to affect Complex I genes more often than other tumor types and more often then expected by chance alone. Additionally, D-loop mutations appear to be underrepresented in parathyroid adenomas compared to other tumor types. In contrast to many other studies, this study rigorously examines the entire mitochondrial genome in matched tumor/normal pairs from the same individual while many analogous studies of other types of tumors focused primarily on coding regions and/or the D-loop region or used population rather than control DNA from the same patient. Due to the wide range of mutation frequencies in the literature and discrepancies in study design, statistical analyses between our data and those from other related studies cannot generally be performed. The importance of observed differences between parathyroid adenomas and other tumor types remains unclear.
Our results raise the possibility that NADH dehydrogenase subunit genes might play a key role in tumorigenesis, perhaps especially in low-replication tissues. Few studies have directly examined this issue. The contribution of specific mtDNA mutations, affecting the ATP synthase subunit 6 [21, 22], to tumorigenesis has been examined, however little is known about the potential pathogenetic significance of the majority of mtDNA mutations identified to date. The contribution of specific germline sequence variants as predisposition alleles to specific tumor types has been suggested, however the majority of germline mtDNA sequence variants appear unlikely to directly contribute to tumorigenesis because of the frequency of such changes in individuals without tumors. Muller-Hocker, et al. examined respiratory chain function in normal parathyroid glands and hyperplastic parathyroids in uremic patients, finding defects in Complex III and IV in both normal and diseased glands. They also reported that chief cells in 20% of uremia-associated (secondary) hyperplastic parathyroid glands had defects of the respiratory chain, and that heart and muscle cells had more respiratory chain defects than parathyroid oxyphil cells. Together with reports of respiratory chain defects in normal, adenomatous, and hyperplastic parathyroid glands, the preferential clustering within NADH dehydrogenase subunit genes of the somatic mitochondrial mutations in our series of adenomas reinforces the intriguing suspicion that these mutations may be playing an integral role in parathyroid tumorigenesis. Given that mitochondria have recently been implicated in the regulation of intracellular calcium homeostasis [24, 25] and apoptosis , and that calcium sensing and signaling are crucial to parathyroid function and regulated PTH secretion, the possible relationship between calcium homeostasis and mitochondrial complex I enzyme function deserves further scrutiny.
Our data also provide additional support for the hypothesis that mitochondria accumulate in oxyphil cells due to abnormalities of the mitochondrial genome. Previous studies of Hurthle cell  and other oncocytic tumors  have demonstrated an association between a large deletion of mtDNA (Δ4977) and oncocytic/oxyphilic tumors in comparison to non-oncocytic tumors of the same tissue type. However, given its presence in both normal and tumor cells of various tissues, the potential pathogenetic significance of such a finding remains unclear. As such, the presence or absence of this deletion was not addressed by our experimental approach and was considered beyond the scope of this manuscript. Functioning oxyphil cell parathyroid adenomas are more likely to contain somatic mtDNA mutations than chief cell adenomas and the somatic mutations in oxyphil adenomas are more likely to be functionally significant. These differences between oxyphil and chief cell adenomas are statistically significant with a p-value of 0.02 (Fisher's Exact Test). Our results show marked similarity to the recent findings in thyroid tumors by Gasparre et al that disruptive mutations (nonsense and frameshift mutations) in Complex I genes are associated with an oncocytic phenotype . It is therefore quite plausible that these molecular genetic lesions may make key contributions to determining the oxyphilic phenotype, as mitochondria may increase their replication rate as a compensatory mechanism for altered energy metabolism due to mtDNA mutations. Further studies of normal oxyphil cells that could shed additional light on this hypothesis should be encouraged.
Yao K, Singer FR, Roth SI, Sassoon A, Ye C, Giuliano AE: Weight of normal parathyroid glands in patients with parathyroid adenomas. J Clin Endocrinol Metab. 2004, 89 (7): 3208-3213. 10.1210/jc.2003-031184.
Fliss MS, Usadel H, Caballero OL, Wu L, Buta MR, Eleff SM, Jen J, Sidransky D: Facile detection of mitochondrial DNA mutations in tumors and bodily fluids. Science. 2000, 287 (5460): 2017-2019. 10.1126/science.287.5460.2017.
Habano W, Sugai T, Nakamura SI, Uesugi N, Yoshida T, Sasou S: Microsatellite instability and mutation of mitochondrial and nuclear DNA in gastric carcinoma. Gastroenterology. 2000, 118 (5): 835-841. 10.1016/S0016-5085(00)70169-7.
Jones JB, Song JJ, Hempen PM, Parmigiani G, Hruban RH, Kern SE: Detection of mitochondrial DNA mutations in pancreatic cancer offers a "mass"-ive advantage over detection of nuclear DNA mutations. Cancer Res. 2001, 61 (4): 1299-1304.
Liu VW, Shi HH, Cheung AN, Chiu PM, Leung TW, Nagley P, Wong LC, Ngan HY: High incidence of somatic mitochondrial DNA mutations in human ovarian carcinomas. Cancer Res. 2001, 61 (16): 5998-6001.
Parrella P, Xiao Y, Fliss M, Sanchez-Cespedes M, Mazzarelli P, Rinaldi M, Nicol T, Gabrielson E, Cuomo C, Cohen D, Pandit S, Spencer M, Rabitti C, Fazio VM, Sidransky D: Detection of mitochondrial DNA mutations in primary breast cancer and fine-needle aspirates. Cancer Res. 2001, 61 (20): 7623-7626.
Penta JS, Johnson FM, Wachsman JT, Copeland WC: Mitochondrial DNA in human malignancy. Mutat Res. 2001, 488 (2): 119-133. 10.1016/S1383-5742(01)00053-9.
Polyak K, Li Y, Zhu H, Lengauer C, Willson JK, Markowitz SD, Trush MA, Kinzler KW, Vogelstein B: Somatic mutations of the mitochondrial genome in human colorectal tumours. Nat Genet. 1998, 20 (3): 291-293. 10.1038/3108.
Sanchez-Cespedes M, Parrella P, Nomoto S, Cohen D, Xiao Y, Esteller M, Jeronimo C, Jordan RC, Nicol T, Koch WM, Schoenberg M, Mazzarelli P, Fazio VM, Sidransky D: Identification of a mononucleotide repeat as a major target for mitochondrial DNA alterations in human tumors. Cancer Res. 2001, 61 (19): 7015-7019.
Yeh JJ, Lunetta KL, van Orsouw NJ, Moore FD, Mutter GL, Vijg J, Dahia PL, Eng C: Somatic mitochondrial DNA (mtDNA) mutations in papillary thyroid carcinomas and differential mtDNA sequence variants in cases with thyroid tumours. Oncogene. 2000, 19 (16): 2060-2066. 10.1038/sj.onc.1203537.
Parfitt A: Parathyroid Growth: Normal and Abnormal. The Parathyroids. Edited by: Bilezikian JP. 2001, San Diego: Academic Press, 293-330. 2
Bogenhagen D, Clayton DA: Mouse L cell mitochondrial DNA molecules are selected randomly for replication throughout the cell cycle. Cell. 1977, 11 (4): 719-727. 10.1016/0092-8674(77)90286-0.
Hamperl H: [Onkocytes and onkocytoma.]. Virchows Arch Pathol Anat Physiol Klin Med. 1962, 335: 452-483. 10.1007/BF00957034.
Roth S, Belsley N, Abu-Jawdwh G: Histology for Pathologists. Edited by: Mills S. 2006, Philadelphia: Lippincott-WIlliams and Wilkins, 1149-1165. 3
Munger BL, Roth SI: The cytology of the normal parathyroid glands of man and Virginia deer; a light and electron microscopic study with morphologic evidence of secretory activity. J Cell Biol. 1963, 16: 379-400. 10.1083/jcb.16.2.379.
Apel R, Asa S: The Parathyroid Glands. Endocrine Pathology. Edited by: LiVolsi VAaSLA. 2002, Philadelphia: Churchill Livingstone, 103-147.
Tallini G: Oncocytic tumours. Virchows Arch. 1998, 433 (1): 5-12. 10.1007/s004280050209.
Kogelnik AM, Lott MT, Brown MD, Navathe SB, Wallace DC: MITOMAP: a human mitochondrial genome database – 1998 update. Nucleic Acids Res. 1998, 26 (1): 112-115. 10.1093/nar/26.1.112.
Anderson S, Bankier AT, Barrell BG, de Bruijn MH, Coulson AR, Drouin J, Eperon IC, Nierlich DP, Roe BA, Sanger F, Schreier PH, Smith AJ, Staden R, Young IG: Sequence and organization of the human mitochondrial genome. Nature. 1981, 290 (5806): 457-465. 10.1038/290457a0.
Shadel GS, Clayton DA: Mitochondrial DNA maintenance in vertebrates. Annu Rev Biochem. 1997, 66: 409-435. 10.1146/annurev.biochem.66.1.409.
Petros JA, Baumann AK, Ruiz-Pesini E, Amin MB, Sun CQ, Hall J, Lim S, Issa MM, Flanders WD, Hosseini SH, Marshall FF, Wallace DC: mtDNA mutations increase tumorigenicity in prostate cancer. Proc Natl Acad Sci USA. 2005, 102 (3): 719-724. 10.1073/pnas.0408894102.
Shidara Y, Yamagata K, Kanamori T, Nakano K, Kwong JQ, Manfredi G, Oda H, Ohta S: Positive contribution of pathogenic mutations in the mitochondrial genome to the promotion of cancer by prevention from apoptosis. Cancer Res. 2005, 65 (5): 1655-1663. 10.1158/0008-5472.CAN-04-2012.
Muller-Hocker J, Aust D, Napiwotzky J, Munscher C, Link TA, Seibel P, Schneeweiss SG, Kadenbach B: Defects of the respiratory chain in oxyphil and chief cells of the normal parathyroid and in hyperfunction. Hum Pathol. 1996, 27 (6): 532-541. 10.1016/S0046-8177(96)90158-6.
Bernardi P: Mitochondrial transport of cations: channels, exchangers, and permeability transition. Physiol Rev. 1999, 79 (4): 1127-1155.
Kowaltowski AJ: Alternative mitochondrial functions in cell physiopathology: beyond ATP production. Braz J Med Biol Res. 2000, 33 (2): 241-250. 10.1590/S0100-879X2000000200014.
Green DR, Reed JC: Mitochondria and apoptosis. Science. 1998, 281 (5381): 1309-1312. 10.1126/science.281.5381.1309.
Maximo V, Soares P, Lima J, Cameselle-Teijeiro J, Sobrinho-Simoes M: Mitochondrial DNA somatic mutations (point mutations and large deletions) and mitochondrial DNA variants in human thyroid pathology: a study with emphasis on Hurthle cell tumors. Am J Pathol. 2002, 160 (5): 1857-1865.
Lewis PD, Baxter P, Paul Griffiths A, Parry JM, Skibinski DO: Detection of damage to the mitochondrial genome in the oncocytic cells of Warthin's tumour. J Pathol. 2000, 191 (3): 274-281. 10.1002/1096-9896(2000)9999:9999<::AID-PATH634>3.0.CO;2-U.
Gasparre G, Porcelli AM, Bonora E, Pennisi LF, Toller M, Iommarini L, Ghelli A, Moretti M, Betts CM, Martinelli GN, Ceroni AR, Curcio F, Carelli V, Rugolo M, Tallini G, Romeo G: Disruptive mitochondrial DNA mutations in complex I subunits are markers of oncocytic phenotype in thyroid tumors. Proc Natl Acad Sci USA. 2007, 104 (21): 9001-9006. 10.1073/pnas.0703056104.
Roth SI, Munger BL: The cytology of the adenomatous, atrophic, and hyperplastic parathyroid glands of man. A light- and electron-microscopic study. Virchows Arch Pathol Anat Physiol Klin Med. 1962, 335: 389-410. 10.1007/BF00957031.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1472-6823/7/8/prepub
We thank Dr. Stephen J. Walsh for assistance with statistical calculations and Ms Teresa Kim for assistance with experiments.
This work was supported in part by grant DHHS/NIDCR 5T32-DE07302 from the National Institutes of Health and by the Murray-Heilig Fund in Molecular Medicine.
The author(s) declare that they have no competing interests.
JC-G participated in the design of the study, carried out molecular genetic studies and sequence analyses, and drafted the manuscript. TT participated in the design of the study and carried out molecular genetic studies and sequence analyses. SIR performed histological analyses for all tumor samples. AA participated in the design and coordination of the study, data analysis and interpretation. All authors participated in revising drafts of the manuscript, and all authors read and approved the final manuscript.