Mitochondrial DNA mutations in human colonic crypt stem cells (original) (raw)

Article Stem cells Free access | 10.1172/JCI19435

Martin J. Barron,1 Gillian M. Borthwick,1 Amy Gospel,1 Patrick F. Chinnery,1 David C. Samuels,2 Geoffrey A. Taylor,3 Stefan M. Plusa,4 Stephanie J. Needham,5 Laura C. Greaves,1,6 Thomas B.L. Kirkwood,6 and Douglass M. Turnbull1

1Department of Neurology, and2Department of Mathematics, The Medical School, University of Newcastle upon Tyne, Newcastle upon Tyne, United Kingdom3Medical Research Council Development Centre for Brain Ageing, University of Newcastle upon Tyne, Newcastle upon Tyne, United Kingdom4Department of Surgery,5Department of Pathology, and6Department of Gerontology, The Medical School, University of Newcastle upon Tyne, Newcastle upon Tyne, United Kingdom

Address correspondence to: D.M. Turnbull, Mitochondrial Research Group, School of Neurology, Neurobiology and Psychiatry, The Medical School, University of Newcastle upon Tyne, Newcastle upon Tyne NE2 4HH, United Kingdom. Phone: 44-191-2228334; Fax: 44-191-2228553; E-mail: d.m.turnbull@ncl.ac.uk.

Find articles by Taylor, R. in:[JCI](/search/results?q=author.first%5Fname%3A%22Robert W.%22+author.last%5Fname%3A%22Taylor%22&search%5Ftype=advanced) |PubMed |Google Scholar

Find articles by Barron, M. in:[JCI](/search/results?q=author.first%5Fname%3A%22Martin J.%22+author.last%5Fname%3A%22Barron%22&search%5Ftype=advanced) |PubMed |Google Scholar

Find articles by Borthwick, G. in:[JCI](/search/results?q=author.first%5Fname%3A%22Gillian M.%22+author.last%5Fname%3A%22Borthwick%22&search%5Ftype=advanced) |PubMed |Google Scholar

Find articles by Gospel, A. in:JCI |PubMed |Google Scholar

Find articles by Chinnery, P. in:[JCI](/search/results?q=author.first%5Fname%3A%22Patrick F.%22+author.last%5Fname%3A%22Chinnery%22&search%5Ftype=advanced) |PubMed |Google Scholar

Find articles by Samuels, D. in:[JCI](/search/results?q=author.first%5Fname%3A%22David C.%22+author.last%5Fname%3A%22Samuels%22&search%5Ftype=advanced) |PubMed |Google Scholar

Find articles by Taylor, G. in:[JCI](/search/results?q=author.first%5Fname%3A%22Geoffrey A.%22+author.last%5Fname%3A%22Taylor%22&search%5Ftype=advanced) |PubMed |Google Scholar

Find articles by Plusa, S. in:[JCI](/search/results?q=author.first%5Fname%3A%22Stefan M.%22+author.last%5Fname%3A%22Plusa%22&search%5Ftype=advanced) |PubMed |Google Scholar

Find articles by Needham, S. in:[JCI](/search/results?q=author.first%5Fname%3A%22Stephanie J.%22+author.last%5Fname%3A%22Needham%22&search%5Ftype=advanced) |PubMed |Google Scholar

Find articles by Greaves, L. in:[JCI](/search/results?q=author.first%5Fname%3A%22Laura C.%22+author.last%5Fname%3A%22Greaves%22&search%5Ftype=advanced) |PubMed |Google Scholar

Find articles by Kirkwood, T. in:[JCI](/search/results?q=author.first%5Fname%3A%22Thomas B.L.%22+author.last%5Fname%3A%22Kirkwood%22&search%5Ftype=advanced) |PubMed |Google Scholar

Find articles by Turnbull, D. in:[JCI](/search/results?q=author.first%5Fname%3A%22Douglass M.%22+author.last%5Fname%3A%22Turnbull%22&search%5Ftype=advanced) |PubMed |Google Scholar

Published November 1, 2003 -More info

Published November 1, 2003 -Version history
Received: July 9, 2003; Accepted: September 9, 2003

View PDF

Abstract

The mitochondrial genome encodes 13 essential subunits of the respiratory chain and has remarkable genetics based on uniparental inheritance. Within human populations, the mitochondrial genome has a high rate of sequence divergence with multiple polymorphic variants and thus has played a major role in examining the evolutionary history of our species. In recent years it has also become apparent that pathogenic mitochondrial DNA (mtDNA) mutations play an important role in neurological and other diseases. Patients harbor many different mtDNA mutations, some of which are mtDNA mutations, some of which are inherited, but others that seem to be sporadic. It has also been suggested that mtDNA mutations play a role in aging and cancer, but the evidence for a causative role in these conditions is less clear. The accumulated data would suggest, however, that mtDNA mutations occur on a frequent basis. In this article we describe a new phenomenon: the accumulation of mtDNA mutations in human colonic crypt stem cells that result in a significant biochemical defect in their progeny. These studies have important consequences not only for understanding of the finding of mtDNA mutations in aging tissues and tumors, but also for determining the frequency of mtDNA mutations within a cell.

Introduction

Mitochondria are essential intracellular organelles required for oxidative phosphorylation. One of the unique aspects of mitochondria is the dual genetic control of mitochondrial biosynthesis. Mitochondria contain their own genome, mitochondrial DNA(mtDNA), which is a small (16.5 kb), self-replicating DNA molecule present in multiple copies in the mitochondrial matrix. This genome encodes 13 essential proteins of the mitochondrial respiratory chain as well as 2 ribosomal RNA (rRNA) and 22 transfer RNA(tRNA) genes required for intramitochondrial synthesis of proteins. All the other mitochondrial proteins are nuclear encoded and transported into the mitochondria. Inherited abnormalities of the mitochondrial genome are now increasingly recognized as important causes of disease (1). These mutations can affect all copies of the mitochondrial genome within a cell — termed homoplasmy — or there may be a mixture of mutated and wild-type genomes in the same cell — heteroplasmy (2). Since mitochondria contain multiple mtDNA copies and cells may contain hundreds or thousands of mitochondria, in the presence of heteroplasmy high levels of mutated mtDNA usually must be present within an individual cell before a defect of oxidative phosphorylation is apparent (3).

