Twinkle, the Mitochondrial Replicative DNA Helicase, Is Widespread in the Eukaryotic Radiation and May Also Be the Mitochondrial DNA Primase in Most Eukaryotes

Recently, the human protein responsible for replicative mtDNA helicase activity was identified and designated Twinkle. Twinkle has been implicated in autosomal dominant progressive external ophthalmoplegia (adPEO), a mitochondrial disorder characterized by mtDNA deletions. The Twinkle protein appears to have evolved from an ancestor shared with the bifunctional primase-helicase found in the T-odd bacteriophages. However, the question has been raised as to whether human Twinkle possesses primase activity, due to amino acid sequence divergence and absence of a zinc-finger motif thought to play an integral role in DNA binding. To date, a primase protein participating in mtDNA replication has not been identified in any eukaryote. Here we investigate the wider phylogenetic distribution of Twinkle by surveying and analyzing data from ongoing EST and genome sequencing projects. We identify Twinkle homologues in representatives from five of six major eukaryotic assemblages (“supergroups”) and present the sequence of the complete Twinkle gene from two members of Amoebozoa, a supergroup of amoeboid protists at the base of the opisthokont (fungal/metazoan) radiation. Notably, we identify conserved primase motifs including the zinc finger in all Twinkle sequences outside of Metazoa. Accordingly, we propose that Twinkle likely serves as the primase as well as the helicase for mtDNA replication in most eukaryotes whose genome encodes it, with the exception of Metazoa. [Reviewing Editior: Dr. Rüdiger Cerff]     

Association of G-quadruplex forming sequences with human mtDNA deletion breakpoints

Background

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.

Results

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.

Conclusions

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.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-677) contains supplementary material, which is available to authorized users.     

FANCJ promotes DNA synthesis through G‐quadruplex structures

Our genome contains many G-rich sequences, which have the propensity to fold into stable secondary DNA structures called G4 or G-quadruplex structures. These structures have been implicated in cellular processes such as gene regulation and telomere maintenance. However, G4 sequences are prone to mutations particularly upon replication stress or in the absence of specific helicases. To investigate how G-quadruplex structures are resolved during DNA replication, we developed a model system using ssDNA templates and Xenopus egg extracts that recapitulates eukaryotic G4 replication. Here, we show that G-quadruplex structures form a barrier for DNA replication. Nascent strand synthesis is blocked at one or two nucleotides from the G4. After transient stalling, G-quadruplexes are efficiently unwound and replicated. In contrast, depletion of the FANCJ/BRIP1 helicase causes persistent replication stalling at G-quadruplex structures, demonstrating a vital role for this helicase in resolving these structures. FANCJ performs this function independently of the classical Fanconi anemia pathway. These data provide evidence that the G4 sequence instability in FANCJ^−/− cells and Fancj/dog1 deficient C. elegans is caused by replication stalling at G-quadruplexes.     

Replicating animal mitochondrial DNA

The field of mitochondrial DNA (mtDNA) replication has been experiencing incredible progress in recent years, and yet little is certain about the mechanism(s) used by animal cells to replicate this plasmid-like genome. The long-standing strand-displacement model of mammalian mtDNA replication (for which single-stranded DNA intermediates are a hallmark) has been intensively challenged by a new set of data, which suggests that replication proceeds via coupled leading- and lagging-strand synthesis (resembling bacterial genome replication) and/or via long stretches of RNA intermediates laid on the mtDNA lagging-strand (the so called RITOLS). The set of proteins required for mtDNA replication is small and includes the catalytic and accessory subunits of DNA polymerase γ, the mtDNA helicase Twinkle, the mitochondrial single-stranded DNA-binding protein, and the mitochondrial RNA polymerase (which most likely functions as the mtDNA primase). Mutations in the genes coding for the first three proteins are associated with human diseases and premature aging, justifying the research interest in the genetic, biochemical and structural properties of the mtDNA replication machinery. Here we summarize these properties and discuss the current models of mtDNA replication in animal cells.     

Mitochondrial DNA replication and repair defects: Clinical phenotypes and therapeutic interventions

