Association of G-quadruplex forming sequences with human mtDNA deletion breakpoints
© Dong et al.; licensee BioMed Central Ltd. 2014
Received: 11 April 2014
Accepted: 4 August 2014
Published: 13 August 2014
Mitochondrial DNA (mtDNA) deletions cause disease and accumulate during aging, yet our understanding of the molecular mechanisms underlying their formation remains rudimentary. Guanine-quadruplex (GQ) DNA structures are associated with nuclear DNA instability in cancer; recent evidence indicates they can also form in mitochondrial nucleic acids, suggesting that these non-B DNA structures could be associated with mtDNA deletions. Currently, the multiple types of GQ sequences and their association with human mtDNA stability are unknown.
Here, we show an association between human mtDNA deletion breakpoint locations (sites where DNA ends rejoin after deletion of a section) and sequences with G-quadruplex forming potential (QFP), and establish the ability of selected sequences to form GQ in vitro. QFP contain four runs of either two or three consecutive guanines (2G and 3G, respectively), and we identified four types of QFP for subsequent analysis: intrastrand 2G, intrastrand 3G, duplex derived interstrand (ddi) 2G, and ddi 3G QFP sequences. We analyzed the position of each motif set relative to either 5' or 3' unique mtDNA deletion breakpoints, and found that intrastrand QFP sequences, but not ddi QFP sequences, showed significant association with mtDNA deletion breakpoint locations. Moreover, a large proportion of these QFP sequences occur at smaller distances to breakpoints relative to distribution-matched controls. The positive association of 2G QFP sequences persisted when breakpoints were divided into clinical subgroups. We tested in vitro GQ formation of representative mtDNA sequences containing these 2G QFP sequences and detected robust GQ structures by UV–VIS and CD spectroscopy. Notably, the most frequent deletion breakpoints, including those of the "common deletion", are bounded by 2G QFP sequence motifs.
The potential for GQ to influence mitochondrial genome stability supports a high-priority investigation of these structures and their regulation in normal and pathological mitochondrial biology. These findings emphasize the potential importance of helicases that subsequently resolve GQ to maintain the stability of the mitochondrial genome.
KeywordsG-quadruplex mtDNA deletions Mitochondrial disease Mitochondrial genome instability Non-B DNA Nucleic acid structures
The mitochondrial genome (mtDNA) is a multicopy DNA molecule that encodes essential components of the respiratory chain. A broad spectrum of human diseases [1, 2] stem from mtDNA mutations, which cause mitochondrial dysfunction and multisystem disorders [3, 4]. The most severe mtDNA mutations are deletions, which cause respiratory dysfunction at a lower mutation load than point mutations (reviewed in ). Based on epidemiological studies, single mtDNA deletions cause ~13-30% of primary mitochondrial disease [3, 5] and can arise from matrilineal transmission or sporadically during oogenesis. Multiple mtDNA deletions may also arise in a small proportion of patients afflicted with mitochondrial disease , and frequently occur due to mutation in one of several nuclear genes that influence mtDNA synthesis . Multiple mtDNA deletions also occur in somatic tissue of all individuals and accumulate with age [8–11].
MtDNA deletion breakpoints show non-random distributions that presumably reflect the mechanisms underlying the formation of deletions. Several studies have attempted to identify sequence-related mechanisms for mtDNA deletion formation. The first sequence features identified were short direct repeats that contain or flank the deletion junction [12–14]. MtDNA deletion breakpoints can contain direct repeats of 3–13 bp, but these motifs are not universal features of deletions . Direct repeats remain the most commonly hypothesized sequence-based mechanism for forming mtDNA deletions (ex. ), but leave unexplained a significant proportion of these mutations.
MtDNA contains numerous sequences with the potential to form non-B DNA structures, which may interfere with efficient DNA replication and repair. Supporting this notion, mtDNA fragments that contain deletion hot spots show altered electrophoretic mobility , possibly due to formation of stable non-B form secondary structures. Sequences predicted to form stable stem-loop or cruciform structures in single-stranded DNA overlap only a subset of the total documented 5′ and 3′ deletion breakpoints [17, 18], indicating the potential involvement of other non-B form structures in mtDNA instability.
Among non-B forms of DNA, G-quadruplex (GQ) structures exhibit characteristics that could implicate them in mitochondrial genome instability. GQs are stable secondary structures composed of two or more stacked planar guanine tetrads stabilized by cations such as potassium. These structures can occur in a single strand of DNA or RNA (intramolecular) or among two or more single-strands (intermolecular) . A subset of intermolecular GQs has been speculated to form between runs of guanines derived from each strand of duplex DNA (duplex-derived intermolecular) . In the nucleus, GQ structures can hinder DNA replication  and cause genome instability [21, 22], and sequences with quadruplex forming potential (QFP) have been associated with deletions and duplications in the genomic DNA of cancer cells . Sequences with QFP are prevalent in yeast mtDNA  and the promoters of prokaryotic genes , but little is known about QFP sequences in the mammalian mitochondrial genome. In mammals, mitochondrial RNA containing the conserved sequence block 2 (CSB2) forms a GQ in vitro . This GQ structure is predicted to affect the formation of a stable RNA primer [27, 28] used both in the formation of D-loop structures thought to associate mtDNA with the inner mitochondrial membrane [29, 30], and in priming first strand mtDNA replication initiation . To date, only three-guanine (3G) QFP have been described in mtDNA [24, 26, 32].
The abundance of QFP sequences in yeast mtDNA  and potentially in mammals (this study), coupled with described effects of GQ on nuclear genome stability, make GQ structures of particular interest as potential agents in human mtDNA deletion formation. The potential for GQ structures to associate with mtDNA deletions has previously been suggested for a limited number of 3G QFP . Here, we performed an examination of multiple categories of QFP sequences in human mtDNA. We analyzed enrichment across small distances between QFP and 5′ or 3′ breakpoints, and developed a reciprocal proximity method to analyze QFP and other sequence motifs for systematic association with 5′ and 3′ mtDNA breakpoints. Our results demonstrate that mtDNA QFP sequences show significant association with deletion breakpoints and are enriched with breakpoints over short distances at least as well as direct repeat sequences. All other structures tested failed to show such association. Not previously reported, the two-guanine (2G) QFP motif set is the only QFP motif that is systematically associated with both 5′ and 3′ breakpoints. We further show the capacity of selected 2G QFP to form GQ structures in vitro.
Identification of QFP sequences
Identification of stem-loop/cruciform structures and repeat sequences
The prediction of stem-loop, cruciform and other hairpin-containing features across human mtDNA was carried out using the hybrid-ss-min core program of the UNAFold 3.8 software package . For the sake of simplicity, the abbreviation SC will be used to define any of these different structural elements. The folding of single-stranded DNA was simulated at 37°C in 1 M sodium and no magnesium, using 100 nucleotides as the maximum distance between paired bases in each structure. The circular nature of the mtDNA was considered in the folding prediction. Python scripts were created for the identification of direct repeats (e.g., ATC-ATC), inverted repeats (e.g., ATC-CTA), complementary repeats (e.g., ATC-TAG) and inverted complementary repeats (e.g., ATC-GAT) in human mtDNA. Midpoints were used in all calculations. The circular representation of mtDNA was made using the Circos software package, version 0.62 .
Minimal Distance Analysis (MDA)
Individual enrichment analysis
Thus, a small p-value signifies an enrichment. Similarly, we can evaluate the enrichment of the number of structural motifs at close distances to mtDNA deletion breakpoints.
Enrichment of motifs with both mtDNA deletion breakpoints
Simultaneously for both breakpoints, mtDNA deletions could be associated with 2G QFP or direct repeat sequences. To test this, the distance of the nearest 2G QFP sequence to each end of a deletion is measured. If both of the distances are less than or equal to 10 nt, the deletion is considered to be close to 2G QFP motifs. The actual number of deletions close to 2G QFP motifs is compared to the control, which is calculated by rotating the 2G QFP set relative to the deletion set for all controls. The p-value is calculated from binomial distribution, as is the case for single-end enrichment. For direct repeat sequences, this is performed by two types of calculations. The first type is identical to that done for 2G QFP sequences, by considering all direct repeat sequences as equivalent entries. For the second type, each deletion is compared with each direct repeat pair. The distance from the 5′ end of the deletion to the 5′ end of the direct repeat pair, and the distance from the 3′ end of the same deletion to the 3′ end of the same direct repeat pair, are calculated. If both distances are less than or equal to 10 nt, the deletion is considered to be close to the direct repeat pair.
All reagents were molecular biology grade or better from Sigma (Sigma-Aldrich, St. Louis, MO) and oligonucleotides were synthesized by IDT (Coralville, IA) with standard purification. All spectroscopic analyses were performed as previously reported [39, 40]. In brief, ~5.0 μM DNA samples were annealed by heating to 95°C for 5 min followed by slow cooling to room temperature. Thermal difference spectra (TDS) were collected on a Cary 300 Varian Spectrophotometer equipped with a Peltier-thermostated cuvette holder. Circular dichroism (CD), thermal melting and wavelength scans were performed on an AVIV 410 spectropolarimeter equipped with a Peltier heating unit. While monitoring the signal at 264 nm, DNA samples were heated from 5 to 95°C at a rate of 10°C/min with an equilibration time of 1 min, an averaging time of 10 s, and a bandwidth of 1 nm in 10 mM lithium cacodylate buffer, pH 7.2, supplemented with 150 mM KCl. The data were analyzed using a two-state Van’t Hoff equation assuming constant enthalpy, ΔH. CD wavelength scans were collected on samples both before and after melting. Three scans, over a wavelength range of 220 to 350 nm, were recorded at 4°C with a bandwidth of 1 nm, an averaging time of 1 s, and a wavelength step of 1 nm. These scans were averaged, baseline corrected, smoothed and converted to molar ellipticity. All data analysis was done in Origin 8.1 (OriginLab, Northampton, MA).
Identification of four types of QFP sequences in human mtDNA
We began our investigation by generating a list of human mtDNA sequences with the potential to form G-quadruplexes in single-stranded (intrastrand) or double-stranded (duplex-derived interstrand, or ddi) DNA. Algorithms designed to predict QFP sequences have been developed for both single-stranded (reviewed in ) and double-stranded sequences . The two mtDNA strands have asymmetric nucleotide distribution, with adenine and guanine purines enriched on the heavy strand. As described in Methods, we identified 2G and 3G ddi QFP sequences, and intrastrand 2G QFP and 3G QFP sequences on the heavy strand of human mtDNA. We found only 2G QFP on the light strand due to the scarcity of guanine residues there. Strikingly, the number of QFP sequences in human mtDNA was significantly higher than the number found in random genomes with the identical nucleotide content and asymmetry (not shown). The positions of the non-overlapping intrastrand QFP sequences are shown in Figure 2 relative to a reference map for human mtDNA (NC_012920).
2G QFP and 3G QFP motifs, but not ddi QFP, are associated with human mtDNA deletion breakpoints
We next tested the null hypothesis that 5′ or 3′ mtDNA deletion breakpoints (b-p) were randomly positioned relative to QFP sequence sets (motif). We selected a set of 730 human mtDNA deletions that included associated clinical data , and identified all unique 5′ and 3′ deletion breakpoints in this collection (Figure 1). The positional enrichment of 3G QFP and 2G QFP sequences relative to unique mtDNA deletion breakpoints is shown in Figure 2A. Because QFP sequences could lie inside or outside the deleted sequence, we calculated the minimal distance in either direction from each breakpoint in a set (i.e. 5′ or 3′ b-p set) to the nearest QFP sequence midpoint (b-p to motif). To establish distribution-matched controls, the midpoint positions for a QFP sequence set were rotated around the circular mtDNA genome in 1 nt increments to generate 16,568 control data sets. In this minimal distance analysis (MDA), p-values were calculated by determining the frequency at which a control set’s average minimal distance from b-p to motif was smaller than or equal to the value for the actual data set (examples shown in Figure 2B,C). As such, an MDA p-value of 0.001 (b-p to motif) indicates that the average distance from a breakpoint to nearest QFP sequence motif was smaller than controls for all but 16 positions in the 16 kb genome. Furthermore, our controls are more restrictive than other randomized control sets , where changes in randomized distribution of motifs or breakpoints is sufficient to reject the null hypothesis.
Quadruplex forming potential (QFP) and other structural motifs analyzed for minimum distance relationship to 5’ or 3’ human mtDNA deletion breakpoints
Actual/control average distance
b-p to motif
Motif to b-p
b-p to motif
Motif to b-p
inverted comp. rpts.
Direct repeat motifs, but not other DNA motifs, are associated with human mtDNA breakpoints
In addition to QFP motifs, we investigated the association of other DNA motifs with deletions by performing MDA on several types of sequence features. Stem-loop or cruciform (SC) structures predicted by the UNAFold algorithm are thought to associate with a limited number of the most common mtDNA deletion breakpoints . Using UNAFold, we generated a list of predicted SC sequences and found that these sequences occur frequently (Table 1). However, MDA showed no association of SC structures with either 5′ or 3′ breakpoints. We conclude that these sequences identified by UNAFold, treated equally, do not explain the diverse population of deletion breakpoints, which is consistent with findings of others . SC sequences have been shown, however, to strongly associate with specific deletion breakpoints among the most common deletions [17, 18], which supports the notion that multiple mechanisms may contribute to mtDNA deletion formation.
Additional DNA sequence motifs have been suggested to contribute to mtDNA instability. Recently, association studies have suggested a negative correlation for inverted complementary repeat sequences in mtDNA in maximum life span across species, possibly through increased deletion mutagenesis . Furthermore, direct repeats of 11–15 bp have been frequently suggested to contribute to deletion breakpoint sites [12–14, 17, 43]. Direct repeat sequence pairs, inverted complementary repeat pairs, and individual stem-loop/cruciform structures are visualized (Additional file 3: Figure S2) to illustrate potential associations of these sequences with both 5′ and 3′ breakpoints in mtDNA deletions.
We next tested direct and inverted complementary repeat sets for association with mtDNA deletion breakpoints, followed by complementary and inverted repeats, which are not expected to basepair elsewhere in the genome. We selected an 11 nt repeat length because it was long enough to limit the number of triplet repeats but short enough to give the largest number of repeats within the range of perfect repeats associated with disease [12–14, 17, 43]. For comparison by MDA, each component of a pair was analyzed as a separate unit. This repeat pair uncoupling allowed us to test the null hypothesis that sequences capable of basepairing elsewhere in the genome are unassociated with 5′ or 3′ deletion breakpoints. We tested this hypothesis with b-p to motif and motif to b-p MDA (differences between MDA are illustrated in Figure 1B,C) and found significant association between direct repeat sequences and both 5′ and 3′ breakpoints in the b-p to motif MDA (Table 1). Although the p-values (p = 0.04, 0.03 for 5′ and 3′, respectively) are not as small as those for 2G QFP, we do believe that the association is significant, given that the p-values are much smaller in the enrichment studies when 5′ and 3′ ends are considered individually (see next section) and simultaneously (see Discussion). Inverted complementary repeats showed no significant breakpoint association, nor did complementary and inverted repeat sequences. The location of 2G QFP, 3G QFP, and direct repeat sequences relative to paired breakpoints in the mtDNA deletion set are shown in Figure 2A. The complete table of all p-values for all motifs and repeats is found in Additional file 2: Table S1. It should be noted that the p-values represent the significance of enrichment between motif and breakpoint relative to control data sets, but they do not allow qualitative comparison of two motifs with statistically significant enrichment.
Human mtDNA deletion breakpoints are enriched near 2G QFP, 3G QFP, and direct repeat motifs
Association of 2G QFP with the "common deletion"
Association of motifs with clinical categories of human mtDNA deletions
To establish how motifs associate with mtDNA deletions from specific clinical categories, we examined all QFP and repeat sequences for association with seven groups of recently described mtDNA deletions . The so-called "common deletion" was the only deletion found in all seven groups. For each category, the distribution of 5’ and 3’ deletion breakpoints across mtDNA is shown by frequency histogram in Additional file 5: Figure S3. By MDA, 2G QFP sequences showed significant association with the 3′ breakpoints of the majority of clinical categories (six of seven), while 5’ breakpoints were associated with 2G QFP sequences in two of the categories (Additional file 6: Table S3A,B). Both 3G ddi QFP and 2G ddi QFP sequences showed association with individual categories, but never with both 5′ and 3′ breakpoints in the same category. Direct repeat sequences showed significant association with 3′ breakpoints from four of seven groups, but only one category of 5′ breakpoints. SC sequences and other control sequences failed to show broad association with any individual clinical category of deletions (Additional file 2: Table S1).
We noted that 2G QFP sequences showed nearly significant association with 5′ breakpoints in patients with single mtDNA deletions. To better understand this, we identified disease-specific mtDNA deletions from within this clinical group and determined the MDA values for 2G QFP sequence association with 5′ and 3′ breakpoints (Additional file 6: Table S3C,D). The distribution of mtDNA deletion breakpoints for Kearns-Sayre syndrome (KSS), progressive external opthalmoplegia (PEO), and Pearson’s syndrome (PS) is shown in a frequency histogram in Additional file 7: Figure S4. MDA showed significant association of 2G QFP sequences with KSS and PEO deletion breakpoints but not those found in PS. Comparisons of mtDNA deletion locations among these populations showed limited overlap between PS and both KSS and PEO sets, but extensive overlap between KSS and PEO sets (Additional file 8: Figure S5). These results suggest that 2G QFP sequences are in close association with at least a subset of single mtDNA deletion breakpoints.
2G QFP sequences form bona fide G-quadruplexes in vitro
Summary of DNA biophysical characterization of representative 2G QFP and corresponding control sequences
GGA TGG GGT GGG GAG G
AGA TGG AGT GGA GAG G
GGG GGA TGC GGG GG
AGG AGA TGC AGG AG
GGA GGG TGG ATG G
AGA GGG TAG ATG G
This study is the first to establish 2G QFP sequence prevalence in mammalian mtDNA, demonstrate their association with mtDNA breakpoints, and illustrate the capacity of 2G QFP mtDNA sequences to form quadruplex secondary structures. Because 3G QFP sequences have also been found within a few nucleotides of nuclear DNA deletion breakpoints , we suggest that both 2G and 3G GQ could contribute broadly to deletion formation in humans through a mechanism common to both nuclear and mitochondrial genomes.
Our work also significantly improves upon previous findings that 3G QFP are associated with mtDNA deletions . The Oliveira study utilized a different QFP-determining protocol to identify five 3G QFP sequences (we report nine) and examined the relative enrichment of 5′ + 3′ pooled breakpoints identified from the Mitomap database (http://www.mitomap.org) to these motif sequences. Two of the 3G QFP were located in the regions with prevalent 5′ or 3′ breakpoints and thus showed a high density of breakpoints within 100 bp, but not with random or partially randomized breakpoint sets. The authors concluded that these two 3G QFP were significantly associated with breakpoints; however, they did not exclude the effects of the non-random distribution of the breakpoints in their random controls. Our study improves upon these findings by expanding the 3G QFP sequences tested, using non-redundant separated 5′ and 3′ breakpoints, examining association from the motif and from the breakpoint, and using distribution matched controls.
One attractive potential mechanism for GQ contribution to deletion formation is that GQ structures contribute to mtDNA deletions during replication. Strong evidence supports replication-dependent deletion formation in the nucleus of cancer cells and yeast. The separation of DNA strands during replication is thought to leave the lagging strand vulnerable to GQ structure formation . In yeast, when GQ structures are encountered by the replisome, DNA elongation is slowed . The hallmarks of replication pausing were reversed by expression of PIF1 helicase, which has strong G-quadruplex resolving activity . Repeat tandem arrays containing 3G QFP sequences are inherently unstable in yeast, but this sequence instability is prevented by expression of PIF1 . The absence of PIF1, or exposure to G-quadruplex stabilizing compounds (ligands), slows replication and increases sequence instability at GQ repeat sequences . Finally, the sites of nuclear genome instability in cancer cells frequently contain QFP, with a strong presence of 2G QFP sequences .
In mtDNA, the displaced heavy strand contains the majority of 2G QFP sequences (this study), which could be exposed by structural transitions induced by transcription, replication, or topological stress. It is worth noting that heavy strand DNA can be found at low levels as a single stranded DNA-RNA hybrid [50, 51] or bound by mtSSB , whose displacement or processing would allow exposed single-stranded DNA to fold into intramolecular GQ structures. RNA-DNA hybrids containing a GQ structure in the DNA have been suggested as mechanism of genome instability . However they might arise, unresolved GQ structures in mtDNA may lead to collapsed replication forks, similar to those observed in the nucleus near stabilized GQ sequences . If failed replication intermediates are repaired rather than degraded, the opportunity may arise to create recombinant mtDNA molecules . We suggest that the increased density of predicted GQ structures in a region such as the boundaries of the common deletion or the deletion "hot spot" near the D-loop 3’ end [17, 18] (Figure 4) would increase the incidence of replication stalling, leading to the formation of single-stranded DNA or double-strand breaks, the latter of which are known to cause mtDNA deletion formation [55, 56] presumably through a non-homologous end joining (NHEJ) repair mechanism. SC structures have also been postulated to cause replication pausing, and have been associated with specific deletion breakpoints .
Circumstances that enhance the likelihood of replication fork collapse could increase mtDNA deletion formation rates. Besides structural interference by folded GQ DNA, nucleotide imbalances, defective polymerase, or oxidative damage could contribute to premature replication termination, the repair of which could lead to mtDNA deletions. Moreover, gene mutation or changes in expression levels during aging could also contribute to mtDNA deletion formation. Consistent with a role for replication in mtDNA deletion formation, the few known genes whose mutation is known to cause mtDNA deletion can be grouped into core replication subunits (POLG, POLG2, and Twinkle), nucleotide metabolism (TP and RRM2B), and nucleotide transport (ANT1, DNC, and potentially MPV17 ) (reviewed in [58, 59]). Altogether, these findings support replication pausing as a potential first step in mtDNA deletion formation.
Although the significant association of QFP sequences with breakpoints inspires the mechanisms proposed above, our results do not exclude others. The QFP sequences might also stimulate deletion formation during DNA repair through the formation of distant trans interactions. In addition, QFP might work in parallel with direct repeats or SC sequences, which have long been proposed to facilitate deletion formation during DNA repair [18, 54]. Interestingly, among all of the DNA motifs tested in this study, direct repeat is the only other sequence type significantly associated with breakpoints. To further demonstrate the significance of the association, we asked what the chances are that both ends of a deletion are in close vicinity to QFP or direct repeat sequences. Both p-values are small, p < 10-14 (QFP) and p < 10-12 (direct repeat) for the number of deletions with both ends within a 10 nt distance to the midpoint of the motifs (Additional file 9: Figure S6). For direct repeat sequences, we further determined the association of a deletion with a direct repeat pair, which decreases the p-value to <10-14 for the 10 nt distance. These findings strongly support the significant association of both QFP and direct repeat sequences with individual, or in the case of direct repeats, paired, deletion breakpoints.
GQ structures are surprisingly stable, and their resolution requires the action of specific proteins. In the nucleus, ATP-dependent helicases are dedicated to resolving these structures, and replication through QFP sequences is decreased in their absence [22, 48]. The presence of dedicated activities to manage the resolution of these structures suggests that GQ cause replication interference. In the nucleus, stabilization of GQ triggers the ataxia telangiectasia and Rad3-related protein (ATR) DNA damage response, consistent with single-stranded DNA formation at stalled replication forks  and exposure of single-stranded DNA gaps behind the replication fork . As multiple different nuclear helicases are required to maintain nuclear genome stability and thus prevent disease, the maintenance of mtDNA stability likely requires multiple helicase activities in mitochondria, including Twinkle primase/helicase (PEO1) [62, 63], the 3′-5′ helicase RECQL4 , and the 5′-3′ helicase PIF1 [65–67]. Future studies describing interplay among replication, DNA repair, helicases, and GQ-forming sequences are likely to provide important insights into mitochondrial genome stability.
In summary, our findings support the hypothesis that sequences predicted to form GQ structures in mtDNA are enriched with mtDNA deletion breakpoints. Among these multiple types of QFP, the 2G QFP motif set are in closest proximity to the largest number of breakpoints and overlap the most common breakpoints. The prediction that sequences of this motif set have the ability to form GQ was validated on a subset of sequences using in vitro structural analysis. Together, these findings suggest that GQ structures could be a major contributor to human mtDNA instability. Further experiments are needed to confirm these structures in vivo and establish their contribution to normal and pathological mitochondrial biology.
We would like to thank Sion Williams for his thoughtful comments on this manuscript. DWD is supported by a grant from the University Research Foundation, University of Pennsylvania. FBJ was supported by NIH R01-AG021521 and FBJ and KC were supported by NIH P01-AG031862. BAK and JEK were supported by NIH P30-DK050306. FP and JD are supported by the Portuguese Foundation for Science and Technology (FCT), Fundo Social Europeu and Programa Operacional Potencial Humano [PTDC/CVT/100881/2008 and Investigator FCT program]; LAY was supported by funds from Swarthmore College.
- Greaves LC, Reeve AK, Taylor RW, Turnbull DM: Mitochondrial DNA and disease. J Pathol. 2011, 226: 274-286.PubMedView ArticleGoogle Scholar
- Schon EA, DiMauro S, Hirano M: Human mitochondrial DNA: roles of inherited and somatic mutations. Nat Rev Genet. 2012, 13: 878-890.PubMed CentralPubMedView ArticleGoogle Scholar
- Schaefer AM, McFarland R, Blakely EL, He L, Whittaker RG, Taylor RW, Chinnery PF, Turnbull DM: Prevalence of mitochondrial DNA disease in adults. Ann Neurol. 2008, 63: 35-39.PubMedView ArticleGoogle Scholar
- Schaefer AM, Taylor RW, Turnbull DM, Chinnery PF: The epidemiology of mitochondrial disorders–past, present and future. Biochim Biophys Acta. 2004, 1659: 115-120.PubMedView ArticleGoogle Scholar
- Chinnery PF, Johnson MA, Wardell TM, Singh-Kler R, Hayes C, Brown DT, Taylor RW, Bindoff LA, Turnbull DM: The epidemiology of pathogenic mitochondrial DNA mutations. Ann Neurol. 2000, 48: 188-193.PubMedView ArticleGoogle Scholar
- Hudson G, Deschauer M, Taylor RW, Hanna MG, Fialho D, Schaefer AM, He L-P, Blakely E, Turnbull DM, Chinnery PF: POLG1, C10ORF2, and ANT1 mutations are uncommon in sporadic progressive external ophthalmoplegia with multiple mitochondrial DNA deletions. Neurology. 2006, 66: 1439-1441.PubMedView ArticleGoogle Scholar
- Spinazzola A, Zeviani M: Disorders from perturbations of nuclear-mitochondrial intergenomic cross-talk. J Intern Med. 2009, 265: 174-192.PubMedView ArticleGoogle Scholar
- Bender A, Krishnan KJ, Morris CM, Taylor GA, Reeve AK, Perry RH, Jaros E, Hersheson JS, Betts J, Klopstock T, Taylor RW, Turnbull DM: High levels of mitochondrial DNA deletions in substantia nigra neurons in aging and Parkinson disease. Nat Genet. 2006, 38: 515-517.PubMedView ArticleGoogle Scholar
- Linnane AW, Marzuki S, Ozawa T, Tanaka M: Mitochondrial DNA mutations as an important contributor to ageing and degenerative diseases. Lancet. 1989, 1: 642-645.PubMedView ArticleGoogle Scholar
- Herbst A, Pak JW, McKenzie D, Bua E, Bassiouni M, Aiken JM: Accumulation of mitochondrial DNA deletion mutations in aged muscle fibers: evidence for a causal role in muscle fiber loss. J Gerontol A Biol Sci Med Sci. 2007, 62: 235-245.PubMed CentralPubMedView ArticleGoogle Scholar
- Kraytsberg Y, Simon DK, Turnbull DM, Khrapko K: Do mtDNA deletions drive premature aging in mtDNA mutator mice?. Aging Cell. 2009, 8: 502-506.PubMed CentralPubMedView ArticleGoogle Scholar
- Schon EA, Rizzuto R, Moraes CT, Nakase H, Zeviani M, DiMauro S: A direct repeat is a hotspot for large-scale deletion of human mitochondrial DNA. Science. 1989, 244: 346-349.PubMedView ArticleGoogle Scholar
- Mita S, Rizzuto R, Moraes CT, Shanske S, Arnaudo E, Fabrizi GM, Koga Y, DiMauro S, Schon EA: Recombination via flanking direct repeats is a major cause of large-scale deletions of human mitochondrial DNA. Nucleic Acids Res. 1990, 18: 561-567.PubMed CentralPubMedView ArticleGoogle Scholar
- Samuels DC, Schon EA, Chinnery PF: Two direct repeats cause most human mtDNA deletions. Trends Genet. 2004, 20: 393-398.PubMedView ArticleGoogle Scholar
- Lakshmanan LN, Gruber J, Halliwell B, Gunawan R: Role of direct repeat and stem-loop motifs in mtDNA deletions: cause or coincidence?. PLoS One. 2012, 7: e35271-PubMed CentralPubMedView ArticleGoogle Scholar
- Hou JH, Wei YH: The unusual structures of the hot-regions flanking large-scale deletions in human mitochondrial DNA. Biochem J. 1996, 318 (Pt 3): 1065-1070.PubMed CentralPubMedView ArticleGoogle Scholar
- Zeviani M, Servidei S, Gellera C, Bertini E, DiMauro S, DiDonato S: An autosomal dominant disorder with multiple deletions of mitochondrial DNA starting at the D-loop region. Nature. 1989, 339: 309-311.PubMedView ArticleGoogle Scholar
- Damas J, Carneiro J, Gonçalves J, Stewart JB, Samuels DC, Amorim A, Pereira F: Mitochondrial DNA deletions are associated with non-B DNA conformations. Nucleic Acids Res. 2012, 40: 7606-7621.PubMed CentralPubMedView ArticleGoogle Scholar
- Bochman ML, Paeschke K, Zakian VA: DNA secondary structures:stability and function ofG-quadruplex structures. Nat Rev Genet. 2012, 13: 770-780.PubMed CentralPubMedView ArticleGoogle Scholar
- Cao K, Ryvkin P, Johnson FB: Computational detection and analysis of sequences with duplex-derived interstrand G-quadruplex forming potential. Methods. 2012, 57: 3-10.PubMed CentralPubMedView ArticleGoogle Scholar
- Lopes J, Piazza A, Bermejo R, Kriegsman B, Colosio A, Teulade-Fichou M-P, Foiani M, Nicolas A: G-quadruplex-induced instability during leading-strand replication. EMBO J. 2011, 30: 4033-4046.PubMed CentralPubMedView ArticleGoogle Scholar
- Ribeyre C, Lopes J, Boulé J-B, Piazza A, Guédin A, Zakian VA, Mergny J-L, Nicolas A: The yeast Pif1 helicase prevents genomic instability caused by G-quadruplex-forming CEB1 sequences in vivo. PLoS Genet. 2009, 5: e1000475-PubMed CentralPubMedView ArticleGoogle Scholar
- De S, Michor F: DNA secondary structures and epigenetic determinants of cancer genome evolution. Nat Struct Mol Biol. 2011, 18: 950-955.PubMed CentralPubMedView ArticleGoogle Scholar
- Capra JA, Paeschke K, Singh M, Zakian VA: G-quadruplex DNA sequences are evolutionarily conserved and associated with distinct genomic features in Saccharomyces cerevisiae. PLoS Comput Biol. 2010, 6: e1000861-PubMed CentralPubMedView ArticleGoogle Scholar
- Rawal P, Kummarasetti VBR, Ravindran J, Kumar N, Halder K, Sharma R, Mukerji M, Das SK, Chowdhury S: Genome-wide prediction of G4 DNA as regulatory motifs: role in Escherichia coli global regulation. Genome Res. 2006, 16: 644-655.PubMed CentralPubMedView ArticleGoogle Scholar
- Wanrooij PH, Uhler JP, Simonsson T, Falkenberg M, Gustafsson CM: G-quadruplex structures in RNA stimulate mitochondrial transcription termination and primer formation. Proc Natl Acad Sci U S A. 2010, 107: 16072-16077.PubMed CentralPubMedView ArticleGoogle Scholar
- Lee DY, Clayton DA: Properties of a primer RNA-DNA hybrid at the mouse mitochondrial DNA leading-strand origin of replication. J Biol Chem. 1996, 271: 24262-24269.PubMedView ArticleGoogle Scholar
- Wanrooij PH, Uhler JP, Shi Y, Westerlund F, Falkenberg M, Gustafsson CM: A hybrid G-quadruplex structure formed between RNA and DNA explains the extraordinary stability of the mitochondrial R-loop. Nucleic Acids Res. 2012, 40: 10334-10344.PubMed CentralPubMedView ArticleGoogle Scholar
- Albring M, Attardi G: Association of a protein structure of probable membrane derivation with HeLa cell mitochondrial DNA near its origin of replication. Proc Natl Acad Sci U S A. 1977, 74: 1348-1352.PubMed CentralPubMedView ArticleGoogle Scholar
- Kopek BG, Shtengel G, Xu CS, Clayton DA, Hess HF: Correlative 3D superresolution fluorescence and electron microscopy reveal the relationship of mitochondrial nucleoids to membranes. Proc Natl Acad Sci U S A. 2012, 109: 6136-6141.PubMed CentralPubMedView ArticleGoogle Scholar
- Chang DD, Clayton DA: Priming of human mitochondrial DNA replication occurs at the light-strand promoter. Proc Natl Acad Sci U S A. 1985, 82: 351-355.PubMed CentralPubMedView ArticleGoogle Scholar
- Oliveira PH, da Silva CL, Cabral JMS: An appraisal of human mitochondrial DNA instability: new insights into the role of non-canonical DNA structures and sequence motifs. PLoS One. 2013, 8: e59907-PubMed CentralPubMedView ArticleGoogle Scholar
- Hazel P, Huppert J, Balasubramanian S, Neidle S: Loop-length-dependent folding of G-quadruplexes. J Am Chem Soc. 2004, 126: 16405-16415.PubMedView ArticleGoogle Scholar
- Petraccone L, Erra E, Duro I, Esposito V, Randazzo A, Mayol L, Mattia CA, Barone G, Giancola C: Relative stability of quadruplexes containing different number of G-tetrads. Nucleosides Nucleotides Nucleic Acids. 2005, 24: 757-760.PubMedView ArticleGoogle Scholar
- Risitano A, Fox KR: Influence of loop size on the stability of intramolecular DNA quadruplexes. Nucleic Acids Res. 2004, 32: 2598-2606.PubMed CentralPubMedView ArticleGoogle Scholar
- Kikin O, D’antonio L, Bagga PS: QGRS Mapper: a web-based server for predicting G-quadruplexes in nucleotide sequences. Nucleic Acids Res. 2006, 34 (Web Server issue): W676-W682.PubMed CentralPubMedView ArticleGoogle Scholar
- Markham NR, Zuker M: UNAFold: software for nucleic acid folding and hybridization. Methods Mol Biol. 2008, 453: 3-31.PubMedView ArticleGoogle Scholar
- Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA: Circos: an information aesthetic for comparative genomics. Genome Res. 2009, 19: 1639-1645.PubMed CentralPubMedView ArticleGoogle Scholar
- Bhattacharjee AJ, Ahluwalia K, Taylor S, Jin O, Nicoludis JM, Buscaglia R, Brad Chaires J, Kornfilt DJP, Marquardt DGS, Yatsunyk LA: Induction of G-quadruplex DNA structure by Zn(II) 5,10,15,20-tetrakis(N-methyl-4-pyridyl)porphyrin. Biochimie. 2011, 93: 1297-1309.PubMedView ArticleGoogle Scholar
- Smith JS, Chen Q, Yatsunyk LA, Nicoludis JM, Garcia MS, Kranaster R, Balasubramanian S, Monchaud D, Teulade-Fichou M-P, Abramowitz L, Schultz DC, Johnson FB: Rudimentary G-quadruplex-based telomere capping in Saccharomyces cerevisiae. Nat Struct Mol Biol. 2011, 18: 478-485.PubMed CentralPubMedView ArticleGoogle Scholar
- Huppert JL: Hunting G-quadruplexes. Biochimie. 2008, 90: 1140-1148.PubMedView ArticleGoogle Scholar
- Yang J-N, Seluanov A, Gorbunova V: Mitochondrial inverted repeats strongly correlate with lifespan: mtDNA inversions and aging. PLoS One. 2013, 8: e73318-PubMed CentralPubMedView ArticleGoogle Scholar
- Moraes CT, DiMauro S, Zeviani M, Lombes A, Shanske S, Miranda AF, Nakase H, Bonilla E, Werneck LC, Servidei S: Mitochondrial DNA deletions in progressive external ophthalmoplegia and Kearns-Sayre syndrome. N Engl J Medi. 1989, 320: 1293-1299.View ArticleGoogle Scholar
- Guo X, Popadin KY, Markuzon N, Orlov YL, Kraytsberg Y, Krishnan KJ, Zsurka G, Turnbull DM, Kunz WS, Khrapko K: Repeats, longevity and the sources of mtDNA deletions: evidence from ‘deletional spectra’. Trends Genet. 2010, 26: 340-343.PubMed CentralPubMedView ArticleGoogle Scholar
- Damas J, Samuels DC, Carneiro J, Amorim A, Pereira F: Mitochondrial DNA rearrangements in health and disease–a comprehensive study. Hum Mutat. 2014, 35: 1-14.PubMedView ArticleGoogle Scholar
- Myers S, Freeman C, Auton A, Donnelly P, McVean G: A common sequence motif associated with recombination hot spots and genome instability in humans. Nat Genet. 2008, 40: 1124-1129.PubMedView ArticleGoogle Scholar
- Mergny J-L, Li J, Lacroix L, Amrane S, Chaires JB: Thermal difference spectra: a specific signature for nucleic acid structures. Nucleic Acids Res. 2005, 33: e138-PubMed CentralPubMedView ArticleGoogle Scholar
- Paeschke K, Capra JA, Zakian VA: DNA replication through G-quadruplex motifs is promoted by the Saccharomyces cerevisiae Pif1 DNA helicase. Cell. 2011, 145: 678-691.PubMed CentralPubMedView ArticleGoogle Scholar
- Paeschke K, Bochman ML, Garcia PD, Cejka P, Friedman KL, Kowalczykowski SC, Zakian VA: Pif1 family helicases suppress genome instability at G-quadruplex motifs. Nature. 2013, 497: 458-462.PubMed CentralPubMedView ArticleGoogle Scholar
- Kolesar JE, Wang CY, Taguchi YV, Chou S-H, Kaufman BA: Two-dimensional intact mitochondrial DNA agarose electrophoresis reveals the structural complexity of the mammalian mitochondrial genome. Nucleic Acids Res. 2013, 41: e58-e58.PubMed CentralPubMedView ArticleGoogle Scholar
- Pohjoismäki JLO, Holmes JB, Wood SR, Yang M-Y, Yasukawa T, Reyes A, Bailey LJ, Cluett TJ, Goffart S, Rigby RE, Jackson AP, Spelbrink JN, Griffith JD, Crouch RJ, Jacobs HT, Holt IJ: Mammalian mitochondrial DNA replication intermediates are essentially duplex but contain extensive tracts of RNA/DNA hybrid. J Mol Biol. 2010, 397: 1144-1155.PubMed CentralPubMedView ArticleGoogle Scholar
- Van Tuyle GC, Pavco PA: The rat liver mitochondrial DNA-protein complex: displaced single strands of replicative intermediates are protein coated. J Cell Biol. 1985, 100: 251-257.PubMedView ArticleGoogle Scholar
- Kim N, Jinks-Robertson S: Transcription as a source of genome instability. Nat Rev Genet. 2012, 13: 204-214.PubMed CentralPubMedGoogle Scholar
- Krishnan KJ, Reeve AK, Samuels DC, Chinnery PF, Blackwood JK, Taylor RW, Wanrooij S, Spelbrink JN, Lightowlers RN, Turnbull DM: What causes mitochondrial DNA deletions in human cells?. Nat Genet. 2008, 40: 275-279.PubMedView ArticleGoogle Scholar
- Srivastava S, Moraes CT: Double-strand breaks of mouse muscle mtDNA promote large deletions similar to multiple mtDNA deletions in humans. Hum Mol Genet. 2005, 14: 893-902.PubMed CentralPubMedView ArticleGoogle Scholar
- Fukui H, Moraes CT: Mechanisms of formation and accumulation of mitochondrial DNA deletions in aging neurons. Hum Mol Genet. 2009, 18: 1028-1036.PubMed CentralPubMedView ArticleGoogle Scholar
- Garone C, Rubio JC, Calvo SE, Naini A, Tanji K, DiMauro S, Mootha VK, Hirano M: MPV17 mutations causing adult-onset multisystemic disorder with multiple mitochondrial DNA DeletionsMPV17 mutations. Arch Neurol. 2012, 69: 1648-1651.PubMed CentralPubMedView ArticleGoogle Scholar
- Copeland WC: Defects in mitochondrial DNA replication and human disease. Crit Rev Biochem Mol Biol. 2012, 47: 64-74.PubMed CentralPubMedView ArticleGoogle Scholar
- Copeland WC: Inherited mitochondrial diseases of DNA replication. Annu Rev Med. 2008, 59: 131-146.PubMed CentralPubMedView ArticleGoogle Scholar
- Douarre C, Mergui X, Sidibe A, Gomez D, Alberti P, Mailliet P, Trentesaux C, Riou J-F: DNA damage signaling induced by the G-quadruplex ligand 12459 is modulated by PPM1D/WIP1 phosphatase. Nucleic Acids Res. 2013, 41: 3588-3599.PubMed CentralPubMedView ArticleGoogle Scholar
- Schwab RA, Nieminuszczy J, Shin-ya K, Niedzwiedz W: FANCJ couples replication past natural fork barriers with maintenance of chromatin structure. J Cell Biol. 2013, 201: 33-48.PubMed CentralPubMedView ArticleGoogle Scholar
- Spelbrink JN, Li FY, Tiranti V, Nikali K, Yuan QP, Tariq M, Wanrooij S, Garrido N, Comi G, Morandi L, Santoro L, Toscano A, Fabrizi GM, Somer H, Croxen R, Beeson D, Poulton J, Suomalainen A, Jacobs HT, Zeviani M, Larsson C: Human mitochondrial DNA deletions associated with mutations in the gene encoding Twinkle, a phage T7 gene 4-like protein localized in mitochondria. Nat Genet. 2001, 28: 223-231.PubMedView ArticleGoogle Scholar
- Jemt E, Farge G, Backstrom S, Holmlund T, Gustafsson CM, Falkenberg M: The mitochondrial DNA helicase TWINKLE can assemble on a closed circular template and support initiation of DNA synthesis. Nucleic Acids Res. 2011, 39: 9238-9249.PubMed CentralPubMedView ArticleGoogle Scholar
- Croteau DL, Rossi ML, Canugovi C, Tian J, Sykora P, Ramamoorthy M, Wang ZM, Singh DK, Akbari M, Kasiviswanathan R, Copeland WC, Bohr VA: RECQL4 localizes to mitochondria and preserves mitochondrial DNA integrity. Aging Cell. 2012, 11: 456-466.PubMed CentralPubMedView ArticleGoogle Scholar
- Kazak L, Reyes A, Duncan AL, Rorbach J, Wood SR, Brea-Calvo G, Gammage PA, Robinson AJ, Minczuk M, Holt IJ: Alternative translation initiation augments the human mitochondrial proteome. Nucleic Acids Res. 2013, 41: 2354-2369.PubMed CentralPubMedView ArticleGoogle Scholar
- George T, Wen Q, Griffiths R, Ganesh A, Meuth M, Sanders CM: Human Pif1 helicase unwinds synthetic DNA structures resembling stalled DNA replication forks. Nucleic Acids Res. 2009, 37: 6491-6502.PubMed CentralPubMedView ArticleGoogle Scholar
- Sanders CM: Human Pif1 helicase is a G-quadruplex DNA-binding protein with G-quadruplex DNA-unwinding activity. Biochem J. 2010, 430: 119-128.PubMedView 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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.