In addition to inherited mtDNA mutations, sporadic mtDNA mutations causing disease are common, and it has been noted for some time that mtDNA mutations occur with aging in postmitotic tissues (4–7). Both these events suggest that mtDNA mutation is a fairly frequent event, a view supported by the high sequence divergence within the human population. It is suspected that the high mutation rate is due to an environment in which there is a high production of free radicals by the respiratory chain, compounded by limited repair of this genome compared with nuclear DNA. Even if mutations occur at high frequency, the mechanism by which these mutations accumulate to high levels in individual cells remains uncertain, but mathematical models would favor clonal expansion by random genetic drift (8, 9).

Previous work on somatic mtDNA mutations has concentrated on the observations in postmitotic tissues, such as muscle and the CNS, since these tissues are predominantly involved in patients with mtDNA disease. Even in these tissues, systematic studies of the nature and incidence of mutations is limited. A small number of reports cite the accumulation of mtDNA mutations in replicative cells (6, 10). An important question in mitochondrial genetics is whether there is evidence of mtDNA mutations in stem cells, and if so, do these mutations accumulate to levels sufficient to result in a biochemical defect. If so, this could have profound implications for our understanding of the potential importance of these mutations in aging and cancer. To investigate the existence of mtDNA mutations in stem cells, we decided to study colonic crypt stem cells. Colonic crypt stem cells are present at the crypt base, and thus all the cells within the crypt must be derived from these stem cells. Thus, we would not have to isolate individual stem cells to look for mtDNA mutations, but could in theory study their progeny. While there is uncertainty about the number of stem cells within human colonic crypts, studies in mice suggest that crypts are maintained by small numbers of stem cells (11). For example, in chimeric mouse strains a single crypt often shows evidence of clonality in adult mice, suggesting that the majority of cells within an individual crypt are the progeny of a single stem cell (12). Studies in humans are much more difficult because strategies that directly mark cell fates, such as aggregation chimeras and mutagenesis, are impractical in humans (13); however, current evidence would support that the situation is similar to mice.

Methods

Patients and colonic samples. Colons were collected from 16 patients with no family history of colon cancer (age range: 50–80 years) undergoing resection for colonic tumors. Each sample of normal colonic mucosa was collected at a distance of at least 12 cm from the edge of the tumor. We also obtained colonoscopy biopsies from 12 individuals undergoing investigation of altered bowel function in which no pathology was found (age range: 36–78 years).

Histochemistry. Colon samples were mounted for sectioning and frozen in isopentane previously cooled to –160°C in liquid nitrogen. Standard H&E histology was performed on 12-μm sections. A dual histochemical assay was used to determine the respiratory activity of the crypts. Cryostat sections (12 μm) were first incubated in cytochrome c oxidase medium (100 μM cytochrome c, 4 mM diaminobenzidine tetrahydrochloride, and 20 μg/ml catalase in 0.2 M phosphate buffer, pH 7.0) at 37°C for 50 minutes. Sections were then washed in PBS (3 × 5 minutes) and incubated in succinate dehydrogenase (SDH) medium (130 mM sodium succinate, 200 μM phenazine methosulphate, 1 mM sodium azide, 1.5 mM nitroblue tetrazolium in 0.2 M phosphate buffer, pH 7.0) at 37°C for 40 minutes. Sections were washed in PBS, pH 7.4 (3 × 5 minutes), dehydrated in a graded ethanol series (70%, 95%, 2 × 100%), cleared in Histoclear (National Diagnostics, Atlanta, Georgia, USA), and mounted in DPX (BDH Laboratory Supplies, Poole, United Kingdom).

The percentage of cytochrome c oxidase–deficient crypts was calculated from transverse sections counted at several different levels, and approximately 100 crypts were evaluated at each level. Only crypts that were more than 50% cytochrome c oxidase deficient were included.

Immunocytochemistry. Immunocytochemistry was performed on 10-μm cryostat sections air-dried at room temperature for 1 hour. The sections were fixed in 4% paraformaldehyde in 0.1 M phosphate buffer for 10 minutes, followed by rinsing in distilled water and then PBS containing 0.1% Triton X-100 for 10 minutes. The sections were further permeabilized in a graded (70%, 95%, 100% vol/vol) methanol series over a period of 1 hour. Endogenous peroxidase activity was quenched at this stage by the addition of hydrogen peroxide (0.3% vol/vol) to the 95% (vol/vol) methanol. The sections were rinsed in PBS, then incubated with anti–cytochrome c oxidase subunit Ab (anti-subunit I, 2 μg/ml, anti-subunit II and IV, 5 μg/ml; Molecular Probes Inc., Eugene, Oregon, USA) diluted in 4% BSA in PBS for 1 hour at room temperature. Following three washes in PBS, the sections were incubated with a peroxidase-conjugated Ab (rabbit anti-mouse Ig’s, 50 μg/ml; DAKO Ltd., Ely, United Kingdom) for 1 hour at room temperature. Following washing in PBS, peroxidase activity was detected by incubating in 1.4 mM 3,3′-diaminobenzadine tetrahydrochloride and 0.01% (vol/vol) hydrogen peroxide in 0.1 M phosphate buffer, pH 7.4, for 5 minutes at room temperature.

