Article Text

PDF

Heterogeneous tissue distribution of a mitochondrial DNA polymorphism in heteroplasmic subjects without mitochondrial disorders
  1. Elmar Kirchesa,
  2. Matthias Michaela,
  3. Michaela Warich-Kirchesa,
  4. Thomas Schneiderb,
  5. Serge Weisa,
  6. Guido Krausea,
  7. Christian Mawrina,
  8. Knut Dietzmanna
  1. aInstitute of Neuropathology, Otto-von-Guericke- University, Leipziger Strasse 44, 39120 Magdeburg, Germany, bDivision of Neurosurgery, Otto-von-Guericke- University, Magdeburg, Germany
  1. Dr Kirches, elmar.kirches{at}medizin.uni-magdeburg.de

Abstract

CONTEXT Several maternally inherited point mutations of the mitochondrial genome cause mitochondrial disorders, but the correlation between genotype and phenotype remains obscure in many cases. The same mutation may cause various diseases, probably because of a different tissue distribution.

OBJECTIVE To assess the role of random somatic segregation in generating interperson differences by analysis of an apparently neutral polymorphism.

DESIGN Screening of 81 brain samples from subjects without mitochondrial disorders and selection of five necropsy cases showing a high level of heteroplasmy for the polymorphism.

MAIN OUTCOME MEASURES A proportion of various distinct genotypes in the mtDNA pool of the tissues, identified by fluorescent PCR products, representing a short polycytosine tract of variable length in the mitochondrial displacement loop.

RESULTS Differences were found between organs or groups of organs within subjects, pointing towards somatic segregation of mtDNA. In addition, marked differences of this organ distribution occurred between subjects, which cannot be explained by tissue specific selection.

CONCLUSIONS The observed interperson differences can be explained by somatic segregation, which occurs randomly at various developmental stages. Besides tissue specific selection, this process might participate in the distribution of pathogenic mtDNA mutations.

  • mtDNA
  • polymorphism
  • HVR2
  • heteroplasmy

Statistics from Altmetric.com

The mitochondrial genome (mtDNA) is a double stranded circular DNA 16.6 kb in length, which encodes 22 tRNAs, two rRNAs, and 13 polypeptides of the oxidative phosphorylation (OXPHOS) complexes.1 The displacement loop (D loop) is the only major non-coding segment of mtDNA. Besides important elements for replication and transcription, it contains two hypervariable regions, designated HVR1 and HVR2, which have been widely used in evolutionary studies and population genetics2 and, recently, have received considerable attention in the field of forensic medicine.3 4

Since 1988, an increasing number of mtDNA mutations have been described, which are associated with either maternally inherited or sporadic disorders.5 6 Some classical mitochondrial diseases, like myoclonic epilepsy and ragged red fibres (MERRF)7 or mitochondrial encephalomyopathy with lactic acidosis and stroke-like episodes (MELAS),8 are caused by point mutations in tRNAs. All these defects are maternally inherited, since only mtDNA derived from the oocytes is transmitted to the embryo. Large deletions or duplications in the range of several kb are typically observed in chronic progressive external ophthalmoplegia (CPEO), in Kearns-Sayre syndrome (KSS), and Pearson's bone marrow syndrome, which are usually sporadic diseases.5 Most of these mutations are heteroplasmic, that is, both mutant and wild type DNA coexist in the same person.

