Heteroplasmy and Evidence for Recombination in the Mitochondrial Control Region of the Flatfish Platichthys flesus

Galice Hoarau*,2, Suzanne Holla*, Rachel Lescasse{dagger}, Wytze T. Stam* and Jeanine L. Olsen*

*Department of Marine Biology, Centre for Ecological and Evolutionary Studies, University of Groningen;
{dagger}Laboratoire Biologie Cellulaire 4, Université Paris-Sud


    Abstract
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Acknowledgements
 References
 
The general assumption that mitochondrial DNA (mtDNA) does not undergo recombination has been challenged recently in invertebrates. Here we present the first direct evidence for recombination in the mtDNA of a vertebrate, the flounder Platichthys flesus. The control region in the mtDNA of this flatfish is characterized by the presence of a variable number of tandem repeats and a high level of heteroplasmy. Two types of repeats were recognized, differing by two C-T point mutations. Most individuals carry a pure "C" or a pure "T" array, but one individual showed a compound "CT" array. Such a compound array is evidence for recombination in the mtDNA control region from the flounder.


    Introduction
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Acknowledgements
 References
 
It has generally been assumed that animal mitochondrial DNA (mtDNA) does not undergo recombination. This assumption has come from indirect evidence such as failure to detect excision repair activity and crossover products (Clayton, Doda, and Friedberg 1975Citation ), failure to detect mtDNA recombination in interspecific somatic cell hybrids (Zuckerman et al. 1984Citation ), clustering of mtDNA molecules preventing physical contact (Satoh and Kuriowa 1991Citation ), and high mutation rate (Howell 1997Citation ). However, two recent studies have directly demonstrated mtDNA recombination in two invertebrates, a nematode Meloidogyne javanica (Lunt and Hyman 1997Citation ) and the mussel Mytilus galloprovincialis (Ladoukakis and Zouros 2001Citation ). With respect to vertebrates, only indirect evidence for mtDNA recombination has been found in the form of high frequencies of homoplasies (Eyre-Walker, Smith, and Maynard Smith 1999Citation ) and lower levels of linkage disequilibrium between more distance polymorphic sites (Awadalla, Eyre-Walker, and Maynard Smith 1999Citation ).

Here, we provide the first direct evidence of recombination in the mtDNA of a vertebrate, the flounder Platichthys flesus (Teleostei: Pleuronectiformes [flatfish]). In flounder, as in many other fishes (Nesbo, Arab, and Jakobsen 1998Citation , Ludwig et al. 2000Citation ), the 5'-end of the mitochondrial control region is characterized by the presence of a variable number of tandem repeats (VNTRs) and a high level of length heteroplasmy for this region. The cause of variation in the repeat number is not fully understood but is variously attributed to slipped strand mispairing, illegitimate elongation, and termination-associated sequence (TAS)-based replication (Ludwig et al. 2000Citation ). It may also be caused by mtDNA recombination, although this is more difficult to demonstrate when VNTRs consist of perfect repeats.


    Materials and Methods
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Acknowledgements
 References
 
A total of 168 individuals of flounder collected from five locations in the northern Atlantic and the Mediterranean were used to analyze the variation of the VNTRs and the extent of heteroplasmy (table 1 ). Total genomic DNA was extracted using a Cetyltrimethylammonium bromide protocol (Hoarau et al. 2002Citation ). Polymerase chain reaction (PCR) amplification of the 5'-end of the mitochondrial control region was carried out in 10 µl volumes containing 1 µl of 1/10 diluted DNA, 1x reaction buffer, 0.2 mM each of dNTP, 2 mM MgCl2, 0.25 U Taq DNA polymerase (Promega), and 0.5 µM each of the primers, DLF (5'-CCA CCT CTA ACT CCC AAA GC-3') (modified from primer A, Lee et al. 1995Citation ) and DLR3 (5'-GGT TAT TAT ACA TGA ATG TCC T-3'). For fragment analysis, primer DLF was fluorescently labeled with the 6-FAM (carboxyfluorescein) dye label (Applied Biosystem). Amplifications were performed in either a PTC100 (MJ Research) or a Mastercycler gradient cycler (Eppendorf). The reaction profile was as follows: initial denaturation at 94°C for 3 min; 30 cycles of denaturation at 94°C for 30 s, annealing at 50°C for 30 s, and extension at 72°C for 1 min; and a final extension at 72°C for 10 min. Length variation was visualized on an ABI Prism-377 automatic sequencer (Applied Biosystems) using a 4.5% denaturing polyacrylamide gel with an internal lane standard and GenescanTM software for determination of the fragment size. PCR products from 18 individuals were ligated into pGEM-T Easy plasmid (Promega) and cloned in Escherichia coli JM109 (Promega) according to the manufacturer's instructions. PCR amplification of 6–48 clones per individual (288 clones in total) was performed as described above and 3–6 clones per individual (84 clones in total) were sequenced using a BigDye Cycle Sequencing kit (Applied Biosystem) and an ABI Prism-377 automatic sequencer (Applied Biosystems). Negative controls were used at every step of the experiment.