Three-dimensional reconstruction of colonic crypts. Three-dimensional reconstruction of crypts was performed using KS300 3D software (Imaging Associates Ltd., Thame, United Kingdom) on digital images of 59 serial sections (8 μm) reacted for cytochrome c oxidase and SDH histochemistry.

Isolation of total DNA from individual and partial crypts. Fresh frozen sections (20 μm) of colon were mounted on polyethylenenaphthalate slides (Leica Microsystems UK Ltd., Milton Keynes, United Kingdom). Samples were subjected to dual histochemistry as described above and air-dried after dehydration. Single and partial crypts were cut into sterile 0.5-ml PCR tubes using a Leica laser microdissection system. Following centrifugation (7,000 g for 10 minutes), the crypt was lysed in 10 μl of cell lysis buffer (50 mM Tris-HCl, pH 8.5, 1 mM EDTA, 0.5% Tween-20, 200 ng/ml proteinase K) at 55°C for 2 hours and then 95°C for 10 minutes to denature the proteinase K.

Real-time PCR analysis of mtDNA. A quantitative real-time PCR approach employing fluorogenic probes was used to quantify the amount of mtDNA in two different regions of the mitochondrial genome, one that is rarely deleted in patients and aging (ND1 gene) and one that is frequently deleted (ND4 gene). PCR primers and fluorogenic probes (Applied Biosystems, Warrington, United Kingdom) for regions of ND1 (forward primer, L3485–3504; reverse primer, H3553–3532; probe, L3506–3529), and ND4 (forward primer, L12087–12109; reverse primer, H12170–12140; probe, L12111–12138) were synthesized, and 5 μl of DNA lysate from a single crypt was amplified separately with each of the ND1 and ND4 primer/probe combinations as described previously (14).

Sequencing of mtDNA from individual colonic crypts. The entire sequence of the mitochondrial genome from microdissected crypts was determined using the single-cell lysate (see above) as the DNA template and a two-stage amplification protocol essentially as described previously (15). The primary PCR reactions involved amplification of the mitochondrial genome in 11 fragments of approximately 2 kb using a series of overlapping primer pairs. These initial large PCR reaction products decrease the risk of amplifying pseudogenes when extracting DNA from small quantities of DNA (16). All PCR amplifications were performed in a 50 μl volume containing 1× PCR buffer (10 mM Tris-HCl, pH 8.3, 1.5 mM MgCl2, 50 mM KCl, 0.001% wt/vol gelatin), 0.2 mM dNTPs, 0.6 μM primers, 1 U AmpliTaq Gold DNA polymerase (Applied Biosystems), and 1 μl single-cell lysate. Reaction conditions were 94°C for 12 minutes and 38 cycles of 94°C for 1 minute, 60°C for 1 minute, and 72°C for 2 minutes. The final extension was for 8 minutes.

The secondary PCR reactions involved amplification of the primary PCR products with 28 primer pairs specifically to generate overlapping fragments of between 600–700 bp that span the entire sequence of the human mitochondrial genome (15). To facilitate the direct sequencing of PCR-amplified products, all primer pairs are tagged with M13 sequence and designed to anneal optimally at 58°C. All reactions proceeded for 30 cycles and used 2 μl of primary reaction product as DNA template. Samples were subsequently purified to remove unincorporated primers and were directly sequenced using BigDye terminator cycle-sequencing chemistries (Applied Biosystems) on an ABI 377 automated DNA sequencer. The sequences obtained were compared with the revised Cambridge reference sequence using Sequence Navigator and Factura software (Applied Biosystems).

Individual mutations were confirmed in all instances by repeating the second-round PCR and resequencing the product. For some products, the mutation was confirmed by restriction fragment–length polymorphism analysis using restriction enzymes that specifically recognized the mutation. Finally, to ensure mutations were not introduced during PCR, we microdissected two separate segments of the same cytochrome c oxidase–deficient crypt. Sequencing of the entire mitochondrial genome from both partial crypts revealed not only the presence of an identical mtDNA mutation, but no other detectable sequence changes.

Results

Respiratory chain deficiency in colonic crypts and crypt stem cells. Cells deficient in cytochrome c oxidase but with normal SDH activity are excellent biomarkers of mtDNA defects (17). We wished to determine if such cells were present in colonic epithelium and prepared frozen sections of normal colonic epithelium. The mucosal anatomy was preserved with our procedure for collection and mounting of sections, as shown by the crypt morphology using H&E (Figure 1, a and b). Some crypts from every patient we examined were cytochrome c oxidase deficient (Figure 1, c–h). To ensure that our techniques were completely specific, we performed control experiments on colonic epithelium using specific inhibitors of mitochondrial enzyme activity (sodium azide for cytochrome c oxidase and sodium malonate for SDH) (Figure 1, i–l). These experiments confirm that the reaction is specific and that our observations reflect a genuine decrease in cytochrome c oxidase activity. We found that these cytochrome c oxidase–deficient crypts were distributed throughout the colon because we examined samples obtained from both right and left hemicolectomy.

Figure 1

Respiratory chain deficiency in normal human colonic mucosa. (a) H&E preparations showing normal mucosal structure. Scale bar: 100 μm. (b) Higher magnification of a, showing normal crypt structure. Scale bar: 50 μm. (c) Normal cytochrome c oxidase activity (brown) in colonic crypts following dual cytochrome c oxidase and succinate dehydrogenase histochemistry. Scale bar: 100 μm. (d) Absence of histochemically detectable cytochrome c oxidase activity in colonic crypts (blue). Note delineation of crypt territories at the luminal surface. Scale bar: 100 μm. (e) Single cytochrome c oxidase–deficient colonic crypt. Scale bar: 100 μm. (f) Higher magnification of e, showing crypt base. Scale bar: 20 μm. (g) Presence of multiple, adjacent cytochrome c oxidase–deficient crypts. Scale bar: 100 μm. (h) Transverse section showing a similar cluster of cytochrome c oxidase–deficient crypts. Scale bar: 100 μm. The specificity of the dual cytochrome c oxidase and succinate dehydrogenase histochemical assay was established. (i) Colonic mucosa incubated in standard cytochrome c oxidase medium. (j) Colonic mucosa incubated as i, but in the presence of 2.5 mM sodium azide, a specific inhibitor of cytochrome c oxidase. (k) Colonic mucosa incubated in standard succinate dehydrogenase medium. (l) Colonic mucosa incubated as k, but in the presence of 50 mM sodium malonate, a competitive inhibitor of succinate dehydrogenase. Scale bars: 100 μm (i–k).