Mitochondria is a unique cellular organelle involved in multiple cellular processes and is critical for maintaining cellular homeostasis. This semi-autonomous organelle contains its circular genome – mtDNA (mitochondrial DNA), that undergoes continuous cycles of replication and repair to maintain the mitochondrial genome integrity. The majority of the mitochondrial genes, including mitochondrial replisome and repair genes, are nuclear-encoded. Although the repair machinery of mitochondria is quite efficient, the mitochondrial genome is highly susceptible to oxidative damage and other types of exogenous and endogenous agent-induced DNA damage, due to the absence of protective histones and their proximity to the main ROS production sites. Mutations in replication and repair genes of mitochondria can result in mtDNA depletion and deletions subsequently leading to mitochondrial genome instability. The combined action of mutations and deletions can result in compromised mitochondrial genome maintenance and lead to various mitochondrial disorders. Here, we review the mechanism of mitochondrial DNA replication and repair process, key proteins involved, and their altered function in mitochondrial disorders. The focus of this review will be on the key genes of mitochondrial DNA replication and repair machinery and the clinical phenotypes associated with mutations in these genes.     

Mitochondrial DNA deletions are associated with non-B DNA conformations

Mitochondrial DNA (mtDNA) deletions are a primary cause of mitochondrial disease and are believed to contribute to the aging process and to various neurodegenerative diseases. Despite strong observational and experimental evidence, the molecular basis of the deletion process remains obscure. In this study, we test the hypothesis that the primary cause of mtDNA vulnerability to breakage resides in the formation of non-B DNA conformations, namely hairpin, cruciform and cloverleaf-like elements. Using the largest database of human mtDNA deletions built thus far (753 different cases), we show that site-specific breakage hotspots exist in the mtDNA. Furthermore, we discover that the most frequent deletion breakpoints occur within or near predicted structures, a result that is supported by data from transgenic mice with mitochondrial disease. There is also a significant association between the folding energy of an mtDNA region and the number of breakpoints that it harbours. In particular, two clusters of hairpins (near the D-loop 3'-terminus and the L-strand origin of replication) are hotspots for mtDNA breakage. Consistent with our hypothesis, the highest number of 5'- and 3'-breakpoints per base is found in the highly structured tRNA genes. Overall, the data presented in this study suggest that non-B DNA conformations are a key element of the mtDNA deletion process.     

Nanopore sequencing identifies a higher frequency and expanded spectrum of mitochondrial DNA deletion mutations in human aging

Mitochondrial DNA (mtDNA) deletion mutations cause many human diseases and are linked to age-induced mitochondrial dysfunction. Mapping the mutation spectrum and quantifying mtDNA deletion mutation frequency is challenging with next-generation sequencing methods. We hypothesized that long-read sequencing of human mtDNA across the lifespan would detect a broader spectrum of mtDNA rearrangements and provide a more accurate measurement of their frequency. We employed nanopore Cas9-targeted sequencing (nCATS) to map and quantitate mtDNA deletion mutations and develop analyses that are fit-for-purpose. We analyzed total DNA from vastus lateralis muscle in 15 males ranging from 20 to 81 years of age and substantia nigra from three 20-year-old and three 79-year-old men. We found that mtDNA deletion mutations detected by nCATS increased exponentially with age and mapped to a wider region of the mitochondrial genome than previously reported. Using simulated data, we observed that large deletions are often reported as chimeric alignments. To address this, we developed two algorithms for deletion identification which yield consistent deletion mapping and identify both previously reported and novel mtDNA deletion breakpoints. The identified mtDNA deletion frequency measured by nCATS correlates strongly with chronological age and predicts the deletion frequency as measured by digital PCR approaches. In substantia nigra, we observed a similar frequency of age-related mtDNA deletions to those observed in muscle samples, but noted a distinct spectrum of deletion breakpoints. NCATS-mtDNA sequencing allows the identification of mtDNA deletions on a single-molecule level, characterizing the strong relationship between mtDNA deletion frequency and chronological aging.     

Twinkle and POLG defects enhance age-dependent accumulation of mutations in the control region of mtDNA

Autosomal dominant and/or recessive progressive external ophthalmoplegia (ad/arPEO) is associated with mtDNA mutagenesis. It can be caused by mutations in three nuclear genes, encoding the adenine nucleotide translocator 1, the mitochondrial helicase Twinkle or DNA polymerase γ (POLG). How mutations in these genes result in progressive accumulation of multiple mtDNA deletions in post- mitotic tissues is still unclear. A recent hypothesis suggested that mtDNA replication infidelity could promote slipped mispairing, thereby stimulating deletion formation. This hypothesis predicts that mtDNA of ad/arPEO patients will contain frequent mutations throughout; in fact, our analysis of muscle from ad/arPEO patients revealed an age-dependent, enhanced accumulation of point mutations in addition to deletions, but specifically in the mtDNA control region. Both deleted and non-deleted mtDNA molecules showed increased point mutation levels, as did mtDNAs of patients with a single mtDNA deletion, suggesting that point mutations do not cause multiple deletions. Deletion breakpoint analysis showed frequent breakpoints around homopolymeric runs, which could be a signature of replication stalling. Therefore, we propose replication stalling as the principal cause of deletion formation.     

Werner syndrome protein suppresses the formation of large deletions during the replication of human telomeric sequences

Werner syndrome (WS) is a disorder characterized by features of premature aging and increased cancer that is caused by loss of the RecQ helicase WRN. Telomeres consisting of duplex TTAGGG repeats in humans protect chromosome ends and sustain cellular proliferation. WRN prevents the loss of telomeres replicated from the G-rich strand, which can form secondary G-quadruplex (G4) structures. Here, we dissected WRN roles in the replication of telomeric sequences by examining factors inherent to telomeric repeats, such as G4 DNA, independently from other factors at chromosome ends that can also impede replication. For this we used the supF shuttle vector (SV) mutagenesis assay. We demonstrate that SVs with [TTAGGG]₆ sequences are stably replicated in human cells, and that the repeats suppress the frequency of large deletions despite G4 folding potential. WRN depletion increased the supF mutant frequency for both the telomeric and non-telomeric SVs, compared with the control cells, but this increase was much greater (27-fold) for telomeric SVs. The higher SV mutant frequencies in WRN-deficient cells were primarily due to an increase in large sequence deletions and rearrangements. However, WRN depletion caused a more dramatic increase in deletions and rearrangements arising within the telomeric SV (70-fold), compared with non-telomeric SV (8-fold). Our results indicate that WRN prevents large deletions and rearrangements during replication, and that this role is particularly important in templates with telomeric sequence. This provides a possible explanation for increased telomere loss in WS cells.     

Werner syndrome protein suppresses the formation of large deletions during the replication of human telomeric sequences

Werner syndrome (WS) is a disorder characterized by features of premature aging and increased cancer that is caused by loss of the RecQ helicase WRN. Telomeres consisting of duplex TTAGGG repeats in humans protect chromosome ends and sustain cellular proliferation. WRN prevents the loss of telomeres replicated from the G-rich strand, which can form secondary G-quadruplex (G4) structures. Here, we dissected WRN roles in the replication of telomeric sequences by examining factors inherent to telomeric repeats, such as G4 DNA, independently from other factors at chromosome ends that can also impede replication. For this we used the supF shuttle vector (SV) mutagenesis assay. We demonstrate that SVs with [TTAGGG]₆ sequences are stably replicated in human cells, and that the repeats suppress the frequency of large deletions despite G4 folding potential. WRN depletion increased the supF mutant frequency for both the telomeric and non-telomeric SVs, compared with the control cells, but this increase was much greater (27-fold) for telomeric SVs. The higher SV mutant frequencies in WRN-deficient cells were primarily due to an increase in large sequence deletions and rearrangements. However, WRN depletion caused a more dramatic increase in deletions and rearrangements arising within the telomeric SV (70-fold), compared with non-telomeric SV (8-fold). Our results indicate that WRN prevents large deletions and rearrangements during replication, and that this role is particularly important in templates with telomeric sequence. This provides a possible explanation for increased telomere loss in WS cells.     

Mutagenic capacity of endogenous G4 DNA underlies genome instability in FANCJ-defective C. elegans

To safeguard genetic integrity, cells have evolved an accurate but not failsafe mechanism of DNA replication. Not all DNA sequences tolerate DNA replication equally well [1]. Also, genomic regions that impose structural barriers to the DNA replication fork are a potential source of genetic instability [2, 3]. Here, we demonstrate that G4 DNA-a sequence motif that folds into quadruplex structures in vitro [4, 5]-is highly mutagenic in vivo and is removed from genomes that lack dog-1, the C. elegans ortholog of mammalian FANCJ [6, 7], which is mutated in Fanconi anemia patients [8-11]. We show that sequences that match the G4 DNA signature G3-5N1-3G3-5N1-3G3-5N1-3G3-5 are deleted in germ and somatic tissues of dog-1 animals. Unbiased aCGH analyses of dog-1 genomes that were allowed to accumulate mutations in >100 replication cycles indicate that deletions are found exclusively at G4 DNA; deletion frequencies can reach 4% per site per animal generation. We found that deletion sizes fall short of Okazaki fragment lengths [12], and no significant microhomology was observed at deletion junctions. The existence of 376,000 potentially mutagenic G4 DNA sites in the human genome could have major implications for the etiology of hereditary FancJ and nonhereditary cancers.     

Mitochondrial DNA deletions and differential mitochondrial DNA content in Rhesus monkeys: Implications for aging

The purpose of this study was to determine the relationship between mitochondrial DNA (mtDNA) deletions, mtDNA content and aging in rhesus monkeys. Using 2 sets of specific primers, we amplified an 8 kb mtDNA fragment covering a common 5.7 kb deletion and the entire 16.5 kb mitochondrial genome in the brain and buffy-coats of young and aged monkeys. We studied a total of 66 DNA samples: 39 were prepared from a buffy-coat and 27 were prepared from occipital cortex tissues. The mtDNA data were assessed using a permutation test to identify differences in mtDNA, in the different monkey groups. Using real-time RT-PCR strategy, we also assessed both mtDNA and nuclear DNA levels for young, aged and male and female monkeys. We found a 5.7 kb mtDNA deletion in 81.8% (54 of 66) of the total tested samples. In the young group of buffy-coat DNA, we found 5.7 kb deletions in 7 of 17 (41%), and in the aged group, we found 5.7 kb deletions in 12 of 22 (54%), suggesting that the prevalence of mtDNA deletions is related to age. We found decreased mRNA levels of mtDNA in aged monkeys relative to young monkeys. The increases in mtDNA deletions and mtDNA levels in aged rhesus monkeys suggest that damaged DNA accumulates as rhesus monkeys age and these altered mtDNA changes may have physiological relevance to compensate decreased mitochondrial function.     

The <Emphasis Type="Italic">Dictyostelium discoideum</Emphasis> homologue of Twinkle,Twm1, is a mitochondrial DNA helicase,an active primase and promotes mitochondrial DNA replication

Background

DNA replication requires contributions from various proteins, such as DNA helicases; in mitochondria Twinkle is important for maintaining and replicating mitochondrial DNA. Twinkle helicases are predicted to also possess primase activity, as has been shown in plants; however this activity appears to have been lost in metazoans. Given this, the study of Twinkle in other organisms is required to better understand the evolution of this family and the roles it performs within mitochondria.

Results

Here we describe the characterization of a Twinkle homologue, Twm1, in the amoeba Dictyostelium discoideum, a model organism for mitochondrial genetics and disease. We show that Twm1 is important for mitochondrial function as it maintains mitochondrial DNA copy number in vivo. Twm1 is a helicase which unwinds DNA resembling open forks, although it can act upon substrates with a single 3′ overhang, albeit less efficiently. Furthermore, unlike human Twinkle, Twm1 has primase activity in vitro. Finally, using a novel in bacterio approach, we demonstrated that Twm1 promotes DNA replication.

Conclusions

We conclude that Twm1 is a replicative mitochondrial DNA helicase which is capable of priming DNA for replication. Our results also suggest that non-metazoan Twinkle could function in the initiation of mitochondrial DNA replication. While further work is required, this study has illuminated several alternative processes of mitochondrial DNA maintenance which might also be performed by the Twinkle family of helicases.

     

Two direct repeats cause most human mtDNA deletions

Mitochondrial DNA (mtDNA) deletions are a common cause of human mitochondrial disease and also occur as part of normal aging. However, it is unknown how the deletions actually occur. To gain further insight, we studied the sequences that flank 263 different human mtDNA deletions. The distribution of deletion breakpoints did not correspond to the basic parameters of wild-type mtDNA that are thought to predispose to deletion formation. But there was a striking correspondence to the position of two 13-bp direct repeats beginning at nucleotides 8470 and 13 447. The vast majority of different mtDNA deletions appear to be related to these two repeats, suggesting a common mechanism related to mtDNA replication.     

The Yeast Pif1 Helicase Prevents Genomic Instability Caused by G-Quadruplex-Forming CEB1 Sequences In Vivo

In budding yeast, the Pif1 DNA helicase is involved in the maintenance of both nuclear and mitochondrial genomes, but its role in these processes is still poorly understood. Here, we provide evidence for a new Pif1 function by demonstrating that its absence promotes genetic instability of alleles of the G-rich human minisatellite CEB1 inserted in the Saccharomyces cerevisiae genome, but not of other tandem repeats. Inactivation of other DNA helicases, including Sgs1, had no effect on CEB1 stability. In vitro, we show that CEB1 repeats formed stable G-quadruplex (G4) secondary structures and the Pif1 protein unwinds these structures more efficiently than regular B-DNA. Finally, synthetic CEB1 arrays in which we mutated the potential G4-forming sequences were no longer destabilized in pif1Δ cells. Hence, we conclude that CEB1 instability in pif1Δ cells depends on the potential to form G-quadruplex structures, suggesting that Pif1 could play a role in the metabolism of G4-forming sequences.     

A257T linker region mutant of T7 helicase-primase protein is defective in DNA loading and rescued by T7 DNA polymerase

The helicase and primase activities of the hexameric ring-shaped T7 gp4 protein reside in two separate domains connected by a linker region. This linker region is part of the subunit interface between monomers, and point mutations in this region have deleterious effects on the helicase functions. One such linker region mutant, A257T, is analogous to the A359T mutant of the homologous human mitochondrial DNA helicase Twinkle, which is linked to diseases such as progressive external opthalmoplegia. Electron microscopy studies show that A257T gp4 is normal in forming rings with dTTP, but the rings do not assemble efficiently on the DNA. Therefore, A257T, unlike the WT gp4, does not preassemble on the unwinding DNA substrate with dTTP without Mg(II), and its DNA unwinding activity in ensemble assays is slow and limited by the DNA loading rate. Single molecule assays measured a 45 times slower rate of A257T loading on DNA compared with WT gp4. Interestingly, once loaded, A257T has almost WT-like translocation and DNA unwinding activities. Strikingly, A257T preassembles stably on the DNA in the presence of T7 DNA polymerase, which restores the ensemble unwinding activity of A257T to ～75% of WT, and the rescue does not require DNA synthesis. The DNA loading rate of A257T, however, remains slow even in the presence of the polymerase, which explains why A257T does not support T7 phage growth. Similar types of defects in the related human mitochondrial DNA helicase may be responsible for inefficient DNA replication leading to the disease states.     

Somatic mtDNA Mutation Spectra in the Aging Human Putamen

The accumulation of heteroplasmic mitochondrial DNA (mtDNA) deletions and single nucleotide variants (SNVs) is a well-accepted facet of the biology of aging, yet comprehensive mutation spectra have not been described. To address this, we have used next generation sequencing of mtDNA-enriched libraries (Mito-Seq) to investigate mtDNA mutation spectra of putamen from young and aged donors. Frequencies of the “common” deletion and other “major arc” deletions were significantly increased in the aged cohort with the fold increase in the frequency of the common deletion exceeding that of major arc deletions. SNVs also increased with age with the highest rate of accumulation in the non-coding control region which contains elements necessary for translation and replication. Examination of predicted amino acid changes revealed a skew towards pathogenic SNVs in the coding region driven by mutation bias. Levels of the pathogenic m.3243A>G tRNA mutation were also found to increase with age. Novel multimeric tandem duplications that resemble murine control region multimers and yeast ρ⁻ mtDNAs, were identified in both young and aged specimens. Clonal ∼50 bp deletions in the control region were found at high frequencies in aged specimens. Our results reveal the complex manner in which the mitochondrial genome alters with age and provides a foundation for studies of other tissues and disease states.     

The mitochondrial DNA helicase TWINKLE can assemble on a closed circular template and support initiation of DNA synthesis

Mitochondrial DNA replication is performed by a simple machinery, containing the TWINKLE DNA helicase, a single-stranded DNA-binding protein, and the mitochondrial DNA polymerase γ. In addition, mitochondrial RNA polymerase is required for primer formation at the origins of DNA replication. TWINKLE adopts a hexameric ring-shaped structure that must load on the closed circular mtDNA genome. In other systems, a specialized helicase loader often facilitates helicase loading. We here demonstrate that TWINKLE can function without a specialized loader. We also show that the mitochondrial replication machinery can assemble on a closed circular DNA template and efficiently elongate a DNA primer in a manner that closely resembles initiation of mtDNA synthesis in vivo.