Pause-melting misalignment: a novel model for the birth and motif indel of tandem repeats in the mitochondrial genome
© Shi et al.; licensee BioMed Central Ltd. 2013
Received: 2 September 2012
Accepted: 11 February 2013
Published: 15 February 2013
Tandem repeats (TRs) in the mitochondrial (mt) genome control region have been documented in a wide variety of vertebrate species. The mechanism by which repeated tracts originate and undergo duplication and deletion, however, remains unclear.
We analyzed DNA sequences of mt genome TRs (mtTRs) in the ridged-eye flounder (Pleuronichthys cornutus), and characterized DNA sequences of mtTRs from other vertebrates using the data available in GenBank. Tandem repeats are concentrated in the control regions; however, we found approximately 16.6% of the TRs elsewhere in the mt genome. The flounder mtTRs possess three motif types with hypervariable characteristics at the 3′ end of the control region (CR).
Based on our analysis of this larger dataset of mtTR sequences, we propose a novel model of Pause Melting Misalignment (PMM) to describe the birth and motif indel of tandem repeats. PMM is activated during a pause event in mitochondrial replication in which a dynamic competition between the nascent (N) heavy strand and the displaced (D) heavy strand may lead to the melting of the N-strand from the template (T) light strand. When mispairing occurs during rebinding of the N-strand, one or several motifs can be inserted or deleted in both strands during the next round of mt-replication or repair. This model can explain the characteristics of TRs in available vertebrate mt genomes.
The mitochondrial (mt) genome appears to evolve faster than the single-copy fraction of the nuclear genome by approximately a factor of ten . The portion of the mt genome that evolves most rapidly is the control region (CR), or the D-loop region, presumably due to lower selection pressure [1–3]. Within the control region of vertebrate mt genomes, tandem repeats (TRs) located at the CR termini (the 3′- or 5′- end) evolve most rapidly, with extensive variations in motif size and copy number resulting from motif insertions and deletions (indels) in TR tracts [4–7].
To date, several hypotheses have been proposed to explain the birth and indel of mt genome tandem repeats (mtTRs). Two related theories explaining the mechanism by which motif indels are created in nuclear genomes have been applied to similar events in mt genomes. In one such theory, recombination, which has been confirmed in mt genomes [5, 8–11], is responsible. In this case, double strand breakage and rejoining of DNA lead to the gain or loss of mtTR motifs. An alternative hypothesis, the Slipped-Strand Mispairing (SSM) model , not only accounts for microsatellite indels [13–15], but also for motif indels in mtTRs [6, 16, 17]. In its simplest form, the SSM model involves local denaturation and displacement of the strands of a DNA duplex, followed by mispairing of complementary bases at the site of an existing similar sequences or short tandem repeat. When followed by replication or repair, misalignment could lead to indels of one or several simple repeats.
A few models have been proposed to explain particular indel events in mt genomes [18–21], and subsequently have been adopted in some other cases [5, 22–24]. For example, the Illegitimate Elongation Model (IEM) suggested by Buroker  involves a competitive equilibrium between the nascent (N) heavy strand and the displaced (D) heavy strand while the N-strand is arrested at a termination-associated sequence (TAS) in the 5′ end of the CR. This model suggests that frequent misalignment occurs prior to elongation in the repeat region and is facilitated by a stable secondary structure in the D-strand. This mechanism could result in both length and sequence heteroplasmy of TRs. Another influential model, the Improper Initiation Model (IIM), was proposed by Broughton and Dowling . In this process, mt-replication begins at the upstream tRNA Phe gene, an improper origin of mt-replication. If mt-replication is completed at the regular termination site, the sequence between tRNA Phe and the normal origin of replication will be repeated on the nascent heavy strand. Under these circumstances, the repeated sequence will always be identical in length and position. Other mechanisms used to explain the origin of special motif indels include that of Taylor and Breden , who present a model for mini-satellite birth from non-contiguous repeats in the mt genome of the guppy, Poecilia reticulate, and that of Hayasaka , used to explain the motif indels of mtTRs involving conserved sequence blocks 2 and 3 (CSB2, CSB3) in the CR of the Japanese monkey, Macaca fuscata.
Each of these models possesses certain limitations. Previous studies have provided some evidence of recombination in mt genomes, but the incidence of mitochondrial DNA recombination is thought to be rare [7, 25–27]. In fact, indels of repeat regions in mtDNA are too abundant to be explained solely by the recombination model. A prediction of the SSM model, due to its intra-helical nature, is that there should be an appreciable bias toward the duplication of short simple repeat units; therefore, the likelihood of a tandem repeat forming would decrease rapidly as the length of the repeat motif increases . Consequently, this model is unlikely to explain the higher proportion of long motif repeats. Both IEM and IIM best explain the formation of motif indels only in a particular position of the CR.
How do tandem repeats form, and subsequently insert and delete in mt genomes, and why do they have high sequence and length heteroplasmy even within an individual? In order to address these questions, we analyzed the sequences of highly variable mtTRs in the ridged-eye flounder (Pleuronichthys cornutus), and explored the characteristics of vertebrate mtTRs using data from GenBank. We then developed a novel model for the birth and motif indel of TRs, in the hope of further understanding the evolutionary mechanism of mtTRs variability, and the evolution of mt genomes more generally.
The amplification of P. cornutusCRs
We used 20 individuals of P. cornutus, thirteen from Dongshan, East China Sea and seven from Qingdao, Yellow Sea. A portion of the epaxial musculature was excised from each fresh specimen, and immediately stored at −70°C. Total genomic DNA was extracted using an SQ Tissue DNA Kit (OMEGA, China) following the manufacturer’s protocol. Two primer pairs (L-DL: ACTCCCAAAGCCAGGATTCT, H-DL: GAGGGTGAGGTTTAACGGGGG; L-AP1: GCCTGTAGCTTTTTAGGTAT, H-AP2: AAGCATAACACTGAAGATG) were designed for amplification of P. cornutus control region sequences. The PCR was performed in a 25μl reaction volume containing 2.0 mM MgCl2, 0.4 mM dNTP, 0.5 μM of each primer, 1.0 U Taq polymerase (Takara, Dalian, China), 2.5 μl 10x Taq buffer, and approximately 50 ng DNA template. PCR cycling conditions included an initial denaturation at 95°C for 3 min, 35 cycles of a denaturation at 94°C for 45 s, an annealing step at 48°C for 45 s, and elongation at 72°C for 3 min with a final extension at 72°C for 5 min. The PCR products were detected in 1.0% agarose gels and purified using a Takara Agarose Gel DNA Purification Kit (Takara, China). Purified products were inserted into the PMD18-T vector (Takara, China), and then transformed into E. coli competent cells. The final cloned products were sequenced in both directions using the ABI 3730 genetic Analyzer. Primers used in sequencing were the same as those used for PCR, and new primers were designed for walking sequencing. The sequenced fragments were assembled using DAMBE v 5  and BioEdit v 7.0.1 .
GenBank database and TRs analysis
The sequences of 1726 vertebrate mt genomes were retrieved from the NCBI genome database in March, 2012 (http://www.ncbi.nlm.nih.gov/genome). Repeated sequences were identified using Tandem Repeats Finder v4.03  with the following parameters: match +2, mismatch −7, indel −7. Only those repeats whose scores exceeded 50 were reported. Looser or stricter parameters were also tested in the same software packages, but these resulted in more paradoxical repeats and the loss of some repeats from the dataset.
The location of TRs in mt genomes from GenBank was manually checked, and TRs located in the CRs were classified by their location relative to the 3′ or 5′ end. Those TRs without a clear location in the CR were not used in the statistical analysis. Utilizing SPSS software, the distributions of motif lengths in the 3′ and 5′ ends and in protein coding genes were analyzed.
Results and discussion
The configurations of P. cornutusmtTRs
The sequences of M-TR motif in 20 individuals of P. cornutus
Hypervariable characteristics of 3′-terminal repeats
Comparisons of 109 3′-TRs showed that the length and the arrangement of 3′-TR and 3′-TR-v differ considerably not only between individuals, but also in different clones from single individuals (Figure 2). Extensive length polymorphism in these tracts is caused by differing motif copy number, ranging from four (DS-04, 32 bp) to 51 copies (QD-07, 408 bp) between individuals, and from four (DS-04, 32 bp) to 26 copies (DS-04, 208 bp) within individuals (Figure 2). A variable number of 3′-TR-v motifs are scattered throughout the 3′-TR tracts. The intra-individual differences are as follows: (i) One motif (3′-TR-v1) in DS-05(4), and two motifs (3′-TR-v1 and 3′-TR) in DS-11(3) are inserted. (ii) From one to six motifs are deleted at DS-13(6), QD-06(100), QD-07. (iii) While we can be sure that the 39th copy of the motif is deleted in QD-07(8, 16, and 174), site of deletion within other clones cannot be easily identified. For example, a deletion of the 3′-TR-v2 in QD-07(6) could have occurred at any site from the 5th to the 29th copy of the motif (see Figure 2).
The variation of copy number and motif organization between individuals QD-02, QD-03, and QD-04 is negligible as the difference of copy number is no more than one copy. TR diversity is seen between individuals, particularly in comparison to QD-07, as the motif organization and repeat length in this individual are entirely different, and the sequence of the 3′-TR-v1 motif TGTAATGT, is changed to 3′-TR-v2 (TATAATGT). The motif organization within a TR in individuals from different locations is not correlated with geographic location, indicating that TR organization is not geographically partitioned. For example, the TR organization of individuals QD-02, 03, 04 and DS-06, and 10 differ by a few 3′-TR motifs, but the TR organization of most other individuals vary to a great extent (Figure 2)
We suggest that: (i) There is remarkable sequence heterogeneity of 3′-TR tracts not only between individuals but also within individuals. (ii) Indel events occur randomly rather than at fixed sites within TR tracts. (iii) The indel unit of mtTRs could be one or multiple copies, suggesting that the indel tempo of mtTRs can be greatly increased by indel events of multiple motifs. (iv) There are two single base variants of 3′-TR:3′-TR-v1 and 3′-TR-v2, and the arrangement of 3′-TR and 3′-TR-v is highly diverse.
Hypervariable features cannot be explained by available models
These data suggest that the low frequency of recombination events could not lead to such high levels of mtTR variation in both organization and length, thus it is impossible to explain the hypervariable features of P. cornutus TRs by the recombination model alone. Although the SSM model could explain such frequent variation, this model more easily explains the formation of shorter simple motifs (microsatellites) than the formation of longer repeated motifs . If SSM is responsible for the mtTR variability in P. cornutus, simple shorter motifs would appear in higher abundance. For example, the repeated motif of 3′-TR-v2 (TA)2ATGT would include simple repeats, such as (TA)3, 4,…ATGT rather than (TATAATGT)2, 3,… alone. However, the fact is that no such simple repeats were found in any of the 109 TR fragments suggests that the SSM is not the most suitable model to account for the hypervariable features in P. cornutus mtTRs. Similarly, it seems the IEM model cannot explains our data as this mechanism only takes effect at the 5′ end of the CR, while we observed repeated tracts at the 3′ end of the CR in this flatfish. Another mitochondrial model, IIM, would yield identical repeated arrays in length and location from tRNA Phe to the origin of heavy strand replication, and this feature is not seen in P. cornutus mtTRs. Furthermore, the probability of false initiation of mt-replication from tRNA Phe , a condition for the IIM model, is not sufficiently high to explain the extremely hypervariable features in P. cornutus mtTRs. The model described by Hayasaka et al.  is suitable for the formation of TRs surrounding the CSB region, not a feature of the mtTRs in P. cornutus.
Repeated tracts in the mt genomes of other species such as dogs, wolves, and halibut [5, 31] share similar hypervariable features to those of P. cornutus. Thus far, no mechanistic model can explain the origin of the observed variation in these species.
The characteristics of tandem repeats in other vertebrate mt genomes
Characteristics of repeated tracts based in vertebrate mt genomes from NCBI
Species with TRs
The TR motifs at the 3′ end are typically of two sizes: 2–4 bp and 6–12 bp (Figure 3b); those at 5′ end are typically much longer (17 to 19 bp; Figure 3c). Further analyses show that most motifs at the 5′ end (>60%) contain TAS blocks, which are not observed in TRs at the 3′ end. The motif length in coding sequences is usually 12–18 bp, and tends to vary in multiples of three bases (12, 15, 18; Figure 3d). It is assumed that a leading reason for the distribution of 17–19 bp motifs at the 5′ end is due to the selection pressure on TASs, because the TAS length is approximately 15 bp , while selection pressure to maintain reading frame leads to a three-base variation in coding sequences based on the triplet code.
It is generally believed that the indel of short motifs (microsatellites) can be explained by the SSM model [13, 14]. The 1–4 bp repeats located at the CR 3′end (Figure 3b) may also be explained by this mechanism. There are, however, far more mtTRs with 6–22 bp motifs than those with shorter ones (Figure 3a). It is difficult to explain the generation of these larger motifs by the SSM model alone because the incidence of SSM decreases rapidly with increasing motif length. Additionally if long-motif mtTRs are generated by recombination, as happens in nuclear TRs , the low probability of mt recombination is not likely to generate the large number of mtTRs observed. Therefore, a new model of TR motif generation is needed.
A novel model of Pause-Melting Misalignment (PMM)
Currently, no model can satisfactorily explain the variability of motif indels in mtTRs, particularly outside the CR region; the uneven distribution of motif sizes; and the concentration of mtTRs at both termini in the CR. A novel model of Pause Melting Misalignment (PMM) is proposed.
The pause event of mitochondrial replication
During mt-replication, the pause event is the trigger for PMM. A pause event of mt-replication at the 5′ end in the CRs is likely since elongation of N-strands is arrested at TAS blocks at the 5′ end [33, 34]. Similarly, we predict frequent pause events at the 3′ end in the CRs, because the replication initiation and closure domain of N-strands exist at this domain, many replication events such as segregation of N-strand molecules and superhelical turns  occur while the replication forks are traversing this section. These cases are likely to create frequent pause events in the replication process.
In addition to frequent replication pauses at the 5′ and 3′ ends, Koike and Wolstenholme  point out that mtDNA replication is a discontinuous process. Based on the observation of different molecular forms of mtDNA in electron micrographs, they found that at least 44% of the mt genome contains discrete positions at which DNA synthesis tends to be arrested. By using a combination of one- and two-dimensional agarose gel electrophoresis and solution hybridization to strand-specific probes, Mayhook et al.  proved that replication pauses exist at ATP6/COIII gene positions. These mt-replication pauses could occur not only at the CR, but also at other positions including protein coding sequences, tRNAs, and rRNAs. These pause events could trigger the PMM process, and lead to the formation and motif indel of mtTRs.
Dynamic competition and melting
Folding and mispairing
Due to the competition of the D-strand, a portion of the N-strand is displaced from the T-strand, thus creating a short single-stranded status for both the N-strand and T-strand. These regions may then fold to form a stable intra-molecular secondary structure (Figure 4c), producing the shortened single strand of both strands, and then the reannealing of these regions could lead to mispairing. Induced stable secondary structures are thought to be a leading reason for mispairing during mtDNA replication [18, 19], however, in our vertebrate mtTR data, the lengths of many motifs are 6, 8 bp or even 2–5 bp. Motifs of these lengths are not known to form stable secondary structures. Even if some motifs are known to form hairpin structures, such hairpins are usually unstable due to their high free energy (exceeding the criterion value, -0.6 kcal/mol). Thus, secondary structure formation does not appear to be necessary for the TR indel events.
The generation of M-TR motifs in P. cornutus
The PMM model can readily explain key features of mtTRs in vertebrates
The PMM model differs from existing models of mtTR generation in several key aspects. Unlike the SSM model, the PMM model is specifically proposed for the origination and indel of motifs in mtTRs. It emphasizes the pause events as the trigger, followed by the melting of the N-strand from the T-strand, and subsequent folding and the mispairing. In contrast to the IEM and IIM models, the PMM model can explain the formation and indel of mtTRs in any locations of the mt genome. Furthermore, the PMM model can explain the origin of many characteristics of mtTRs in vertebrate mt genomes that the other models cannot.
In summary, the following characteristics of mtTRs can be explained by the PMM model: (i) The formation and motif indels of mtTR at any location in the mt genome can be produced by PMM mechanisms, even those outside the CR. Previous models do not address the production of indels at these locations. (ii) The concentration of large mtTRs at the 3′ and 5′ ends in the CRs is explained by the PMM model. Because there are stable pause events at both ends of the CR during mt-replication and the selection pressure upon the control region is less than that on other functional regions , additional TRs can be retained at these locations. (iii) Shorter motif lengths at the 3′ end could be due to reduced competition between the N-strand and the T-strand compared to that at the 5′ CR end. (iv) The generation of the M-TR motif in P. cornutus provides evidence for the origin of mtTRs by PMM mechanisms. (v) The PMM model can explain the hypervariable features of 3′-TRs in P. cornutus. Because mt-replication frequently pauses at the 3′ end of the CR, PMM would be triggered to generate a new copy or motif indels in existing mtTRs. The newborn motifs would subsequently increase the incidence of mispairing, which in turn would increase the frequency of motif indels. Furthermore, through the PMM model, more than one motif could be inserted or deleted in a single event, hence the generation tempo of motif indels greatly increases.
MtTRs of P. cornutus exhibit hypervariable characteristics, thus the 3′-TR-v motifs could act as molecular markers for the exploration of the newly generated of mtTRs as well as future indels. Based on analyses of the mtTRs in P. cornutus and vertebrate mt genomes, we developed a novel model of Pause Melting Misalignment for the formation and motif indel of tandem repeats in the mt genome. The model can explain the characters of TRs in available mt genomes, and we hope this model will provide an important basis for the explanation of repeated region evolution in mt genomes.
This work was supported by the Natural Science Foundation of China (30870283, 31071890). We thank Prof. Ziniu Yu and Dr. Xiangyun Wu for critical review and advice.
- Brown WM, George M, Wilson AC: Rapid evolution of animal mitochondrial DNA. Proc Natl Acad Sci U S A. 1979, 76 (4): 1967-1971. 10.1073/pnas.76.4.1967.PubMed CentralView ArticlePubMedGoogle Scholar
- McMillan WO, Palumbi SR: Rapid rate of control-region evolution in Pacific butterflyfishes (Chaetodontidae). J Mol Evol. 1997, 45 (5): 473-484. 10.1007/PL00006252.View ArticlePubMedGoogle Scholar
- White MM, Martin HR: Structure and conservation of tandem repeats in the mitochondrial DNA control region of the Least Brook lamprey (Lampetra aepyptera). J Mol Evol. 2009, 68 (6): 715-723. 10.1007/s00239-009-9246-1.View ArticlePubMedGoogle Scholar
- Mundy NI, Winchell CS, Woodruff DS: Tandem Repeats and Heteroplasmy in the Mitochondrial DNA Control Region of the Loggerhead Shrike (Lanius ludovicianus). J Hered. 1996, 87 (1): 21-26. 10.1093/oxfordjournals.jhered.a022948.View ArticlePubMedGoogle Scholar
- Mjelle KA, Karlsen BO, Jorgensen TE, Moum T, Johansen SD: Halibut mitochondrial genomes contain extensive heteroplasmic tandem repeat arrays involved in DNA recombination. BMC Genomics. 2008, 9: 10-10.1186/1471-2164-9-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Okada K, Yamazaki Y, Yokobori S, Wada H: Repetitive sequences in the lamprey mitochondrial DNA control region and speciation of Lethenteron. Gene. 2010, 465 (1–2): 45-52.View ArticlePubMedGoogle Scholar
- Ray DA, Densmore LD: Repetitive sequences in the crocodilian mitochondrial control region: Poly-A sequences and heteroplasmic tandem repeats. Mol Biol Evol. 2003, 20 (6): 1006-1013. 10.1093/molbev/msg117.View ArticlePubMedGoogle Scholar
- Piganeau G, Gardner M, Eyre-Walker A: A broad survey of recombination in animal mitochondria. Mol Biol Evol. 2004, 21 (12): 2319-2325. 10.1093/molbev/msh244.View ArticlePubMedGoogle Scholar
- Hoarau G, Holla S, Lescasse R, Stam WT, Olsen JL: Heteroplasmy and evidence for recombination in the mitochondrial control region of the flatfish Platichthys flesus. Mol Biol Evol. 2002, 19 (12): 2261-2264. 10.1093/oxfordjournals.molbev.a004049.View ArticlePubMedGoogle Scholar
- Lunt DH, Hyman BC: Animal mitochondrial DNA recombination. Nature. 1997, 387 (6630): 247-10.1038/387247a0.View ArticlePubMedGoogle Scholar
- Sammler S, Bleidorn C, Tiedemann R: Full mitochondrial genome sequences of two endemic Philippine hornbill species (Aves: Bucerotidae) provide evidence for pervasive mitochondrial DNA recombination. BMC Genomics. 2011, 12: 35-10.1186/1471-2164-12-35.PubMed CentralView ArticlePubMedGoogle Scholar
- Levinson G, Gutman GA: Slipped-strand mispairing: a major mechanism for DNA sequence evolution. Mol Biol Evol. 1987, 4 (3): 203-221.PubMedGoogle Scholar
- Torres-Cruz J, van der Woude MW: Slipped-strand mispairing can function as a phase variation mechanism in Escherichia coli. J Bacteriol. 2003, 185 (23): 6990-6994. 10.1128/JB.185.23.6990-6994.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Ellegren H: Microsatellite mutations in the germline: implications for evolutionary inference. Trends Genet. 2000, 16 (12): 551-558. 10.1016/S0168-9525(00)02139-9.View ArticlePubMedGoogle Scholar
- Mayer C, Leese F, Tollrian R: Genome-wide analysis of tandem repeats in Daphnia pulex - a comparative approach. BMC Genomics. 2010, 11 (227): 1-228.Google Scholar
- Dieringer D, Schlotterer C: Two distinct modes of microsatellite mutation processes: evidence from the complete genomic sequences of nine species. Genome Res. 2003, 13 (10): 2242-2251. 10.1101/gr.1416703.PubMed CentralView ArticlePubMedGoogle Scholar
- Leclercq S, Rivals E, Jarne P: DNA slippage occurs at microsatellite loci without minimal threshold length in humans: a comparative genomic approach. Genome Biol Evol. 2010, 2: 325-335. 10.1093/gbe/evq023.PubMed CentralView ArticlePubMedGoogle Scholar
- Buroker NE, Brown JR, Gilbert TA, O’Hara PJ, Beckenbach AT, Thomas WK, Smith MJ: Length heteroplasmy of sturgeon mitochondrial DNA: an illegitimate elongation model. Genetics. 1990, 124 (1): 157-163.PubMed CentralPubMedGoogle Scholar
- Broughton RE, Dowling TE: Length variation in mitochondrial DNA of the minnow Cyprinella spiloptera. Genetics. 1994, 138 (1): 179-190.PubMed CentralPubMedGoogle Scholar
- Taylor JS, Breden F: Slipped-strand mispairing at noncontiguous repeats in Poecilia reticulata: A model for minisatellite birth. Genetics. 2000, 155 (3): 1313-1320.PubMed CentralPubMedGoogle Scholar
- Hayasaka K, Ishida T, Horai S: Heteroplasmy and polymorphism in the major noncoding region of mitochondrial DNA in Japanese monkeys: association with tandemly repeated sequences. Mol Biol Evol. 1991, 8 (4): 399-415.PubMedGoogle Scholar
- Purdue JR, Oleksyk TK, Smith MH: Independent occurrences of multiple repeats in the control region of mitochondrial DNA of white-tailed deer. J Hered. 2006, 97 (3): 235-243.View ArticlePubMedGoogle Scholar
- Mancini E, De Biase A, Mariottini P, Bellini A, Audisio P: Structure and evolution of the mitochondrial control region of the pollen beetle Meligethes thalassophilus (Coleoptera: Nitidulidae). Genome. 2008, 51 (3): 196-207. 10.1139/G07-116.View ArticlePubMedGoogle Scholar
- Kuhn K, Streit B, Schwenk K: Conservation of structural elements in the mitochondrial control region of Daphnia. Gene. 2008, 420 (2): 107-112. 10.1016/j.gene.2008.05.020.View ArticlePubMedGoogle Scholar
- Eyre-Walker A, Awadalla P: Does human mtDNA recombine?. J Mol Evol. 2001, 53 (4–5): 430-435.View ArticlePubMedGoogle Scholar
- Ladoukakis ED, Zouros E: Direct evidence for homologous recombination in mussel (Mytilus galloprovincialis) mitochondrial DNA. Mol Biol Evol. 2001, 18 (7): 1168-1175. 10.1093/oxfordjournals.molbev.a003904.View ArticlePubMedGoogle Scholar
- Ciborowski KL, Consuegra S, Garcia de Leaniz C, Beaumont MA, Wang J, Jordan WC: Rare and fleeting: an example of interspecific recombination in animal mitochondrial DNA. Biol Lett. 2007, 3 (5): 554-557. 10.1098/rsbl.2007.0290.PubMed CentralView ArticlePubMedGoogle Scholar
- Xia X, Xie Z: DAMBE: software package for data analysis in molecular biology and evolution. J Hered. 2001, 92 (4): 371-373. 10.1093/jhered/92.4.371.View ArticlePubMedGoogle Scholar
- Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999, 41: 95-98.Google Scholar
- Benson G: Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999, 27 (2): 573-580. 10.1093/nar/27.2.573.PubMed CentralView ArticlePubMedGoogle Scholar
- Savolainen P, Arvestad L, Lundeberg J: mtDNA tandem repeats in domestic dogs and wolves: mutation mechanism studied by analysis of the sequence of imperfect repeats. Mol Biol Evol. 2000, 17 (4): 474-488. 10.1093/oxfordjournals.molbev.a026328.View ArticlePubMedGoogle Scholar
- Toth G, Gaspari Z, Jurka J: Microsatellites in different eukaryotic genomes: survey and analysis. Genome Res. 2000, 10 (7): 967-981. 10.1101/gr.10.7.967.PubMed CentralView ArticlePubMedGoogle Scholar
- Doda JN, Wright CT, Clayton DA: Elongation of displacement-loop strands in human and mouse mitochondrial DNA is arrested near specific template sequences. P Natl Acad Sci USA. 1981, 78 (10): 6116-6120. 10.1073/pnas.78.10.6116.View ArticleGoogle Scholar
- Shadel GS, Clayton DA: Mitochondrial DNA maintenance in vertebrates. Annu Rev Biochem. 1997, 66: 409-435. 10.1146/annurev.biochem.66.1.409.View ArticlePubMedGoogle Scholar
- Clayton DA: Replication of animal mitochondrial DNA. Cell. 1982, 28 (4): 693-705. 10.1016/0092-8674(82)90049-6.View ArticlePubMedGoogle Scholar
- Koike K, Wolstenholme DR: Evidence for discontinuous replication of circular mitochondrial DNA molecules from Novikoff rat ascites hepatoma cells. J Cell Biol. 1974, 61 (1): 14-25. 10.1083/jcb.61.1.14.PubMed CentralView ArticlePubMedGoogle Scholar
- Mayhook AG, Rinaldi AM, Jacobs HT: Replication origins and pause sites in sea urchin mitochondrial DNA. Proc Biol Sci. 1992, 248 (1321): 85-94. 10.1098/rspb.1992.0046.View ArticlePubMedGoogle Scholar
- Haasnoot CAG, Debruin SH, Hilbers CW, Vandermarel GA, Vanboom JH: Loopstructures in synthetic oligonucleotides - hairpin stability and structure studied as a function of loop elongation. J Biosci (Bangalore). 1985, 8 (3–4): 767-780.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.