In addition to detecting whole crypts showing cytochrome c oxidase deficiency, we observed crypts showing a deficiency of this enzyme in only part of the crypt. To investigate the phenomenon of partial cytochrome c oxidase deficiency further, we serially sectioned (8 μm) colonic mucosa transversely and performed 3D reconstruction of the crypts following dual cytochrome c oxidase/SDH histochemistry (Figure 2a). In these crypts there is a ribbon of cytochrome c oxidase–deficient cells within an otherwise normal crypt that is entirely compatible with the view that there are multiple stem cells in some crypts and that what we are observing reflects involvement of one stem cell and its progeny.

Figure 2

Partial respiratory chain deficiency in some crypts and immunocytochemistry of cytochrome c oxidase deficiency. (a) Reconstructed images were produced from 59 transverse sections (8 μm) covering the full thickness of the colonic lamina propria. Representative images of the sections reacted for cytochrome c oxidase and succinate dehydrogenase are shown in the panel on the extreme right, with the reconstructed crypt arrowed. Scale bar: 50 μm. The position of each image corresponds to the approximate position in the adjacent reconstructed images. The images of the reconstructed colonic crypt are each rotated clockwise through 45°C to show a continuous ribbon of cytochrome c oxidase–deficient cells originating at the crypt base. (b–d) Representative serial transverse sections of colonic crypts showing immunoreactivity of cytochrome c oxidase subunit I (b), immunoreactivity of cytochrome c oxidase subunit IV (c), and dual cytochrome c oxidase/succinate dehydrogenase histochemistry (d). Scale bar: 50 μm (b–d).

We also studied the expression of cytochrome c oxidase subunits I (Figure 2b) and II (both mtDNA encoded) and subunit IV (nuclear encoded) (Figure 2c) in nine patients. In all patients there was an absence of immunoreactivity against subunits I (Figure 2b) and II in cytochrome c oxidase–deficient crypts. Immunoreactivity against subunit IV was reduced in all cytochrome c oxidase–deficient crypts (Figure 2c), but not decreased to the level seen with subunits I and II.

Mitochondrial genome analysis. To establish if the cytochrome c oxidase deficiency, observed in the colonic crypts was due to mtDNA mutations we identified individual crypts on longitudinal sections and isolated the crypt by laser microdissection. We searched for deletions of the major arc of the mitochondrial genome by real-time PCR and for point mutations of mtDNA by sequencing using techniques that have been modified to enable the investigation of single cells (14, 15).

Our studies have focused on five patients (age range 69–75 years; cytochrome c oxidase deficient crypts 2.8–27.9%). Using real-time PCR to screen for mtDNA deletions, we failed to detect any mtDNA deletions in either 40 cytochrome c oxidase–positive or 40 cytochrome c oxidase–deficient crypts. From each patient we then completely sequenced eight cytochrome c oxidase–deficient crypts and four cytochrome c oxidase–positive crypts. For each crypt we compared the sequence to the revised Cambridge reference sequence (18) and the sequence obtained from the other crypts from that patient. We detected numerous mtDNA point mutations in the crypts from these patients in both the heteroplasmic and homoplasmic state (Table 1). In some crypts, we detected only a single mtDNA point mutation, while in others there were multiple mutations (up to four per crypt). Mutations were seen in the cytochrome c oxidase–normal crypts as well as the deficient crypts. In the cytochrome c oxidase–normal crypts, the mutations seen are neutral polymorphisms or mutations that are not predicted to disrupt cytochrome c oxidase deficiency. Such changes include (i) pathogenic mutations (e.g., G15084A; ref. 19) involving genes encoding structural subunits of other respiratory chain complexes (I and III); (ii) mutations that occur at sites unlikely to induce an enzyme defect (e.g., G9985A affects subunit III of cytochrome c oxidase but alters the penultimate amino acid of this subunit); and (iii) mutations that are present at levels of heteroplasmy that are unlikely to affect mitochondrial protein synthesis (e.g., T10020C within the tRNAGly gene).

Table 1

Mitochondrial DNA sequence analysis of human colonic crypts

In contrast, in the cytochrome c oxidase–deficient crypts we detected pathogenic mutations, many involving the cytochrome c oxidase subunits. For example one of the mutations (C9537ins) in COIII (Figure 3b) has been reported previously to be pathogenic (20). In another cytochrome c oxidase–deficient crypt from a different patient we also detected a single nucleotide deletion in this same short, poly-C tract (C9537del) that induces a stop codon (Figure 3c). Other mutations include frameshift mutations, stop mutations, or highly conserved amino acid changes within the CO genes (Table 1). In some of the crypts in which we did not detect mutations involving the cytochrome c oxidase subunits, we detected high levels of heteroplasmy or even homoplasmic mutations of tRNA or rRNA genes (Table 1).

Figure 3