View this table:
[in this window]
[in a new window]
 
Table 1 Heteroplasmy in Platichthys flesus

 

    Results and Discussion
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Acknowledgements
 References
 
Nearly all (165 of 168) individuals were heteroplasmic (table 1 ), containing 3–9 length variants per individual. The basic repeat unit is 19 bp long and the number of repeats ranges from 1 to 10. From the total set of 168 individuals, 18 individuals were chosen for further characterization. First, PCR products were cloned from these individuals and between 6 and 48 clones per individual (total of 288) were then reamplified and run on a gel to confirm that only one band was present as expected from individual clones. This was an important check because PCR artifacts (resulting in multiple bands) associated with these types of repeat motifs have been reported by Campbell, Sturm, and Barker (2001)Citation . We found single sharp bands in 274 clones, whereas 14 amplifications produced double bands. Products involving more than two bands were never observed. Thus, although artifacts accounted for 4.8% of the observations, the high level of heteroplasmy observed is real though perhaps slightly overestimated.

Again from the same set of 18 individuals, 3–6 clones per individual were sequenced (total of 84) to characterize the array motifs themselves and their number. Two types of core repeats were detected, a "C" type and a "T" type, which differed by two point mutations (fig. 1 ). Among the 18 individuals, 13 contained only the pure "C" type array, 4 only the pure "T" type array, and 1 individual a compound array of "C" and "T" (fig. 1A, B, and C , respectively). The compound array is indicative of mtDNA recombination. It is highly unlikely that the compound array arose by mutation as the differences between the "C" type and the "T" type involve two independent mutations and neither seems to be associated with secondary structure such as hairpin (no palindrome in the sequence of the repeat unit). Moreover these two positions are highly conserved in the closely related species, Pleuronectes platessa (unpublished data).



View larger version (11K):
[in this window]
[in a new window]
 
Fig. 1.—Alignment of the VNTR of mtDNA control region in Platichthys flesus. {block} C type repeat, T type repeat. The two point mutations are boldface. Ind1 shows a pure "C" type array, Ind2 a pure "T" type array, and Ind3 a "CT" compound array. GenBank AY143162 to AY143164

 
Once again, however, the issue of possible PCR artifacts arises because Taq polymerase can produce chimeric DNA molecules (i.e., PCR jumping sensu Paabo, Irvin, and Wilson 1990Citation ) when damaged fragments of a mixed template of DNAs are joined during the amplification process. Mixed templates of undamaged fragments do not result in PCR jumping (Ladoukakis and Zouros 2001Citation ). To check for this possibility, equal amounts of DNA from pure "C" and pure "T" individuals were mixed in the same PCR reaction in an effort to create a chimeric product. The PCR product was cloned and 16 clones were sequenced. Ten were pure "C" and six were pure "T". There were no "CT" arrays detected. Thus a PCR artifact is unlikely.

Further support for true recombination comes from the fact that the recombinant "CT" array was found in each of three independent DNA extractions, PCR, and sequencing reactions from the recombinant individual. In addition, the entire DNA extraction, amplification, cloning, and sequencing procedure was conducted in an independent laboratory where no studies on fish have ever been conducted. Twelve clones were sequenced and recombinant "CT" arrays were found in seven clones and a parental "T" was found in the five other clones (table 2 ). Interestingly, only one of the "parental" types ("T") is found together with recombinant arrays. This suggests that recombination took place at least a generation ago and that the other "parental" type ("C") was lost by random drift during ontogenesis or gametogenesis (Chinnery et al. 2000Citation ). Variation in the number of "C" or "T" repeats in recombinant arrays further suggests that other mechanisms, such as slipped strand mispairing, besides (or in addition to) recombination are involved in the evolution of the number of repeats.