Mitochondrial disorders can show considerable variation in the severity of symptoms, and the same mutation may cause a large variety of diseases. The MELAS syndrome may be taken as an example. The severity of the syndrome can vary greatly within a maternal lineage and only severely affected family members may develop full blown symptoms of the central nervous system. Furthermore, the most common mutation, the base substitution A3243G, can sometimes cause different and less severe disorders like CPEO or maternally inherited diabetes with deafness (MIDD). It may cause MELAS/KSS overlap syndromes or neurogenic muscle weakness, ataxia, and retinitis pigmentosa (NARP). These and other examples of phenotypic diversity can be found in the MITOMAP database (http://infinity.gen.emory.edu/mitomap.html) and in the review by Zeviani et al.6

It is widely accepted that the percentage of mutant mtDNA in the affected tissues controls the severity of the phenotype in many disorders, but it is more difficult to explain why the same mutation can cause a variety of distinct diseases. It is believed that the heterogeneous tissue distribution might cause this variability. Differences in the percentage of mutant mtDNA distributed between tissues of a subject could be explained by tissue specific selective pressure. However, the strong interperson differences of the tissue distribution of inherited point mutations cannot be explained without the participation of random segregation of heteroplasmic mtDNA genotypes during embryogenesis and/or effects of nuclear background. Although random segregation is widely hypothesised, to our knowledge it has never been shown in human necropsy material using neutral polymorphisms in non-diseased subjects. In this study, we show a high degree of heteroplasmy in a polycytosine tract (c-tract) of the mitochondrial HVR2 region (fig 1) in postmitotic tissues (81 brain samples). In a second step, we used this marker to analyse the tissue distribution of c-tract variants in various organs of five necropsy cases with a high level of heteroplasmy.

Figure 1

Diagram showing the mitochondrial control region between the tRNAs proline (Pro) and phenylalanine (Phe) of the mitochondrial genome. The numbers refer to the position of primers allowing the amplification of the complete hypervariable regions 1 and 2 (HVR1, HVR2) and to the positions of the thymine residues within the two incomplete c-tracts. The variable first part of the HVR2 c-tract investigated in this study (underlined) extends from np 303 of the human mtDNA to np 309. The tract of cytosine residues following the thymine in position 310 usually contains six cytosines in contrast to the originally published sequence.1

Materials and methods

SAMPLING AND STANDARD DNA ISOLATION

Fifty five brain samples were derived by laser microdissection from the tumour margin of paraffin embedded surgical astroglioma samples. The samples had originally been obtained during screening for the detection of somatic mutations in these tumours. The cases were selected only when tumour free normal brain parenchyma could be isolated by laser microdissection from the paraffin slides. Fifty two cases were glioblastomas, while three patients suffered from low grade astrocytomas. Laser microdissection of 0.03 mm2 areas from 5 μm thick slides was performed by an experienced neuropathologist to exclude contamination of the brain tissue samples with tumour tissue. Furthermore, areas with necrosis or other microscopically visible alterations as well as blood vessels were excluded by this technique. Twenty six necropsied brains and other organs of seven of these 26 patients were derived from subjects who had suffered from various diseases, but not mitochondrial disorders. In the necropsy cases, the tissue was taken from a 1.5 cm2, dewaxed, 10 μm thick paraffin section followed by proteinase K digestion, standard phenol/chloroform extraction, and ethanol precipitation of DNA.9. The dried DNA pellet was dissolved in 30 μl TE buffer. Five μl of this solution were used for subsequent PCR reactions. A serial control section had been stained in each case with haematoxylin and eosin and was studied by an experienced neuropathologist. The organs examined are listed in fig 2. Both sexes were represented in the study (53% male and 47% female, n=81) and the ages ranged from 8 to 89 years.

Figure 2

Electropherograms obtained during fragment length analysis of the mitochondrial HVR2 c tract in various organs of five patients with a high level of heteroplasmy. In three subjects (A), the quantitative peak pattern was well conserved among various organs, while it was not in two others (B). Patient 1: female, 65 years. Patient 2: male, 39 years. Patient 3: female, 66 years. Patient 4: male, 34 years. Patient 5: female, 61 years.

LASER MICRODISSECTION

Four cm2 pieces of a 1.3 μm polyethylene membrane (LPC membrane, PALM, Bernried, Germany) were mounted on cover slides. Five μm thick paraffin sections were placed on the poly-L-lysine treated membranes and incubated overnight at 37°C, followed by H&E staining (two minutes Mayer's haemalaun, one minute 2% eosin). After drying at room temperature overnight, small circles of 0.03 mm2 were microdissected by laser pressure capturing (LPC) using a Robot Microbeam System (PALM, Bernried, Germany). The dissected material was treated with 30 μl of 1 × PCR buffer containing 200 μg/ml proteinase K, digested for two hours at 55°C, and incubated for 10 minutes at 94°C. Ten μl of this solution were used for subsequent PCR reactions.

PCR ANALYSIS OF FRAGMENT LENGTH

Ten μl of the crude DNA extracts from microdissected samples or 5 μl of conventionally prepared DNA from single paraffin sections diluted with 5 μl TE buffer were used for subsequent PCR reactions. Goldstar-polymerase (Eurogentec, Belgium) (0.5 units) and all other PCR components were added in a total volume of 40 μl 1 × PCR buffer. PCR cycle conditions were one minute at 94°C, one minute at 58°C, and one minute at 72°C (35 cycles). The 6-FAM 5′ labelled forward primer was a 20 bp sequence starting at np 266 of mtDNA1and the 20 bp reverse primer started at np 427. After controlling the products on silver stained PAGE gels, the samples were digested withHaeIII, which creates a blunt end at np 333 immediately behind the c-tract of interest. Using cloned and sequenced HVR2 fragments from an earlier study10 and known HVR2 blood sequences, we could clearly determine the absolute c-tract length on a 373A sequencer (Applied Biosystems, Foster City, CA) by comparison with known sequences. Mixing experiments of known variants showed that a certain proportion of two variants always resulted in the same proportion of both peaks in the sequencing gels. Furthermore, a fraction of 10% of a minor variant with a length difference of 1 bp could be detected as a distinct peak. For peak recognition and determination of peak amplitudes, the GeneScan system from Applied Biosystems was used. The procedure for fragment length analysis is a simplified non-radioactive version of the trimmed PCR method introduced by Marchington et al 11 and is less sensitive in the detection of low levels of heteroplasmy. However, the sensitivity was sufficient for the purposes of this study.

Results

The highly sensitive analysis of the c-tract length allowed the detection of heteroplasmy in 54 of the 81 brain samples analysed (67%). In five selected necropsy cases with a high level of heteroplasmy in the brain, several additional organs were examined. Fig2 shows the electrophoretic peak patterns obtained from these tissues. Each peak in these electropherograms represents an mtDNA clone with a certain c-tract length. Two neighbouring peaks differ by 1 bp. While the absolute amplitudes depend on various experimental parameters (such as the amount of DNA in the sample), the relative amplitudes (peak ratios) remain unchanged in a constant proportion of different mtDNAs. This fact had been checked by mixing experiments with cloned DNA. Although the amplitude ratios are not numerically identical to the proportions of the corresponding mtDNA clones, they were used to define the peak ratios in table 1, since amplitudes can be measured more precisely than the peak areas. According to the electrophoretic mobility of cloned and sequenced reference DNA in the gels, the peaks in the diagrams could be labelled with their absolute c-tract length. It is known12 and is supported by our own sequencing data10 that usually only the first half of the incomplete HVR2 c-tract exhibits length variability (fig 1). Therefore, the peak designations in fig 2 reflect only the number of cytosine residues within this region.

Table 1

Peak ratios in various organs of patients 1, 2, and 3

In patient 1, the same five clones could be identified in all organs examined, but their ratios differed strikingly between the brain and the four body organs (fig 2A). Table 1 shows a detailed analysis of peak ratios of the four major genotypes of this patient, which were significantly above background in all organ samples (C7, C8, C9, and C10). It can be seen that all three ratios differed by less than 13% within the two endodermal derivatives, that is, lung and liver, and by less than 15% between the two mesodermal derivatives, that is, heart and kidney. On the other hand, some of the peak ratios differed up to 1.9-fold between endodermal and mesodermal derivatives, up to 6.8-fold between brain (neuroectoderm) and endodermal derivatives and up to 12.5-fold between brain and mesodermal derivatives.

When the peak ratio was calculated for the two dominant peaks in patient 2 (C8, C9), it differed by less than 11% in the liver, kidney, heart, and brain. The peak ratios between these organs and spleen differed by the factor 1.5. In addition, the high amplitude of the C7 peak in the spleen separated this organ from the others, in which this clone was barely detectable (fig 2A).

In patient 3, the peak ratios C7/C8 and C7/C9 differed by 23% and 3%, respectively, between the two endodermal derivatives, lung and liver, while the brain (neuroectoderm) and kidney (mesoderm) showed qualitatively different patterns (fig 2A). All tissues from patient 4 contained the same four c-tract genotypes, but in highly varying proportions, while in patient 5 qualitative/quantitative differences occurred between all tissues (fig 2B). In the deltoid muscle of patient 5, one peak could not be attributed with sufficient precision to a certain length variant in repeated measurements. In fig 2B, this peak was designated C10, although it seemed to consist of two independent variants with similar electrophoretic mobility.

To elucidate the stability of the HVR2 c tract further, we analysed in brain tissue of two cases without heteroplasmy (homoplasmic controls) the cerebral cortex (temporal and occipital lobes), cerebellar cortex, kidney, liver, lung, heart, skeletal muscle, spleen, colon, and oesophagus. In all the tissues of these two cases (a 63 year old male and a 89 year old female) only one identical c-tract genotype (C7) was detected.

Discussion

RELEVANCE OF THE HVR2 C-TRACT IN THE ANALYSIS OF mtDNA SEGREGATION IN OOCYTES AND TISSUES

A genetic bottleneck occurring during mtDNA germline transmission in mammals is an important concept in mitochondrial genetics for explaining rapid fixation of mutations. A complete switch to another homoplasmic genotype between mother and offspring has been observed in Holstein cows.13 Furthermore, estimates for the number of segregating units constituting the mtDNA pool of the offspring showed a strong variation depending on the species investigated and the mathematical models used.14-18 The possibility of a narrow bottleneck in humans would provide an explanation for the rapid acquisition of high mutational loads, which can sometimes be observed during transmission of mitochondrial disorders. In a female KSS patient with multiple mtDNA rearrangements including duplications, a postmortem analysis of the ovary was performed by Marchingtonet al.19 Patients with duplications may have a propensity for germline transmission.20 Rearranged molecules were detected in the germline of the KSS case, but the percentage of mutant DNA varied markedly among 15 oocytes.

Moreover, Marchington et al 11showed that the length of the HVR2 c-tract is a polymorphic marker, which can be used to determine differences in the mtDNA pools in mature human oocytes in any female. This method allowed a direct analysis of heteroplasmic mtDNA genotypes in the germline of normal subjects without mitochondrial disorders, and supported the occurrence of an early genetic bottleneck during oocyte maturation. In that study, some mature oocytes were found to be heteroplasmic in their HVR2 c-tract. This suggests that heteroplasmy in this region should sometimes be transmitted to the offspring. Therefore, the tissue distribution of this marker in some selected subjects with a high level of heteroplasmy provides a source of information about the fate of transmitted mtDNA during embryogenesis and further development. A prerequisite would be sufficient stability of the marker (low level of somatic mutations). This stability was tested in our study by analysing the two elderly homoplasmic subjects (63 and 89 years), showing the same homoplasmic genotype in all 11 tissues examined, which is in accordance with the results of Marchington et al 11obtained from two to five organs from five necropsy cases.11

TISSUE DISTRIBUTION OF MUTATIONAL LOAD AND SELECTIVE PRESSURE

The segregation of mtDNA within the mother's germline can be one factor responsible for the unpredictable level of heteroplasmy transmitted to the offspring. It can therefore partly explain the marked differences in phenotypic severity sometimes observed in a maternal lineage. The specific tissue distribution of mutations found in a certain disorder may be explained by the varying energy demands, differences between proliferating and postmitotic cells, and other factors, which might cause a tissue specific selective pressure. The accumulation of mtDNA deletions in muscle tissue of CPEO or KSS patients might be explained by a replicative advantage of shorter molecules in postmitotic tissues, where the disturbed mitochondria lead to no essential disadvantage for cell survival. On the other hand, deletions might be eliminated from rapidly dividing haematopoetic stem cells by negative selective pressure at the cellular level, since mitochondrial dysfunction may hamper cell division. Somatic segregation might lead to high levels of mutant mtDNA in some stem cells and low levels in others. Cells with a high mutational load may stop to divide. The positive selection of deleted molecules in postmitotic tissues might explain the accumulation of deletions in skeletal muscle over time associated with progression of disease.21

The pathogenic transition A3243G decreases over time in the blood,22 probably because the mitochondrial dysfunction hampers cell division. At first sight, a transition should confer no replicative advantage to the mtDNA molecules carrying it. Therefore, an increase over time in the mutational load of a point mutation in skeletal muscle23 seems more difficult to explain than an accumulation of deleted molecules. However, cell culture experiments with cybrids carrying the transition A3243G showed a marked replicative advantage.24 Moreover, not all pathogenic point mutations show similar behaviour. The transversion T8993G is not elimininated from blood. Generally, substitutions at nucleotide position 8993 were found at uniformly high heteroplasmic levels in all tissues examined so far.25

Tissue specific selection cannot explain how marked differences in the tissue distribution of a base substitution might occur, leading to a variety of different diseases. In addition, not all pathogenic mutations are prone to selective pressure.25

RANDOM SEGREGATION VERSUS TISSUE SPECIFIC SELECTION

In our study, nearly identical quantitative patterns of the heteroplasmic HVR2 c-tract were detected within the organ pairs lung/liver (endoderm) and heart/kidney (mesoderm) in a 65 year old woman (patient 1 in fig 2A, table 1), while the pattern differed markedly in all three germinal layers (endoderm, mesoderm, exoderm). The differences between the organ groups can be explained by random segregation of mtDNA. It is most likely that the process occurred at early stages of embryogenesis. This could explain the unequal distribution of the different genotypes to the three germinal layers in this case. In patient 2, four organs had a similar quantitative peak pattern, while the spleen had a different one. This result and the complex patterns observed in the other patients are compatible with random somatic segregation of mtDNA, which is not coupled to distinct stages of embryogenesis.

The variability of the HVR2 c-tract is a polymorphism, which is thought to confer no selective advantage or disadvantage to the cells carrying different genotypes. Nevertheless, a possible role of selection for the tissue distribution of this marker has to be discussed.

Surprisingly, some mtDNA polymorphisms in mammals seem to be not truly neutral and a tissue specific selection for some polymorphisms seems to occur. Jenuth et al 26 derived heteroplasmic mice with a hybrid mitochondrial genome containing mtDNA from the mouse strains BALB/cBy1 (BALB) and NZB/BINJ (NZB). An analysis of the tissue distribution of both genomes in the hybrid mice showed a decrease of the NZB genotype over one year postpartum in the blood and spleen and an increase in kidney and liver, while in several other tissues no significant changes over time could be observed. Furthermore, the nuclear background (BALB or BALB/NZB) did not influence the distribution of the mitochondrial genomes.

The occurrence of a tissue specific selection for polymorphisms was clearly established by these experiments, although the mechanism remains unknown. Differences in replicative properties caused by sequence variation in the D loop were unlikely, since the genomes did not differ in the conserved sequence elements. Another speculative explanation, suggested by the authors, is a dependence of mitochondrial turnover on mtDNA genotype. At least a selection at the cellular level could occur, based on mitochondrial function. Most mitochondrial proteins are encoded by nuclear genes. Since tissue specific differences in such proteins might exist, some of these isoforms might cooperate more effectively with certain isoforms of mtDNA encoded proteins.

In our study, a highly variable tissue distribution was found for a polymorphic marker in five people. Although compatible with random somatic segregation, this situation cannot be explained by tissue specific selection. In patient 1, the peak ratio C7/C8 was much higher in kidney, when compared with liver or lung, while the situation was reversed in patients 3 and 4 (table 1, fig 2B). In patient 2, the ratio C7/C8 was quite similar in kidney and liver (table 1). Nevertheless, the occurrence of selection is not ruled out by these results. Regarding the study of Jenuth et al,26 it seems possible that in some patients additional polymorphisms in mtDNA coding regions are coupled with certain c-tract genotypes, thereby providing a basis for selection in some organs.

Taken together, our results suggest that different mtDNA clones, which have passed the bottleneck in the female germline can be distributed in varying proportions to various organs during embryogenesis leading to a strong imbalance in the composition of mtDNA pools among organs. Marked interperson differences will result from the fact that this process occurs randomly during embryogenesis without being associated with certain developmental stages or tissues.

Acknowledgments

We gratefully acknowledge the excellent technical assistance of Mrs I Schellhase, G Küchler, I Otto, and S Hartmann. We thank Dr R Schneider-Stock and Prof Dr A Roessner (Magdeburg) for the use of the microdissection device. This work was supported by the State of Sachsen Anhalt, Contract Grant Sponsor 2526A/0086H and by the Deutsche Forschungsgemeinschaft, Contract Grant Sponsor KI 530/2-1.

References

View Abstract

Request permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.