Examples of frameshift mutations within cytochrome c oxidase structural genes associated with cytochrome c oxidase deficiency. (a) Sequencing electropherogram detailing a region of COIII containing a short, poly-C tract of six residues (nucleotides 9532–9537) amplified from a cytochrome c oxidase positive–reacting colonic crypt. (b) Sequencing of this region in an alternative cytochrome c oxidase–deficient crypt revealed a homoplasmic single nucleotide insertion (C9537ins). The frameshift induced by this insertion creates a premature stop codon predicting the synthesis of a truncated cytochrome c oxidase subunit III polypeptide of 110 amino acids and as such is likely to be the cause of the enzyme defect (20). (c) In another cytochrome c oxidase–deficient colonic crypt, we identified a single C nucleotide deletion within this poly-C tract (C9537del) present at near homoplasmic levels. Similar to the C9537ins mutation, this change also leads to the premature truncation of the cytochrome c oxidase III subunit, resulting in cytochrome c oxidase deficiency. Abbreviations for the amino acid residues are as follows:Asn, asparagine; Gln, glutamine; Gly, glycine; Ile, isoleucine; Leu, leucine; Pro, proline; Thr, threonine. Ter represents termination codon.

The majority of mutations we detected were transitions involving G residues. This is similar to the observation made by Polyak et al. (21) when they performed mtDNA sequencing of human colonic tumors and probably results from G residues being the preferred target for oxidative damage. This would certainly be compatible with the concept that the mutations are acquired and secondary to the reactive oxygen species generated in the mitochondria (22).

Finally, patient 5 harbored two heteroplasmic mutations (G9025A and G9612A) that were present in several crypts (both cytochrome c oxidase positive and deficient). These mutations segregated together with levels varying from 0% to 100% within different crypts. It is likely that these are inherited mutations that are segregating through random drift, similar to a recent mouse model of mtDNA heteroplasmy (23).

Mathematical modeling of mtDNA mutation and replication in stem cells. Our observation of mtDNA mutations in individual crypts raises questions as to the mechanism by which these mutations develop in the stem cells. Consequently, we developed a computer model of mtDNA replication and segregation in human colonic crypt stem cells based upon experimentally derived parameters (9, 24–26). We studied the effects of a continuous somatic mutation rate on the mtDNA population within 1,000 simulated stem cells over an 80-year period. The simulation results are well described by simple probability functions. If mutations occur at a rate R per day per genome, then in a total cellular population size of N genomes the number of mutation events (E) that occur by time t is E = RNt. The probability (P) of a new mutation undergoing a clonal expansion by random drift is P = 1/N. Through simple probability arguments the theoretical curves are given by _P_0 = (1 – P)E, _P_1 = P_0_EP/(1 – P), _P_2 = (1/2)_P_1(E – 1)P/(1 – P) and _P_>2 = 1 –_P_0 – P1 – _P_2. With these simple parameters we saw the accumulation of colonic crypt cells containing high levels of mutant mtDNA through random genetic drift due to the random partitioning of mitochondrial genotypes during cell division (Figure 4a). The number of cells containing homoplasmic or near-homoplasmic mutants was dependent upon the chosen mutation rate.

Figure 4

Mathematical modeling of mitochondrial DNA mutation and replication in colonic crypt stem cells. Simulated cells contained approximately 10,000 mtDNA and divided every 24 hours. Each symbol is the result of 300 (at high mutation rates) to 3,000 (at low mutation rates) independent simulations. (a) The percentage of simulated crypt stem cells containing detectable (>30%) mutant mtDNA during a human lifespan. Each curve corresponds to a different mutation rate indicated on the figure. The model predicts that the random partitioning of mitochondrial genomes during crypt cell division will cause random genetic drift and lead to the clonal expansion of somatic mtDNA mutations during human life. (b) The percentage of simulated crypt stem cells containing a detectable amount of mtDNA mutations as a function of the mtDNA mutation rate (circles, no mutations, > 30%; triangles, one mutation, >30%; inverted triangles, two mutations, >30%; asterisk, three or more mutations. The curves represent the theoretical simple probability functions). The model predicts that the mutation rate must be approximately 5 × 10–5 per genome per day in order to simultaneously generate crypt stem cells that contain no mutations, and some that have two or more detectable mutations, at age 80 years. Previous simulations demonstrate that the random partitioning of individual genomes during stem cell division will cause random genetic drift and clonal expansion of somatic mtDNA mutations, and the speed of the random drift is dependent upon the frequency of the crypt stem cell divisions. Changing the rate of stem cell division to once per 48 hours did not alter the simulation results.

By giving each new mtDNA mutation a specific label, we were able to study the number of different mtDNA mutations within each simulated stem cell. Previous work has shown that our single-cell sequencing techniques detect levels of heteroplasmy of more than 30% for the lesser species (15). We therefore measured the number of different mtDNA mutations at more than 30% within each colonic crypt stem cell at age 80 years over a range of different mutation rates (Figure 4b). At lower mutation rates (10–6 mutations per genome per day) we saw only single clonally expanded mutations in a small percentage (<10%) of cells and a high proportion of stem cells containing only wild-type mtDNA. By contrast, at a higher mutation rate (2 × 10–4 mutations per genome per day) the vast majority (>95%) of stem cells contained three or more clonally expanded mutant species, and very few contained only wild-type mtDNA. Only at intermediate mutation rates (5 × 10–5 mutations per genome per day) did we see both single clonally expanded mtDNA mutations and multiple clonally expanded species affecting 10–30% of the simulated crypt stem cells and a similar number of cells containing only wild-type mtDNA. This has important implications for our understanding of somatic mtDNA mutation. Our data provide the first estimate of the somatic mtDNA mutation rate in humans in vivo of approximately 5 × 10–5 mutations per genome per day, corresponding to a rate of 1.10 mutations per site per million years. This is strikingly similar to the pedigree-derived mtDNA substitution rates (27) that are consistently higher than mutation rates derived from phylogenetic analysis (28). Moreover, the rate we have derived from the study of colonic crypts is similar to the rate measured for mtDNA (2 × 10–8 mutations per site per cell division) in cultured cells in vitro (29) and substantially greater than that of nuclear DNA (6 × 10–11 mutations per site per cell division).