View this table:
[in this window]
[in a new window]
 
Table 2 Type of Repeat Arrays Found in the Recombinant Individual

 
Recombination can explain the evolution of mtDNA VNTRs, but it is usually difficult to trace when VNTRs exist as perfect arrays. Because flounder VNTRs are imperfect and the two different types of repeat unit are found, the demonstration of mtDNA recombination is straightforward.

In principle, an alternative hypothesis to recombination could be that the compound array is the ancestral form and that the derived type of array we have observed evolved by a series of duplication-and-loss events. In this scenario, the two mutations (i.e., switching between the C to T type and vice versa) occurred in a duplicated region, such that both it and the ancestral form would be retained in the same individual. The different number of basic repeats could then be accounted for by a history of duplication and loss. Given the range in the number of repeat regions found (1–10), such duplications and losses might occur fairly regularly. However, we do not favor this alternative hypothesis because (1) if the "CT" arrays are indeed the ancestral form, then they should have been found more frequently in P. flesus and in closely related flatfish such as Pleuronectes platessa (plaice)—plaice contains only one repeat, suggesting that the ancestral form was one repeat and (2) this hypothesis requires that the two mutations be retained in the same individual, whereas the recombination scenario requires only that the two types be present in the same population.

The mtDNA recombination we found between "C" and "T" type arrays in flounder has several implications for vertebrate mtDNA in general. Before such recombination is possible, the "C" and "T" type mitochondria must be present in the same cell and their DNA must be coupled. This implies paternal leakage followed by fusion of the mitochondria. Paternal leakage has been reported for mice (Gyllensten et al. 1991Citation ) and anchovies (Magoulas and Zouros 1993Citation ), whereas fusion of mitochondria has been demonstrated in Drosophila (Yaffe 1999Citation ), and the enzymes necessary for recombination have been found in human mitochondria (Thyagarajan, Padua, and Campbell 1996Citation ). It appears therefore that all these properties are present in flounder and in vertebrates more generally.

The importance of recombination in vertebrate mitochondria has broad implications across several fields, ranging from human mitochondrial diseases (Schon 2000Citation ) to the compromise of phylogenetic and population studies that assume strict clonal inheritance of mtDNA (Schierup and Hein 2000Citation ). In the case of human mitochondrial diseases, mtDNA recombination will greatly change the mode and patterns of inheritance, which in turn may affect current diagnostic methods. Recombination can also affect the accuracy of phylogenetic reconstruction (Posada and Crandall 2002Citation ), inferences related to demographic history, and the application of molecular clocks (Schierup and Hein 2000Citation ).


    Acknowledgements
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Acknowledgements
 References
 
We thank J. Coyer and M-P. Oudot-Le Secq and two anonymous reviewers for comments and corrections on the earlier versions of the manuscript—especially with respect to the alternative hypothesis to recombination. We also thank H. Van der Veer and R. Gesser for providing samples, S. Boele-Bos for laboratory assistance, and L. Van de Zande for laboratory facilities. This research was supported by NWO-Prioriteit Programma "Sustainable Use of Marine Natural Resources" Project No. 885-10.101-P.


    Footnotes
 
Edward Holmes, Reviewing Editor

Keywords: heteroplasmy recombination mitochondria DLoop VNTR Platichthys flesus Back