Discussion

We hypothesized that mtDNA mutations may accumulate in stem cells in view of a high mutation rate in the mitochondrial genome and the possibility of accumulation of these mutations by a combination of relaxed replication and random segregation between daughter cells. Our investigations concentrated on colonic crypt stem cells because the progeny of these stem cells can be clearly identified and thus detecting the mutations and the resulting biochemical deficiency would be easier than for many other stem cells. Our investigations show that these stem cells do indeed accumulate mtDNA mutations to high levels and if the mutations involve an important base in either a structural or RNA gene then they can result in a biochemical defect. We have no reason to believe that these findings will be limited to colonic crypt stem cells and thus believe that our observations have implications for other stem cells.

Interestingly, in the cytochrome c oxidase–deficient crypts we detected predominantly mutations in the protein-encoding genes, which are very different from the mutation spectrum we see in patients with mtDNA disease. In our experience of investigating patients with inherited mtDNA defects, the number of mutations involving RNA genes is much greater than the mutations involving structural genes. This may indicate that the more severe phenotype associated with subunit mutations is less compatible with transmission through the germline, an observation in keeping with the relatively rare transmission of mtDNA deletions (30).

In some cytochrome c oxidase–deficient crypts we did not detect mtDNA point mutations, or, if they were present, they were likely to be polymorphic variants rather than pathogenic mutations. This was especially true of patients 4 and 5, the patients who had the highest number of cytochrome c oxidase–deficient crypts. There are several possibilities to explain the cytochrome c oxidase deficiency in these crypts, including mtDNA deletions that were not detected by our real-time assay or multiple mtDNA point mutations present below the limit of detection by automated sequencing. A further possibility would be mutations of nuclear genes. Our immunocytochemical studies revealed an absence of cytochrome c oxidase subunits I and II (encoded by the mitochondrial genome) and partial deficiency of subunit IV (encoded by the nuclear genome). There was no difference in the immunocytochemical reactions between those patients with high levels of mtDNA mutations and those (patient 4 and 5) in which the mtDNA mutations were not detected. These results could imply involvement of either the mitochondrial or nuclear genome. Nuclear gene mutations causing systemic cytochrome c oxidase deficiency are increasingly being recognized and may involve a nuclear-encoded cytochrome c oxidase structural subunit (although no mutations have been described to date) (31), a cytochrome c oxidase assembly factor, or a protein involved in mtDNA maintenance (32–37). If, indeed, the cytochrome c oxidase deficiency is due to a nuclear gene defect, it is uncertain as to whether this involves a single gene or several different genes, with each deficient crypt due to a new mutation. It is recognized that nuclear gene defects can cause a mosaic picture of cytochrome c oxidase deficiency and therefore the involvement of different crypts does not differentiate between the possibilities of a mitochondrial and nuclear gene involvement (38). We suspect that a nuclear gene defect causing either mtDNA depletion or multiple deletions would have been detected by our real-time assay, which not only detects deletions but also quantifies mtDNA copy number (14). We failed, however, to detect any difference in copy number between cytochrome c oxidase–deficient and cytochrome c oxidase–normal crypts in any of our patients. Another possibility is that the observation of low cytochrome c oxidase activity relates to a downregulation of the respiratory chain activity in response to the colonic tumor. This we think is unlikely because the cytochrome c oxidase–deficient crypts were also observed in patients without tumors, and they occur in the presence of normal SDH activity (Figure 1). Finally, if nuclear genetic defects (or altered regulation) are the cause of the cytochrome c oxidase deficiency in patient 4 and 5, what then could be the contribution of the nuclear genome to the cytochrome c oxidase–deficient crypts in our other patients? At this stage we can only speculate that we believe that the mtDNA defect is pathogenic in some crypts, but further studies are required to determine if the mutations are sufficient alone to generate the biochemical defect.

What then is the significance of our observations for mitochondrial genetics and stem cell biology? Our observations highlight how frequently these mutations occur and, when present, that they can accumulate to high levels, sufficient, we believe, to cause a biochemical defect. The observed mutation load in these stem cells is high; we sequenced the complete mitochondrial genome in 60 crypts and detected a total of 59 different mtDNA point mutations. In addition, there is likely to have been a significant number of mtDNA mutations present within these stem cells that were not detected. Automated fluorescent-sequencing techniques and real-time PCR will only detect moderate levels of a heteroplasmic change (>25% and >40%, respectively) (14, 39), and low-level heteroplasmy would be missed. If a significant number of stem cells accumulate mutations, it becomes clear that within an individual there will be a mixture of mutations, which may influence tissue function.

Our studies have shown that mtDNA mutations accumulate in cytochrome c oxidase–deficient crypts, and these crypts show an age-related increase in their number (Figure 5 and Table 2). These studies also represent the first description of pathogenic mtDNA mutations in stem cells with human aging. These observations support a role for mtDNA mutations in the aging process, which until now has been based mainly on an age-related accumulation of mtDNA mutations in postmitotic cells. Current understanding indicates that aging is caused by accumulation of stochastic molecular and cellular damage (40, 41), and there is evidence from studies in mice that functional impairment of intestinal stem cells occurs with age and is associated with an altered DNA damage response (42). The heterogeneous nature of such stochastic damage at the cellular level makes its detection more difficult and is fully consistent with the diversity we have observed. While our observations are in the colon, we believe they raise important questions about the wider role of mtDNA mutations in stem cell aging.

Figure 5

Incidence of cytochrome c oxidase–deficient colonic crypts per decade. Only crypts in which at least 50% of all cells were cytochrome c oxidase deficient were included. The crypts were counted on at least ten different transverse sections for 28 individual patients. The average value for each patient was included in the mean ± SEM value given for each decade. The number of patients per decade was 35–44 years, four patients; 45–54 years, seven patients; 55–64 years, five patients; 65–74 years, six patients; 75–84 years, six patients. (See also Table 2.) _R_2 is a measure of the conformation to an exponential increase where 0 represents no association and 1 is a perfect fit.