Address for correspondence and reprints: Galice Hoarau, Department of Marine Biology, Centre for Ecological and Evolutionary Studies, University of Groningen, P.O. Box 14, 9750 AA Haren, The Netherlands. E-mail: g.hoarau{at}biol.rug.nl Back


    References
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Acknowledgements
 References
 

    Awadalla P., A. Eyre-Walker, J. Maynard Smith, 1999 Linkage disequilibrium and recombination in hominid mitochondrial DNA Science 284:2090-2091[Medline]

    Campbell N. J. H., R. A. Sturm, S. C. Barker, 2001 Large mitochondrial repeats multiplied during the polymerase chain reaction Mol. Ecol. Notes 1:336.[ISI]

    Chinnery P. F., D. R. Thorburn, D. C. Samuels, S. L. White, H.-H. M. Dahl, D. M. Turnbull, R. N. Lightowlers, N. Howell, 2000 The inheritance of mitochondrial DNA heteroplasmy: random drift, selection or both? Trends Genet 16:500-505[ISI][Medline]

    Clayton D. A., J. N. Doda, E. C. Friedberg, 1975 The absence of a pyrimidine dimer repair mechanism in mammalian mitochondria Proc. Natl. Acad. Sci. USA 71:2777-2781[ISI]

    Eyre-Walker A., N. H. Smith, J. Maynard Smith, 1999 How clonal are human mitochondria? Proc. R. Soc. Lond. B Biol. Sci 266:477-483[ISI][Medline]

    Gyllensten U., D. Wharton, A. Josefsson, A. C. Wilson, 1991 Paternal inheritance of mitochondrial DNA in mice Nature 352:255-257[ISI][Medline]

    Hoarau G., A. D. Rijnsdorp, H. W. Van der Veer, W. T. Stam, J. L. Olsen, 2002 Population structure of plaice (Pleuronectes platessa L.) in northern Europe: microsatellites revealed large scale spatial and temporal homogeneity Mol. Ecol 11:1165-1176[ISI][Medline]

    Howell N., 1997 MtDNA recombination: what do in vitro data mean? Am. J. Hum. Genet 61:19-22[ISI][Medline]

    Ladoukakis E. D., E. Zouros, 2001 Direct evidence for homologous recombination in mussel (Mytilus galloprovincialis) mitochondrial DNA E. Mol. Biol. Evol 18:1168-1175

    Lee W. J., J. Conroy, W. H. Howell, T. D. Kocher, 1995 Structure and evolution of teleost mitochondrial control regions J. Mol. Evol 41:54-66[ISI][Medline]

    Ludwig A., B. May, L. Debus, I. Jenneckens, 2000 Heteroplasmy in the mtDNA control region of sturgeon (Acipenser, Huso, and Scaphirhynchus) Genetics 156:1933-1947[Abstract/Free Full Text]

    Lunt D. H., B. C. Hyman, 1997 Animal mitochondrial DNA recombination Nature 387:247.[ISI][Medline]

    Magoulas A., E. Zouros, 1993 Restriction-site heteroplasmy in anchovy (Engraulis encrasicholus) indicates incidental biparental inheritance of mitochondrial DNA Mol. Biol. Evol 10:319-325[Free Full Text]

    Nesbo C. L., M. O. Arab, K. S. Jakobsen, 1998 Heteroplasmy, length and sequence variation in the mtDNA control regions of three percid fish species (Perca fluviatilis, Acerina cernua, Stizostedion lucioperca) Genetics 148:1907-1919[Abstract/Free Full Text]

    Paabo S., D. M. Irwin, A. C. Wilson, 1990 DNA damage promotes jumping between templates during enzymatic amplification J. Biol. Chem 265:4718-4721[Abstract/Free Full Text]

    Posada D., K. A. Crandall, 2002 The effect of recombination on the accuracy of phylogeny estimation J. Mol. Evol 54:396-402[ISI][Medline]

    Satoh M., T. Kuriowa, 1991 Organization of multiple nucleoids and DNA molecules in mitochondria of a human cell Exp. Cell Res 196:137-140[ISI][Medline]

    Schierup M. H., J. Hein, 2000 Consequences of recombination on traditional phylogenetic analysis Genetics 156:879-891[Abstract/Free Full Text]

    Schon E. A., 2000 Mitochondrial genetics and disease Trends Biochem. Sci 25:555-560[ISI][Medline]

    Thyagarajan B., R. A. Padua, C. Campbell, 1996 Mammalian mitochondria possess homologous DNA recombination activity J. Biol. Chem 271:27536-27543[Abstract/Free Full Text]

    Yaffe P. M., 1999 The machinery of mitochondrial inheritance and behavior Science 283:1493-1497[Abstract/Free Full Text]

    Zuckerman S. H., J. F. Solus, F. P. Gillespie, J. M. Eisenstadt, 1984 Retention of both parental mitochondrial DNA species in mouse-Chinese hamster somatic cell hybrids Somat. Cell Mol. Genet 10:85-91[ISI][Medline]

Accepted for publication August 22, 2002.