Table 2

Details of patients studied for the presence of cytochrome c oxidase–deficient crypts

Our findings are also relevant to reports of mtDNA mutations in colon (21) and other forms of cancer (43). The presence of mtDNA mutations within the colonic crypt stem cells suggests that the mutations are present prior to the development of the colon cancer. Thus, mtDNA mutations in other cancers suggest that these mutations might also originate in the stem cells in these tissues. It is as yet uncertain whether the mtDNA mutation and its associated biochemical defect are involved in the development of malignancy. There is a precedent, however, in that defects of oxidative metabolism caused by nuclear mutations may lead to the development of tumors (44, 45).

Acknowledgments

We are grateful to N. Howell for access to the MitoKor database and J. Elson for help with interpretation of sequence analysis. We would also like to thank R.N. Lightowlers, J. Mathers, E. Williams, and J. Burn for helpful discussions and provision of normal colonic material. These studies were supported by the Wellcome Trust, the Medical Research Council, and Cancer Research UK. A. Gospel was in receipt of a bursary from the PPP Foundation. Finally, we are grateful for the insightful comments of the anonymous reviewers of this manuscript, some of which have contributed to the discussion.

Footnotes

See the related Commentary beginning on page 1312.

David C. Samuels’ present address is: Virginia Bioinformatics Institute, Virginia Polytechnical Institute and State University, Blacksburg, Virginia 24061, USA.

Robert W. Taylor and Martin J. Barron contributed equally to this work.

Conflict of interest: The authors have declared that no conflict of interest exists.

Nonstandard abbreviations used: mitochondrial DNA (mtDNA); ribosomal RNA (rRNA); transfer RNA (tRNA); succinate dehydrogenase (SDH).

References

Chinnery, PF, et al. The epidemiology of pathogenic mitochondrial DNA mutations. _Ann. Neurol. 2000. 48:188-193.
Lightowlers, RN, Chinnery, PF, Turnbull, DM, Howell, N. Mammalian mitochondrial genetics: heredity, heteroplasmy and disease. _Trends Genet. 1997. 13:450-455.
Sciacco, M, Bonilla, E, Schon, EA, DiMauro, S, Moraes, CT. Distribution of wild-type and common deletion forms of mtDNA in normal and respiration-deficient muscle fibers from patients with mitochondrial myopathy. _Hum. Mol. Genet. 1994. 3:13-19.
Brierley, EJ, Johnson, MA, Lightowlers, RN, James, OF, Turnbull, DM. Role of mitochondrial DNA mutations in human aging: implications for the central nervous system and muscle. _Ann. Neurol. 1998. 43:217-223.
Cortopassi, G, Arnheim, N. Detection of specific mitochondrial DNA deletion in tissues of older humans. _Nucleic Acids Res. 1990. 275:169-180.
Michikawa, Y, Mazzucchelli, F, Bresolin, N, Scarlato, G, Attardi, G. Aging-dependent large accumulation of point mutations in the human mtDNA control region for replication. _Science. 1999. 286:774-779.
Murdock, DG, Christacos, NC, Wallace, DC. The age-related accumulation of a mitochondrial DNA control region mutation in muscle, but not brain, detected by a sensitive PNA-directed PCR clamping based method. _Nucleic Acids Res. 2000. 28:4350-4355.
Coller, HA, et al. High frequency of homoplasmic mitochondrial DNA mutations in human tumors can be explained without selection. _Nat. Genet. 2001. 28:147-150.
Chinnery, PF, Samuels, DC. Relaxed replication of mtDNA: a model with implications for the expression of disease. _Am. J. Hum. Genet. 1999. 64:1158-1165.
Nekhaeva, E, et al. Clonally expanded mtDNA point mutations are abundant in individual cells of human tissues. _Proc. Natl. Acad. Sci. U. S. A. 2002. 99:5521-5526.
Bach, SP, Renehan, AG, Potten, CS. Stem cells: the intestinal stem cell as a paradigm. _Carcinogenesis. 2000. 21:469-476.
Schmidt, GH, Winton, DJ, Ponder, BA. Development of the pattern of cell renewal in the crypt-villus unit of chimaeric mouse small intestine. _Development. 1988. 103:785-790.
Yatabe, Y, Tavare, S, Shibata, D. Investigating stem cells in human colon by using methylation patterns. _Proc. Natl. Acad. Sci. U. S. A. 2001. 98:10839-10844.
He, L, et al. Detection and quantification of mitochondrial DNA deletions in individual cells by real-time PCR. _Nucleic Acids Res. [serial online] 2002. 30:e68. http://nar.oupjournals.org/cgi/content/full/30/14/e68
Taylor, RW, Taylor, GA, Durham, SE, Turnbull, DM. The determination of complete human mitochondrial DNA sequences in single cells: implications for the study of somatic mitochondrial DNA point mutations. _Nucleic Acids Res. [serial online] 2001. 29:e74. http://nar.oupjournals.org/cgi/content/full/29/15/e74
Hirano, M, et al. Apparent mtDNA hetzeroplasmy in Alzheimer’s disease patients and normals due to PCR amplification of nuclear-embedded mtDNA pseudogenes. _Proc. Natl. Acad. Sci. U. S. A. 1997. 94:14894-14899.
Johnson, MA, Bindoff, LA, Turnbull, DM. Cytochrome c oxidase activity in single muscle fibers: assay techniques and diagnostic applications. _Ann. Neurol. 1993. 33:28-35.
Andrews, RM, et al. Reanalysis and revision of the Cambridge Reference Sequence for human mitochondrial DNA. _Nat. Genet. 1999. 23:147. (Letter).
Andreu, AL, et al. Exercise intolerance due to mutations in the cytochrome b gene of mitochondrial DNA. _N. Engl. J. Med. 1999. 341:1037-1044.
Tiranti, V, et al. A novel frameshift mutation of the mtDNA COIII gene leads to impaired assembly of cytochrome c oxidase in a patient affected by Leigh-like syndrome. _Hum. Mol. Genet. 2000. 9:2733-2742.
Polyak, K, et al. Somatic mutations of the mitochondrial genome in human colorectal tumours. _Nat. Genet. 1998. 20:291-293.
Cadenas, E, Davies, KJ. Mitochondrial free radical generation, oxidative stress, and aging(1). _Free Radic. Biol. Med. 2000. 29:222-230.
Jenuth, JP, Peterson, AC, Shoubridge, EA. Tissue-specific selection for different mtDNA genotypes in heteroplasmic mice. _Nat. Genet. 1997. 16:93-95.
Shadel, GS, Clayton, DA. Mitochondrial DNA maintenance in vertebrates. _Annu. Rev. Biochem. 1997. 66:409-435.
Gross, NJ, Getz, GS, Rabinowitz, M. Apparent turnover of mitochondrial deoxyribonucleic acid and mitochondrial phospholipids in the tissues of the rat. _J. Biol. Chem. 1969. 244:1552-1562.
Elson, JL, Samuels, DC, Turnbull, DM, Chinnery, PF. Random intracellular drift explains the clonal expansion of mitochondrial DNA mutations with age. _Am. J. Hum. Genet. 2001. 68:802-806.
Howell, N, Kubacka, I, Mackey, DA. How rapidly does the human mitochondrial genome evolve? _Am. J. Hum. Genet. 1996. 59:501-509.
Heyer, E, et al. Phylogenetic and familial estimates of mitochondrial substitution rates: study of control region mutations in deep-rooting pedigrees. _Am. J. Hum. Genet. 2001. 69:1113-1126.
Khrapko, K, et al. Mitochondrial mutational spectra in human cells and tissues. _Proc. Natl. Acad. Sci. U. S. A. 1997. 94:13798-13803.
Shanske, S, et al. Identical mitochondrial DNA deletion in a woman with ocular myopathy and in her son with Pearson syndrome. _Am. J. Hum. Genet. 2002. 71:679-683.
Adams, PL, Lightowlers, RN, Turnbull, DM. Molecular analysis of cytochrome c oxidase deficiency in Leigh’s syndrome. _Ann. Neurol. 1997. 41:268-270.
Papadopoulou, LC, et al. Fatal infantile cardioencephalomyopathy with COX deficiency and mutations in SCO2, a COX assembly gene. _Nat. Genet. 1999. 23:333-337.
Tiranti, V, et al. Mutations of SURF-1 in Leigh disease associated with cytochrome c oxidase deficiency. _Am. J. Hum. Genet. 1998. 63:1609-1621.
Zhu, Z, et al. _SURF_1, encoding a factor involved in the biogenesis of cytochrome c oxidase, is mutated in Leigh syndrome. _Nat. Genet. 1998. 20:337-343.
Mootha, VK, et al. Identification of a gene causing human cytochrome c oxidase deficiency by integrative genomics. _Proc. Natl. Acad. Sci. U. S. A. 2003. 100:605-610.
Antonicka, H, et al. Mutations in COX15 produce a defect in the mitochondrial heme biosynthetic pathway, causing early-onset fatal hypertrophic cardiomyopathy. _Am. J. Hum. Genet. 2003. 72:101-114.
Valnot, I, et al. Mutations of the SCO1 gene in mitochondrial cytochrome c oxidase deficiency with neonatal-onset hepatic failure and encephalopathy. _Am. J. Hum. Genet. 2000. 67:1104-1109.
Rahman, S, et al. Cytochrome oxidase immunohistochemistry: clues for genetic mechanisms. _Brain. 2000. 123:591-600.
Herrnstadt, C, et al. A high frequency of mtDNA polymorphisms in HeLa cell sublines. _Mutat. Res. 2002. 501:19-28.
Kirkwood, TB, Austad, SN. Why do we age? _Nature. 2000. 408:233-238.
Herndon, LA, et al. Stochastic and genetic factors influence tissue-specific decline in ageing C. elegans. _Nature. 2002. 419:808-814.
Martin, K, Potten, CS, Kirkwood, TB. Age-related changes in irradiation-induced apoptosis and expression of p21 and p53 in crypt stem cells of murine intestine. _Ann. N. Y. Acad. Sci. 2000. 908:315-318.
Fliss, MS, et al. Facile detection of mitochondrial DNA mutations in tumors and bodily fluids. _Science. 2000. 287:2017-2019.
Neumann, HP, et al. Germ-line mutations in nonsyndromic pheochromocytoma. _N. Engl. J. Med. 2002. 346:1459-1466.
Tomlinson, IP, et al. Germline mutations in FH predispose to dominantly inherited uterine fibroids, skin leiomyomata and papillary renal cell cancer. _Nat. Genet. 2002. 30:406-410.
Wu, CH, et al. The Protein Information Resource: an integrated public resource of functional annotation of proteins. _Nucleic Acids Res. 2002. 30:35-37.
MITOMAP: A human mitochondrial genome database. http://www.mitomap.org.
Herrnstadt, C, et al. Reduced-median-network analysis of complete mitochondrial DNA coding-region sequences for the major African, Asian, and European haplogroups. _Am. J. Hum. Genet. 2002. 70:1152-1171.
Taylor, R.W., et al. 2003. mtDNA polymorphisms in patient colonic crypts. http://www.ncl.ac.uk/nnp/nnp/research/mnrg/mtdnaage/colon.htm.

Version history

Version 1 (November 1, 2003): No description